Luogu P1344 追查坏牛奶

Luogu P1344 追查坏牛奶

---

题目

题目描述

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

输入输出格式

输入格式：
第一行: 两个整数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

题目大意
给定一个n个点m条边的有向图，要求砍断一些边，使得从1点走到n点的行为变得不可能，求和砍的边的最小的权值和砍的边的数量。

思路

考虑把图中的点分成2个集合S，T且起点和终点不在同一个集合中,如果把“起点在S中，终点在T中”的边全部删除，就无法从s到达t了。这样的集合划分
(S,T)称为一个s-t割他的容量$C(S,T)$定义为：

$$C(S,T)=\sum_{u\in S 且 v\in T} c(u,v)$$
即所有起点在S集合内，终点在T集合内的边的容量之和。
显然，可以找到一种集合划分方式使得 $C(S,T)$最小。我把最小的 $C(S,T)$记做 $mC(S,T)$(最小割)
那么我们有结论：

1. 对于任意的s-t流f和任意的s-t割(S,T),有$f \leq C(S,T)$.
2. 存在情况，使得$f=C(S,T)$（思考一下，我们怎么深搜求最大流？）

解释：当搜索没有找到任何s-t道路时，把已访问结点集合看成S，令$T=V-S$(此时割$(S,T)$就是图的最小割。)则在残量网络中S和T分离，因此在原图中跨越S和T的所有边均满载（这样的边才不会存在于残量网络中），且没有从T回到S的流量，因此成立$f≤(S,T)$成立。

所以：最大流等于最小割。
知道这个结论，第一问代码便容易实现。
对于第二问，只需要对每一个容量都乘以2000(或任意较大数)再+1，那么答案就等于求出来的最小割整除2000，砍的边数就等于求出来的最大流mod2000了。

代码

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
const long long maxn=1000*2+10;
inline void read(long long &misaka) {
    long long ans=0,flag=1;
    char ch=getchar();
    while(ch<'0' || ch>'9') {
        if(ch=='-')
            flag=-1;
        ch=getchar();
    }
    while(ch>='0' && ch<='9') {
        ans=ans*10+(long long)(ch-'0');
        ch=getchar();
    }
    misaka=ans*flag;
}


long long n,m;
struct Edge{
    long long nxt,to;
    long long c;
}edge[maxn];
long long head[maxn],cnt=0;
inline void add(long long u,long long v,long long w){
    edge[cnt].c=w;
    edge[cnt].to=v;
    edge[cnt].nxt=head[u];
    head[u]=cnt++;
}

long long pf=0,ans;
bool vis[maxn];
inline bool dfs(long long x,long long f){
    if(x>=n){
        ans+=f;
        pf=f;
        return true;
    }
    vis[x]=true;
    for(long long i=head[x];~i;i=edge[i].nxt){
        if(edge[i].c>0 && !vis[edge[i].to]){
            if(dfs(edge[i].to,min(f,edge[i].c))){
                edge[i].c=edge[i].c-pf;
                edge[i^1].c=edge[i^1].c-(-pf);
                return true;
            }
        }
    }
    return false;
}

int main(){
    ios::sync_with_stdio(false);
    memset(head,-1,sizeof(head));
    read(n);read(m);
    for(long long i=1;i<=m;i++){
        long long a,b,c;
        read(a);read(b);read(c);
        add(a,b,c*2000+1);
        add(b,a,0);
    }
    pf=0;memset(vis,0,sizeof(vis));
    while(dfs(1,100000000)){
        pf=0;
        memset(vis,0,sizeof(vis));
    }
    cout<<ans/2000<<" "<<ans%2000;
    return 0;
}