• Time limit 5000 ms
• Memory limit 524288 kB
• OS Windows

解题思路

因为升级只能一级一级地升，所以所求期望满足了区间加的性质，可以一级一级地算，然后求前缀和、输出（状态不好，临博客涕零，不知所言）

接着扔链接（留坑）

• 这篇博客推公式的时候好像有些地方下标有点错，而且推公式的过程省了不少（雾）
• 这篇推公式的方向对我来说好新颖，进行非整数次的氪金，这样就不用死磕无穷级数了……（我的那篇也可以计算快点了）

源代码

参考别人博客写的，感觉像抄的一样

#include
    
     const int mod=1e9+7;const int MAXN=5e5+5;int T;int n,m;long long inv(long long x)//快速幂求逆元 inv[a]=a^(p-2){    long long res=1LL;    long long b=mod-2;    while(b){        if(b&1) res=res*x%mod;        x*=x;        x%=mod;        b>>=1;    }    return res%mod;}long long sum[MAXN];int main(){    scanf("%d",&T);    while(T--)    {        scanf("%d%d",&n,&m);        sum[0] = 0;        for(int i = 1; i <= n; ++i)        {            long long r, s, x, a;            scanf("%lld%lld%lld%lld", &r, &s, &x, &a);            long long p =  r * inv(s) % mod;            sum[i] = (sum[i - 1] + (a + (1 - p + mod) % mod * (a + sum[i - 1] - sum[x - 1] + mod) % mod * inv(p) % mod) % mod + mod) % mod;//关键就是这个了        }        while(m--)        {            int l, r;            scanf("%d%d", &l, &r);            long long ans = (sum[r - 1] - sum[l - 1] + mod) % mod;//这里为啥来着……            printf("%lld\n", ans);        }    }    return 0;}