第一次写最小流真的有点方..
先膜一发Oxer神犇的二分答案做法传送门
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<vector>
#include<cmath>
#include<queue>
#include<set>
#include<map>
#define N 100006
#define mx 111111111
using namespace std;
int sc()
{
int i=0; char c=getchar();
while(c>'9'||c<'0')c=getchar();
while(c>='0'&&c<='9')i=i*10+c-'0',c=getchar();
return i;
}
int dis[N],q[N];
int head[N],nxt[N],lst[N],c[N];
int a[211][211],pos[211],d[211];
int n,m,tot,s,t,S,T;
void insert(int x,int y,int z)
{
lst[++tot]=y;nxt[tot]=head[x];head[x]=tot;c[tot]=z;
lst[++tot]=x;nxt[tot]=head[y];head[y]=tot;c[tot]=0;
}
bool BFS()
{
memset(dis,0,sizeof(dis));
dis[q[1]=S]=1;
for(int l=1,r=2;l<r;l++)
{
int x=q[l];
for(int i=head[x];i;i=nxt[i])
if(c[i]&&!dis[lst[i]])
dis[q[r++]=lst[i]]=dis[x]+1;
}
return dis[T];
}
int dfs(int x,int f)
{
if(x==T)return f;
int ww=0,w;
for(int i=head[x];i;i=nxt[i])
if(c[i]&&dis[lst[i]]==dis[x]+1)
{
w=dfs(lst[i],min(c[i],f-ww));
c[i]-=w,c[i^1]+=w,ww+=w;
if(ww==f)return ww;
}
if(!ww)dis[x]=0;
return ww;
}
bool jud(int x)
{
memset(head,0,sizeof(head));
memset(d,0,sizeof(d));
int ans=0,sum=0;
tot=1;s=n+1,t=n+2,S=n+3,T=n+4;
insert(t,s,x);
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
if(a[i][j])
d[i]--,d[j]++,insert(i,j,mx);
for(int i=1;i<=n;i++)
insert(s,i,mx),insert(i,t,mx);
for(int i=1;i<=n;i++)
if(d[i]>0) insert(S,i,d[i]),sum+=d[i];
else if(d[i]<0) insert(i,T,-d[i]);
while(BFS())ans+=dfs(S,mx);
return sum==ans;
}
int main()
{
n=sc();
for(int i=1;i<=n;i++)
{
int m=sc();
while(m--)a[i][sc()]=1;
}
int l=1,r=n*n,ans=n*n;
while(l<=r)
{
int mid=l+r>>1;
if(jud(mid))ans=mid,r=mid-1;
else l=mid+1;
}
cout<<ans;
return 0;
}
其实也可以不用二分答案,直接搞一个类似于退流的东西
PoPoQQQ大爷的第二种做法传送门
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<vector>
#include<queue>
#include<cmath>
#include<set>
#include<map>
#define N 20006
#define mx 1111111111
using namespace std;
int sc()
{
int i=0;char c=getchar();
while(c>'9'||c<'0')c=getchar();
while(c>='0'&&c<='9')i=i*10+c-'0',c=getchar();
return i;
}
int dis[N],du[N],q[N];
int head[N],nxt[N],lst[N],c[N];
int n,tot=1,s,t,S,T,ans;
void insert(int x,int y,int z)
{
lst[++tot]=y;nxt[tot]=head[x];head[x]=tot;c[tot]=z;
lst[++tot]=x;nxt[tot]=head[y];head[y]=tot;c[tot]=0;
}
bool BFS()
{
for(int i=1;i<=T;i++)dis[i]=0;
dis[q[1]=S]=1;
for(int l=1,r=2;l<r;l++)
{
int x=q[l];
for(int i=head[x];i;i=nxt[i])
if(c[i]&&!dis[lst[i]])
dis[q[r++]=lst[i]]=dis[x]+1;
}
return dis[T];
}
int dfs(int x,int f)
{
if(x==T)return f;
int ww=0,w;
for(int i=head[x];i;i=nxt[i])
if(c[i]&&dis[lst[i]]==dis[x]+1)
{
w=dfs(lst[i],min(f-ww,c[i]));
c[i]-=w,c[i^1]+=w,ww+=w;
if(f==ww)return ww;
}
if(!ww)dis[x]=0;
return ww;
}
int main()
{
n=sc();s=n+1,t=n+2,S=n+3,T=n+4;
for(int i=1;i<=n;i++)
{
int m=sc(),x;
du[i]-=m;
while(m--)
du[x=sc()]++,
insert(i,x,mx);
insert(s,i,mx),insert(i,t,mx);
}
for(int i=1;i<=n;i++)
if(du[i]>0) insert(S,i,du[i]);
else if(du[i]<0) insert(i,T,-du[i]);
insert(t,s,mx);
while(BFS())dfs(S,mx);
ans=c[tot];c[tot]=c[tot-1]=0;
for(int i=1;i<=n;i++)
for(int j=lst[i];j;j=nxt[j])
if(lst[j]==S||lst[j]==T)
c[j]=0,c[j^1]=0;
insert(S,t,mx),insert(s,T,mx);
while(BFS())ans-=dfs(S,mx);
cout<<ans;
return 0;
}
做法三:其实完全不用像算法二那样搞一个退流的的操作,可以先不连
t−>s
的边跑一边网络流,然后再连起来跑一边网络流,然后第二编的最大流即为答案。
在算法二的基础上稍作变动即可
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<vector>
#include<queue>
#include<cmath>
#include<set>
#include<map>
#define N 20006
#define mx 1111111111
using namespace std;
int sc()
{
int i=0;char c=getchar();
while(c>'9'||c<'0')c=getchar();
while(c>='0'&&c<='9')i=i*10+c-'0',c=getchar();
return i;
}
int dis[N],du[N],q[N];
int head[N],nxt[N],lst[N],c[N];
int n,tot=1,s,t,S,T,ans;
void insert(int x,int y,int z)
{
lst[++tot]=y;nxt[tot]=head[x];head[x]=tot;c[tot]=z;
lst[++tot]=x;nxt[tot]=head[y];head[y]=tot;c[tot]=0;
}
bool BFS()
{
for(int i=1;i<=T;i++)dis[i]=0;
dis[q[1]=S]=1;
for(int l=1,r=2;l<r;l++)
{
int x=q[l];
for(int i=head[x];i;i=nxt[i])
if(c[i]&&!dis[lst[i]])
dis[q[r++]=lst[i]]=dis[x]+1;
}
return dis[T];
}
int dfs(int x,int f)
{
if(x==T)return f;
int ww=0,w;
for(int i=head[x];i;i=nxt[i])
if(c[i]&&dis[lst[i]]==dis[x]+1)
{
w=dfs(lst[i],min(f-ww,c[i]));
c[i]-=w,c[i^1]+=w,ww+=w;
if(f==ww)return ww;
}
if(!ww)dis[x]=0;
return ww;
}
int main()
{
n=sc();s=n+1,t=n+2,S=n+3,T=n+4;
for(int i=1;i<=n;i++)
{
int m=sc(),x;
du[i]-=m;
while(m--)
du[x=sc()]++,
insert(i,x,mx);
insert(s,i,mx),insert(i,t,mx);
}
for(int i=1;i<=n;i++)
if(du[i]>0) insert(S,i,du[i]);
else if(du[i]<0) insert(i,T,-du[i]);
//insert(t,s,mx);
while(BFS())dfs(S,mx);
insert(t,s,mx);
//ans=c[tot];c[tot]=c[tot-1]=0;
//for(int i=1;i<=n;i++)
//for(int j=lst[i];j;j=nxt[j])
//if(lst[j]==S||lst[j]==T)
//c[j]=0,c[j^1]=0;
//insert(S,t,mx),insert(s,T,mx);
while(BFS())ans+=dfs(S,mx);
cout<<ans;
return 0;
}