[BZOJ 4562] 食物链 拓扑排序+递推

题目传送门:【BZOJ 4562】

题目大意:现在给你 n 个物种和 m 条能量流动关系(1 ≤ N ≤ 100000,0 ≤ m ≤ 200000),求其中的食物链条数。物种的名称从 1 到 n 编号,M条能量流动关系形如(a1 b1),(a2 b2),(a3 b3)……(am-1 bm-1),(am bm)
其中(ai bi)表示能量从物种 ai流向物种 bi。注意单独的一种孤立生物不算一条食物链。数据保证符合生物学特点,且不会有重复的能量流动关系出现。

题目分析:

由题,根据题目性质,我们可以发现,整个图就是一个 DAG。所以,我们就根据这个性质,从所有的起始点开始走,然后递推求出每个点对应的食物链条数;最后统计所有终止点的总数,全部加起来即为答案。

起始点:指入度为 0,出度不为 0 的点。
终止点:指出度为 0,入度不为 0 的点。
因为入度/出度都为 0 的点不能被统计,所以我们要在开始时单独判断一次。

下面附上代码:

[cpp] view plain copy
print ?
  1. #include  
  2. #include  
  3. #include  
  4. #include  
  5. using namespace std;  
  6. const int MX=100005;  
  7.   
  8. struct Edge{  
  9.     int to,next;  
  10. }edge[MX*4];  
  11. int n,m,head[MX],now=0,ind[MX],outd[MX];  
  12. int dp[MX],ans=0;  
  13. queue<int> q,st,ed;  
  14.   
  15. inline void adde(int u,int v){  
  16.     edge[++now].to=v;  
  17.     edge[now].next=head[u];  
  18.     head[u]=now;  
  19. }  
  20.   
  21. void bfs(int s){                            //名为 bfs,实际上是拓扑排序   
  22.     dp[s]=1;  
  23.     q.push(s);  
  24.     while (!q.empty()){  
  25.         int u=q.front();  
  26.         q.pop();  
  27.         for (int i=head[u];i;i=edge[i].next){  
  28.             int v=edge[i].to;  
  29.             dp[v]+=dp[u];  
  30.             ind[v]–;  
  31.             if (ind[v]==0) q.push(v);  
  32.         }  
  33.     }  
  34. }  
  35.   
  36. int main(){  
  37.     int a,b;  
  38.     scanf(”%d%d”,&n,&m);  
  39.     for (int i=1;i<=m;i++){  
  40.         scanf(”%d%d”,&a,&b);  
  41.         adde(a,b);  
  42.         ind[b]++,outd[a]++;  
  43.     }  
  44.     for (int i=1;i<=n;i++){                      //判断哪些点是起始点,哪些点是终止点,哪些点不能被统计   
  45.         if (ind[i]==0 && outd[i]==0) continue;  
  46.         if (ind[i]==0){  
  47.             st.push(i);  
  48.         }  
  49.         if (outd[i]==0){  
  50.             ed.push(i);  
  51.         }  
  52.     }  
  53.     while (!st.empty()){  
  54.         int i=st.front();  
  55.         bfs(i);  
  56.         st.pop();  
  57.     }  
  58.     while (!ed.empty()){  
  59.         int i=ed.front();  
  60.         ans+=dp[i];  
  61.         ed.pop();  
  62.     }  
  63.     printf(”%d”,ans);  
  64.     return 0;  
  65. }  

你可能感兴趣的:(各大OJ专题(POJ,BZOJ,hdu等),拓扑排序)