LG1344 「USACO4.4」Pollutant Control 最小割

问题描述

LG1344

---

题解

我太菜了，我一开始竟然没有看出这是个最小割裸题。。。

两个询问。

第一个询问，直接跑最小割就好了。

第二个询问，建图的时候边权建 (1) ，代表割掉这条边需要 (1) 的代价。

---

(mathrm{Code})

#include<bits/stdc++.h>
using namespace std;

template <typename Tp>
void read(Tp &x){
	x=0;char ch=1;int fh;
	while(ch!='-'&&(ch>'9'||ch<'0')) ch=getchar();
	if(ch=='-') ch=getchar(),fh=-1;
	else fh=1;
	while(ch>='0'&&ch<='9') x=(x<<1)+(x<<3)+ch-'0',ch=getchar();
	x*=fh;
}

const int maxn=37;
const int maxm=2007;

int n,m,S,T;

struct Graph{
	int Head[maxn],to[maxm],Next[maxm],w[maxm],tot=1;
	Graph(){
		memset(Head,0,sizeof(Head));
		memset(Next,0,sizeof(Next));
	}
	void add(int x,int y,int z){
		to[++tot]=y,Next[tot]=Head[x],Head[x]=tot,w[tot]=z;
	}
}G[2];

int d[maxn];

bool bfs(int gp){
	memset(d,0,sizeof(d));
	queue<int>q;q.push(S);d[S]=1;
	while(!q.empty()){
		int x=q.front();q.pop();
		for(int i=G[gp].Head[x];i;i=G[gp].Next[i]){
			int y=G[gp].to[i];
			if(d[y]||!G[gp].w[i]) continue;
			q.push(y);d[y]=d[x]+1;
			if(y==T) return true;
		}
	}
	return false;
}

int dfs(int x,int flow,int gp){
	if(x==T) return flow;
	int rest=flow;
	for(int i=G[gp].Head[x];i&&rest;i=G[gp].Next[i]){
		int y=G[gp].to[i];
		if(d[y]!=d[x]+1||!G[gp].w[i]) continue;
		int k=dfs(y,min(rest,G[gp].w[i]),gp);
		if(!k) d[y]=0;
		else{
			G[gp].w[i]-=k,G[gp].w[i xor 1]+=k;
			rest-=k;
		}
	}
	return flow-rest;
}

int ans;

int main(){
	read(n);read(m);
	for(int i=1,x,y,z;i<=m;i++){
		read(x);read(y);read(z);
		G[1].add(x,y,z);G[1].add(y,x,0);
		G[0].add(x,y,1);G[0].add(y,x,0);
	}
	S=1,T=n;
	while(bfs(1)){
		int t;
		while(t=dfs(S,0x3f3f3f3f,1)) ans+=t;
	}
	printf("%d ",ans);ans=0;
	while(bfs(0)){
		int t;
		while(t=dfs(S,0x3f3f3f3f,0)) ans+=t;
	}
	printf("%d
",ans);
	return 0;
}