Havel-Hakimi定理及其应用
Havel-Hakimi定理:由非负整数组成的非增序列s:A[0], A[1],.....,A[n](n>=2,A[0]>=1)是可图的,当且仅当序列s1:A[1]-1,A[2]-1, A[A[0]]-1,A[A[0]+1],.....A[n]是可图的。
Havel-Hakimi定理实际上给出了一个序列s构图方法:把序列s按照非递增顺序排好序后,从中选出度数最大的顶点(设为v1),将它与度数次大的前A[0]个顶点连边,然后这个顶点就可以不管了,即在序列中删除首项A[0]后,并把后面 的A[0]个度数减1;再把剩下的序列重新按非增序列排序。按照上述过程连边。。。直到建出完整的图或者出现负度数等明显不合理的情况为止。
应用:http://poj.org/problem?id=1659
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int M = 15;
struct Point {
int val, no;
};
Point point[M];
int edge[M][M];
void print(int t, int n) {
int i;
for(i=t; i<n; ++i) printf("%d\t", point[i].val);
printf("\n");
}
bool cmp(Point pa, Point pb) {
return pa.val > pb.val;
}
bool havelHakimi(int n) {
bool flag = true;
sort(point, point+n, cmp);
//print(0, n);
int i, j;
for(i=0; i<n; ++i) {
int t = point[i].val;
for(j=i+1; j<=i+t; ++j) {
point[j].val--;
if(point[j].val < 0) {flag = false; break;}
edge[point[i].no][point[j].no] = edge[point[j].no][point[i].no] = 1;
}
sort(point+i+1, point+n, cmp);
//print(i+1, n);
if(false == flag) break;
}
return flag;
}
int main() {
int t, n;
int i, j;
//freopen("/home/lsy/Desktop/2.txt", "r", stdin);
scanf("%d", &t);
while(t--) {
scanf("%d", &n);
for(i=0; i<n; ++i) {
scanf("%d", &point[i].val);
point[i].no = i;
}
memset(edge, 0, sizeof(edge));
if(false == havelHakimi(n)) printf("NO\n");
else {
printf("YES\n");
for(i=0; i<n; ++i) {
for(j=0; j<n; ++j) {
if(0 != j) printf(" ");
printf("%d", edge[i][j]);
}
printf("\n");
}
}
if(t) printf("\n");
}
return 0;
}