COGS [USACO Mar07] 奶牛交通

                                              [USACO Mar07] 奶牛交通

农场中,由于奶牛数量的迅速增长,通往奶牛宿舍的道路也出现了严重的交通拥堵问题.FJ打算找出最忙碌的道路来重点整治.

这个牧区包括一个由M (1 ≤ M ≤ 50,000)条单行道路(有向)组成的网络,以及 N (1 ≤ N ≤ 5,000)个交叉路口(编号为1..N),每一条道路连接两个不同的交叉路口.奶牛宿舍位于第N个路口.每一条道路都由编号较小的路口通向编号较大的路 口.这样就可以避免网络中出现环.显而易见,所有道路都通向奶牛宿舍.而两个交叉路口可能由不止一条边连接.

在准备睡觉的时候,所有奶牛都从他们各自所在的交叉路口走向奶牛宿舍,奶牛只会在入度为0的路口,且所有入度为0的路口都会有奶牛.

帮助FJ找出最忙碌的道路,即计算所有路径中通过某条道路的最大次数.答案保证可以用32位整数存储.

输入格式:

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

输出格式:

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

样例输入:

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

样例输出

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


通过某一条路的次数就是通过某一条路两端点路径条数的乘积
f[i]代表从一个入度为0的点到f[i]的路径条数
g[i]代表从n到i点的路径条数

f[i]=Sum{ f[j] } 存在边(j,i)
g[i]=Sum{ g[j] } 存在边(i,j) 

f[k]=1 k为入度为0的点
g[n]=1;
ans=Max{ f[a]*g[b] } 存在边(a,b) 

#include <cctype>
#include <cstdio>

const int MAXN=5010;

int n,m;

int f[MAXN],g[MAXN];

struct node {
    int to;
    int next;
    node() {}
    node(int to,int next):to(to),next(next) {}
};
node e[MAXN*10],r[MAXN*10];

int head[MAXN],Head[MAXN],tot;

inline void read(int&x) {
    int f=1;register char c=getchar();
    for(x=0;!isdigit(c);c=='-'&&(f=-1),c=getchar());
    for(;isdigit(c);x=x*10+c-48,c=getchar());
    x=x*f;
}

inline int max(int a,int b) {return a<b?b:a;}

inline void add(int x,int y) {
    e[++tot]=node(y,head[x]);
    head[x]=tot;
    r[tot]=node(x,Head[y]);
    Head[y]=tot;
}

int hh() {
    freopen("cowtraffic.in","r",stdin);
    freopen("cowtraffic.out","w",stdout);
    read(n);read(m);
    int x,y;
    for(register int i=1;i<=m;++i) {
        read(x);read(y);
        add(x,y);
        if(x>y) x^=y^=x^=y;//其实没必要 
    }
    g[n]=1;
    for(register int i=1;i<=n;++i) 
      if(!Head[i]) f[i]=1;
      
    for(register int i=1;i<=n;++i)
      for(register int j=head[i];j;j=e[j].next)
        f[e[j].to]+=f[i];
        
    for(register int i=n;i>=1;--i) 
      for(register int j=Head[i];j;j=r[j].next)
        g[r[j].to]+=g[i]; 
        
    int ans=-1;
    for(register int i=1;i<=n;++i) 
      for(register int j=head[i];j;j=e[j].next)
        ans=max(ans,f[i]*g[e[j].to]);
    printf("%d
",ans);
    return 0;
}

int sb=hh();
int main(int argc,char**argv) {;}