BZOJ 1497 最大权闭合子图

#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<string>
#include<cstring>
#include<iostream>
#include<iomanip>
#include<algorithm>
#define INF 1000000000
using namespace std;
struct bian
{
int r,f;
}a[2000000];
int tot=1;
int d[1000000];
int fir[2000000];
int nex[2000000];
int s,t;

void add_edge(int l,int r,int f)
{
a[++tot].r=r;
a[tot].f=f;
nex[tot]=fir[l];
fir[l]=tot;
}
bool bfs()
{
static int dui[1000000];
int top=1,my_final=2;
memset(d,-1,sizeof(d));
d[s]=0;
dui[top]=s;
while(top<my_final)
{
int u=dui[top];
for(int o=fir[u];o!=0;o=nex[o])
{
if(d[a[o].r]!=-1 || !a[o].f) continue;
dui[my_final++]=a[o].r;
d[a[o].r]=d[u]+1;
if(a[o].r==t) return true;
}
top++;
}
return false;
}
int dinic(int u,int flow)
{
if(u==t) return flow;
int left=flow;
for(int o=fir[u];o!=0&&left;o=nex[o])
{
if(d[a[o].r]==d[u]+1 && a[o].f)
{
int temp=dinic(a[o].r,min(left,a[o].f));
if(temp==0) d[a[o].r]=-1;
left-=temp;
a[o].f-=temp;
a[o^1].f+=temp;
}
}
return flow-left;
}
int main()
{
int n,m;
scanf("%d%d",&m,&n);
s=0;t=999999;
for(int i=1;i<=m;i++)
{
int x;
scanf("%d",&x);
add_edge(n+i,t,x);
add_edge(t,n+i,0);
}
int ans=0;
for(int i=1;i<=n;i++)
{
int x,y,z;
scanf("%d%d%d",&x,&y,&z);
ans+=z;
add_edge(s,i,z);
add_edge(i,s,0);
add_edge(i,n+x,INF);
add_edge(n+x,i,0);
add_edge(i,n+y,INF);
add_edge(n+y,i,0);
}
while(bfs()) ans-=dinic(s,INF);
cout<<ans;
return 0;
}