[NOI 2008]志愿者招募（费用流）

【题目描述】：

一共有n天，每天需要a[i]个人，一共有m类志愿者，每类有无限个，可以从s[i]服务到t[i]，需要花费c[i]的费用。问要满足要求，最小的费用。

【题目分析】：

这题有一个经典的也很奇妙的构图法：利用不等式构图~这个我不懂，想看的话请到byvoid神牛的blog里面搜~~

我的想法很朴素，流的话就表示人。用容量来限制需要的人数。这样我有n+3个点，一个源点（0号），一个汇点（n+2号）。然后其他的点编号为1~n+1。从源点到1~n连一条容量为a[i]，费用为0的边。再从2~n+1向汇点连一条容量为a[i-1]，费用为0的边。从2~n+1向前一个点连一条容量正无穷，费用为0的边。对于每类志愿者(s,t,c)，连一条从s[i]到t[i]+1容量为正无穷，费用为c[i]的边。这样我们再来看看这个图，对于每天的要求必须通过从前往后走的志愿者的边和回向边才能到达当天对应的流的出口，因为一定有解，所以这个网络流的最小费用最大流的费用就是结果~~

这题的模型我是参照的经典的餐巾问题想出来的~据说当时在考场上难倒了一大批神牛，但是确实很变态~费用流貌似用SPFA的那个会超时，所以要用zkw~~zkw的话参见之前的一片BLOG

【代码】：

program employee; const inf=1000000; type ntype=record data,w,c,op,next:longint; end; var a,g,di,dis:array[0..2000] of longint; v:array[0..2000] of boolean; node:array[1..50000] of ntype; ans,tot,n,m,i,s,t,x,y,w:longint; function min(x,y:longint):longint; begin if x<y then exit(x) else exit(y); end; procedure link(x,y,c,w:longint); var p,t:longint; begin inc(tot); p:=tot; node[p].data:=y; node[p].c:=c; node[p].w:=w; node[p].next:=g[x]; g[x]:=p; inc(tot); t:=tot; node[t].data:=x; node[t].c:=0; node[t].w:=-w; node[t].next:=g[y]; g[y]:=t; node[p].op:=t; node[t].op:=p; end; function aug(x,flow:longint):longint; var p,tmp:longint; begin if x=t then begin inc(ans,flow*dis[s]); exit(flow); end; v[x]:=true; p:=di[x]; while p<>0 do begin if (not v[node[p].data]) and (node[p].c>0) and (dis[node[p].data]+node[p].w=dis[x]) then begin tmp:=aug(node[p].data,min(node[p].c,flow)); if tmp>0 then begin di[x]:=p; dec(node[p].c,tmp); inc(node[node[p].op].c,tmp); exit(tmp); end; end; p:=node[p].next; end; di[x]:=0; exit(0); end; function modlabel:boolean; var i,p:longint; delta:longint; begin delta:=inf; for i:=s to t do if v[i] then begin p:=g[i]; while p<>0 do begin if not v[node[p].data] and (node[p].c>0) then delta:=min(delta,dis[node[p].data]+node[p].w-dis[i]); p:=node[p].next; end; end; if delta=inf then exit(true); for i:=s to t do if v[i] then inc(dis[i],delta); fillchar(v,sizeof(v),0); exit(false); end; begin assign(input,'employee.in'); assign(output,'employee.out'); reset(input); rewrite(output); readln(n,m); for i:=1 to n do begin read(a[i]); link(0,i,a[i],0); end; for i:=2 to n+1 do begin link(i,i-1,inf,0); link(i,n+2,a[i-1],0); end; for i:=1 to m do begin readln(x,y,w); link(x,y+1,inf,w); end; s:=0; t:=n+2; repeat di:=g; while aug(s,inf)>0 do fillchar(v,sizeof(v),0); until modlabel; writeln(ans); close(input); close(output); end.