用计算机怎么算lg1344,LG1344 「USACO4.4」Pollutant Control 最小割(示例代码)

问题描述

题解

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

两个询问。

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

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

(mathrm{Code})

#include

using namespace std;

template

void read(Tp &x){

x=0;char ch=1;int fh;

while(ch!='-'&&(ch>'9'||ch

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));

queueq;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;

}