最小割

最小割定义

• 这篇文章解释的挺好!

---

1. 割 : 点的划分方式, 将图中的所有点划分成2个集合,源点s,汇点t, $s\in S,t\in T$.
2. 割的容量 ： 表示所有从S到T的边的容量之和,$c(S,T)={\Sigma }_{u\in S,v\in T}c(u,v)$
3. 最小割问题: 求一个割的方法使得割的容量最小

最大流 = 最小割

• 求最小割和最大流可以看做是一个问题

---

最小割例题 luogu P1344

问题

1. 最小割的值
2. 最小割边的个数

直接求最小割

• 第一问直接求最大流
• 对于第二问重新建图 将边权改为1,那么最大流=最小割=最小割的边数

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

const int maxn = 50;
const int maxm = 1e4+5;
const ll inf = 1e9+1;

struct Edge{
    int v,next;
    ll c;
} edge[maxm];

int n,m,tot,source,sink;
int head[maxn],dep[maxn];

void add(int u,int v,ll c) {
    edge[tot].v=v; edge[tot].c=c;
    edge[tot].next=head[u]; head[u]=tot++;
}

void addedge(int u,int v,ll c) {
    add(u,v,c);
    add(v,u,0);
}

bool bfs() {
    queue<int> que;
    memset(dep,-1,sizeof dep);
    dep[source] = 0; // 深度标记数组
    que.push(source);
    while (!que.empty()) {
        int cur = que.front(); que.pop();
        for (int i=head[cur]; i!=-1; i=edge[i].next) {
            int to = edge[i].v;
            if (edge[i].c>0 && dep[to]==-1) {
                dep[to] = dep[cur] + 1;
                que.push(to);
            }
        }
    }
    return dep[sink] != -1;
}

ll dfs(int u,ll delta) {
    if (u == sink)
        return delta;
    ll flow = 0;
    for (int i=head[u]; i!=-1; i=edge[i].next) {
        int to = edge[i].v;
        if (edge[i].c>0 && dep[to]==dep[u]+1) {
            ll temp = dfs(to,min(delta-flow,edge[i].c));
            edge[i].c -= temp;
            edge[i^1].c += temp;
            flow += temp;
        }
    }
    if (!flow) dep[u] = -1;
    return flow;
}

ll Dinic() {
    ll ans = 0;
    while(bfs()) {
        while (1) {
            ll temp = dfs(source,inf);
            if (!temp) break;
            ans += temp;
        }
    }
    return ans;
}

struct Storage {
    int u,v;
} storage[maxm];


int main() {
    memset(head,-1,sizeof head);
    scanf("%d%d",&n,&m);
    source = 1;
    sink = n;
    for (int i=1; i<=m; ++i) {
        int u,v;
        ll c;
        scanf("%d %d %lld",&u,&v,&c);
        addedge(u,v,c);
        storage[i].u = u;
        storage[i].v = v;
    }
    ll maxflow = Dinic();
    // 重新建图求 最小割的数量

    memset(head,-1,sizeof head);
    memset(edge,0,sizeof edge);
    tot = 0;

    for (int i=1; i<=m; ++i)
        addedge(storage[i].u,storage[i].v,1);

    ll cut_count = Dinic();

    cout << maxflow << " " << cut_count << endl;
    return 0;
}

最小割转最大流

• 优秀思路!!!

1. 建图时边权$w=k\ast w+1$, 若$S$为最小割边集 则 $w_i\in S$

$w_1+w_2+w_3+...+w_n=ans$

$w_1\ast k+w_2\ast k+w_3\ast k+...+w_n\ast k=ans\ast k$

$(w_1\ast k+1)+(w_2\ast k+1)+(w_3\ast k+1)+...+(w_n\ast k+1)=ans\ast k+n=an{s}_1$

• 最小割: $ans=an{s}_1/k$
• 最大流: $cut=an{s}_1\%k$

---

1. 注意long long

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

const int maxn = 50;
const int maxm = 1e4+5;
const ll mod = 2000;
const ll inf = 1e9+1;

struct Edge {
    int v,next;
    ll c;
}edge[maxm];

int head[maxn],dep[maxn];
int n,m,tot,source,sink;

void add(int u,int v,ll c){
    edge[tot].v=v; edge[tot].c=c;
    edge[tot].next=head[u]; head[u]=tot++;
}
void addedge(int u,int v,ll c) {
    add(u,v,c);
    add(v,u,0);
}

bool bfs() {
    queue<int> que;
    memset(dep,-1,sizeof dep);
    dep[source] = 0;
    que.push(source);
    while(!que.empty()) {
        int cur = que.front(); que.pop();
        for (int i=head[cur]; i!=-1; i=edge[i].next) {
            int to = edge[i].v;
            if (edge[i].c>0 && dep[to] == -1) {
                dep[to] = dep[cur] + 1;
                que.push(to);
            }
        }
    }
    return dep[sink] != -1;
}

ll dfs(int u,int delta) {
    if (u == sink)
        return delta;
    ll flow = 0;
    for (int i=head[u]; i!=-1; i=edge[i].next) {
        int to = edge[i].v;
        if (edge[i].c>0 && dep[to]==dep[u]+1) {
            ll temp = dfs(to,min(delta-flow,edge[i].c));
            edge[i].c -= temp;
            edge[i^1].c += temp;
            flow += temp;
        }
    }
    if (!flow) dep[u] = -1;
    return flow;
}

ll Dinic() {
    ll ans=0;
    while (bfs()) {
        while(1) {
            ll temp = dfs(source,inf);
            if (!temp) break;
            ans += temp;
        }
    }
    return ans;
}

int main() {
    memset(head,-1,sizeof head);
    scanf("%d%d",&n,&m);
    for (int i=1; i<=m; ++i) {
        int u,v;
        ll c;
        scanf("%d%d%lld",&u,&v,&c);
        addedge(u,v,c*mod+1);
    }
    source = 1;
    sink = n;
    ll ans = Dinic();
    cout << ans/mod << " " << ans%mod << endl;
    return 0;
}

• 2019/11/18 0:36:23 真晚