zoj3165二分图点权最大独立集+dinic匹配时的匹配数

来源:互联网 发布:java小游戏程序代码 编辑:程序博客网 时间:2024/06/06 01:47

题意:要想将一些学生分组,其中有一些学生不能跟另外一些学生在一起,问应该如何分组,在满足要求的前提下,使得人数最多。

分析:就是一个方格取数的模型。二分图点权最大独立集。。。

特点是要求出二分图两边匹配的数目。

#include <iostream>#include <stdio.h>using namespace std;const int oo=1e9;const int mm=111111;const int mn=999;int node,src,dest,edge;int ver[mm],flow[mm],next[mm],end[mm],net[mm];int head[mn],work[mn],dis[mn],q[mn],vis[mn],from[mn],to[mn];void prepare(int _node,int _src,int _dest){    node=_node,src=_src,dest=_dest;    for(int i=0;i<node;i++) head[i]=-1,vis[i]=0;    edge=0;}void addedge(int u,int v,int c){    ver[edge]=v;end[edge]=u;flow[edge]=c;next[edge]=head[u];net[edge]=head[u];head[u]=edge++;    ver[edge]=u;end[edge]=v;flow[edge]=0;next[edge]=head[v];net[edge]=head[v];head[v]=edge++;}bool Dinic_bfs(){    int i,u,v,l,r=0;    for(i=0;i<node;i++)        dis[i]=-1;    dis[q[r++]=src]=0;    for(l=0;l<r;l++)        for(i=head[u=q[l]];i>=0;i=next[i])            if(flow[i]&&dis[v=ver[i]]<0)            {                dis[q[r++]=v]=dis[u]+1;                if(v==dest) return 1;            }    return 0;}int Dinic_dfs(int u,int exp){    if(u==dest) return exp;    for(int &i=work[u],v,tmp;i>=0;i=next[i])    {        if(flow[i]&&dis[v=ver[i]]==dis[u]+1&&(tmp=Dinic_dfs(v,min(exp,flow[i])))>0)        {            flow[i]-=tmp;            flow[i^1]+=tmp;            return tmp;        }    }    return 0;}int Dinic_flow(){    int i,ret=0,delta;    while(Dinic_bfs())    {        for(i=0;i<node;i++)            work[i]=head[i];        while(delta=Dinic_dfs(src,oo))  ret+=delta;    }    return ret;}void dfs(int u){     vis[u]=1;     for(int i=head[u],v;i!=-1;i=next[i])             if(flow[i]&&!vis[v=ver[i]]) dfs(v);}int main(){    int m,n,s,i,j,a,b,sum,u,v,k;    while(~scanf("%d%d%d",&m,&n,&s))    {        prepare(n+m+2,0,n+m+1);        sum=0;        for(i=1;i<=m;i++)        {            scanf("%d",&a);            sum+=a;            addedge(src,i,a);        }        for(i=1;i<=n;i++)        {            scanf("%d",&a);            sum+=a;            addedge(i+n,dest,a);        }        for(i=1;i<=s;i++)        {            scanf("%d%d",&a,&b);            addedge(a,b+n,oo);        }        int tmp=Dinic_flow();        printf("%d ",sum-tmp);        dfs(0);        u=v=0;        for(i=1;i<=m;i++)            if(vis[i])  {from[u++]=i;}        for(i=1;i<=n;i++)            if(!vis[i+m])            {                to[v++]=i;            }        printf("%d %d\n",u,v);        k=0;        for(i=0;i<u;i++)//这种输出方式值得注意。。。        {            if(!k)            {                k=1;printf("%d",from[i]);            }            else    printf(" %d",from[i]);        }        puts("");        k=0;        for(i=0;i<v;i++)        {            if(!k)            {                k=1;printf("%d",to[i]);            }            else printf(" %d",to[i]);        }        puts("");    }    return 0;}