[CTSC2008]祭祀(构造方案)

前面的话

这道题显然就是最长反链

根据 (Dilworth) 定理：最小链覆盖数 = 最长反链长度

然后传递闭包跑匹配即可

(luogu)交了一下，(WA) 了

(QAQ)

本来各种 (OJ) 上都是只要求最长反链，不需要构造方案

虽然原题要构造

然后 (luogu) 上的同志写了个 (SPJ)， 然后 (luogu) 就要输出方案了

切不掉很难受

luogu传送门

Sol

先放两个博客:

r_64

某神仙

首先建图后发现最大独立集和最小点覆盖互为补集

而这个图中的最大独立集就是最长反链可能，猜的...(因为两两不可到达)

然后只要知道怎么构造最小点覆盖就好了

第一步

对于一个二分图，这样来做：

先最大匹配
每次从左边找到一个未匹配点增广(假的增广，显然增广不了，因为已经是最大匹配)
然后标记点

最后左边没有标记过的点和右边标记过的点就是最小点覆盖

伪证：因为一条假的增广路一定是左边的点作为开头和结尾的，所以选右边的就能覆盖这个假的增广路

去掉这些点就是要找的最长反链

第二步

找到所有可以出现在最长反链上的点

枚举每个点，删掉它以及可以到达它和它可以到达的点(邻居)

再求最长反链，如果大小减小了 (1)，这个点就是可以的(显然选这个点不会和当前的集合冲突)

代码

# include <bits/stdc++.h>
# define IL inline
# define RG register
# define Fill(a, b) memset(a, b, sizeof(a))
using namespace std;
typedef long long ll;

IL int Input(){
    RG char c = getchar(); RG int x = 0, z = 1;
    for(; c < '0' || c > '9'; c = getchar()) z = c == '-' ? -1 : 1;
    for(; c >= '0' && c <= '9'; c = getchar()) x = (x << 1) + (x << 3) + (c ^ 48);
    return x * z;
}

const int maxn(205);
const int maxm(1005);

int n, m, ans, match[maxn], to[maxn], vis[maxn], idx, f[maxn][maxn], g[maxn][maxn];
int ban[maxn], s[maxn], t[maxn];

IL int Dfs(RG int u){
	if(ban[u]) return 0;
    for(RG int i = 1; i <= n; ++i)
        if(f[u][i] && vis[i] != idx && !ban[i]){
            vis[i] = idx;
            if(!match[i] || Dfs(match[i])){
                to[u] = i, match[i] = u;
                return 1;
            }
        }
    return 0;
}

IL void Calc(RG int u){
    if(s[u]) return;
    s[u] = 1;
    for(RG int i = 1; i <= n; ++i)
        if(f[u][i] && !t[i]) t[i] = 1, Calc(match[i]);
}

int main(RG int argc, RG char* argv[]){
    n = Input(), m = Input();
    for(RG int i = 1; i <= m; ++i) g[Input()][Input()] = 1;
    for(RG int i = 1; i <= n; ++i)
        for(RG int j = 1; j <= n; ++j)
            for(RG int k = 1; k <= n; ++k) g[j][k] |= g[j][i] & g[i][k];
	memcpy(f, g, sizeof(f)), ans = n;
    for(RG int i = 1; i <= n; ++i) ++idx, ans -= Dfs(i);
    printf("%d
", ans);
    for(RG int i = 1; i <= n; ++i) if(!to[i]) Calc(i);
    for(RG int i = 1; i <= n; ++i) printf("%d", s[i] && !t[i]);
	puts("");
	for(RG int nw = 1; nw <= n; ++nw){
		for(RG int i = 1; i <= n; ++i) match[i] = to[i] = 0;
		RG int ret = 0, nn = 0;
		Fill(f, 0), Fill(ban, 0);
		for(RG int i = 1; i <= n; ++i)
			if(g[i][nw] || g[nw][i] || i == nw) ban[i] = 1;
			else ++nn;
		ret = nn;
		for(RG int i = 1; i <= n; ++i)
			for(RG int j = 1; j <= n; ++j)
				if(!ban[i] && !ban[j]) f[i][j] = g[i][j];
		for(RG int i = 1; i <= n; ++i) if(!ban[i]) ++idx, ret -= Dfs(i);
		printf("%d", ret == ans - 1);
	}
    return 0;
}