码迷,mamicode.com
首页 > 其他好文 > 详细

BZOJ 1877:[SDOI2009]晨跑(最小费用最大流)

时间:2016-04-14 19:20:04      阅读:226      评论:0      收藏:0      [点我收藏+]

标签:

                                                                    晨跑
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.
View Code

 

BZOJ 1877:[SDOI2009]晨跑(最小费用最大流)

标签:

原文地址:http://www.cnblogs.com/qtyytq/p/5392245.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!