No.24 - POJ3469 邻接表网络流

最新推荐文章于 2019-11-01 22:42:46 发布

ShellDawn

最新推荐文章于 2019-11-01 22:42:46 发布

阅读量173

点赞数

分类专栏： POJ 文章标签： POJ

本文链接：https://blog.csdn.net/ShellDawn/article/details/96728710

版权

POJ 专栏收录该内容

37 篇文章 0 订阅

订阅专栏

经典网络流模版题：

首先，求最小割就是求最大流。

稀疏图，采用数组模拟二维动态链表。

$\red{小trick}$ ：在会产生回边的地方，直接加入双向边，互相索引，提高DFS速度

$\red{大trick}$ ：stdio.h 比 cstdio 快一倍

// ShellDawn
// POJ3469
// No.24

#include<stdio.h>
#include<string.h>
#include<queue>
#define MM(x,y) memset(x,y,sizeof(x))
#define INF 0x3f3f3f3f
using namespace std;

#define maxn 20010
#define maxm 200010
int n;
struct Edge{
    int to;
    int v;
    int rvs; // 指向反向边
    int pre;
};
// maxm每一条边都可能加回边，maxm*2
// maxn每个点都链接两个CPU，maxn*2
Edge E[maxm*2 + maxn*2]; 
int cnt = 1; // 0节点为所有根
int pre[maxn]; // 指向队尾
int V[maxn];

void add(int a,int b,int v,int r){
    E[cnt].to = b;
    E[cnt].v = v;
    E[cnt].pre = pre[a];
    E[cnt].rvs = r;
    pre[a]  = cnt;
    cnt++;
}

bool BFS(int s){
    MM(V,0);
    queue<int> q;
    q.push(s);
    V[s] = 1;
    while(!q.empty()){
        int now = q.front();
        q.pop();
        // 剪枝
        if(V[now] == V[n+1]) continue;
        for(int i=pre[now];i!=0;i=E[i].pre){
            if(E[i].v > 0 && V[E[i].to] == 0){
                V[E[i].to] = V[now] + 1;
                q.push(E[i].to);
            }
        }
    }
    if(V[n+1] > 0) return true;
    return false;
}

int DFS(int now,int minflow){
    if(now == n+1) return minflow;
    int flow = 0;
    for(int i=pre[now];i!=0 && flow < minflow;i=E[i].pre){
        if(E[i].to != n+1 && V[E[i].to] == V[n+1]) continue;
        if(V[E[i].to] == V[now] + 1 && E[i].v > 0){
            int f = DFS(E[i].to,min(E[i].v,minflow-flow));
            flow += f;
            E[i].v -= f;
            E[E[i].rvs].v += f;
        }
    }
    if(flow == 0) V[now] = 0;
    return flow;
}

int main(){
    int m;
    scanf("%d%d",&n,&m);
    cnt = 1;
    MM(pre,0);
    for(int i=1;i<=n;i++){
        int a,b;
        scanf("%d%d",&a,&b);
        add(0,i,a,0);
        add(i,n+1,b,0);
    }
    while(m--){
        int a,b,v;
        scanf("%d%d%d",&a,&b,&v);
        // 构造回边
        int t = cnt;
        add(a,b,v,t+1);
        add(b,a,v,t);
    }
    int ans = 0;
    while(BFS(0)) ans+=DFS(0,INF);
    printf("%d\n",ans);
    return 0;
}