BZOJ 1877：[SDOI2009]晨跑（最小费用最大流）

                                                                    晨跑
Description
Elaxia最近迷恋上了空手道，他为自己设定了一套健身计划，比如俯卧撑、仰卧起坐等 等，不过到目前为止，他坚持下来的只有晨跑。 现在给出一张学校附近的地图，这张地图中包含N个十字路口和M条街道，Elaxia只能从 一个十字路口跑向另外一个十字路口，街道之间只在十字路口处相交。Elaxia每天从寝室出发 跑到学校，保证寝室编号为1，学校编号为N。 Elaxia的晨跑计划是按周期（包含若干天）进行的，由于他不喜欢走重复的路线，所以 在一个周期内，每天的晨跑路线都不会相交（在十字路口处），寝室和学校不算十字路 口。Elaxia耐力不太好，他希望在一个周期内跑的路程尽量短，但是又希望训练周期包含的天 数尽量长。 除了练空手道，Elaxia其他时间都花在了学习和找MM上面，所有他想请你帮忙为他设计 一套满足他要求的晨跑计划。
Input
第一行：两个数N,M。表示十字路口数和街道数。 接下来M行，每行3个数a,b,c，表示路口a和路口b之间有条长度为c的街道（单向）。
Output
两个数，第一个数为最长周期的天数，第二个数为满足最长天数的条件下最短的路程长 度。
Sample Input
7 10
1 2 1
1 3 1
2 4 1
3 4 1
4 5 1
4 6 1
2 5 5
3 6 6
5 7 1
6 7 1
Sample Output
2 11
HINT
对于30%的数据，N ≤ 20，M ≤ 120。
对于100%的数据，N ≤ 200，M ≤ 20000

分析:

  拆点费用流，每个点拆成两个并连一条边，容量为1（起始点终点为m），费用为0，其它各边容量为1，费用为其长度，跑最小费用最大流即可。

代码：

program bzoj1877;
const
  nn=1000;
type
  point=record
     next,f,t,cap,v:longint;
  end;
var
  e:array[0..400001]of point;
  head,p,d:array[0..401]of longint;
  q:array[0..nn+1]of longint;
  g:array[0..401]of boolean;
  n,i,m,x,y,v,tot,s,t:longint;
function min(x,y:longint):longint;
begin
  if x<y then min:=x else min:=y;
end;
procedure add(x,y,c,v:longint);
begin
  inc(tot);e[tot].next:=head[x];head[x]:=tot; e[tot].f:=x;e[tot].t:=y;e[tot].cap:=c; e[tot].v:=v;
  inc(tot);e[tot].next:=head[y];head[y]:=tot; e[tot].f:=y;e[tot].t:=x; e[tot].cap:=0; e[tot].v:=-v;
end;
function spfa(s,t:longint):boolean;
var i,x,y,h,tail:longint;
begin
  for i:=1 to n do begin d[i]:=maxlongint div 3; g[i]:=false; end;
  g[s]:=true; d[s]:=0; h:=0; tail:=1;  q[tail]:=s;
  while h<>tail do
   begin
     inc(h);x:=q[h];g[x]:=false; i:=head[x]; if h=nn then h:=0;
     while i<>0 do
      begin
        y:=e[i].t;
        if (e[i].cap>0)and(d[y]>d[x]+e[i].v) then
         begin
           d[y]:=d[x]+e[i].v; p[y]:=i;
           if not g[y] then
             begin
               g[y]:=true; inc(tail); q[tail]:=y; if tail=nn then tail:=0;
             end;
          end;
        i:=e[i].next;
      end;
   end;
   exit(d[t]<maxlongint div 3);
end;
procedure mcf(s,t:longint);
var ans,f,x,flow:longint;
begin
  ans:=0; f:=0; flow:=0; x:=0;
  while spfa(s,t) do
   begin
     f:=maxlongint; x:=t;
     while x<>s do begin f:=min(e[p[x]].cap,f);x:=e[p[x]].f; end;
     x:=t;
     while x<>s do
      begin dec(e[p[x]].cap,f); inc(e[p[x]-ord(p[x] mod 2=0)+ord(p[x] mod 2>0)].cap,f);x:=e[p[x]].f; end;
     inc(ans,d[t]*f); inc(flow,f);
   end;
  writeln(flow,' ',ans);
end;
begin 
  readln(n,m);
  fillchar(head,sizeof(head),0);
  s:=1; t:=n;     add(s,s+n,m,0); add(t,t+n,m,0);
  for i:=2 to n-1 do add(i,i+n,1,0);
  for i:=1 to m do
   begin
    readln(x,y,v);
    add(x+n,y,1,v);
   end;
  n:=n*2;
  mcf(s,t+n div 2);
end.