虽然是志愿者招募的变形,但是也有一定的新意
给出一棵内向树,每条边有一定的权值,目标是所有的边权小于等于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;
}