HDU6598 Harmonious Army

http://acm.hdu.edu.cn/showproblem.php?pid=6598
题意：n个士兵，有两种职业。有m对士兵待安排职业，对于每一对的两人，若都是魔法师，增加a力量，都是战士，增加c力量，一个魔法师一个战士，增加b力量。现在要你安排这n个人的职业，使得总力量最大。
题解：用最小割建模。对于每一对的两个人，s与两人相连，容量为a/2，两人与t相连，容量为c/2，两人互相连两条边，容量为a/2+c/2-b。设最小割是f，答案就是$({\sum }_{i=1}^{m}a+c)-f$。
( 这道题与bzoj3511类似。)

#include<bits/stdc++.h>
using namespace std;
const int maxn=600;
typedef long long ll;
const ll INF=(1LL<<50);

struct Dinic{
    struct Edge{
        int from,to;
        ll cap,flow;
    };
    int n,m0,m,s,t;
    vector<Edge> edges;
    vector<int> G[maxn];
    bool vis[maxn];
    ll d[maxn];
    int cur[maxn];
    ll sum;

    void init()
    {
        int u,v,x,y,z;
        edges.clear();
        for(int i=0;i<=n+1;i++)G[i].clear();
        s=0,t=n+1;
        sum=0;
        for(int i=1;i<=m0;i++)
        {
            cin>>u>>v>>x>>y>>z;
            sum+=x+z;
            AddEdge(s,u,x);
            AddEdge(s,v,x);
            AddEdge(u,t,z);
            AddEdge(v,t,z);
            AddEdge(u,v,x+z-2*y);
            AddEdge(v,u,x+z-2*y);
        }
    }

    void AddEdge(int f,int t,ll c)
    {
        edges.push_back((Edge){f,t,c,0});
        edges.push_back((Edge){t,f,0,0});
        m=edges.size();
        G[f].push_back(m-2);
        G[t].push_back(m-1);
    }

    bool bfs()
    {    
        memset(vis,0,sizeof(vis));
        queue<int> Q;
        Q.push(s);
        d[s]=0;
        vis[s]=1;
        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];
    }

    ll dfs(int x,ll a)
    {
        if(x==t || a==0)return a;
        ll 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)break;
            }
        }
        return flow;
    }

    ll MaxFlow()
    {
        ll flow=0;
        while(bfs())
        {
            memset(cur,0,sizeof(cur));
            flow+=dfs(s,INF);
        }
        return flow;
    }
}ans;

int main()
{
    //freopen("input.in","r",stdin);
    while(cin>>ans.n>>ans.m0)
    {
        ans.init(); 
        ll flow=ans.MaxFlow();
        cout<<ans.sum-flow/2<<"\n";
    }
    return 0;
}