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

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

 

题解

题目有点怪。dp可过

我们考虑每条路x--y要经过多少次，这取决于都多少种路径可以到x（记为f[x]），有多少种路径可以从y到n（几位g[y]）。这条路经过的次数即为f[x]*g[y]。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cstdlib>
#include<cmath>
using namespace std;
int n,m,z1,z2,head[5002],h[5002];
struct bian {int to,nx;} e[50002],e2[50002];
int x[50002],y[50002],ck[5002];
int f[5002],g[5002],ans;
void insert(int x,int y)
{
	z1++; e[z1].to=y; e[z1].nx=head[x]; head[x]=z1;
	z2++; e2[z2].to=x; e2[z2].nx=h[y]; h[y]=z2;
}
void init()
{
	scanf("%d%d",&n,&m);
	int i;
	for(i=1;i<=m;i++)
	   {scanf("%d%d",&x[i],&y[i]);
	    insert(x[i],y[i]);
	    ck[y[i]]=1;
	   }
}
void dp1(int x)
{
	if(h[x]==0) {f[x]=1; return;}
	int i;
	for(i=h[x];i;i=e2[i].nx)
	   {if(!f[e2[i].to]) dp1(e2[i].to);
	    f[x]+=f[e2[i].to];
	   }
}
void dp2(int x)
{
	int i;
	if(head[x]==0) {g[x]=1; return;}
	for(i=head[x];i;i=e[i].nx)
	   {if(!g[e[i].to]) dp2(e[i].to);
	    g[x]+=g[e[i].to];
	   }
}
int main()
{
	init(); dp1(n);
	int i;
	for(i=1;i<=n;i++)
	   {if(!ck[i]) dp2(i);}
	for(i=1;i<=m;i++)
	   ans=max(ans,f[x[i]]*g[y[i]]);
	printf("%d\n",ans);
	return 0;
}