1638: [Usaco2007 Mar]Cow Traffic 奶牛交通

1638: [Usaco2007 Mar]Cow Traffic 奶牛交通

Time Limit: 5 Sec  Memory Limit: 64 MB
Submit: 618  Solved: 217
[Submit][Status]

Description

农场中,由于奶牛数量的迅速增长,通往奶牛宿舍的道路也出现了严重的交通拥堵问题.FJ打算找出最忙碌的道路来重点整治. 这个牧区包括一个由M (1 ≤ M ≤ 50,000)条单行道路(有向)组成的网络,以及 N (1 ≤ N ≤ 5,000)个交叉路口(编号为1..N),每一条道路连接两个不同的交叉路口.奶牛宿舍位于第N个路口.每一条道路都由编号较小的路口通向编号较大的路口.这样就可以避免网络中出现环.显而易见,所有道路都通向奶牛宿舍.而两个交叉路口可能由不止一条边连接. 在准备睡觉的时候,所有奶牛都从他们各自所在的交叉路口走向奶牛宿舍,奶牛只会在入度为0的路口,且所有入度为0的路口都会有奶牛. 帮助FJ找出最忙碌的道路,即计算所有路径中通过某条道路的最大次数.答案保证可以用32位整数存储.

Input

第一行:两个用空格隔开的整数:N,M.

第二行到第M+1行:每行两个用空格隔开的整数ai,bi,表示一条道路从ai到bi.

Output

第一行: 一个整数,表示所有路径中通过某条道路的最大次数.

Sample Input

7 7
1 3
3 4
3 5
4 6
2 3
5 6
6 7

Sample Output

4
样例说明:

1 4
\ / \
3 6 -- 7
/ \ /
2 5
通向奶牛宿舍的所有路径:

1 3 4 6 7
1 3 5 6 7
2 3 4 6 7
2 3 5 6 7

HINT

 

Source

Silver

题解：我想揍死这个出题人——题目中明明说好的32位整数就可以的，但是当我天真的交了个没开int64的程序时，WA！！！然后我换成了int64，别的啥都没改，AC！！！这这这。。。好了步入正题——其实只需要顺着来一遍求出从出发点到各个点的路径数，再反着求一遍从各个点到奶牛宿舍的路径数（不难写的递推，此题嘛，呵呵呵，连拓扑排序都免了。。。）然后每个边的通过次数=出发点到此边的源点的路径数×此边的汇点到奶牛宿舍的路径数，这样子，才O(2N+M)，轻松水过。。。

 type
     point=^node;
     node=record
                g:longint;
                next:point;
     end;
 var
    i,j,k,m,n:longint;
    l:int64;
    a,b:array[..] of point;
    c,d:array[..] of int64;
    p:point;
 procedure add(x,y:longint);
           var p:point;
           begin
                new(p);
                p^.g:=y;
                p^.next:=a[x];
                a[x]:=p;

                new(p);
                p^.g:=x;
                p^.next:=b[y];
                b[y]:=p;
           end;
 begin
      readln(n,m);
      for i:= to m do
          begin
               readln(j,k);
               add(j,k);
          end;
      for i:= to n do
          begin
               p:=b[i];
               if p=nil then
                  c[i]:=
               else
                   begin
                        c[i]:=;
                        while p<>nil do
                              begin
                                   c[i]:=c[i]+c[p^.g];
                                   p:=p^.next;
                              end;
                   end;
          end;
      d[n]:=;
      for i:=n- downto  do
          begin
               p:=a[i];
               if p=nil then
                  d[i]:=
               else
                   begin
                        d[i]:=;
                        while p<>nil do
                              begin
                                   d[i]:=d[i]+d[p^.g];
                                   p:=p^.next;
                              end;
                   end;
          end;
      l:=;
      for i:= to n do
          begin
               p:=a[i];
               while p<>nil do
                     begin
                          k:=c[i]*d[p^.g];
                          //writeln(i,' ',p^.g,' - ',c[i],' ',d[p^.g]);
                          if k>l then l:=k;
                          p:=p^.next;
                     end;
          end;
      writeln(l);
 end.