ACM模版

弦图判断

/* * 弦图判断 * INIT: g[][]置为邻接矩阵; * CALL: mcs(n); peo(n); * 第一步: 给节点编号 mcs(n) * 设已编号的节点集合为A, 未编号的节点集合为B * 开始时A为空, B包含所有节点. * for num=n-1 downto 0 do { * 在B中找节点x, 使与x相邻的在A集合中的节点数最多, * 将x编号为num，并从B移入A。 * } * 第二部：检查peo(n) * for num=0 to n-1 do { * 对编号为num的点x，设所有编号>num且与x相邻的点集为C * 在C中找出编号最小的节点y， * 若C中存在x != y，使得y与x之间无边，则此图不是弦图。 * } * 检查完毕, 则此图是弦图. */
const int V = 10010;
int g[V][V], order[V], inv[V], tag[V];

void mcs(int n)
{
    int i, j, k;
    memset(tag, 0, sizeof(tag));
    memset(order, -1, sizeof(order));
    for (i = n - 1; i >= 0; i--)
    {   // vertex:0~n-1
        for (j = 0; order[j] >= 0; j++);
        for (k = j + 1; k < n; k++)
        {
            if (order[k] < 0 && tag[k] > tag[j])
            {
                j = k;
            }
        }
        order[j] = i, inv[i] = j;
        for (k = 0; k < n; k++)
        {
            if (g[j][k])
            {
                tag[k]++;
            }
        }
    }
    return ;
}

int peo(int n)
{
    int i, j, k, w, min;
    for (i = n - 2; i >= 0; i--)
    {
        j = inv[i], w = -1, min = n;
        for (k = 0; k < n; k++)
        {
            if (g[j][k] && order[k] > order[j] && order[k] < min)
            {
                min = order[k], w=k;
            }
        }
        if (w < 0)
        {
            continue;
        }
        for (k = 0; k < n; k++)
        {
            if (g[j][k] && order[k] > order[w] && !g[k][w])
            {
                return 0;   // no
            }
        }
    }
    return 1;               // yes
}