BZOJ 4011 HNOI2015 落忆枫音

最新推荐文章于 2019-08-21 18:03:00 发布

Fuxey

最新推荐文章于 2019-08-21 18:03:00 发布

阅读量865

点赞数

分类专栏：朱刘算法 dp 文章标签：朱刘算法 dp

本文链接：https://blog.csdn.net/Fuxey/article/details/50864481

版权

dp 同时被 2 个专栏收录

23 篇文章 0 订阅

订阅专栏

朱刘算法

2 篇文章 0 订阅

订阅专栏

故事很有意思啊QAQ 考试的时候应该配乐诗朗诵才对……

提示：
1. 朱刘算法里面一个核心思想被运用啦，我有一篇博文讲这玩意，边栏里可以看到。
2. 此题并不难，搞一搞就应该可以出来的……

详细题解在代码后：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <string>
#include <queue>
#include <algorithm>
#include <vector>
#include <deque>
#include <set>
#include <map>

using namespace std;
typedef long long ll;
const ll maxn = 1e5+1e2;
const ll modu = 1e9 + 7;
ll n , m , t , s , in[maxn] , d[maxn];
vector<ll> g[maxn];

void exGcd(ll a , ll b , ll& d , ll& x , ll& y)
{
    if(!b) x = 1 , y = 0 , d = a;
    else 
    {
        exGcd(b, a%b, d, y, x);
        y -= x*(a/b);
    }
}

ll inv(ll a)
{
    ll x , y , d;
    exGcd(a, modu, d, x, y);

    return x;
}

ll dp(ll u)
{
    if(d[u]!=-1) return d[u];

    d[u] = 0;
    for(ll i=0;i<g[u].size();i++) d[u] = (d[u] + dp(g[u][i]))%modu;

    return d[u] = (d[u] * inv(in[u]))%modu;
}

int main()
{   
    cin>>n>>m>>t>>s;

    while(m--)
    {
        ll a , b;
        scanf("%lld%lld" ,&a , &b);
        if(b==1) continue;
        in[b] ++;
        g[a].push_back(b);
    }
    in[s] ++;

    memset(d , -1 , sizeof d); 
    d[t] = inv(in[t]);
    dp(s);
    ll res = 1;
    for(ll i=2;i<=n;i++) res = (res * in[i])%modu;

    if(s!=1) res = ((res - res*d[s])%modu + modu)%modu;
    cout<<res<<endl;

    return 0;
}