hdu_4289

跟3491基本上一样,最小割,不同的是该题,原点和会点都可以设置边
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
#define MAXV    402
#define MAXE    80001
#define INF        0x7fffffff
using namespace std;

struct Edge{
    int v, val, next;
}edge[MAXE];

int e_cnt;
int head[MAXV], cur[MAXV], pre[MAXV], dis[MAXV], gap[MAXV];

void insert_arc(int u, int v, int w)
{
    edge[e_cnt].v = v; edge[e_cnt].val = w;
    edge[e_cnt].next = head[u];    head[u] = e_cnt ++;
    swap(u, v);
    edge[e_cnt].v = v; edge[e_cnt].val = 0;
    edge[e_cnt].next = head[u];    head[u] = e_cnt ++;
}

void bfs(int des)
{
    int u, v;
    queue<int> q;
    memset(gap, 0, sizeof(gap));
    memset(dis, -1, sizeof(dis));
    dis[des] = 0;    gap[0] = 1;
    q.push(des);
    while( !q.empty() ) {
        u = q.front(); q.pop();
        for(int i = head[u]; i != -1; i = edge[i].next) {
            v = edge[i].v;
            if( edge[i].val || dis[v] != -1 )     continue;
            dis[v] = dis[u]+1;    q.push(v);    gap[dis[v]] ++;
        }
    }
}

int i_sap(int src, int des, int VN)
{
    int u(src), v, min_dis, max_flow(0), path_flow(INF);
    bfs(src);    pre[u] = u;
    for(int i = 0; i <= VN; i ++)
        cur[i] = head[i];
    while( dis[src] < VN ) {
loop:
        for(int &i = cur[u]; i != -1; i = edge[i].next) {
            v = edge[i].v;
            if( !edge[i].val || dis[u] != dis[v]+1 )    continue;
            pre[v] = u;    path_flow = min(path_flow, edge[i].val);
            u = v;
            if( v == des ) {
                for(u = pre[u]; v != src; v = u, u = pre[u]) {
                    edge[cur[u]].val -= path_flow;    edge[cur[u]^1].val += path_flow;
                }
                max_flow += path_flow; path_flow = INF;
            }
            goto loop;
        }
        min_dis = VN;
        for(int i = head[u]; i != -1; i = edge[i].next) {
            v = edge[i].v;
            if( edge[i].val && min_dis > dis[v] ) {
                cur[u] = i; min_dis = dis[v];
            }
        }
        gap[dis[u]] --;
        if( !gap[dis[u]] )    break;
        dis[u] = min_dis+1;
        gap[dis[u]] ++; u = pre[u];
    }
    return max_flow;
}

int main(int argc, char const *argv[])
{
    //freopen("test.in", "r", stdin);
    int vertex, arc, src, des, v, f, t, ss, ee;
    while( ~scanf("%d %d %d %d", &vertex, &arc, &src, &des) ) {
        e_cnt = 0;    memset(head, -1, sizeof(head));
        ss = 0;    ee = 2*vertex+1;
        insert_arc(ss, src, INF);    insert_arc(des+vertex, ee, INF);

        for(int i = 1; i <= vertex; i ++) {
            scanf("%d", &v);
            insert_arc(i, i+vertex, v);    //insert_arc(i+vertex, i, ve_val[i]);
        }

        for(int i = 0; i < arc; i ++) {
            scanf("%d %d", &f, &t);
            insert_arc(f+vertex, t, INF);
            insert_arc(t+vertex, f, INF);
        }

        printf("%d\n", i_sap(ss, ee, ee));
    }
    return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值