HDU3996 Gold Mine 最大闭合权

　　题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3996

　　很无语加蛋碎的题目！题意描述不清，然后还有超int（这里wa，果断被坑）= =！方法很简单，很容易看出是最大闭合权。

//STATUS:G++_AC_109MS_2068KB
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<iostream>
#include<string>
#include<algorithm>
#include<vector>
#include<queue>
#include<stack>
#include<map>
using namespace std;
#define LL __int64
#define pii pair<int,int>
#define Max(a,b) ((a)>(b)?(a):(b))
#define Min(a,b) ((a)<(b)?(a):(b))
#define mem(a,b) memset(a,b,sizeof(a))
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
const int MAX=2510,INF=0x3f3f3f3f;
const LL LLNF=0x3f3f3f3f3f3f3f3fLL;

struct Edge{
    int u,v;
    LL cap;
};

vector<Edge> e;
vector<int> g[MAX];
int d[MAX],cur[MAX],p[MAX],num[MAX],val[MAX],ma[110][25],li[MAX][55][2],cou[MAX];
int n,m,k,T,s,t,mt;

void adde(int a,int b,LL val){
    e.push_back((Edge){a,b,val});
    g[a].push_back(mt++);
    e.push_back((Edge){b,a,0});
    g[b].push_back(mt++);
}

LL augment()
{
    int x=t;
    LL a=LLNF;
    while(x!=s){
        a=Min(a,e[p[x]].cap);
        x=e[p[x]].u;
    }
    x=t;
    while(x!=s){
        e[p[x]].cap-=a;
        e[p[x]^1].cap+=a;
        x=e[p[x]].u;
    }
    return a;
}

LL isap()
{
    int i,x,ok,min;
    LL flow=0;
    mem(d,0);mem(num,0);
    num[0]=t+1;
    mem(cur,0);
    x=s;
    while(d[s]<=t){
        if(x==t){
            flow+=augment();
            x=s;
        }
        ok=0;
        for(i=cur[x];i<g[x].size();i++){
            Edge& et=e[g[x][i]];
            if(et.cap && d[x]==d[et.v]+1){
                ok=1;
                p[et.v]=g[x][i];
                cur[x]=i;
                x=et.v;
                break;
            }
        }
        if(!ok){
            min=t;
            for(i=0;i<g[x].size();i++){
                Edge& et=e[g[x][i]];
                if(et.cap && d[et.v]<min)min=d[et.v];
            }
            if(--num[d[x]]==0)break;
            num[d[x]=min+1]++;
            cur[x]=0;
            if(x!=s)x=e[p[x]].u;
        }
    }
    return flow;
}

void bfs()
{
    int x,i,j;
    mem(d,0);
    queue<int> q;
    q.push(s);
    d[s]=1;
    while(!q.empty()){
        x=q.front();q.pop();
        for(i=0;i<g[x].size();i++){
            Edge& et=e[g[x][i]];
            if(et.cap && !d[et.v]){
                d[et.v]=1;
                q.push(et.v);
            }
        }
    }
}

int main()
{
  //  freopen("in.txt","r",stdin);
    int i,j,q,x,y,ca=0;
    LL ans;
    scanf("%d",&T);
    while(T--)
    {
        k=mt=ans=0;
        e.clear();
        scanf("%d",&n);
        for(i=1;i<=n;i++){
            scanf("%d",&m);
            for(j=1;j<=m;j++){
                ++k;
                scanf("%d%d%d",&x,&y,&cou[k]);
                val[k]=y-x;
                ma[i][j]=k;
                for(q=1;q<=cou[k];q++)
                    scanf("%d%d",&li[k][q][0],&li[k][q][1]);
            }
        }
        s=0,t=k+1;
        for(i=0;i<=t;i++)g[i].clear();
        for(i=1;i<=k;i++){
            if(val[i]>0)adde(s,i,val[i]);
            else if(val[i]<0)adde(i,t,-val[i]);
            for(j=1;j<=cou[i];j++)
                adde(i,ma[li[i][j][0]][li[i][j][1]],LLNF);
        }

        isap();
        bfs();
        for(i=1;i<=k;i++)
            if(d[i]>0)ans+=val[i];
        printf("Case #%d: %I64d\n",++ca,ans);
    }
    return 0;
}