ZOJ 2676 Network Wars

ZOJ_2676

    这个题目可以像最优比率生成树那样用0-1分数规划去做，只不过最优比率生成树每次是求一棵生成树，而这个题目要求一个最小割。

    对于每次二分，在建图的时候可能会出现边权为负的边，由于通过分析后可知这些边一定会选，所以就没必要再将其加入到新建的图中了。

    最后要输出最小割集，求最小割集的时候可以用DFS遍历做完最大流之后的图，且只能沿没有满流的边向下走。遍历完成后，如果某条边有一个端点遍历到了而另外一个端点没有遍历到，那么就说明这条边是属于最小割集中的。

#include<stdio.h>
#include<string.h>
#include<algorithm>
#define zero 1e-8
#define MAXD 110
#define MAXM 1610
#define INF 0x3f3f3f3f
using namespace std;
int N, M, S, T, first[MAXD], work[MAXD], e, next[MAXM], v[MAXM], id[MAXM], use[MAXM], vis[MAXD], q[MAXD], d[MAXD];
double flow[MAXM];
struct Graph
{
    int u[MAXM], v[MAXM], w[MAXM], id[MAXM], e;
    void init()
    {
        e = 0;
    }
    void add(int x, int y, int z, int flag)
    {
        u[e] = x, v[e] = y, w[e] = z, id[e] = flag;
        ++ e;
    }
}g;
void init()
{
    int i, x, y, z;
    g.init();
    for(i = 1; i <= M; i ++)
    {
        scanf("%d%d%d", &x, &y, &z);
        g.add(x, y, z, i);
    }
}
void add(int x, int y, double f, int flag)
{
    v[e] = y, flow[e] = f, id[e] = flag;
    next[e] = first[x], first[x] = e ++;
}
double initgraph(double r)
{
    int i;
    double ans = 0;
    memset(use, 0, sizeof(use));
    memset(first, -1, sizeof(first));
    e = 0;
    for(i = 0; i < g.e; i ++)
    {
        if(g.w[i] - r <= zero)
            use[g.id[i]] = 1, ans += g.w[i] - r;
        else
        {
            add(g.u[i], g.v[i], g.w[i] - r, g.id[i]), add(g.v[i], g.u[i], 0, g.id[i]);
            add(g.v[i], g.u[i], g.w[i] - r, g.id[i]), add(g.u[i], g.v[i], 0, g.id[i]);
        }
    }
    return ans;
}
int bfs()
{
    int i, j, rear = 0;
    memset(d, -1, sizeof(d));
    d[S] = 0, q[rear ++] = S;
    for(i = 0; i < rear; i ++)
        for(j = first[q[i]]; j != -1; j = next[j])
            if(flow[j] > zero && d[v[j]] == -1)
            {
                d[v[j]] = d[q[i]] + 1;
                if(v[j] == T)
                    return 1;
                q[rear ++] = v[j];
            }
    return 0;
}
double dfs(int cur, double a)
{
    if(cur == T)
        return a;
    double t;
    for(int &i = work[cur]; i != -1; i = next[i])
        if(flow[i] > zero && d[v[i]] == d[cur] + 1)
            if((t = dfs(v[i], flow[i] < a ? flow[i] : a)) > zero)
            {
                flow[i] -= t;
                flow[i ^ 1] += t;
                return t;
            }
    return 0;
}
double dinic()
{
    double ans = 0, t;
    while(bfs())
    {
        memcpy(work, first, sizeof(first));
        while((t = dfs(S, INF)) > zero)
            ans += t;
    }
    return ans;
}
void DFS(int cur)
{
    int i;
    vis[cur] = 1;
    for(i = first[cur]; i != -1; i = next[i])
        if(flow[i] > zero && !vis[v[i]])
            DFS(v[i]);
}
void print()
{
    int i, j, n = 0, flag = 0;
    memset(vis, 0, sizeof(vis));
    DFS(S);
    for(i = 1; i < N; i ++)
        for(j = first[i]; j != -1; j = next[j])
            if((vis[i] & !vis[v[j]]) || (!vis[i] && vis[v[j]]))
                use[id[j]] = 1;
    for(i = 1; i <= M; i ++)
        if(use[i])
            ++ n;
    printf("%d\n", n);
    for(i = 1; i <= M; i ++)
        if(use[i])
        {
            flag ? printf(" ") : flag = 1;
            printf("%d", i);
        }
    printf("\n");
}
void solve()
{
    int i;
    double ans, min, max, mid;
    max = 10000010, min = 0;
    S = 1, T = N;
    for(i = 0; i < 50; i ++)
    {
        mid = (min + max) / 2;
        ans = initgraph(mid);
        ans += dinic();
        if(ans > 0)
            min = mid;
        else
            max = mid;
    }
    print();
}
int main()
{
    while(scanf("%d%d", &N, &M) == 2)
    {
        init();
        solve();
    }
    return 0;
}