Codeforces Round #432 (Div. 1, based on IndiaHacks Final Round 2017) D. Tournament Construction（dp + 构造）

题意

一个竞赛图的度数集合是由该竞赛图中每个点的出度所构成的集合。

现给定一个 $m$ 个元素的集合，第 $i$ 个元素是 $a_i$ 。（此处集合已经去重）

判断其是否是一个竞赛图的度数集合，如果是，找到点数最小的满足条件的竞赛图，并构造方案。

$m, a_i \le 30$ ， $a_i$ 互不相同。

题解

首先给出结论:假如给出每个点的出度,那么这些点能形成一个竞赛图当且仅当排序后的序列 $d_1, d_2, d_3, \dots , d_n$ 满足对于所有 $k < n$ 有 $\displaystyle \sum _{k_i=1}^{k} d_i ≥ {k \choose 2}$ 且 $\displaystyle \sum_{i = 1} ^ {n} d_i = {n \choose 2}$ 。

考虑证明，不难发现对于竞赛图的任意点集 $S$ ，都必须满足 $\displaystyle \sum_{i \in S} d_i \ge {|S| \choose 2}$ ，因为点集 $S$ 之间的一对点至少会存在一条边贡献一个出度。

并且由于

(\binom{n}{2}) = \sum_{i = 1}^{n} i - 1 \leq n \times max {a_{i}}

${n \choose 2} = \sum_{i = 1} ^ {n} i - 1 \le n \times \max \{a_i\}$

所以发现最后的点数是 $O(\max \{a_i\})$ 级别的，其实最多就是 $2 \max \{a_i\} + 1$ 。

那么就可以 $dp$ 了，先将给定集合中的元素从小到大排序。

设 $f[i][j][sum]$ 表示考虑完前 $i$ 个元素，当前图中已经有 $j$ 个点了，度数之和为 $sum$ 是否可行。

对于这个 $dp$ 转移时枚举下一个元素出现了几次就行了。

然后我们只需要找到一个最小的 $i$ 使得 $\displaystyle f[n][i][{i \choose 2}]$ 满足就行了，这个就是合法点数。

似乎一定可以构造出方案。。不存在无解qwq 至于原因，可以去问 zhou888 。

然后需要构造方案，我们转移 $dp$ 的时候多记一下当前如何转移的就行了，最后逆推回去就得到了度数可重集合。

然后我们用这个度数可重集合来构造答案，每次将度数从小到大排序，然后把第一个点向后面出度个点连边就行了，不难发现这样一定可以构造出一组合法解。

最后复杂度就是 $O(n{a_i}^3)$ 的，可以通过qwq

总结

对于竞赛图可以考虑任意一个子集度数都不少于 $\displaystyle {n \choose 2}$ 的性质。

然后考虑构造图的时候可以贪心的去暴力分配度数，或者用网络流模拟这个过程也行qwq

代码

用了 goto 以及 lambda 函数等骚操作。。

#include <bits/stdc++.h>

#define For(i, l, r) for(register int i = (l), i##end = (int)(r); i <= i##end; ++i)
#define Fordown(i, r, l) for(register int i = (r), i##end = (int)(l); i >= i##end; --i)
#define Set(a, v) memset(a, v, sizeof(a))
#define Cpy(a, b) memcpy(a, b, sizeof(a))
#define debug(x) cout << #x << ": " << (x) << endl
#define DEBUG(...) fprintf(stderr, __VA_ARGS__)

using namespace std;

template<typename T> inline bool chkmin(T &a, T b) {return b < a ? a = b, 1 : 0;}
template<typename T> inline bool chkmax(T &a, T b) {return b > a ? a = b, 1 : 0;}

inline int read() {
    int x(0), sgn(1); char ch(getchar());
    for (; !isdigit(ch); ch = getchar()) if (ch == '-') sgn = -1;
    for (; isdigit(ch); ch = getchar()) x = (x * 10) + (ch ^ 48);
    return x * sgn;
}

void File() {
#ifdef zjp_shadow
	freopen ("D.in", "r", stdin);
	freopen ("D.out", "w", stdout);
#endif
}

const int Maxn = 35, N = 65;

bitset<2510> dp[Maxn][N], pre[Maxn][N];

int n, m, a[N], d[N], sum[N]; bitset<N> vis, G[N];

inline int Comb(int x) {
	return x * (x - 1) / 2;
}

int main () {

	File();

	n = read();
	For (i, 1, n) a[i] = read();
	sort(a + 1, a + n + 1);
	For (i, 1, n) sum[i] = sum[i - 1] + a[i];

	dp[0][0][0] = true;
	int pos = n;
	For(j, 1, 61) For (i, 1, n) {
		For (k, max(sum[i], Comb(j)), j * 30) {
			if (dp[i - 1][j - 1][k - a[i]]) dp[i][j][k] = pre[i][j][k] = true;
			if (dp[i][j - 1][k - a[i]]) dp[i][j][k] = !(pre[i][j][k] = false);
		}
		if (dp[n][j][Comb(j)]) { n = j; goto tag; }
	}
tag: ;

	 int deg = Comb(n);
	 Fordown (i, n, 1)
		 d[i] = a[pos], pos -= pre[pos][i][deg], deg -= d[i];

	 int id[N], len;
	 For (i, 1, n) {
		 len = 0;
		 For (j, 1, n) if (!vis[j]) id[++ len] = j;
		 sort(id + 1, id + len + 1, [&](int lhs, int rhs) { return d[lhs] < d[rhs]; });
		 vis[id[1]] = true;
		 For (j, 2, d[id[1]] + 1) G[id[1]][id[j]] = true;
		 For (j, d[id[1]] + 2, len) G[id[j]][id[1]] = true, -- d[id[j]];
	 }

	 printf ("%d\n", n);

	 For (i, 1, n) 
		 For (j, 1, n + 1)
		 putchar (j == jend ? '\n' : G[i][j] + '0');

	 return 0;

}

__EOF__

本文作者：zjp_shadow
本文链接：https://www.cnblogs.com/zjp-shadow/p/9769240.html
关于博主：评论和私信会在第一时间回复。或者直接私信我。
版权声明：本博客所有文章除特别声明外，均采用 BY-NC-SA 许可协议。转载请注明出处！
声援博主：如果您觉得文章对您有帮助，可以点击文章右下角【推荐】一下。您的鼓励是博主的最大动力！