ZOJ-1015-Fishing net

题意：给定一个图。判断是不是弦图？

思路：（1）神马是弦图？对于一个无向图，若该图的任意一个长度大于3的环中至少存在一条边连接这个环上不相邻的两点，则此图称作弦图。

（2）什么是团？团是原图的一个子图，子图就是包含了原图的某些点，那么就要包含这些点之间的边。并且团不是一般的子图而是一个完全子图，就是这个子图的任意两个顶点之间都有边。

（3）完美消除序列：原图的一个点的序列(每个点出现且恰好出现一次)v1, v2,……, vn满足{vi, vi+1,…,vn}的组成的子图为团。

（4）一个无向图是弦图当且仅当它有一个完美消除序列。

（5）如何计算完美消除序列？最大势算法： 从n到1的顺序依次给点标号(标号为i的点出现在完美消除序列的第i个)。 设label[i]表示第i个点与多少个已标号的点相 邻，每次选择label[i]最大的未标号的点进行标号。注意这里只是计算出了完美消除序列，但是在求出这个之后还没有判定是不是弦图。

（6）如何从完美消除序列判断原图是不是弦 图？最朴素的办法是依次判断 {vi+1,…,vn}中所有与vi相邻的点是否构成了一个团。可以这样优化：设{vi+1,…,vn}中所有与vi相邻的点依次为 vj1,……,vjk。只需判断vj1是否与vj2,……,vjk相邻即可。

最大势算法

维护每个结点的“势”。并维护n个链表，分别表示势为0~n-1的结点各有哪些。
维护一个best值表示当前最大的势。
再开一个bool数组记录每个点是否已被取出放入序列。

初始时每个点势为0，都放在0号链表里。best=0。
①更新操作：某结点u的势从x变为x+1时，只要把u插到x+1链表的头部，无需将其从x链表中删除（所以不用写双向链表）。然后更新best值。
②取最大操作：取最大势结点时从best链表的头开始取，如果取到的是已被放入序列的结点就删除并继续往后取，仍取不到则–best，直到取到一个未被放入序列的点。
操作①每次O(1)，总共执行O(m)次。best++也被执行O(m)次，那么操作②中–best也只有O(m)次。操作①每次插一个结点，则总共有O(n+m)个结点，而操作②中每往后取一个就删一个，所以也是O(n+m)。

时间复杂度为O(n+m)的算法，n是图的点数，m是图的边数。
第一步：给节点编号
设已编号的节点集合为A，未编号的节点集合为B
开始时A为空，B包含所有节点。
for num=n-1 downto 0 do
{
在B中找节点x，使与x相邻的在A集合中的节点数最多，将x编号为num，
并从B移入A
}
第二步：检查
for num=0 to n-1 do
{
对编号为num的节点x，设所有编号大于num且与x相邻的节点集合为C，
在集合C中找出编号最小的节点y，如果集合C中存在不等于y的节点z，
且y与z间没有边，则此图不是弦图，退出。
}
检查完了，则此图是弦图

( scanf("%d %d", &n, &m), n+m) – 这是“逗号”表达式。
第一个 表达式 是 scanf("%d%d", &n, &m),
第二个 表达式 是 n+m

“逗号”表达式 运算顺序从左到右依表达式次序执行，
但是 ，“逗号”表达式 的结果 等于 最后一个 表达式的结果。

#include <cstdio>
#include <cstring>
using namespace std;

/**
初始化所有节点的势为0，然后每次取出势最大的节点并为之标号(第几个取出的)，并给它的所有相邻节点的势 +1s加1，重复直到所有点被取出。 
然而每个图用最大势算法都能得到一个序列。所以，还要加入判断是否为完美消除序列的步骤。 
大概是，对每个点选中所有比它的标号大的点，判断其中标号最小的一个是否与其它选中点都相邻，不是则不是完美消除序列。
*/

int count,visited[1024],label[1024],set[1024],g[1024][1024];
void Relabel(int n)
{
    int i,j,id;
    memset(visited,0,sizeof(visited));
    memset(label,0,sizeof(label));
    count=0;
    for(i=n;i>=1;i--)//从后往前确定序列的点
    {
        id=1;
        for(j=1;j<=n;j++)
        {
            if(!visited[j]&&label[id]<label[j])//每次取点选择标号最大的点之一
            {
                id=j;
            }
        }
        set[i]=id;
        visited[id]=1;
        for(j=1;j<=n;j++)//每取一个点都把还没加入序列的和它相连的点的标号+1
        {
            if(id!=j&&g[id][j]&&!visited[j])
                label[j]++;
        }
    }
}
int check(int n)
{
    int i,j;
    for(i=1;i<=n;i++)//对每个点选中所有比它的标号大的点，判断其中标号最小的一个是否与其它选中点都相邻 ，不是则不是完美消除序列 
    {
        int id;
        for(j=i+1;j<=n;j++)
        {
            if(g[set[i]][set[j]])//与选中的点相邻，标号最小的 
            {
                id=j;
                break;
            }
        }
        for(j=i+1;j<=n;j++)//标号最小的，是否与选中的点的其他相邻点都相邻 
        {
            if(g[set[i]][set[j]]&&id!=j&&g[set[id]][set[j]]==0)
            return 0;
        }
    }
    return 1;
 
}
int main()
{
	//freopen("tmp.txt", "r", stdin); 
    int n,m,i;
    while(scanf("%d%d",&n,&m),m+n)
    {
        memset(g,0,sizeof(g));
        for(i=0;i<m;i++)
        {
            int a,b;
            scanf("%d%d",&a,&b);
            g[a][b]=g[b][a]=1;
        }
        Relabel(n);
        if(check(n))
        {
            printf("Perfect\n\n");
        }
        else
            printf("Imperfect\n\n");
    }
    return 0;
}