[xsy1232]Magic

题意：一个无向图，每个点有$a_i,b_i$，对任意点$i$你都可以花费$b_i$的费用将$a_i$变为$0$，最后你还要付出$\sum\limits_{i=1}^n\max\limits_{(i,j)}a_j$的费用，求最小费用

把每个点拆成两个点，对每条边新建一个点，将一个点的出边按终点的$a$从大到小排序后建图如下

只有让某个前缀$a_i$全部变为$0$，一个点产生的费用才会改变，这个建图保证了如果要割$a_i$就必须要割掉$b_{1\cdots i-1}$，并且不会割掉两个$a$

听myh说是套路但我不会...

#include<stdio.h>
#include<algorithm>
#include<vector>
#include<string.h>
using namespace std;
const int inf=2147483647;
int h[52010],nex[204010],to[204010],cap[204010],M=1,S,T;
void ins(int a,int b,int c){
	M++;
	to[M]=b;
	cap[M]=c;
	nex[M]=h[a];
	h[a]=M;
}
void add(int a,int b,int c){
	ins(a,b,c);
	ins(b,a,0);
}
int dis[52010],q[52010];
bool bfs(){
	int head,tail,x,i;
	memset(dis,-1,sizeof(dis));
	head=tail=1;
	q[1]=S;
	dis[S]=0;
	while(head<=tail){
		x=q[head++];
		for(i=h[x];i;i=nex[i]){
			if(cap[i]&&dis[to[i]]==-1){
				dis[to[i]]=dis[x]+1;
				if(to[i]==T)return 1;
				q[++tail]=to[i];
			}
		}
	}
	return 0;
}
int cur[52010];
int dfs(int x,int flow){
	if(x==T)return flow;
	int us=0,i,t;
	for(i=cur[x];i&&flow;i=nex[i]){
		if(cap[i]&&dis[to[i]]==dis[x]+1){
			t=dfs(to[i],min(flow,cap[i]));
			cap[i]-=t;
			cap[i^1]+=t;
			us+=t;
			flow-=t;
			if(cap[i])cur[x]=i;
		}
	}
	if(us==0)dis[x]=-1;
	return us;
}
int dicnic(){
	int ans=0;
	while(bfs()){
		memcpy(cur,h,sizeof(h));
		ans+=dfs(S,inf);
	}
	return ans;
}
vector<int>g[1010];
int a[1010],b[1010];
bool cmp(int x,int y){return a[x]>a[y];}
int main(){
	int n,m,i,x,y,N;
	scanf("%d%d",&n,&m);
	for(i=1;i<=n;i++)scanf("%d",a+i);
	for(i=1;i<=n;i++)scanf("%d",b+i);
	for(i=1;i<=m;i++){
		scanf("%d%d",&x,&y);
		g[x].push_back(y);
	}
	N=2*n;
	S=N+m+1;
	T=S+1;
	for(i=1;i<=n;i++)add(S,i,inf);
	for(x=1;x<=n;x++){
		sort(g[x].begin(),g[x].end(),cmp);
		for(i=0;i<(int)g[x].size();i++){
			N++;
			add(i?N-1:x,N,a[g[x][i]]);
			add(N,g[x][i]+n,inf);
		}
	}
	for(i=1;i<=n;i++)add(i+n,T,b[i]);
	printf("%d",dicnic());
}

转载于:https://www.cnblogs.com/jefflyy/p/9552771.html