poj3281 拆点+网络流

来源：互联网发布：网络红猫编辑：程序博客网时间：2024/05/16 15:38
和上题一样，水过。
代码：
#include<iostream>#include<cstring>#include<cstdio>#define MAXN 815#define INF 1e8#define min(a,b) (a<b?a:b)#define max(a,b) (a>b?a:b)using namespace std;struct edge{    int u,v,w,next;}E[200000];int head[MAXN],ecnt;int gap[MAXN],cur[MAXN],pre[MAXN],dis[MAXN];int l,r,mid;int N,M,scr,sink,vn,num;int k,c1,m;void Insert(int u,int v,int w){    E[ecnt].u=u;    E[ecnt].v=v;    E[ecnt].w=w;    E[ecnt].next=head[u];    head[u]=ecnt++;    E[ecnt].u=v;    E[ecnt].v=u;    E[ecnt].w=0;    E[ecnt].next=head[v];    head[v]=ecnt++;}int Sap(int s,int t,int n)//核心代码(模版){    int ans=0,aug=INF;//aug表示增广路的流量    int i,v,u=pre[s]=s;    for(i=0;i<=n;i++)    {        cur[i]=head[i];        dis[i]=gap[i]=0;    }    gap[s]=n;    bool flag;    while(dis[s]<n)    {        flag=false;        for(int &j=cur[u];j!=-1;j=E[j].next)//一定要定义成int &j,why        {            v=E[j].v;            if(E[j].w>0&&dis[u]==dis[v]+1)            {                flag=true;//找到容许边                aug=min(aug,E[j].w);                pre[v]=u;                u=v;                if(u==t)                {                    ans+=aug;                    while(u!=s)                    {                        u=pre[u];                        E[cur[u]].w-=aug;                        E[cur[u]^1].w+=aug;//注意                    }                    aug=INF;                }                break;//找到一条就退出            }        }        if(flag) continue;        int mindis=n;        for(i=head[u];i!=-1;i=E[i].next)        {            v=E[i].v;            if(E[i].w>0&&dis[v]<mindis)            {                mindis=dis[v];                cur[u]=i;            }        }        if((--gap[dis[u]])==0) break;        gap[dis[u]=mindis+1]++;        u=pre[u];    }    return ans;}int n,f,d;int main(){   while(scanf("%d%d%d",&n,&f,&d)!=EOF)   {       memset(head,-1,sizeof(head));ecnt=0;       scr=0;sink=f+d+2*n+1;vn=sink+1;       for(int i=1;i<=f;i++)        Insert(scr,i,1);       for(int i=1;i<=d;i++)        Insert(i+f,sink,1);       for(int i=1;i<=n;i++)       {           Insert(i+d+f,i+d+f+n,1);           int x,y;int v;           scanf("%d%d",&x,&y);           for(int j=1;j<=x;j++)           {               scanf("%d",&v);               Insert(v,i+f+d,1);           }           for(int j=1;j<=y;j++)           {               scanf("%d",&v);               Insert(i+f+d+n,v+f,1);           }       }    printf("%d\n",Sap(scr,sink,vn));   }   return 0;}