Time Limit: 5000MS | Memory Limit: 10000K | |||
Total Submissions: 6076 | Accepted: 2636 | Special Judge |
Description
未名湖附近共有N个大小湖泊L1, L2, ..., Ln(其中包括未名湖),每个湖泊Li里住着一只青蛙Fi(1 ≤i ≤ N)。如果湖泊Li和Lj之间有水路相连,则青蛙Fi和Fj互称为邻居。现在已知每只青蛙的邻居数目x1,x2, ..., xn,请你给出每两个湖泊之间的相连关系。
Input
第一行是测试数据的组数T(0 ≤ T ≤ 20)。每组数据包括两行,第一行是整数N(2 < N < 10),第二行是N个整数,x1,x2,..., xn(0 ≤ xi ≤ N)。
Output
对输入的每组测试数据,如果不存在可能的相连关系,输出"NO"。否则输出"YES",并用N×N的矩阵表示湖泊间的相邻关系,即如果湖泊i与湖泊j之间有水路相连,则第i行的第j个数字为1,否则为0。每两个数字之间输出一个空格。如果存在多种可能,只需给出一种符合条件的情形。相邻两组测试数据之间输出一个空行。
Sample Input
3
7
4 3 1 5 4 2 1
6
4 3 1 4 2 0
6
2 3 1 1 2 1
Sample Output
YES
0 1 0 1 1 0 1
1 0 0 1 1 0 0
0 0 0 1 0 0 0
1 1 1 0 1 1 0
1 1 0 1 0 1 0
0 0 0 1 1 0 0
1 0 0 0 0 0 0 NO YES
0 1 0 0 1 0
1 0 0 1 1 0
0 0 0 0 0 1
0 1 0 0 0 0
1 1 0 0 0 0
0 0 1 0 0 0
Source
havel定理:
给定一个非负整数序列{dn},若存在一个无向图使得图中各点的度与此序列一一对应,则称此序列可图化。
进一步,若图为简单图,则称此序列可简单图化。
判定过程:(1)对当前数列排序,使其呈递减,(2)从S【2】开始对其后S【1】个数字-1,(3)
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,到这一步出现了负数,因此
该序列是不可图的。
havel定理的应用:
对于一个给定的度序列,看能不能形成一个简单无向图。
Havel算法的思想简单的说如下:
(1)对序列从大到小进行排序。
(2)设最大的度数为t,把最大的度数置0,然后把最大度数后(不包括自己)的t个度数分别减1(意思就是
把度数最大的点与后几个点进行连接)
(3)如果序列中出现了负数,证明无法构成。如果序列全部变为0,证明能构成,跳出循环.前两点不出现,
就跳回第一步!
感想:
看到越来越多转化为图连边来分析问题的做法,好奇妙啊好奇妙~
代码:
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std; int f[12][12];
struct node
{
int degree,flag;
}a[12];
bool cmp(node x,node y)
{
return x.degree>y.degree;
}
int main()
{
int T,i,j,n;
scanf("%d",&T);
while(T--)
{
memset(f,0,sizeof(f));
scanf("%d",&n);
for(i=1;i<=n;i++)
{
scanf("%d",&a[i].degree);
a[i].flag=i;
}
bool succ=1;
while(1)
{
sort(a+1,a+n+1,cmp);
if(a[1].degree==0)
break;
for(i=2;i<=a[1].degree+1;i++)
{
a[i].degree--;
f[a[1].flag][a[i].flag]=1;
f[a[i].flag][a[1].flag]=1;
if(a[i].degree<0)
{
succ=0;
break;
}
}
if(!succ) break;
a[1].degree=0;
}
if(!succ) puts("NO\n");
else
{
puts("YES");
for(i=1;i<=n;i++)
{
for(j=1;j<=n;j++)
printf("%d%c",f[i][j],j==n?'\n':' ');
}
puts("");
}
}
return 0;
}
12014030 |
Accepted |
164K |
0MS |
1306B |
2013-08-20 14:14:16 |