题目

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

题解

Havel-Hakimi定理

给定一个非负整数序列{dn},若存在一个无向图使得图中各点的度与此序列一一对应,则称此序列可图化
进一步,若图为简单图,则称此序列可简单图化
可图化的判定:d1+d2+……dn=0(mod 2)
关于具体图的构造,我们可以简单地把奇数度的点配对,剩下的全部搞成自环

可简单图化的判定(Havel定理):把序列排成不增序,即d1>=d2>=……>=dn,则d可简单图化当且仅当d’={d2-1,d3-1,……d(d1+1)-1, d(d1+2),d(d1+3),……dn}可简单图化
简单的说,把d排序后,找出度最大的点(设度为d1),把它与度次大的d1个点之间连边,然后这个点就可以不管了,一直继续这个过程,直到建出完整的图,或出现负度等明显不合理的情况

简单来说就是从多到少开始尽可能多的连接
用一个 优先队列 维护当前状态下还能连接的个数最多的点

虽然当时不知道这个定理,但是感觉这样是对的……

代码

/*
By:OhYee
Github:OhYee
Blog:http://www.oyohyee.com/
Email:oyohyee@oyohyee.com

かしこいかわいい?
エリーチカ!
要写出来Хорошо的代码哦~
*/
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <string>
#include <iostream>
#include <vector>
#include <list>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <functional>
#include <bitset>
using namespace std;

const int maxn = 15;

struct naive {
    int n;
    int f;
    bool operator < (const naive &rhs)const {
        return f < rhs.f;
    }
};

bool Map[maxn][maxn];

priority_queue <naive> Q1;
queue<naive> Q2;

void Do() {
    int n;
    cin >> n;

    while(!Q1.empty()) Q1.pop();
    while(!Q2.empty()) Q2.pop();
    bool flag = false;
    memset(Map,false,sizeof(Map));

    for(int i = 0;i < n;i++) {
        naive t;
        cin >> t.f;
        t.n = i;
        Q1.push(t);
    }

    while(!Q1.empty()) {
        naive qu = Q1.top();
        Q1.pop();
        while(!Q1.empty() && qu.f > 0) {
            naive qv = Q1.top();
            Q1.pop();
            if(qv.f == 0)
                continue;

            int u = qu.n;
            int v = qv.n;
            if(Map[u][v] == false) {
                Map[u][v] = Map[v][u] = true;
                qu.f--;
                qv.f--;
            }
            if(qv.f > 0)
                Q2.push(qv);
        }
        if(qu.f > 0) {
            flag = true;
            break;
        }
        while(!Q2.empty()) {
            Q1.push(Q2.front());
            Q2.pop();
        }
    }

    if(flag) {
        cout << "NO" << endl;
    } else {
        cout << "YES" << endl;
        for(int i = 0;i < n;i++) {
            for(int j = 0;j < n;j++) {
                if(j != 0)
                    cout << " ";
                cout << Map[i][j] ? 1 : 0;

            }
            cout << endl;
        }
    }
}

int main(){
    int T;
    cin >> T;
    while(T--) {
        Do();
        if(T)
            cout << endl;
    }

    return 0;
}