【題目描述】：

一共有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.