[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<cstdio>  
2. #include<cstring>  
3. #include<algorithm>  
4. #include<queue>  
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. }