题意:就是删掉所有边,对于每个点,可以删掉他的出边和入边(都有花费的),要求花费最小。
求最小点权覆盖集,不懂可以看07年胡波涛的论文。。。
把每个点拆开,源点连到对应入边的点,容量就是这个点删掉入边的花费,出边连到汇点,流量同样,然后从入边的那个点向出边的点连边就是了。。
跑一遍最大流之后,沿残余网络遍历,出边遍历到了就要选,入边没遍历到就选。
/*************************************************************************
> File Name: poj2125.cpp
> Author: tjw
> Mail:
> Created Time: 2014年11月02日 星期日 10时28分50秒
************************************************************************/
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#include<vector>
#include<stack>
#include<map>
#define ll long long
#define ls k<<1
#define rs k<<1|1
using namespace std;
const int MAXN=210;
const int MAXE=12000;
const int INF=1<<30;
struct EDGE
{
int v,next;
int cap;
}edge[MAXE];
int head[MAXN],size;
void init()
{
memset(head,-1,sizeof(head));
size=0;
}
void add_edge(int u,int v,int cap)
{
edge[size].v=v;
edge[size].cap=cap;
edge[size].next=head[u];
head[u]=size++;
edge[size].v=u;
edge[size].cap=0;
edge[size].next=head[v];
head[v]=size++;
}
int gap[MAXN],dist[MAXN],pre[MAXN],pe[MAXN];
int ISAP(int s,int t)
{
int cur_flow,u,neck,i,temp;
int max_flow;
memset(gap,0,sizeof(gap));
memset(dist,0,sizeof(dist));
memset(pre,-1,sizeof(pre));
for(i=0;i<=t+1;i++)
pe[i]=head[i];
max_flow=0;
u=s;
while(dist[s]<t+1)
{
if(u==t)
{
cur_flow=INF+1;
for(i=s;i!=t;i=edge[pe[i]].v)
{
if(cur_flow>edge[pe[i]].cap)
{
cur_flow=edge[pe[i]].cap;
neck=i;
}
}
for(i=s;i!=t;i=edge[pe[i]].v)
{
temp=pe[i];
edge[temp].cap-=cur_flow;
edge[temp^1].cap+=cur_flow;
}
max_flow+=cur_flow;
u=neck;
}
for(i=pe[u];i!=-1;i=edge[i].next)
{
if(edge[i].cap&&dist[u]==dist[edge[i].v]+1)
break;
}
if(i!=-1)
{
pe[u]=i;
pre[edge[i].v]=u;
u=edge[i].v;
}
else
{
if(0==--gap[dist[u]])
break;
pe[u]=head[u];
temp=t+1;
for(i=head[u];i!=-1;i=edge[i].next)
{
if(edge[i].cap)
temp=min(temp,dist[edge[i].v]);
}
dist[u]=temp+1;
++gap[dist[u]];
if(u!=s)
u=pre[u];
}
}
return max_flow;
}
bool vis[MAXN];
void dfs(int u)
{
if(vis[u])
return;
vis[u]=1;
for(int i=head[u];i!=-1;i=edge[i].next)
{
if(edge[i].cap)
dfs(edge[i].v);
}
}
int main()
{
int n,m,i,j;
while(scanf("%d%d",&n,&m)==2)
{
init();
int s=0,t=2*n+1;
int w;
for(i=1;i<=n;i++)
{
scanf("%d",&w);
add_edge(i+n,t,w);
}
for(i=1;i<=n;i++)
{
scanf("%d",&w);
add_edge(s,i,w);
}
int u,v;
while(m--)
{
scanf("%d%d",&u,&v);
add_edge(u,v+n,INF);
}
int ans=ISAP(s,t);
printf("%d\n",ans);
memset(vis,0,sizeof(vis));
dfs(s);
int num=0;
for(i=1;i<=n;i++)
{
if(vis[i+n])
num++;
if(!vis[i])
num++;
}
printf("%d\n",num);
for(i=1;i<=n;i++)
{
if(vis[i+n])
{
printf("%d +\n",i);
}
if(!vis[i])
printf("%d -\n",i);
}
}
return 0;
}