CSU 1804 有向无环图 (拓扑排序+乱搞)

有向无环图
Time Limit: 5000MS Memory Limit: 131072KB 64bit IO Format: %lld & %llu
Submit

Status

Description
Bobo 有一个 n 个点，m 条边的有向无环图（即对于任意点 v，不存在从点 v 开始、点 v 结束的路径）。
为了方便，点用 1,2,…,n 编号。 设 count(x,y) 表示点 x 到点 y 不同的路径数量（规定 count(x,x)=0），Bobo 想知道

除以 (10 9+7) 的余数。
其中，a i,b j 是给定的数列。
Input
输入包含不超过 15 组数据。
每组数据的第一行包含两个整数 n,m (1≤n,m≤10 5).
接下来 n 行的第 i 行包含两个整数 a i,b i (0≤a i,b i≤10 9).
最后 m 行的第 i 行包含两个整数 u i,v i，代表一条从点 u i 到 v i 的边 (1≤u i,v i≤n)。
Output
对于每组数据，输出一个整数表示要求的值。
Sample Input
3 3
1 1
1 1
1 1
1 2
1 3
2 3
2 2
1 0
0 2
1 2
1 2
2 1
500000000 0
0 500000000
1 2
Sample Output
4
4
250000014

先拓扑排序一下，根据拓扑序我们维护一个链条（单一路径），对每个节点做合并即可（因为可能有多条路径，每条都要算）。

#include "cstring"#include "iostream"#include "string.h"#include "cstdio"#include "queue"#define MOD 1000000007#define MAX 100005using namespace std;int main(){    int n,m;    while(~scanf("%d%d",&n,&m))    {        vector<int> list[MAX];        long long a[MAX],b[MAX];        int degree[MAX];        long long sum[MAX];        memset(degree,0,sizeof(degree));        memset(sum,0,sizeof(sum));        for(int i=1;i<=n;i++)        {            scanf("%lld%lld",&a[i],&b[i]);            sum[i]=b[i];        }        for(int i=1;i<=m;i++)        {            int fa,son;            scanf("%d%d",&fa,&son);            degree[fa]++;            list[son].push_back(fa);        }        queue<int> original;        queue<int> sorted;        for(int i=1;i<=n;i++)        {            if(degree[i]==0)                original.push(i);        }        //topology sort        while(!original.empty())        {            int son=original.front();            original.pop();            for(int i=0;i<list[son].size();i++)            {                int fa=list[son][i];                degree[fa]--;                if(degree[fa]==0)                    original.push(fa);            }            sorted.push(son);        }        long long ans=0;        while(!sorted.empty())        {            int son=sorted.front();            sorted.pop();            for(int i=0;i<list[son].size();i++)            {                int fa=list[son][i];                ans=(ans+(a[fa]*sum[son])%MOD)%MOD;                sum[fa]=(sum[fa]+sum[son])%MOD;            }        }        printf("%lld\n",ans%MOD);    }}