zoukankan      html  css  js  c++  java
  • hdu 3452 Bonsai 最小割

    #include <cstdio>
    #include <vector>
    #include <queue>
    #include <cstring>
    using namespace std;
    #define maxn 1500
    #define INF 100000
    struct Edge
    {
        int from, to, cap, flow;
    };
    
    int n, m, s, t;
    vector<Edge> edges;    // 边数的两倍
    vector<int> G[maxn];   // 邻接表,G[i][j]表示结点i的第j条边在e数组中的序号
    bool vis[maxn];        // BFS使用
    int d[maxn];           // 从起点到i的距离
    int cur[maxn];         // 当前弧指针
    int deg[maxn];
    int min(int a,int b)
    {
        if(a<b) return a;
        else return b;
    }
    void AddEdge(int from, int to, int cap)
    {
        int len;
        Edge temp;
        temp.from=from;temp.to=to;temp.cap=cap;temp.flow=0;
        edges.push_back(temp);
        temp.from=to;temp.to=from;temp.cap=0;temp.flow=0;
        edges.push_back(temp);
        len = edges.size();
        G[from].push_back(len-2);
        G[to].push_back(len-1);
    }
    bool BFS()
    {
        memset(vis, 0, sizeof(vis));
        queue<int> Q;
        Q.push(s);
        vis[s] = 1;
        d[s] = 0;
        while(!Q.empty())
        {
            int x = Q.front(); Q.pop();
            for(int i = 0; i < G[x].size(); i++)
            {
                Edge& e = edges[G[x][i]];
                if(!vis[e.to] && e.cap > e.flow)
                {
                    vis[e.to] = 1;
                    d[e.to] = d[x] + 1;
                    Q.push(e.to);
                }
            }
        }
        return vis[t];
    }
    int DFS(int x, int a)
    {
        if(x == t || a == 0) return a;
        int flow = 0, f;
        for(int& i = cur[x]; i < G[x].size(); i++)
        {
            Edge& e = edges[G[x][i]];
            if(d[x] + 1 == d[e.to] && (f = DFS(e.to, min(a, e.cap-e.flow))) > 0)
            {
                e.flow += f;
                edges[G[x][i]^1].flow -= f;
                flow += f;
                a -= f;
                if(a == 0) break;
            }
        }
        return flow;
    }
    int Maxflow()
    {
        int flow = 0;
        while(BFS())
        {
            memset(cur, 0, sizeof(cur));
            flow += DFS(s, INF);
        }
        return flow;
    }
    int main()
    {
        int i,j,x,y,z;
        while(1)
        {
            scanf("%d%d",&n,&s);
            if(n==0&&s==0) break;
            memset(deg,0,sizeof(deg));
            t=n+1;
            for(i=1;i<=n+10;i++)
                G[i].clear();
            for(i=1;i<=n-1;i++)
            {
                scanf("%d%d%d",&x,&y,&z);
                deg[x]++;
                deg[y]++;
                AddEdge(x,y,z);
                AddEdge(y,x,z);
            }
            for(i=1;i<=n;i++)
            {
                if(i==s) continue;
                if(deg[i]==1)
                    AddEdge(i,t,INF);
            }
            printf("%d
    ",Maxflow());
        }
        return 0;
    }
    


     

  • 相关阅读:
    UVA 10480 Sabotage (最大流最小割)
    bzoj2002 [Hnoi2010]Bounce 弹飞绵羊 (分块)
    poj3580 SuperMemo (Splay+区间内向一个方向移动)
    bzoj1500: [NOI2005]维修数列 (Splay+变态题)
    hdu3436 Queue-jumpers(Splay)
    hdu4710 Balls Rearrangement(数学公式+取模)
    hdu1890 Robotic Sort (splay+区间翻转单点更新)
    zoj2112 Dynamic Rankings (主席树 || 树套树)
    poj3581 Sequence (后缀数组)
    notepa++ Emmet的安装方法
  • 原文地址:https://www.cnblogs.com/vermouth/p/3710148.html
Copyright © 2011-2022 走看看