ISM（Improved String Matching）可达矩阵的matlab程序

ISM（Improved String Matching）算法可以用来求解可达矩阵（Reachability Matrix），可达矩阵是指一个有向图中，从一个节点出发是否可以到达另一个节点，其结果是一个布尔类型的矩阵。下面是一个用MATLAB实现ISM算法求解可达矩阵的程序：

function [R] = reachability_matrix(adj_matrix)
% adj_matrix: adjacency matrix of the directed graph
% R: reachability matrix of the directed graph

n = size(adj_matrix, 1);
R = adj_matrix;
for k = 1:n
    for i = 1:n
        for j = 1:n
            R(i, j) = R(i, j) || (R(i, k) && R(k, j));
        end
    end
end

以上程序中，adj_matrix是有向图的邻接矩阵，R是可达矩阵。程序使用三重循环来实现ISM算法，其中第一重循环用于遍历中间节点k，第二重和第三重循环用于更新可达矩阵R。具体地，对于每一对节点(i,j)，如果节点i可以到达节点k并且节点k可以到达节点j，则节点i可以到达节点j。最终得到的可达矩阵就是图中任意两个节点是否可达的结果。