【HDOJ6118】度度熊的交易计划（费用流）

题意：

度度熊参与了喵哈哈村的商业大会，但是这次商业大会遇到了一个难题：

喵哈哈村以及周围的村庄可以看做是一共由n个片区，m条公路组成的地区。

由于生产能力的区别，第i个片区能够花费a[i]元生产1个商品，但是最多生产b[i]个。

同样的，由于每个片区的购买能力的区别，第i个片区也能够以c[i]的价格出售最多d[i]个物品。

由于这些因素，度度熊觉得只有合理的调动物品，才能获得最大的利益。

据测算，每一个商品运输1公里，将会花费1元。

那么喵哈哈村最多能够实现多少盈利呢？

1<=n<=500,
1<=m<=1000,
1<=a[i],b[i],c[i],d[i],k[i]<=1000,
1<=u[i],v[i]<=n

思路：明显的一个最小费用流模型 写的是最大费用流，注意没有要求流量最大

刚开始TooNaive先跑了一波最短路然后建的图

而后发现直接照原图建边，费用流过程中就会自动调整到最短路上

还是要学习一个

SPFA找的时候如果当前继续跑流会使利润减小就停止，即dis[src]<0

 const inf=;
 var q:array[..]of longint;
     f:array[..,..]of longint;
     a,b,c,d,dis:array[..]of longint;
     head,vet,len1,len2,next:array[..]of longint;
     pre:array[..,..]of longint;
     inq:array[..]of boolean;
     fan:array[..]of longint;
     n,m,i,tot,t1,s,source,src,x,y,z,j:longint;
     ans1,ans2:int64;

 procedure add(a,b,c,d:longint);
 begin
  inc(tot);
  next[tot]:=head[a];
  vet[tot]:=b;
  len1[tot]:=c;
  len2[tot]:=d;
  head[a]:=tot;
 end;

 function min(x,y:int64):int64;
 begin
  if x<y then min:=x
   else min:=y;
 end;

 function spfa:boolean;
 var u,e,v,t,w,i:longint;
 begin
  for i:= to s do
  begin
   dis[i]:=-(maxlongint>>);
   inq[i]:=false;
  end;
  t:=; w:=; q[]:=source; dis[source]:=; inq[source]:=true;
  while t<w do
  begin
   inc(t); u:=q[t mod ]; inq[u]:=false;
   e:=head[u];
   while e<> do
   begin
    v:=vet[e];
    if (len1[e]>)and(dis[u]+len2[e]>dis[v]) then
    begin
     dis[v]:=dis[u]+len2[e];
     pre[v,]:=u;
     pre[v,]:=e;
     if not inq[v] then
     begin
      inc(w); q[w mod ]:=v; inq[v]:=true;
     end;
    end;
    e:=next[e];
   end;
  end;
  //if dis[src]=-(maxlongint>>) then spfa:=false
  // else spfa:=true;
  if dis[src]< then spfa:=false
   else spfa:=true;
 end;

 procedure mcf;
 var k,e:longint;
     t:int64;
 begin
  k:=src; t:=<<;
  while k<>source do
  begin
   t:=min(t,len1[pre[k,]]);
   k:=pre[k,];
  end;
  k:=src;
  while k<>source do
  begin
   e:=pre[k,];
   len1[e]:=len1[e]-t;
   len1[fan[e]]:=len1[fan[e]]+t;
   ans2:=ans2+t*len2[e];
   k:=pre[k,];
  end;
 end;

 begin
  assign(input,'1005.in'); reset(input);
  assign(output,'1005.out'); rewrite(output);
  for i:= to  do
   if i and = then fan[i]:=i+
    else fan[i]:=i-;
  while not eof do
  begin
   read(n,m);
   tot:=;
   if (n=)and(m=) then break;
   for i:= to n do read(a[i],b[i],c[i],d[i]);
   tot:=;
   for i:= to n do head[i]:=;
   source:=n+; src:=n+; s:=n+;
   for i:= to n do
   begin
    add(source,i,b[i],-a[i]);
    add(i,source,,a[i]);
   end;
   for i:= to n do
   begin
    add(i,src,d[i],c[i]);
    add(src,i,,-c[i]);
   end;

   for i:= to m do
   begin
    read(x,y,z);
    add(x,y,maxlongint,-z);
    add(y,x,,z);
    add(y,x,maxlongint,-z);
    add(x,y,,z);
   end;
   ans1:=; ans2:=;
   while spfa do mcf;
   writeln(ans2);
   for i:= to s do head[i]:=;
   for i:= to s do dis[i]:=;
   for i:= to s do
    for j:= to  do pre[i,j]:=;
  end;
  close(input);
  close(output);

 end.