POJ 2125 Destroying The Graph

更详细地思路请参见Amber大牛的论文《最小割模型在信息学竞赛中的应用》

 

对于这道题，如果不考虑权值，就是一个二分图的最小点覆盖问题，加入权值后，变成了二分图的最小点权覆盖问题

 

X和Y分别表示出点和入点的集合，源S向每个出点i连一条边，容量为删除所有从点i出去的边的代价，每个入点i向汇点T连一条边，容量为删除所有进入点i的边的代价，对于每条边i->j，连边二分图中的i->j+n,容量为无穷

 

然后求最大流即可，最后找割边

 

在残留网络中从源点开始DFS，将源点能到达的点和不能到达的点分开标记

 

代码：

#include<iostream> using namespace std; const int MAX=205; const int inf=1<<30; struct node { int u,v,c,next; }g[100000]; int dis[MAX],num[MAX],adj[MAX],cur[MAX],pre[MAX]; int res[MAX],mat[MAX][MAX],vis[MAX]; int s,t,e,vn,n,m,cnt; void add(int u,int v,int c) { g[e].u=u; g[e].v=v; g[e].c=c; g[e].next=adj[u]; adj[u]=e++; g[e].u=v; g[e].v=u; g[e].c=0; g[e].next=adj[v]; adj[v]=e++; } void dfs(int u) { int i,v; vis[u]=1; for(i=adj[u];i!=-1;i=g[i].next) { v=g[i].v; if(!vis[v]&&g[i].c>0) { dfs(v); } } } int sap() { int i,u,v,flow=0,aug=inf,flag; for(i=0;i<=vn;i++) { dis[i]=num[i]=0; cur[i]=adj[i]; } num[0]=vn; u=s; while(dis[s]<vn) { flag=0; for(i=cur[u];i!=-1;i=g[i].next) { v=g[i].v; if(g[i].c&&dis[u]==dis[v]+1) { flag=1; cur[u]=i; aug=min(aug,g[i].c); pre[v]=u; u=v; if(u==t) { flow+=aug; while(u!=s) { u=pre[u]; g[cur[u]].c-=aug; g[cur[u]^1].c+=aug; } aug=inf; } break; } } if(flag) continue; if(--num[dis[u]]==0) return flow; for(dis[u]=vn,i=adj[u];i!=-1;i=g[i].next) { v=g[i].v; if(g[i].c&&dis[v]<dis[u]) { dis[u]=dis[v]; cur[u]=i; } } dis[u]++; num[dis[u]]++; if(u!=s) u=pre[u]; } return flow; } int main() { int i,j,a,b,cost; memset(adj,-1,sizeof(adj)); //memset(mat,0,sizeof(mat)); e=0; scanf("%d%d",&n,&m); s=0; t=n*2+1; vn=t+1; for(i=1;i<=n;i++) { scanf("%d",&cost); add(i+n,t,cost);//进入i } for(j=1;j<=n;j++) { scanf("%d",&cost); add(s,j,cost);//从j出去 } for(i=1;i<=m;i++) { scanf("%d%d",&a,&b); //mat[a][b+n]++; add(a,b+n,inf); } /* for(i=1;i<=n;i++) { for(j=n+1;j<=2*n;j++) { if(mat[i][j]) { add(i,j,inf); } } }*/ printf("%d/n",sap()); memset(vis,0,sizeof(vis)); cnt=0; dfs(s); for(i=0;i<e;i+=2) { if(vis[g[i].u]==1&&vis[g[i].v]==0&&g[i].c==0) { res[cnt++]=i; } } cout<<cnt<<endl; for(i=0;i<cnt;i++) { if(g[res[i]].u==s) { printf("%d -/n",g[res[i]].v); } else { printf("%d +/n",g[res[i]].u-n); } } return 0; }