[NOI2006]最大获利 BZOJ1497 - 最小割

题目链接：右转进入题目

可能我不适合做NOI的题

迄今为止做NOI的题有一半都是MLE 23333

M的范围看小了，不停地RE

建边：

(S,中转站,p),(中转站A&B,用户群,inf),(用户群,T,C)

跑一边最小割。

所有用户群的获利之和减去最小割的值就是答案。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<climits>
#define INF INT_MAX
#define MAXN 56000
#define MAXM 311000
#define build_edge(u,v,c) add_edge(u,v,c),add_edge(v,u,0);
using namespace std;
int lev[MAXN],h[MAXN],cur[MAXN],etop;
bool vis[MAXN];queue<int> q;
struct edges{
int to,resf,pre;
}e[MAXM];
int add_edge(int u,int v,int resc)
{
etop++;
e[etop].to=v;
e[etop].resf=resc;
e[etop].pre=h[u];
h[u]=etop;return 0;
}
bool bfs(int s,int t)
{
memset(vis,false,sizeof(vis));
memset(lev,0,sizeof(lev));
while(!q.empty()) q.
4000
pop();
q.push(s);lev[s]=0;vis[s]=true;
while(!q.empty())
{
int x=q.front();q.pop();
for(int i=h[x];i;i=e[i].pre)
if(!vis[e[i].to]&&e[i].resf)
{
vis[e[i].to]=true;
lev[e[i].to]=lev[x]+1;
q.push(e[i].to);
}
}
return vis[t];
}
int dfs(int s,int t,int a)
{
if(s==t||!a) return a;
int flow=0,f;
for(int &i=cur[s];i;i=e[i].pre)
if(lev[e[i].to]==lev[s]+1&&(f=dfs(e[i].to,t,min(a,e[i].resf)))>0)
{
a-=f;flow+=f;e[i].resf-=f;
e[((i-1)^1)+1].resf+=f;
if(!a) break;
}
return flow;
}
int main()
{
int n,m;scanf("%d%d",&n,&m);
int s=0,t=n+m+1,totans=0;
for(int i=1;i<=n;i++)
{
int p;scanf("%d",&p);
build_edge(s,i,p);
}
for(int i=n+1;i<=m+n;i++)
{
int a,b,c;scanf("%d%d%d",&a,&b,&c);
totans+=c;
build_edge(i,t,c);
build_edge(a,i,INF);
build_edge(b,i,INF);
}
while(bfs(s,t))
{
for(int i=s;i<=t;i++)
cur[i]=h[i];
totans-=dfs(s,t,INF);
}
printf("%d\n",totans);
return 0;
}