【JZOJ6229】【20190621】san

题目

\(n\)个点\(m\)条边的有向图，每个点有点权

你可以选择拓扑序的一个区间的

最大化点权和

$n \le 50 \ , m \le \frac{n*(n-1)}{2} , 0 \le |a_i| \le 200 $

题解

• 一条路径上的点一定是：不选-选-不选

• 把所有正权加起来为sum，然后建最小割模型，割哪段表示选在

• 连边：

  $ (w_i \gt 0) (S,i,w_i) \ (i',T,w_i) $

  $ (w_i \lt 0) (i,i',-w_i) $

• 边有向边\((a,b)\)相当于限制了a割的段不能在b后面

  那连边：
  \((a,b,inf) \ , \ (a',b',inf)\)

• 答案是sum-flow

#include<bits/stdc++.h>
#define inf 0x3f3f3f3f

using namespace std;

const int N=110;
int n,m,vis[N],cur[N],hd[N],o,S,T,d[N];
struct Edge{int v,nt,f;}E[N*N];	
int q[N],head,tail;
void adde(int u,int v,int f){
	E[o]=(Edge){v,hd[u],f};hd[u]=o++;
	E[o]=(Edge){u,hd[v],0};hd[v]=o++;
}

bool bfs(){
	for(int i=S;i<=T;++i)vis[i]=0,d[i]=0;
	head=tail=0;vis[q[++tail]=S]=1;d[S]=1;
	while(head<tail){
		int u=q[++head];
		for(int i=hd[u];~i;i=E[i].nt)if(E[i].f){
			int v=E[i].v;
			if(vis[v])continue;
			vis[v]=1;
			d[v]=d[u]+1;
			q[++tail]=v;
			if(v==T)return true;
		}
	}return false;
}

int dfs(int u,int F){
	if(u==T||!F)return F;
	int flow=0,f;
	for(int i=cur[u];~i;i=E[i].nt){
		int v=E[cur[u]=i].v;
		if(d[v]==d[u]+1&&(f=dfs(v,min(F,E[i].f)))){
			flow+=f,F-=f;
			E[i].f-=f,E[i^1].f+=f;
			if(!F)break;
		}
	}
	if(!flow)d[u]=0;//记得加上这条优化
	return flow;
}

int main(){
	freopen("san.in","r",stdin);
	freopen("san.out","w",stdout);
	scanf("%d%d",&n,&m);
	int ans=0;S=0;T=n*2+1;
	for(int i=S;i<=T;++i)hd[i]=-1;
	for(int i=1,w;i<=n;++i){
		scanf("%d",&w);
		if(w>0)ans+=w,adde(S,i,w),adde(i+n,T,w);
		else adde(i,i+n,-w);
	}
	for(int i=1,u,v;i<=m;++i){
		scanf("%d%d",&u,&v);
		adde(u,v,inf);
		adde(u+n,v+n,inf);
	}
	while(bfs()){
		for(int i=S;i<=T;++i)cur[i]=hd[i];
		ans-=dfs(S,inf);
	}
	cout<<ans<<endl;
	return 0;
}