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

bzoj1638[Usaco2007 Mar]Cow Traffic 奶牛交通

题意：

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

题解：

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

代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
#define inc(i,j,k) for(int i=j;i<=k;i++)
#define dec(i,j,k) for(int i=j;i>=k;i--)
#define maxn 5010
using namespace std;

inline int read(){
    char ch=getchar(); int f=1,x=0;
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}
    while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();
    return f*x;
}
int f1[maxn],f2[maxn],n,m,f[maxn*10],t[maxn*10];
struct e{int t,n;}; e es[maxn*10]; int g[maxn];
int dfs1(int x){
    if(f1[x])return f1[x];
    if(!g[x])f1[x]=1;else{
        for(int i=g[x];i;i=es[i].n){f1[x]+=dfs1(es[i].t);}
    }
    return f1[x];
}
int dfs2(int x){
    if(f2[x])return f2[x]; 
    if(!g[x])f2[x]=1;else{
        for(int i=g[x];i;i=es[i].n){f2[x]+=dfs2(es[i].t);}
    }
    return f2[x];
}
int main(){
    n=read(); m=read(); inc(i,1,m)f[i]=read(),t[i]=read();
    memset(g,0,sizeof(g)); inc(i,1,m)es[i]=(e){f[i],g[t[i]]},g[t[i]]=i; dfs1(n);
    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);
    inc(i,1,m)f1[0]=max(f1[0],f1[f[i]]*f2[t[i]]);
    printf("%d",f1[0]); return 0;
}

 

20160808