Havel-Hakimi定理及其应用

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;
}