省选专练 [NOI2006]最大获利（profit）

模板

最大权闭合图。

对于边权拆除一个点和两个连上的无穷大边

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<queue>
using namespace std;
const int INF=1e9+7;
const int N=1e6+100;
struct Front_star{
	int u,v,w,nxt;
}e[N*4];
int cnt=1;
int first[N]={0};
void addedge(int u,int v,int w){
	cnt++;
	e[cnt].u=u;
	e[cnt].v=v;
	e[cnt].w=w;
	e[cnt].nxt=first[u];
	first[u]=cnt;
} 
void add(int u,int v,int w){
//	cout<<"add "<<u<<" "<<v<<" "<<w<<'\n';
	addedge(u,v,w);
	addedge(v,u,0);
}
int n,m;
int pos=0;
int S=0;
int T=60000;
int d[N]={0};
bool bfs(){
	queue<int> q;
	memset(d,-1,sizeof(d));
	d[S]=1;
	q.push(S);
	while(!q.empty()){
		int x=q.front();
		q.pop();
		for(int i=first[x];i;i=e[i].nxt){
			int v=e[i].v;
			if(e[i].w&&d[v]==-1){
				d[v]=d[x]+1;
				q.push(v);
			}
		}
	}
	return d[T]!=-1;
}
int dfs(int st,int ed,int nowdat){
	int dat=0;
	if(st==ed){
		return nowdat;
	}
	for(int i=first[st];i;i=e[i].nxt){
		int v=e[i].v;
		if(d[v]==d[st]+1&&e[i].w){
			int now=min(nowdat-dat,e[i].w);
			now=dfs(v,ed,now);
			dat+=now;
			e[i].w-=now;
			e[i^1].w+=now;
			if(nowdat==dat){
				return dat;
			}
		}
	}
	if(!dat){
		d[st]=-2;
	}
	return dat;
}
int Dinic(){
	int ans=0;
	while(bfs()){
		ans+=dfs(S,T,INF);
	}
	return ans;
}
int main(){
	scanf("%d%d",&n,&m);
//	cout<<n<<" "<<m<<'\n';
	for(int i=1;i<=n;i++){
		int p;
		scanf("%d",&p);
//		int u=i;
		add(i,T,p);
	}
//	cout<<n<<" "<<m<<endl;
//	cout<<"ready"<<'\n';
	for(int i=1;i<=m;i++){
//		cout<<i<<'\n';
		int a,b,c;
//		cout<<"hello world"<<'\n';
		scanf("%d%d%d",&a,&b,&c);
		pos+=c;
		add(i+n,a,INF);
		add(i+n,b,INF);
		add(S,i+n,c);
	}
	cout<<pos-Dinic();
}

转载于:https://www.cnblogs.com/Leo-JAM/p/10079278.html