poj 3272Cow Traffic解题报告

求一个有向图中，从所有起点（入度为0）到终点（唯一，出度为0）经过的次数最多的边经过的次数，两次拓扑排序，一次向前求出所有起点到某个点的所有路径总数f[i]，一次向后，求出某个点到终点的所有路径总数g[i]，则对一条边(u,v)，f[u]*g[v]即为边所需要经过的次数。

#include<stdio.h>
#include<string.h>
#define N 1005//数据范围没有题中描述的那么大……
#define M 10005
int map[N][N];
int edge[M][2];
int in[N],out[N],f[N],g[N];
int que[N];
void init()
{
    memset(in,0,sizeof(in));
    memset(out,0,sizeof(out));
    memset(f,0,sizeof(f));
    memset(g,0,sizeof(g));
    memset(map,0,sizeof(map));
}
int main()
{
    int n,m,i,j,fr,r,u,v,max;
    scanf("%d%d",&n,&m);
    init();
    for(i=0;i<m;i++)
    {
        scanf("%d%d",&u,&v);
        edge[i][0]=u,edge[i][1]=v;
        map[u][v]++;
        in[v]++;
        out[u]++;
    }
    fr=r=0;
    for(i=1;i<=n;i++)
    {
        if(in[i]==0)
        {
            que[r++]=i;
            f[i]=1;
        }
    }
    while(fr<r)
    {
        u=que[fr++];
        for(i=1;i<=n;i++)
        {
            if(map[u][i])
            {
                f[i]+=(map[u][i]*f[u]);
                in[i]-=map[u][i];
                if(in[i]==0)
                que[r++]=i;
            }
        }
    }
    fr=r=0;
    que[r++]=n;
    g[n]=1;
    while(fr<r)
    {
        u=que[fr++];
        for(i=1;i<=n;i++)
        {
            if(map[i][u])
            {
                g[i]+=(map[i][u]*g[u]);
                out[i]-=map[i][u];
                if(out[i]==0)
                que[r++]=i;
            }
        }
    }
    max=0;
    for(i=1;i<=m;i++)
    {
        if(max<f[edge[i][0]]*g[edge[i][1]])
        max=f[edge[i][0]]*g[edge[i][1]];
    }
    printf("%d\n",max);
    return 0;
}