弦图判断

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}