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

2016.11.12

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

Time Limit: 5 Sec  Memory Limit: 64 MB
Submit: 877  Solved: 332
[Submit][Status][Discuss]

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

考虑求出从入度为0的点到图上每个点的i方案数f[i]，和图上每个点到n点的方案数g[i]，那么对于某条边(a,b),f[i]*g[i]即为通过改道路的路径数。

如何求？dp即可。。。。

dp方程:对于有向边(a=>b)f[b]+=f[a];

dp边界：入度为零的点f[i]=1

附代码：

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<climits>
#include<string>
#include<cstring>
#include<cmath>
#define N 5001
#define M 100010
using namespace std;
int n,m,ans;
int pos,head[N],head1[N];
int a[M],b[M];
int f[N],g[N];
struct edge{int to,next;}e[M],ed[M];
void add(int u,int v)
{
	pos++;e[pos].to=v,e[pos].next=head[u];head[u]=pos;
	ed[pos].to=u,ed[pos].next=head1[v];head1[v]=pos;
}
void dp(int u)
{
	if(!head[u]){f[u]=1;return ;}
	for(int i=head[u];i;i=e[i].next)
		{
			int v=e[i].to;
			if(!f[v])dp(v);
			f[u]+=f[v];
		}
}
void dp1(int u)
{
	if(!head1[u]){g[u]=1;return;}
	for(int i=head1[u];i;i=ed[i].next)
		{
			int v=ed[i].to;
			if(!g[v])dp1(v);
			g[u]+=g[v];
		}
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;i++)
		{
			scanf("%d%d",&a[i],&b[i]);
			add(a[i],b[i]);
		}
	for(int i=1;i<=n;i++)if(!f[i])dp(i);
	dp1(n);
	for(int i=1;i<=m;i++)
		ans=max(ans,g[a[i]]*f[b[i]]);
	printf("%d\n",ans);
}