Kejin Player【2019 HDU多校】【DP 推理过程】

HDU-6656 题目链接

---

  这道题，看到查询的这个关系，所以，我们首先要想到的就是要去证明这个关系是线性关系，才有可差性，这个，在比赛的时候，我通过列图的方式来证明它（一开始想的方向是差分约束），然后竟然无独有偶的竟然歪打误撞给证明了它的线性关系，到目前点u的关系，是可以通过“s -> u”来推理得到的，然后呢，我们怎么去求解“u -> v”呢？

  那么，我们可以看到u -> v上的花费，其实跟"s -> u"的贡献是不影响的，并且是可以直接由"s -> v"减去"s -> u"来得到的，这个手玩吧……

  然后就是这个dp方程了，我们求一个从s起点到该点 i 的期望花费，那么就是dp[ i ]，然后又因为是线性的，所以当我们求到 i 的时候，其实i-1之前都是求出来的，所以，我们假设它有tmp次尝试升级，然后1次成功，其余（tmp - 1）次都失败的情况，那么，dp[i + 1] = dp[i]（前面i项的花费）+ (dp[i] - dp[x[i]])（从失败的跳转点回到i的花费）* (tmp - 1)（失败次数） + (tmp - 1) * a[i]（失败花费） + a[i]（一次成功）。

  化简之后呢，dp[i + 1] = dp[i] + ( dp[i] - dp[x[i]] ) * (tmp - 1) + tmp * a[i]；

  那么，我们就要去把tmp算出来，我们知道ri和si，我们可以用这个来求tmp，(tmp - 1) / tmp = 1 - ri / si；化简过后，tmp = si / ri；然后这上面的运算该取模的取模，该求逆元的求逆元，OK了。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <limits>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#define lowbit(x) ( x&(-x) )
#define pi 3.141592653589793
#define e 2.718281828459045
#define INF 0x3f3f3f3f
#define HalF (l + r)>>1
#define lsn rt<<1
#define rsn rt<<1|1
#define Lson lsn, l, mid
#define Rson rsn, mid+1, r
#define QL Lson, ql, qr
#define QR Rson, ql, qr
#define myself rt, l, r
using namespace std;
typedef unsigned long long ull;
typedef long long ll;
const int maxN = 5e5 + 7;
const ll mod = 1e9 + 7;
int N, Q, x[maxN];
ll r[maxN], s[maxN], a[maxN], dp[maxN];
inline ll fast_mi(ll a, ll b = mod - 2)
{
    ll ans = 1;
    while(b)
    {
        if(b & 1) ans = ans * a % mod;
        a = a * a % mod;
        b >>= 1;
    }
    return ans;
}
int main()
{
    int Cas; scanf("%d", &Cas);
    while(Cas--)
    {
        scanf("%d%d", &N, &Q);
        for(int i=1; i<=N; i++)
        {
            scanf("%lld%lld%d%lld", &r[i], &s[i], &x[i], &a[i]);
            ll tmp = fast_mi(r[i]) * s[i] % mod;
            dp[i + 1] = ( dp[i] + ( dp[i] - dp[x[i]] ) * (tmp - 1) % mod + tmp * a[i] % mod + mod ) % mod;
        }
        int l, r;
        while (Q--)
        {
            scanf("%d%d", &l, &r);
            printf("%lld\n", (dp[r] - dp[l] + mod) % mod);
        }
    }
    return 0;
}