2015年ACM/ICPC北京赛区 D题（最大流）

题目链接：https://icpcarchive.ecs.baylor.edu/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=5276

参考网上题解：

题意：给定一颗类似于dnf、lol点天赋时候的天赋（技能）有向树，技能获得的前提是他的前置技能都获得了，作为一个RMB玩家，你有特权：
1.直接花费一定数量的钱获得某个技能。
2.花费一定数量的钱将一个技能的某一个前置关系取消，即将前置技能到该技能的边消除（但不会获得该前置技能）。
如果正常学习技能的话每一个技能都要花费一定量的时间，问获得指定的技能的最少的花费是多少。

思路：读完题就会知道是图论相关，第一反应是最短路什么的，但是这些复杂的优先级关系并非最短路能处理的。后来看题解发现是拆点+最小割，其实就是最大流模板题，随便套个网络流算法的模板就能过，但是此题重点和难点在于如何建图。
假设要获得的技能为S。
1.将所有点拆成i和i'两个点，在建好的图中就用i和i+n分别表示。
2.若i->j有边（i为j的前置技能），则将i'到j建边，权值为用钱将该边消除的花费。
3.将源点和i建边，边权为正常学习该技能花费的时间（满足前置技能以后再学的花费）。
4.将i与i'建边，权值为用钱直接获得该点的花费（氪金，不用管前置技能）。
5.将S’与汇点建边，权值为inf。
这样最大流跑完了以后，S’与汇点之间正向边的流量就是源点与汇点的最小割，也是要求的最小花费。自己画个简单的图模拟一下很容易就明白为什么这样是对的，但是这种将题中不同条件转化为对应的边的思想是很难一天两天练成的，建图也是图论相关题的奥义所在。多刷题，少看题解，争取早日领略其中奥义，共勉。
 

网络流的题目，还是要尽量少看题解。不过这个输入真的是有问题，看了老半天才发现少了空格。。。迷，最大流的思维还是需要多做题目、仔细思考才能有所领悟。（测试了我的ISAP模板，还不错）

代码：

#include<bits/stdc++.h>
#define ll long long
#define inf 0x3f3f3f3f
#define clear(A,X) memset(A,X,sizeof A)
#define copy(A,B) memcpy(A,B,sizeof A)
using namespace std;
const int maxe=1e6+10;
const int maxq=1e6+10;
const int maxn=1e5+10;

struct Edge{
    int v;
    int c;
    int n;
}edge[maxe];

int adj[maxn],cnte;
int Q[maxq],head,tail;
int d[maxn],cur[maxn],pre[maxn],num[maxn];
int sourse,sink,nv;//sourse：源点，sink：汇点，nv：编号修改的上限
int n,m;
void add(int u,int v,int c)
{
    edge[cnte].v=v;
    edge[cnte].c=c;
    edge[cnte].n=adj[u];
    adj[u]=cnte++;

    edge[cnte].v=u;
    edge[cnte].c=0;
    edge[cnte].n=adj[v];
    adj[v]=cnte++;
}

void rev_bfs()
{
    clear(num,0);
    clear(d,-1);
    d[sink]=0;
    num[0]=1;
    head=tail=0;
    Q[tail++]=sink;
    while(head!=tail)
    {
        int u=Q[head++];
        for(int i=adj[u];~i;i=edge[i].n)
        {
            int v=edge[i].v;
            if(~d[v]) continue;
            d[v]=d[u]+1;
            Q[tail++]=v;
            num[d[v]]++;
        }
    }
}

int ISAP()
{
    copy(cur,adj);
    rev_bfs();
    int flow=0,u=pre[sourse]=sourse,i;
    while(d[sink]<nv)
    {
        if(u==sink)
        {
            int f=inf,neck;
            for(i=sourse;i!=sink;i=edge[cur[i]].v)
            {
                if(f>edge[cur[i]].c)
                {
                    f=edge[cur[i]].c;
                    neck=i;
                }
            }
            for(i=sourse;i!=sink;i=edge[cur[i]].v)
            {
                edge[cur[i]].c-=f;
                edge[cur[i]^1].c+=f;
            }
            flow+=f;
            u=neck;
        }
        for(i=cur[u];~i;i=edge[i].n)
        if(d[edge[i].v]+1==d[u]&&edge[i].c) break;
        if(~i)
        {
            cur[u]=i;
            pre[edge[i].v]=u;
            u=edge[i].v;
        }
        else
        {
            if(0==(--num[d[u]])) break;
            int mind=nv;
            for(i=adj[u];~i;i=edge[i].n)
            {
                if(edge[i].c&&mind>d[edge[i].v])
                {
                    cur[u]=i;
                    mind=d[edge[i].v];
                }
            }
            d[u]=mind+1;
            num[d[u]]++;
            u=pre[u];
        }
    }
    return flow;
}

void init()
{
    clear(adj,-1);
    cnte=0;
}
void work()
{
    int S;
    scanf("%d%d%d",&n,&m,&S);
    int u,v,c;
    init();
    for(int i=0;i<m;++i)
    scanf("%d%d%d",&u,&v,&c),add(u+n,v,c);
    for(int i=1;i<=n;i++)
    scanf("%d",&c),add(0,i,c);
    for(int i=1;i<=n;i++)
    scanf("%d",&c),add(i,i+n,c);
    add(S+n,2*n+1,inf);
    sourse=0;sink=2*n+1;nv=sink+1;
    printf("%d\n",ISAP());
}

int main()
{
    int T,cas=1;
    scanf("%d",&T);
    while(T--) work();
    return 0;
}
/*
2
5 5 5
1 2 5
1 3 5
2 4 8
4 5 10
3 5 15
3 5 7 9 11
100 100 100 200 200
5 5 5
1 2 5
1 3 5
2 4 8
4 5 10
3 5 15
3 5 7 9 11
5 5 5 50 50
*/