题目链接

戳我

\(Solution\)

第一问

这道题要知道一个叫做\(Dilworth\)的定理

最长反链\(=\)最小链覆盖

证明(\(from\ r\_64\)):

所以我们只要求一个最小链覆盖即可

这个很好求

对于每个点拆点,拆成\((x,x')\),\(s->x\)流量为\(1\),\(x'->t\),流量为\(1\)

对于每个相连通的边\((x,y)\),将\(x->y'\)流量为\(1\)

最后用\(n-Dinic()\)即可

但是注意不能根据输入的边连边,需要用传递闭包判一下连通性,因为如下图这种情况,就不对

这样子如果直接连边跑出来的结果是\(3\)但是实际上应该是\(2\)

本来应该是没有后面两问的,但是因为有了\(spj\)导致后面两问又出来了

第二问

它要求一个最长反链,也就是求一个最大独立集

\(from\ zsy\)

第三问

我们删掉一个点和与他相关的点与边,再求一次最长反链,如果答案减少了\(1\)则这个可能在最长反链中否则不可能,具体见代码

\(Code\)

#include<bits/stdc++.h>
using namespace std;
const int inf=1e9;
typedef long long ll;
int read(){
int x=0,f=1;char c=getchar();
while(c<'0'||c>'9') f=(c=='-')?-1:1,c=getchar();
while(c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
return x*f;
}
struct node{
int to,next,v;
}a[2000001];
int head[1000001],cnt=1,tot,ans,n,m,s,t,x,y,z,dep[1000001],f[1001][1001],vis[1001],bj[1001];
void add(int x,int y,int c){
a[++cnt].to=y,a[cnt].next=head[x],a[cnt].v=c,head[x]=cnt;
a[++cnt].to=x,a[cnt].next=head[y],a[cnt].v=0,head[y]=cnt;
}
queue<int> q;
int bfs(){
memset(dep,0,sizeof(dep));
q.push(s);
dep[s]=1;
while(!q.empty()){
int now=q.front();
q.pop();
for(int i=head[now];i;i=a[i].next){
int v=a[i].to;
if(!dep[v]&&a[i].v>0)
dep[v]=dep[now]+1,q.push(v);
}
}
if(dep[t])
return 1;
return 0;
}
int dfs(int k,int list){
if(k==t||!list)
return list;
for(int i=head[k];i;i=a[i].next){
int v=a[i].to;
if(dep[v]==dep[k]+1&&a[i].v>0){
int p=dfs(v,min(list,a[i].v));
if(p){
a[i].v-=p;
a[i^1].v+=p;
return p;
}
}
}
return dep[k]=0;
}
int Dinic(){
int ans=0,k;
while(bfs())
while((k=dfs(s,inf)))
ans+=k;
return ans;
}
void floyd(){
for(int k=1;k<=n;k++)
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
f[i][j]=max(f[i][j],f[i][k]&f[k][j]);
}
void dfs1(int x){
if(!vis[x]&&x<=n) return ;
if(vis[x]&&x>n) return ;
if(x<=n){
vis[x]=0;
for(int i=head[x];i;i=a[i].next)
if(a[i].v)
dfs1(a[i].to);
}
else {
vis[x]=1;
for(int i=head[x];i;i=a[i].next)
if(!a[i^1].v)
dfs1(a[i].to);
}
}
int main(){
n=read(),m=read(),s=0,t=2*n+1;
for(int i=1;i<=m;i++)
x=read(),y=read(),f[x][y]=1;
for(int i=1;i<=n;i++)
add(s,i,1),bj[i]=cnt-1,add(i+n,t,1),f[i][i]=1,vis[i]=1;
floyd();
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
if(i!=j&&f[i][j])
add(i,j+n,1);
printf("%d\n",ans=(n-Dinic()));
for(int i=1;i<=n;i++)
if(a[bj[i]].v)
dfs1(i);
for(int i=1;i<=n;i++)
printf("%d",(vis[i]|vis[i+n])^1);
cout<<endl;
for(int k=1;k<=n;k++){
memset(head,0,sizeof(head)),cnt=1,tot=0;
for(int i=1;i<=n;i++)
if(i!=k&&!f[i][k]&&!f[k][i])
add(s,i,1),add(i+n,t,1),tot++;
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
if(i!=j&&f[i][j]&&i!=k)
add(i,j+n,1);
printf("%d",ans-1==tot-Dinic());
}
}

「CTSC 2008」祭祀的更多相关文章

  1. 「CTSC 2011」排列

    「CTSC 2011」排列 要求不存在公差为 A 或者公比为 B 的子列,那么实际上可以把该问题转化为求一个图的最优拓朴序. 任意差为 A 或者比为 B 的两个数连一条边. 求一个合法序列的答案可以用 ...

  2. 「CTSC 2011」幸福路径

    [「CTSC 2011」幸福路径 蚂蚁是可以无限走下去的,但是题目对于精度是有限定的,只要满足精度就行了. \({(1-1e-6)}^{2^{25}}=2.6e-15\) 考虑使用倍增的思想. 定义\ ...

  3. 「BZOJ 1831」「AHOI 2008」逆序对「贪心」

    题意 给定一个长度为\(n\),值域为\([1,k]\),某些位置不确定的数组,求最小的逆序对.\(n\leq 10^4, k \leq 100\) 题解 这题有人用前缀和优化\(dp\)过了,但是这 ...

  4. 「BZOJ 1791」「IOI 2008」Island「基环树」

    题意 求基环树森林所有基环树的直径之和 题解 考虑的一个基环树的直径,只会有两种情况,第一种是某个环上结点子树的直径,第二种是从两个环上结点子树内的最深路径,加上环上这两个结点之间的较长路径. 那就找 ...

  5. 「CTSC 2013」组合子逻辑

    Tag 堆,贪心 Description 给出一个数列 \(n\) 个数,一开始有一个括号包含 \([1,n]\),你需要加一些括号,使得每个括号(包括一开始的)所包含的元素个数 \(\leq\) 这 ...

  6. Solution -「CTSC 2018」「洛谷 P4602」混合果汁

    \(\mathcal{Description}\)   Link.   \(n\) 种果汁,第 \(i\) 种美味度为 \(d_i\),每升价格 \(p_i\),一共 \(l_i\) 升.\(m\) ...

  7. Solution -「NOI 2008」「洛谷 P3980」志愿者招募

    \(\mathcal{Description}\)   Link.   一项持续 \(n\) 天的任务,第 \(i\) 天需要至少 \(a_i\) 人工作.还有 \(m\) 种雇佣方式,第 \(i\) ...

  8. 「POJ Challenge」生日礼物

    Tag 堆,贪心,链表 Solution 把连续的符号相同的数缩成一个数,去掉两端的非正数,得到一个正负交替的序列,把该序列中所有数的绝对值扔进堆中,用所有正数的和减去一个最小值,这个最小值的求法与「 ...

  9. 前端构建工具之gulp(一)「图片压缩」

    前端构建工具之gulp(一)「图片压缩」 已经很久没有写过博客了,现下终于事情少了,开始写博吧 今天网站要做一些优化:图片压缩,资源合并等 以前一直使用百度的FIS工具,但是FIS还没有提供图片压缩的 ...

随机推荐

  1. 查询SQL Server版本号

    一.查看SQL Server 2005版本号 SELECT @@VERSION 显示以下信息: Microsoft SQL Server 2005 - 9.00.3042.00 (Intel X86) ...

  2. 获取url传来的参数

    //根据传递过来的参数name获取对应的值 function getParameter(name) { var reg = new RegExp("(^|&)" + nam ...

  3. go之环境安装

    1)Linux安装 https://studygolang.com/dl # 在 ~ 下创建 go 文件夹,并进入 go 文件夹 mkdir ~/go && cd ~/go 下载的 g ...

  4. JS和JS是IE上JavaScript或JScript的缩写。

    JS和JS是IE上JavaScript或JScript的缩写.javascript是所有浏览器的开放式标准脚本语言JScript是微软自己的开放式脚本语言标准,只有微软的IE浏览器遵循.JScript ...

  5. Scrapy框架之高级 转

    一.CrawlSpider模板 创建项目 scrapy startproject 项目名称 查看模板 scrapy genspider -l 创建crawl模板 scrapy genspider -t ...

  6. mint-ui下拉加载(demo实例)

    <template> <div class="share"> <div class="header"> <div cl ...

  7. layDate面板出现红色花纹图案

    要使用layDate,有两种方法: 1. 要么在引用layui.js和layui.css,然后通过layui.use('laydate', callback) 加载模块后,调用方法使用. 2. 去la ...

  8. three.js之创建坐标系网格

    <!DOCTYPE html> <html> <head> <meta charset=utf-8> <title>My first thr ...

  9. GitHub的本地与远程

    首先要有一个github账户(这不是废话吗) 在linux中先安装git arch linux : pacman -S git 在终端里输入 ssh-keygen ##一直默认就可以了 将公钥加入到G ...

  10. 重新编译mysqldump,使mysqldump具有进度输出功能

    重新编译mysql,使mysqldump具有进度输出功能 安装编译过程所必须的依赖包和环境 yum install -y gcc cmake boost boost-build boost-devel ...