poj 2125 (最小点权覆盖)

题目链接：http://poj.org/problem?id=2125

题意描述：给你一张有向图，每个点有两个权值（w+，w-），w+表示删掉该点的所有入边所需要的花费，w-表示删掉该点的所有出边所需要的花费，问要删除所有的点需要的最少花费为多少？

分析：因为涉及到了点权和所有的边那么很容易想到最小点权覆盖，但是这里每个点有两个权值怎么处理呢？这样建图： 将点拆成两个点，每个点对应一个权值，有联系的两个点（u,v），u‘到v’‘连一条权值为无穷大的边，v’‘到汇点连一条权值为w+的边，源点到u’连一条权值为w-的边，这样就将原图转换为二分图，这样求的最小个即为原问题的解！最后就是输出问题（找到所有的割边即可）

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;

const int N = 1000;
const int E = 50000;
const int inf = 0xfffffff;

bool vis[N];
int e,head[N];
int dep[N],que[N],cur[N];
struct node
{
	int x,y;
	int nxt;
	int c;
}edge[E];
void addedge(int u,int v,int c)
{
	edge[e].x=u;
	edge[e].y=v;
	edge[e].nxt=head[u];
	edge[e].c=c;
	head[u]=e++;

	edge[e].x=v;
	edge[e].y=u;
	edge[e].nxt=head[v];
	edge[e].c=0;
	head[v]=e++;
}

int maxflow(int s,int t)
{
	int i,j,k,front,rear,top,min,res=0;
	while(1)
	{
		memset(dep,-1,sizeof(dep));
		front=0;
		rear=0;
		que[rear++]=s;
		dep[s]=0;
		while(front!=rear)
		{
			i=que[front++];
			for(j=head[i];j!=-1;j=edge[j].nxt)
				if(edge[j].c&&dep[edge[j].y]==-1)
				{
					dep[edge[j].y]=dep[i]+1;
					que[rear++]=edge[j].y;
				
				}
		}
		if(dep[t]==-1)
			break;
		memcpy(cur,head,sizeof(head)); 
		for(i=s,top=0;;)
		{
			if(i==t)
			{
				min=inf;
				for(k=0;k<top;k++)
					if(min>edge[que[k]].c)
					{
						min=edge[que[k]].c;
						front=k;
					}
				for(k=0;k<top;k++)
				{
					edge[que[k]].c-=min;
					edge[que[k]^1].c+=min;
				}
				res+=min;
				i=edge[que[top=front]].x;
				
			}
			for(j=cur[i];cur[i]!=-1;j=cur[i]=edge[cur[i]].nxt)
				if(dep[edge[j].y]==dep[i]+1&&edge[j].c)
					break;
			if(cur[i]!=-1)
			{
				que[top++]=cur[i];
				i=edge[cur[i]].y;
			}
			else
			{
				if(top==0)
					break;
				dep[i]=-1;
				i=edge[que[--top]].x;
			}
		}
	}
	return res;
}

void dfs(int u)
{
	vis[u] = true;
	for(int i = head[u]; i!=-1;i=edge[i].nxt)
	{
		int y = edge[i].y;
		if(vis[y]==false && edge[i].c)
			dfs(y);
	}
}
int main ()
{
	int n, m;
	int i, x, y;
	while(scanf("%d%d",&n,&m)!=EOF)
	{
		e=0;
		memset(head,-1,sizeof(head));
		for(i=1;i<=n;i++)
		{
			scanf("%d",&x);
			addedge(i+n,2*n+1,x);
		}
		for(i=1;i<=n;i++)
		{
			scanf("%d",&x);
			addedge(0,i,x);
		}
		for(i=0;i<m;i++)
		{
			scanf("%d%d",&x,&y);
			addedge(x,y+n,inf);
		}
		printf("%d\n",maxflow(0,2*n+1));
		memset(vis,0,sizeof(vis));
		dfs(0);
	   int count = 0;

    for( i = 1; i <= 2 * n; i ++ ){

        if( vis[i] ){
            if( i > n )     count ++;        
        }
        else {
            if( i <= n )    count ++;
        }
    }

    printf("%d\n",count);

    for( i = 1; i <= 2 * n; i ++ ){
        if( vis[i] ){            
            if( i > n )     printf("%d +\n",i - n);    
        }
        else {
            if( i <= n )    printf("%d -\n",i);            
        }
    }
	}
	return 0;
}