1804: 有向无环图
Time Limit: 5 Sec Memory Limit: 128 MBSubmit: 189 Solved: 82
[ Submit][ Status][ Web Board]
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
HINT
Source
解题思路:dp[u]表示从u出发所能到达的所有v的count[v]*b[v]之和,不难发现每一个v都对它的前继u有dp[v]+b[v]的贡献,故按照反向拓扑排序进行dp,即从入度为0的节点开始进行dp,建图的时候也反向建图,p[i]表示结点i的入度。最后每个节点乘以a[i]即可。
/* ***********************************************
┆ ┏┓ ┏┓ ┆
┆┏┛┻━━━┛┻┓ ┆
┆┃ ┃ ┆
┆┃ ━ ┃ ┆
┆┃ ┳┛ ┗┳ ┃ ┆
┆┃ ┃ ┆
┆┃ ┻ ┃ ┆
┆┗━┓ 马 ┏━┛ ┆
┆ ┃ 勒 ┃ ┆
┆ ┃ 戈 ┗━━━┓ ┆
┆ ┃ 壁 ┣┓┆
┆ ┃ 的草泥马 ┏┛┆
┆ ┗┓┓┏━┳┓┏┛ ┆
┆ ┃┫┫ ┃┫┫ ┆
┆ ┗┻┛ ┗┻┛ ┆
************************************************ */
#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <stack>
#include <set>
#include <map>
#include <string>
#include <math.h>
#include <stdlib.h>
#include <bitset>
using namespace std;
#define rep(i,a,b) for (int i=(a),_ed=(b);i<=_ed;i++)
#define per(i,a,b) for (int i=(b),_ed=(a);i>=_ed;i--)
#define pb push_back
#define mp make_pair
const int inf_int = 2e9;
const long long inf_ll = 2e18;
#define inf_add 0x3f3f3f3f
#define mod 1000000007
#define LL long long
#define ULL unsigned long long
#define MS0(X) memset((X), 0, sizeof((X)))
#define SelfType int
SelfType Gcd(SelfType p,SelfType q){return q==0?p:Gcd(q,p%q);}
SelfType Pow(SelfType p,SelfType q){SelfType ans=1;while(q){if(q&1)ans=ans*p;p=p*p;q>>=1;}return ans;}
#define Sd(X) int (X); scanf("%d", &X)
#define Sdd(X, Y) int X, Y; scanf("%d%d", &X, &Y)
#define Sddd(X, Y, Z) int X, Y, Z; scanf("%d%d%d", &X, &Y, &Z)
#define reunique(v) v.resize(std::unique(v.begin(), v.end()) - v.begin())
#define all(a) a.begin(), a.end()
typedef pair<int, int> pii;
typedef pair<long long, long long> pll;
typedef vector<int> vi;
typedef vector<long long> vll;
inline int read(){int ra,fh;char rx;rx=getchar(),ra=0,fh=1;while((rx<'0'||rx>'9')&&rx!='-')rx=getchar();if(rx=='-')fh=-1,rx=getchar();while(rx>='0'&&rx<='9')ra*=10,ra+=rx-48,rx=getchar();return ra*fh;}
//#pragma comment(linker, "/STACK:102400000,102400000")
const int maxn = 1e5 + 5;
struct Edge
{
int to,nx;
}edge[maxn*2];
int head[maxn],cnt;
LL dp[maxn],a[maxn],b[maxn],p[maxn];
void addedge(int u,int v)
{
edge[cnt] = Edge{v,head[u]};
head[u] = cnt++;
}
int main()
{
//freopen("in.txt","r",stdin);
//freopen("out.txt","w",stdout);
ios::sync_with_stdio(0);
cin.tie(0);
int n,m;
while(~scanf("%d%d",&n,&m))
{
MS0(dp);
MS0(p);
MS0(head);
cnt = 1;
for(int i=1;i<=n;i++)
scanf("%lld%lld",&a[i],&b[i]);
while(m--)
{
int u,v;
scanf("%d%d",&u,&v);
addedge(v,u);
p[u]++;
}
queue<int>q;
for(int i=1;i<=n;i++)
if(p[i]==0)q.push(i);
while(!q.empty())
{
int u = q.front();
q.pop();
for(int i=head[u];i;i=edge[i].nx)
{
int v = edge[i].to;
if(v==u)continue;
dp[v]=(dp[v]+(dp[u]+b[u])%mod)%mod;
p[v]--;
if(p[v]==0)q.push(v);
}
}
LL ans = 0;
for(int i=1;i<=n;i++)
{
ans = (ans + (dp[i]*a[i])%mod)%mod;
}
printf("%lld\n",ans);
}
return 0;
}