poj 3281 网络流

来源:互联网 发布:java web源码 编辑:程序博客网 时间:2024/06/05 08:34

经典题目,拆点网络流。

#include <iostream>#include <cstdio>#include <cstring>using namespace std;const int maxn=4e2+9,inf=1e9;int n,f,d;int level[maxn],que[maxn];int head[maxn],lon;struct{    int next,to,c;}e[maxn*maxn<<1];void edgeini(){    memset(head,-1,sizeof(head));    lon=-1;}void edgemake(int from,int to,int c){    e[++lon].c=c;    e[lon].to=to;    e[lon].next=head[from];    head[from]=lon;}bool makelevel(int s,int t){    memset(level,0,sizeof(level));    int front=1,end=0;    que[++end]=s;    level[s]=1;    while(front<=end)    {        int u=que[front++];//        cout<<u<<endl;        if(u==t) return true;        for(int k=head[u];k!=-1;k=e[k].next)        {            int v=e[k].to;            if(!level[v]&&e[k].c)            {                que[++end]=v;                level[v]=level[u]+1;            }        }    }    return false;}int dfs(int now,int t,int maxf){    if(t==now) return maxf;    int ret=0;    for(int k=head[now];k!=-1;k=e[k].next)    {        int u=e[k].to;        if(e[k].c&&level[u]==level[now]+1)        {            int f=dfs(u,t,min(maxf-ret,e[k].c));            e[k].c-=f;            e[k^1].c+=f;            ret+=f;            if(ret==maxf) return ret;        }    }    return ret;}int maxflow(int s,int t){    int ret=0;    while(makelevel(s,t))    {        ret+=dfs(s,t,inf);    }    return ret;}int main(){//    freopen("in.txt","r",stdin);    while(scanf("%d%d%d",&n,&f,&d)!=EOF)    {        edgeini();        for(int i=1;i<=n;i++)        {            int tmp1,tmp2,txt;            scanf("%d%d",&tmp1,&tmp2);            for(int j=1;j<=tmp1;j++)            {                scanf("%d",&txt);                edgemake(txt,200+i*2-1,1);                edgemake(200+i*2-1,txt,0);            }            for(int j=1;j<=tmp2;j++)            {                scanf("%d",&txt);                edgemake(200+i*2,txt+100,1);                edgemake(txt+100,200+i*2,0);            }        }        for(int i=1;i<=n;i++)        {            edgemake(200+i*2-1,200+i*2,1);            edgemake(200+i*2,200+i*2-1,0);        }        for(int i=1;i<=f;i++)        {            edgemake(0,i,1);            edgemake(i,0,0);        }        for(int i=1;i<=d;i++)        {            edgemake(i+100,401,1);            edgemake(401,i+100,0);        }        int ans=maxflow(0,401);        cout<<ans<<endl;    }    return 0;}


原创粉丝点击