【流量平衡】hdu3947

虽然是志愿者招募的变形，但是也有一定的新意

给出一棵内向树，每条边有一定的权值，目标是所有的边权小于等于0，给出若干操作ui,vi,li,ci，表示将ui到vi的路径上的边减1的权值需要ci的费用，该操作共可执行li次。

一开始被叉姐的讲解误导，用一种类似点事件的思想，每种操作看做在某个点+1的标记和某个点-1的标记，这个思想或许在志愿者招募适用，但在树上有分叉，就不太可行了（或许是我理解错了叉姐的意思）

言归正传，先对每条父亲边放缩，建出xi+xj+...+xk=wi+yi的等式，然后减去所有儿子边的等式，对于根节点就用0==0去减，就可以得出一系列每个变量等式左右各只出现一次的等式，然后用流量平衡建图即可。

关键是减去所有儿子这一步，我想了半天，如果只是父亲减儿子的话，每个变量出现的次数就不定了。

建图的时候，每个变量分开考虑建起来也相当简单。

超时一次是距离数组开大了，每次memset时间相当可观，差点就改写原始对偶了...

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
const int oo=1073741819;
using namespace std;
int next[2000000],sora[2000000],flow[2000000],cost[2000000],po[2000000],tail[2000000];
int p[2000000],st[2000000],d[250],v[2000000];
int s,t,n,ss,m,ans,sum,lim,T;
bool spfa(int s,int t)
{
    int h,r,ne,na;
    memset(d,127,sizeof(d));
    h=r=0;
    st[r=1]=s,d[s]=0,v[s]=1,p[s]=0;
    for (;h<r;) {
        ne=st[++h];
        for (int i=ne;next[i];) {
            i=next[i],na=sora[i];
            if (flow[i] && (d[ne]+cost[i]<d[na])) {
                d[na]=d[ne]+cost[i],p[na]=i;
                if (!v[na]) {
                    v[na]=1,st[++r]=na;
                }
            }
        }
        v[ne]=0;
    }
    return d[t]<oo;
}
void widen()
{
    int k=oo;
    for (int i=t;i!=s;i=sora[po[p[i]]]) {
        k=min(k,flow[p[i]]);
    }
    sum+=k;
    for (int i=t;i!=s;i=sora[po[p[i]]]) {
        ans+=k*cost[p[i]];
        flow[p[i]]-=k,flow[po[p[i]]]+=k;
    }
}
void origin()
{
    s=0,t=n+1,ss=t,lim=0;
    for (int i=s;i<=t;i++) tail[i]=i,next[i]=0;
}
void link(int x,int y,int z,int c)
{
    ++ss,next[tail[x]]=ss,tail[x]=ss,sora[ss]=y,flow[ss]=z,cost[ss]=c,next[ss]=0;
    ++ss,next[tail[y]]=ss,tail[y]=ss,sora[ss]=x,flow[ss]=0,cost[ss]=-c,next[ss]=0;
    po[ss]=ss-1,po[ss-1]=ss;
    if (x==s) lim+=z;
}
int main()
{
	freopen("input.txt","r",stdin);
	freopen("output.txt","w",stdout);
    cin>>T;
    for (int test=1;T;T--,test++) {
        printf("Case #%d: ",test);
        scanf("%d\n",&n);
        origin();
        for (int i=1;i<=n-1;i++) {
            int x,y,z;
            scanf("%d%d%d\n",&x,&y,&z);
            link(x,y,oo,0);
            link(x,t,z,0);
            link(s,y,z,0);
        }
        scanf("%d\n",&m);
        for (int i=1;i<=m;i++) {
            int x,y,z,c;
            scanf("%d%d%d%d\n",&x,&y,&z,&c);
            link(y,x,z,c);
        }
        sum=0,ans=0;
        for (;spfa(s,t);widen()) ;
//        cout<<sum<<' '<<ans<<endl;
        if (sum==lim) {
            printf("%d\n",ans);
        }
        else {
            printf("-1\n");
        }
    }
	return 0;
}