题意
给一个有n个点m条边的DAG,问最多能选出多少个点使得选出的点两两不能到达。
n≤105,m≤2n
n
≤
10
5
,
m
≤
2
n
分析
根据某定理,最长反链=最小链覆盖,注意这里的链覆盖是允许每个点被重复经过的。
当n比较小的时候,我们就可以通过传递闭包,将原图转为二分图,然后跑最大匹配就好了。
朴素一点的做法就是把每个点拆点,在它们之间连下界为1的边,然后跑上下界最小流,能得到较高的分数但还是比较慢。
考虑在建图的过程中,我们可以直接构造出最大流,也就是把每条s->i->i+n->t都增广一遍。
然后再考虑最小流的退流过程,那么转化后的建图就变为:s向i连流量为1的边,i+n向t连流量为1的边,i+n向i连流量为inf的边,然后对于图上一条x->y的边,从x向y+n连流量为inf的边。
最后n-最大流即为答案。
代码
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<queue>
const int N=200005;
const int inf=1000000000;
int n,m,s,t,cnt,last[N],cur[N],dis[N];
struct edge{int to,c,next;}e[N*10];
std::queue<int> que;
void addedge(int u,int v,int c)
{
e[++cnt].to=v;e[cnt].c=c;e[cnt].next=last[u];last[u]=cnt;
e[++cnt].to=u;e[cnt].c=0;e[cnt].next=last[v];last[v]=cnt;
}
bool bfs()
{
for (int i=s;i<=t;i++) dis[i]=0;
while (!que.empty()) que.pop();
dis[s]=1;que.push(s);
while (!que.empty())
{
int u=que.front();que.pop();
for (int i=last[u];i;i=e[i].next)
if (e[i].c&&!dis[e[i].to])
{
dis[e[i].to]=dis[u]+1;
if (e[i].to==t) return 1;
que.push(e[i].to);
}
}
return 0;
}
int dfs(int x,int maxf)
{
if (x==t||!maxf) return maxf;
int ret=0;
for (int &i=cur[x];i;i=e[i].next)
if (e[i].c&&dis[e[i].to]==dis[x]+1)
{
int f=dfs(e[i].to,std::min(e[i].c,maxf-ret));
e[i].c-=f;
e[i^1].c+=f;
ret+=f;
if (maxf==ret) break;
}
return ret;
}
int dinic()
{
int ans=0;
while (bfs())
{
for (int i=s;i<=t;i++) cur[i]=last[i];
ans+=dfs(s,inf);
}
return ans;
}
int main()
{
scanf("%d%d",&n,&m);
s=0;t=n*2+1;cnt=1;
for (int i=1;i<=n;i++)
{
addedge(s,i,1);
addedge(i+n,t,1);
addedge(i+n,i,inf);
}
for (int i=1;i<=m;i++)
{
int x,y;scanf("%d%d",&x,&y);
addedge(x,y+n,inf);
}
printf("%d",n-dinic());
return 0;
}