ZOJ 3732 Graph Reconstruction Havel_Hakimi定理

链接：http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=5078

题意：有N个点组成无向图，每个点的度为ai，问是否能组成图，并且组成的图方式是否唯一。

思路：Havel_Hakimi定理的应用。

（以下转载）

1，Havel-Hakimi定理主要用来判定一个给定的序列是否是可图的。

2，首先介绍一下度序列：若把图 G 所有顶点的度数排成一个序列 S，则称 S 为图 G 的度序列。

3，一个非负整数组成的有限序列如果是某个无向图的序列，则称该序列是可图的。

4，判定过程：（1）对当前数列排序，使其呈递减，（2）从S【2】开始对其后S【1】个数字-1，（3）一直循环直到当前序列出现负数（即不是可图的情况）或者当前序列全为0 （可图）时退出。

5，举例：序列S：7,7,4,3,3,3,2,1  删除序列S的首项 7 ，对其后的7项每项减1，得到：6,3,2,2,2,1,0，继续删除序列的首项6，对其后的6项每项减1，得到：2,1,1,1,0，-1，到这一步出现了负数，因此该序列是不可图的。

（转载结束）

判断是否成图方式唯一的方法是在每个点寻找连接的边的时候，找到的最后一个点，如果与它之后第一个度不为零的点的度是相同的，则两个点可以交换位置并且不影响建图可能性，而出现至少两种建图方式。

代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <map>
#include <cstdlib>
#include <queue>
#include <stack>
#include <vector>
#include <ctype.h>
#include <algorithm>
#include <string>
#include <set>
#define PI acos(-1.0)
#define maxn 10005
#define INF 0x7fffffff
#define eps 1e-8
typedef long long LL;
typedef unsigned long long ULL;
using namespace std;
struct aa
{
int num;
int degree;
} a[105],b[105];
int way[10005][2];
bool cmp1(aa a,aa b)
{
return a.degree>b.degree;
}
int main()
{
int tot;
while(~scanf("%d",&tot))
{
int top=0,A,B;
bool flag1=0,flag2=0;
for(int i=0; i<tot; i++)
{
a[i].num=b[i].num=i;
scanf("%d",&a[i].degree);
b[i].degree=a[i].degree;
}
for(int i=0; i<tot; i++)
{
sort(b+i,b+tot,cmp1);
if(i+b[i].degree>=tot)
{
flag1=1;
break;
}
int tt=0;
for(int j=i+1; j<tot; j++)
{
if(!b[j].degree)
continue;
tt++;
way[top][0]=b[i].num;
way[top++][1]=b[j].num;
if(tt==b[i].degree)
{
int p=j+1;
while(b[p].degree==0&&p<tot)
p++;
if(p<tot&&b[p].degree==b[j].degree)
flag2=1;
}
b[j].degree--;
if(b[j].degree<0)
{
flag1=1;
break;
}
if(tt==b[i].degree)
break;
}
if(tt<b[i].degree)
flag1=1;
if(flag1==1)
break;
}
if(!flag1&&!flag2)
{
printf("UNIQUE\n");
printf("%d %d\n",tot,top);
for(int i=0; i<top; i++)
{
if(i==0)
printf("%d",way[i][0]+1);
else printf(" %d",way[i][0]+1);
}
printf("\n");
for(int i=0; i<top; i++)
{
if(i==0)
printf("%d",way[i][1]+1);
else printf(" %d",way[i][1]+1);
}
printf("\n");
}
else if(flag1)
printf("IMPOSSIBLE\n");
else
{
printf("MULTIPLE\n");
printf("%d %d\n",tot,top);
for(int i=0; i<top; i++)
{
if(i==0)
printf("%d",way[i][0]+1);
else printf(" %d",way[i][0]+1);
}
printf("\n");
for(int i=0; i<top; i++)
{
if(i==0)
printf("%d",way[i][1]+1);
else printf(" %d",way[i][1]+1);
}
printf("\n");
printf("%d %d\n",tot,top);
top=0;
flag2=0;
for(int i=0; i<tot; i++)
b[i]=a[i];
for(int i=0; i<tot; i++)
{
sort(b+i,b+tot,cmp1);
if(i+b[i].degree>=tot)
{
flag1=1;
break;
}
int tt=0;
for(int j=i+1; j<tot; j++)
{
if(!b[j].degree)
continue;
tt++;
if(tt==b[i].degree&&flag2==0)
{
int p=j+1;
while(b[p].degree==0&&p<tot)
p++;
if(p<tot&&b[p].degree==b[j].degree)
{
flag2=1;
j=p;
}
}
way[top][0]=b[i].num;
way[top++][1]=b[j].num;
b[j].degree--;
if(tt==b[i].degree)
break;
}
}
for(int i=0; i<top; i++)
{
if(i==0)
printf("%d",way[i][0]+1);
else printf(" %d",way[i][0]+1);
}
printf("\n");
for(int i=0; i<top; i++)
{
if(i==0)
printf("%d",way[i][1]+1);
else printf(" %d",way[i][1]+1);
}
printf("\n");
}
}
return 0;
}