pku3680 网络流

额，最近几天没怎么写程序，关键数学初赛考萎了，以后要多做数学了。

这道题我叫了N遍，最后发现是由于他多组数据，l没有初始化。。。。

题目大意：

  一个数轴上，有n条线段，每条的线段有a[i],b[i],w[i]表示起点终点和权值。让你选若干条线段让这些线段的和最大。数轴的某一段区间可以被重复选k次。

简略题解：

  首先膜拜cai0715,看了他的构图方法才做出来的。

  先离散化，就是先把所有线段的端点排序，然后再把原来的a[i],b[i]缩小成p,q，然后，增加源汇s,t，s=0 t=r+1(总区间长度+1)。然后对于数轴，i向i+1连边，费用0流量k。对于每条线段，从p向q连边，费用为-c[i]，流量为1。

  然后就是最小费用最大流了，最后ans取反就可以了。

附程序：

program pku_3680; var pre,other,w,c:array[0..100000] of longint; zt:array[0..500000] of longint; last,d,nb,pro,pra:array[0..402] of longint; a,b,e:array[0..500] of longint; v:array[0..402] of boolean; n,k,l,testcase,i,r,st,ed,ans:longint; {-------------------------} procedure grow(x,y,z,t:longint); begin inc(l); pre[l]:=last[x];last[x]:=l;other[l]:=y;w[l]:=z;c[l]:=t; inc(l); pre[l]:=last[y];last[y]:=l;other[l]:=x;w[l]:=-z;c[l]:=0; end; {-------------------------} procedure qsort(x,y:longint); var g,t,l,r:longint; begin t:=nb[(x+y) div 2]; l:=x;r:=y; repeat while nb[l]<t do inc(l); while nb[r]>t do dec(r); if l<=r then begin g:=nb[l];nb[l]:=nb[r];nb[r]:=g; inc(l);dec(r); end; until l>r; if l<y then qsort(l,y); if x<r then qsort(x,r); end; {-------------------------} function find(x:longint):longint; var s,t,mid:longint; begin s:=1;t:=r; repeat mid:=(s+t) div 2; if nb[mid]>x then t:=mid-1 else s:=mid+1; until s>t; exit(t); end; {-------------------------} procedure init; var i,len:longint; begin readln(n,k); r:=0; for i:=1 to n do begin readln(a[i],b[i],e[i]); inc(r);nb[r]:=a[i]; inc(r);nb[r]:=b[i]; end; qsort(1,r); len:=1; for i:=1 to r do if nb[i]>nb[len] then begin inc(len);nb[len]:=nb[i]; end; r:=len; for i:=1 to n do begin a[i]:=find(a[i]); b[i]:=find(b[i]); end; end; {-------------------------} procedure spfa; var l,r,p,q,i:longint; begin l:=0;r:=1; for i:=1 to ed do d[i]:=maxlongint; d[st]:=0; fillchar(v,sizeof(v),false); v[st]:=true; zt[1]:=st; while l<r do begin inc(l); p:=zt[l];q:=last[p];v[p]:=false; while q<>0 do begin if c[q]>0 then if d[other[q]]>d[p]+w[q] then begin d[other[q]]:=d[p]+w[q]; pro[other[q]]:=p; pra[other[q]]:=q; if not v[other[q]] then begin inc(r);zt[r]:=other[q]; v[other[q]]:=true; end; end; q:=pre[q]; end; end; end; {-------------------------} procedure change; var p,q,tmp:longint; begin p:=ed; tmp:=maxlongint; while p<>st do begin q:=pra[p]; if c[q]<tmp then tmp:=c[q]; p:=pro[p]; end; p:=ed; while p<>st do begin q:=pra[p]; dec(c[q],tmp); if odd(q) then inc(c[q+1],tmp) else inc(c[q-1],tmp); p:=pro[p]; end; dec(ans,d[ed]*tmp); end; {-------------------------} procedure main; var i:longint; begin ans:=0; l:=0; st:=0;ed:=r+1; fillchar(last,sizeof(last),0); for i:=0 to r do grow(i,i+1,0,k); for i:=1 to n do grow(a[i],b[i],-e[i],1); d[ed]:=0; while true do begin spfa; if d[ed]<maxlongint then change else break; end; writeln(ans); end; {-------------------------} begin readln(testcase); for i:=1 to testcase do begin init; main; end; end.