每个点单独放到T集合,那么的话值就相当于∑(Bout+Dout)-∑Bin;
也就是说,出边可以看做一定要出去的流量,入边看做可以进入的流量,那么只要出边的流量都能进入入边也就ok了
所以,对于每条边,设置一个下界是b,上界是b+d.直接跑无源汇有上下界最大流就行了。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>
using namespace std;
struct rec
{
int num,w,next,o;
}e[100010];
int linkk[100010];
int du[100010];
int s[100010];
int h[100010];
int n,m,e_tot;
void insert(int x,int y,int z)
{
e_tot++;e[e_tot]={y,z,linkk[x],e_tot+1};linkk[x]=e_tot;
e_tot++;e[e_tot]={x,0,linkk[y],e_tot-1};linkk[y]=e_tot;
}
int dfs(int x,int total)
{
if (x==n+2) return total;
int leave=total;
int minh=n+1;
int tmp;
for (int p=linkk[x];p&&leave;p=e[p].next)
{
if (e[p].w>0)
{
if (h[x]==h[e[p].num]+1)
{
tmp=dfs(e[p].num,min(e[p].w,leave));
leave-=tmp;
e[p].w-=tmp;
e[e[p].o].w+=tmp;
if (h[n+1]==n+2) return total-leave;
}
minh=min(minh,h[e[p].num]);
}
}
if (leave==total)
{
s[h[x]]--;
if (s[h[x]]==0) h[n+1]=n+2;
h[x]=minh+1;
s[h[x]]++;
}
return total-leave;
}
int main()
{
// freopen("input.txt","r",stdin);
int T;
scanf("%d",&T);
for (int tt=T;T;T--)
{
memset(e,0,sizeof(e));
memset(linkk,0,sizeof(linkk));
memset(du,0,sizeof(du));
e_tot=0;
printf("Case #%d: ",tt-T+1);
scanf("%d%d",&n,&m);
int u,v,D,B;
for (int i=1;i<=m;i++)
{
scanf("%d%d%d%d",&u,&v,&D,&B);
insert(u,v,B);
du[u]-=D;
du[v]+=D;
}
for (int i=1;i<=n;i++)
{
if (du[i]>0) insert(n+1,i,du[i]);
if (du[i]<0) insert(i,n+2,-du[i]);
}
int ans=0;
memset(h,0,sizeof(h));
memset(s,0,sizeof(s));
s[0]=n+2;
while (h[n+1]<n+2)
ans+=dfs(n+1,100000000);
for (int i=1;i<=n;i++)
if (du[i]>0) ans-=du[i];
if (ans==0) printf("happy\n");
else printf("unhappy\n");
}
return 0;
}