[CTSC2008]祭祀(构造方案)
前面的话
这道题显然就是最长反链
根据 \(Dilworth\) 定理:最小链覆盖数 = 最长反链长度
然后传递闭包跑匹配即可
\(luogu\)交了一下,\(WA\) 了
\(QAQ\)
本来各种 \(OJ\) 上都是只要求最长反链,不需要构造方案
虽然原题要构造
然后 \(luogu\) 上的同志写了个 \(SPJ\), 然后 \(luogu\) 就要输出方案了
切不掉很难受
Sol
先放两个博客:
首先建图后发现最大独立集和最小点覆盖互为补集
而这个图中的最大独立集就是最长反链可能,猜的...(因为两两不可到达)
然后只要知道怎么构造最小点覆盖就好了
第一步
对于一个二分图,这样来做:
先最大匹配
每次从左边找到一个未匹配点增广(假的增广,显然增广不了,因为已经是最大匹配)
然后标记点
最后左边没有标记过的点和右边标记过的点就是最小点覆盖
伪证:因为一条假的增广路一定是左边的点作为开头和结尾的,所以选右边的就能覆盖这个假的增广路
去掉这些点就是要找的最长反链
第二步
找到所有可以出现在最长反链上的点
枚举每个点,删掉它以及可以到达它和它可以到达的点(邻居)
再求最长反链,如果大小减小了 \(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\n", 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;
}
[CTSC2008]祭祀(构造方案)的更多相关文章
- BZOJ 1143 1143: [CTSC2008]祭祀river 最长反链
1143: [CTSC2008]祭祀river Description 在遥远的东方,有一个神秘的民族,自称Y族.他们世代居住在水面上,奉龙王为神.每逢重大庆典, Y族都会在水面上举办盛大的祭祀活动. ...
- Bzoj 2718: [Violet 4]毕业旅行 && Bzoj 1143: [CTSC2008]祭祀river 传递闭包,二分图匹配,匈牙利,bitset
1143: [CTSC2008]祭祀river Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 1878 Solved: 937[Submit][St ...
- bzoj 1143: [CTSC2008]祭祀river / 2718: [Violet 4]毕业旅行 -- 二分图匹配
1143: [CTSC2008]祭祀river Time Limit: 10 Sec Memory Limit: 162 MB Description 在遥远的东方,有一个神秘的民族,自称Y族.他们 ...
- BZOJ1143 [CTSC2008]祭祀river 【二分图匹配】
1143: [CTSC2008]祭祀river Time Limit: 10 Sec Memory Limit: 162 MB Submit: 3236 Solved: 1651 [Submit] ...
- bzoj1143(2718)[CTSC2008]祭祀river(最长反链)
1143: [CTSC2008]祭祀river Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 2781 Solved: 1420[Submit][S ...
- BZOJ 1143: [CTSC2008]祭祀river 最长反链
1143: [CTSC2008]祭祀river Time Limit: 1 Sec Memory Limit: 256 MB 题目连接 http://www.lydsy.com/JudgeOnline ...
- 1143: [CTSC2008]祭祀river(最长反链)
1143: [CTSC2008]祭祀river 题目链接:https://www.lydsy.com/JudgeOnline/problem.php?id=1143 Description: 在遥远的 ...
- P4298 [CTSC2008]祭祀
P4298 [CTSC2008]祭祀 传递闭包跑一遍按联通建图 $(1)$最长反链长度=最小链覆盖=n-最大匹配 $(2)$定义作为最大匹配出现在左端点的集合为$S$,作为最大匹配出现在右端点的集合为 ...
- 洛谷 P4298: bzoj 1143: [CTSC2008]祭祀
题目传送门:洛谷 P4298. 题意简述: 给定一个 \(n\) 个点,\(m\) 条边的简单有向无环图(DAG),求出它的最长反链,并构造方案. 最长反链:一张有向无环图的最长反链为一个集合 \(S ...
随机推荐
- 解决pymysql.err.InternalError: (1366, "Incorrect string value: '\\xF0\\x9F\\x8C\\xB8' for column 'headline' at row 1")
解决pymysql.err.InternalError: (1366, "Incorrect string value: '\\xF0\\x9F\\x8C\\xB8' for column ...
- JSP页面开发知识点整理
刚学JSP页面开发,把知识点整理一下. ----------------------------------------------------------------------- JSP语法htt ...
- Vue-cli 2.9 多页配置及多页面之间的跳转问题
vue开发,现在大部分做的都是(SPA)应用,但是,由于,需求不同,我们针对的用户需求变更较为,频繁,如果每次都全量打包更新,给开发的自测,及测试妹子的任务就会多,每次都要重新验证一下才放心.所以,想 ...
- 功能一: 数据库访问DAO层方法定义
功能1: 今天到现在为止 实战课程的访问量 yyyyMMdd courseID 使用数据库来进行存储我们的统计结果 Spark Streaming把统计结果写入到数据库里面 可视化前端根据: yyyy ...
- CF1007D. Ants(树链剖分+线段树+2-SAT及前缀优化建图)
题目链接 https://codeforces.com/problemset/problem/1007/D 题解 其实这道题本身还是挺简单的,这里只是记录一下 2-SAT 的前缀优化建图的相关内容. ...
- springboot(十二)-分布式锁(redis)
什么是分布式锁? 要介绍分布式锁,首先要提到与分布式锁相对应的是线程锁.进程锁. 线程锁:主要用来给方法.代码块加锁.当某个方法或代码使用锁,在同一时刻仅有一个线程执行该方法或该代码段.线程锁只在同一 ...
- DB2 close auto commit
db2 关闭命令行CLP自动提交 --临时关闭自动提交 #db2 "update command options using C off --永久关闭自动提交 ----linux 环境下 # ...
- pull强制覆盖本地
今天我总结的是在项目中经常用到的Git命令,上传和下拉文件. 当然在进行上传和下拉操作之前,你首先要做的就是将本地和Git库连接起来. 连接命令: git remote add origin + 你G ...
- 《LeetBook》leetcode题解(6): ZigZag Conversion[E]
我现在在做一个叫<leetbook>的免费开源书项目,力求提供最易懂的中文思路,目前把解题思路都同步更新到gitbook上了,需要的同学可以去看看 书的地址:https://hk029.g ...
- Callable和Future详解
Java程序员必须掌握的线程知识-Callable和Future Callable和Future出现的原因 创建线程的两种方式:继承Thread类和实现Runnable接口 这两种方式都有一种缺陷,执 ...