这题差不多可以说是一个模板,它能让我们简单,快捷求出最小割以及割的边数。
只要看出题目为最小割。几乎都可以套用此模板。
首先明确你要求的量:最小割与最小割边数。
我们可以建两次图。但那样太过麻烦,也占空间。
那么根据定义最小割的值等于最大流。只需建图时将边权乘上一个大数+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表示在损失最小的前提下,最少要停止的卡车数。
输入输出样例
4 5 1 3 100 3 2 50 2 4 60 1 2 40 2 3 80
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;
}