bzoj1638[Usaco2007 Mar]Cow Traffic 奶牛交通*

bzoj1638[Usaco2007 Mar]Cow Traffic 奶牛交通

题意：

N点M边有向图，每个入度为0的点都有无限只奶牛，现在它们要回宿舍（点1），求通过量最大的路的通过量。N≤5000，M≤50000

题解：

一条路的通过量=到达节点到入度为0节点的方案数*点1到出发节点的方案数（其实我也不知道为什么，这题意完全是模糊的），2次dfs就行了。

代码：

 1 #include <cstdio>
 2 #include <cstring>
 3 #include <algorithm>
 4 #define inc(i,j,k) for(int i=j;i<=k;i++)
 5 #define dec(i,j,k) for(int i=j;i>=k;i--)
 6 #define maxn 5010
 7 using namespace std;
 8 
 9 inline int read(){
10     char ch=getchar(); int f=1,x=0;
11     while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}
12     while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();
13     return f*x;
14 }
15 int f1[maxn],f2[maxn],n,m,f[maxn*10],t[maxn*10];
16 struct e{int t,n;}; e es[maxn*10]; int g[maxn];
17 int dfs1(int x){
18     if(f1[x])return f1[x];
19     if(!g[x])f1[x]=1;else{
20         for(int i=g[x];i;i=es[i].n){f1[x]+=dfs1(es[i].t);}
21     }
22     return f1[x];
23 }
24 int dfs2(int x){
25     if(f2[x])return f2[x]; 
26     if(!g[x])f2[x]=1;else{
27         for(int i=g[x];i;i=es[i].n){f2[x]+=dfs2(es[i].t);}
28     }
29     return f2[x];
30 }
31 int main(){
32     n=read(); m=read(); inc(i,1,m)f[i]=read(),t[i]=read();
33     memset(g,0,sizeof(g)); inc(i,1,m)es[i]=(e){f[i],g[t[i]]},g[t[i]]=i; dfs1(n);
34     memset(g,0,sizeof(g)); inc(i,1,m)es[i]=(e){t[i],g[f[i]]},g[f[i]]=i; inc(i,1,n)if(!f2[i])dfs2(i);
35     inc(i,1,m)f1[0]=max(f1[0],f1[f[i]]*f2[t[i]]);
36     printf("%d",f1[0]); return 0;
37 }

20160808