网络流（最小割）————[USACO4.4]追查坏牛奶Pollutant Control

这题差不多可以说是一个模板，它能让我们简单，快捷求出最小割以及割的边数。

只要看出题目为最小割。几乎都可以套用此模板。

首先明确你要求的量：最小割与最小割边数。

我们可以建两次图。但那样太过麻烦，也占空间。

那么根据定义最小割的值等于最大流。只需建图时将边权乘上一个大数+1，这样我们能求得最大流ans，则最小割的值为ans/a，割的边数为ans%a。

那么为什么呢？

证明：

因为最小割静流cost1+cost2+cost3…+costn=ans，等式左右同乘a在加1，

有wcost1*a+1+cost2*a+1+cost3*a+1…+costn*a+1=ans*a+k，

k不就是我们最小割的边数吗，因为总边数不大于1000，所以k<1000，

这就说明，我们的大数是跟题目定义的，边数有多少，大数必定比边数多。

题目描述

你第一天接手三鹿牛奶公司就发生了一件倒霉的事情：公司不小心发送了一批有三聚氰胺的牛奶。很不幸，你发现这件事的时候，有三聚氰胺的牛奶已经进入了送货网。这个送货网很大，而且关系复杂。你知道这批牛奶要发给哪个零售商，但是要把这批牛奶送到他手中有许多种途径。送货网由一些仓库和运输卡车组成，每辆卡车都在各自固定的两个仓库之间单向运输牛奶。在追查这些有三聚氰胺的牛奶的时候，有必要保证它不被送到零售商手里，所以必须使某些运输卡车停止运输，但是停止每辆卡车都会有一定的经济损失。你的任务是，在保证坏牛奶不送到零售商的前提下，制定出停止卡车运输的方案，使损失最小。

输入输出格式

输入格式：

第一行: 两个整数N(2<=N<=32)、M(0<=M<=1000), N表示仓库的数目，M表示运输卡车的数量。仓库1代 表发货工厂，仓库N代表有三聚氰胺的牛奶要发往的零售商。 第2..M+1行: 每行3个整数Si,Ei,Ci。其中Si,Ei表示这 辆卡车的出发仓库，目的仓库。Ci(0 <= C i <= 2,000,000) 表示让这辆卡车停止运输的损失。

输出格式：

两个整数C、T：C表示最小的损失，T表示在损失最小的前提下，最少要停止的卡车数。

输入输出样例

输入样例#1：

4 5
1 3 100
3 2 50
2 4 60
1 2 40
2 3 80

输出样例#1： 

60 1

上代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<cstdlib>
using namespace std;
const int INF=0x3f3f3f3f;
const int maxn=100001;
long long s,t,n,m;
struct edge{
	long long to,next,cap,op;
}e[maxn*2];
long long head[maxn],cnt,depth[maxn],cur[maxn];
inline void add_edge(long long x,long long y,long long flow){
	e[++cnt].to=y;
	e[cnt].cap=flow;
	e[cnt].next=head[x];
	head[x]=cnt;
	
	e[++cnt].to=x;
	e[cnt].cap=0;
	e[cnt].next=head[y];
	head[y]=cnt;
	
	e[head[x]].op=head[y];
	e[head[y]].op=head[x];
}

inline bool bfs(){
	memset(depth,0,sizeof(depth));
	queue<long long>q;
	q.push(s);
	depth[s]=1;
	while(!q.empty()){
		long long u=q.front();
		q.pop();
		for(long long i=head[u];i;i=e[i].next){
			long long to=e[i].to;
			if(e[i].cap&&!depth[to]){
				depth[to]=depth[u]+1;
				q.push(to);
				if(to==t) 
					return true;
			}
		}
	}
	return false;
}
inline long long dfs(long long now,long long f){
	if(now==t) return f;
	long long rest=f;
	for(long long i=cur[now];i&&rest;i=e[i].next){
		cur[now]=i;
		long long to=e[i].to;
		if(e[i].cap&&depth[to]==depth[now]+1){
			long long flow=dfs(to,min(rest,e[i].cap));
			rest-=flow;
			e[i].cap-=flow;
			e[e[i].op].cap+=flow;
		}
	}
	return f-rest;
}
inline long long dinic(long long start){
	long long ans=0;
	while(bfs()){
		memset(cur,0,sizeof(cur));
		for(long long i=1;i<=n;i++) cur[i]=head[i];
        ans+=dfs(start,INF);		
	}
	return ans;
}
int main(){
	cin>>n>>m;
	s=1,t=n;
    long long x,y,z;
    for(long long i=0;i<m;i++){
        scanf("%lld%lld%lld",&x,&y,&z);
        add_edge(x,y,z*1002+1);
    }
    long long ans=dinic(s);
    cout<<ans/1002<<" "<<ans%1002;
    return 0;
}