#include <iostream>
#include <cstdio>
using namespace std;
const int oo=1e9;
const int mm=1111111;
const int mn=9999;
int node,src,dest,edge;
int ver[mm],flow[mm],next[mm];
int head[mn],work[mn],dis[mn],q[mn],vis[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;flow[edge]=c;next[edge]=head[u];head[u]=edge++;
ver[edge]=u;flow[edge]=0;next[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)//多路增广,从源点开始,用dfs从前一层向后一层反复寻找增广路
{
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;
}
long long Dinic_flow()//求最大流
{
int i,delta;
long long ret=0;
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 n,m,a,b,i,j;
long long ans,sum;
while(~scanf("%d%d",&n,&m))
{
prepare(n+2,0,n+1);
sum=0;ans=0;
for(i=1;i<=n;i++)
{
scanf("%d",&a);
if(a<0) addedge(i,dest,-a);
else {addedge(src,i,a);sum+=a;}
}
for(i=1;i<=m;i++)
{
scanf("%d%d",&a,&b);
addedge(a,b,oo);
}
long long tmp=Dinic_flow();
sum=sum-tmp;
dfs(0);
for(i=1;i<=n;i++) if(vis[i]) ans++;
printf("%I64d %I64d\n",ans,sum);
}
return 0;
}
poj2987
最新推荐文章于 2018-03-09 10:51:15 发布