2017年8月8日提高组T3 题目

Description

小C旅行到了美丽的x市。热爱oi的小C把x市分成了n个建筑物和m条双向街道，每一条街道都有一个通过时间，每个建筑物都有一个观赏值v。接着小C就在想了，如果我从任意一个点出发，经过至少两个点（包括出发点）后回到出发点，那么我能得到的最大平均观赏值是多少。平均观赏值指的是总观赏值/总耗费时间。你可以理解成观赏是不需要耗费时间的，且若多次到达同一建筑物，观赏值只会被计算一次。但小C还要和他的好朋友们开黑打农药，于是他就把这个又简单又无聊的问题交给了准备参加NOIP的你。

Input

第一行两个数分别表示n和m.
第二行n个整数，表示每个建筑物的观赏值v。
接下来m行，每行包含三个整数u,v,t，表示一条道路的两个端点和通过时间。

Output

输出一个实数，保留两位小数，表示最大平均观赏值。

Sample Input

5 7
30
10
10
5
10
1 2 3
2 3 2
3 4 5
3 5 2
4 5 5
5 1 3
5 2 2

Sample Output

6.00

Hint

【样例说明】
如果小C选择1 -> 2 -> 3 -> 5 -> 1的旅行路线，他能得到的总观赏值为60，为此他得花费10单位的时间在走路上。于是小C在这次旅行中的平均观赏值为6。如果他走2 -> 3 -> 5 -> 2的路线，就只能得到30/6 = 5的平均观赏值。并且，任何去参观建筑物4的旅行路线的平均乐趣值都没有超过4。

【数据规模与约定】
对于30%的数据，n<=6,m<=15
对于60%的数据，n<=100,m<=500
对于100%的数据，n<=1000,m<=5000,1<=v<=1000,1<=t<=1000，保证没有自环。

分析：
这题是01分数规划。（大佬说的）
一定要二分答案。
对于二分来的一个数ans,有sum(v[i])/sum(t[i])>=ans
{v指价值，t指时间}
ans才是正确的，可行的，但不是最优的。

化简为
sum(v[i])>=sum(t[i])*ans
sum(v[i])-sum(t[i])*ans>=0

因为路径为环，则点数=边数。
化简有
(v[1]-t[1]*ans)+(v[2]-t[2]**ans)+……+
(v[k]-v[k]*ans)>=0

k为环中点树，且该环为简单环。

在转化一下：
t[i]*ans-v[i]<=0

这就是一个判断是否存在负环的问题了。

用spfa跑最短路，把不等式变为小于等于，其思想和查分约束差不多。

最大的，满足上述条件的即为答案。注意：存在负环才满足。

代码：

const
 maxn=100;
 maxv=5001;
 maxe=1001;
type
 node=record
  x,y,next:longint;
  w:longint;
 end;

var
 list,v,ls,a,sum:array [0..maxe] of longint;
 d:array [0..maxe] of longint;
 g:array [0..maxv] of node;
 i,j,k,p,o,max,n,m:longint;


procedure add(x,y:longint; w:longint);
 begin
  inc(o);
  g[o].x:=x; g[o].y:=y; g[o].w:=w; g[o].next:=ls[x]; ls[x]:=o;
 end;

procedure init;
 var i,j,x,y,w:longint;
begin
 readln(n,m);
 for i:=1 to n do
  begin
   read(a[i]);
   a[i]:=a[i]*1000;
  end;
 for i:=1 to m do
  begin
   readln(x,y,w);
   add(x,y,w);
  end;
end;

function spfa(first:longint;z:longint):boolean;
var
  head,tail,t,i,qe,sw:longint;
begin
  fillchar(sum,sizeof(sum),0);
  fillchar(v,sizeof(v),0);
  tail:=1;
  list[1]:=first;
  for i:=1 to n do
  d[i]:=maxlongint;
  d[first]:=0;
  v[first]:=1;
  while tail<>0 do
    begin
      t:=ls[list[tail]];
      qe:=list[tail];
      tail:=tail-1;
      while t>0 do
        with g[t] do
          begin
            sw:=z*w-a[y];
            if d[x]+sw<d[y] then
              begin
                d[y]:=d[x]+sw;
                if v[y]=0 then
                  begin
                    v[y]:=1;
                    tail:=tail+1;
                    list[tail]:=y;
                    inc(sum[y]);
                    if (sum[y]>=n) then
                     begin
                      exit(true);
                     end;
                  end;
              end;
            t:=next;
          end;
      v[qe]:=0;
    end;
  exit(false);
end;


procedure main;
 var l,r:longint;
     mid:longint;
begin
 max:=0;
 l:=1; r:=1000000;
 while l<=r do
  begin
   mid:=(l+r) shr 1;
   if spfa(1,mid)
    then
     begin
      l:=mid+1;
      if mid>max then max:=mid;
     end
    else r:=mid-1;
  end;
end;

procedure print;
 begin
  writeln(max/1000:0:2);
 end;

begin
 init;
 main;
 print;
end.