您的位置:首页 > 理论基础 > 计算机网络

POJ 2125 Destring the Graph(网络流)

2016-03-05 00:07 507 查看
题意:

跟你一个有向图,每次你可以删除一个点所有出边或者是入边,花费分别是w+,w-,问你怎么删可以把所有的边全部删掉并且花费最小。

啊,网络流的姿势实在太多,这个是最小点权覆盖,你可以想象,每一条边都会被一个点所覆盖,也就是删掉,所以我们把点拆成两个,边都会

被一个点覆盖,出边被出点,入边被入点覆盖,这样想就清楚了,至于最小点权独立集,我们只需要把出点连源点,入点连汇点,其他的按照出

点到汇点相连就可以了,至于为什么,你把图画出来就懂了。最后求一个最大流就是答案,至于输出方案,我们只需要做一次dfs,如果这条边

没有流量了,说明对应的流量是从这里过的,那么就是一种方案。

代码:

[code]//
//  Created by  CQU_CST_WuErli
//  Copyright (c) 2016 CQU_CST_WuErli. All rights reserved.
//
#include <iostream>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <cctype>
#include <cmath>
#include <string>
#include <vector>
#include <list>
#include <map>
#include <queue>
#include <stack>
#include <set>
#include <algorithm>
#include <sstream>
#define CLR(x) memset(x,0,sizeof(x))
#define OFF(x) memset(x,-1,sizeof(x))
#define MEM(x,a) memset((x),(a),sizeof(x))
#define BUG cout << "I am here" << endl
#define lookln(x) cout << #x << "=" << x << endl
#define SI(a) scanf("%d",&a)
#define SII(a,b) scanf("%d%d",&a,&b)
#define SIII(a,b,c) scanf("%d%d%d",&a,&b,&c)
#define rep(flag,start,end) for(int flag=start;flag<=end;flag++)
#define Rep(flag,start,end) for(int flag=start;flag>=end;flag--)
#define Lson l,mid,rt<<1
#define Rson mid+1,r,rt<<1|1
#define Root 1,n,1
const int INF_INT=0x3f3f3f3f;
const long long INF_LL=0x7fffffff;
const int MOD=1e9+7;
const double eps=1e-9;
const double pi=acos(-1);
typedef long long  ll;
using namespace std;

const int N=5e4;
const int M=400;

int n,m;
int in[M],out[M];
int vis[M];

struct Dinic
{
    const static int NE=400;
    const static int NV=5e4;
    int source,sink;
    int pnt
,nxt
,head[M],cap
,cnt;
    int iter[M],level[M];
    void init(int s=0,int t=NE) {
        source=s,sink=t;
        cnt=0;
        OFF(head);
    }
    void add_edge(int u,int v,int c) {
        pnt[cnt]=v;nxt[cnt]=head[u];cap[cnt]=c;head[u]=cnt++;
    }
    bool bfs(int s,int t) {
        OFF(level);
        level[s]=0;
        queue<int> q;
        q.push(s);
        while(q.size()) {
            int x=q.front();q.pop();
            for (int i=head[x];~i;i=nxt[i]) {
                int v=pnt[i];
                if (level[v]==-1 && cap[i]) {
                    level[v]=level[x]+1;
                    q.push(v);
                }
            }
        }
        return level[t]!=-1;
    }
    int dfs(int u,int t,int f) {
        if (u==t || !f) return f;
        int left=f;
        for (int i=iter[u];~i;i=nxt[i]) {
            int v=pnt[i];
            if (level[v]==level[u]+1 && cap[i]) {
                int d=dfs(v,t,min(cap[i],left));
                iter[u]=i;
                cap[i]-=d;
                cap[i^1]+=d;
                left-=d;
                if (!left) return f;
            }
        }
        level[u]=-1;
        return f-left;
    }
    int solve(int s,int t) {
        int Max_flow=0;
        while (bfs(s,t)) {
            for (int i=s;i<=t;i++) iter[i]=head[i];
            Max_flow+=dfs(s,t,INF_INT);
        }
        return Max_flow;
    }
    void dfs(int u) {
        vis[u]=1;
        for (int i=head[u];~i;i=nxt[i]) {
            int v=pnt[i];
            if (cap[i]==0 || vis[v]) continue;
            dfs(v);
        }
    }
}dinic;

int main(int argc, char const *argv[]) {
#ifdef LOCAL
    freopen("C:\\Users\\john\\Desktop\\in.txt","r",stdin);
    // freopen("C:\\Users\\john\\Desktop\\out.txt","w",stdout);
#endif
    while(SII(n,m)==2) {
        rep(i,1,n) SI(in[i]);
        rep(i,1,n) SI(out[i]);
        int source=0,sink=n+n+1;
        dinic.init(source,sink);
        rep(i,1,n) {
            dinic.add_edge(i+n,sink,in[i]);
            dinic.add_edge(i,source,0);
        }
        rep(i,1,n) {
            dinic.add_edge(source,i,out[i]);
            dinic.add_edge(i,source,0);
        }
        rep(i,1,m) {
            int u,v;
            SII(u,v);
            dinic.add_edge(u,v+n,INF_INT);
            dinic.add_edge(v+n,u,0);
        }
        printf("%d\n",dinic.solve(source,sink));
        CLR(vis);
        dinic.dfs(source);
        int ans=0;
        rep(i,1,n) ans+=(!vis[i]+vis[i+n]);
        printf("%d\n",ans);
        rep(i,1,n) {
            if (!vis[i]) printf("%d %c\n",i,'-');
            if (vis[i+n]) printf("%d %c\n",i,'+');
        }
    }
    return 0;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: