2502: 清理雪道|有上下界的网络流之最小流

第一次写最小流真的有点方..
先膜一发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;
}