hdu 6472 艺术台阶·微积分+多项式维护

题解

懒得写了，重点直接看别人的吧

注意了，
由于题目中的h是实数范围内，所以是${\sum }_{i=0}^{b[1]}{\sum }_{j=i}^{b[2}{\sum }_{k=j}^{b[3]}$，而不是${\sum }_{i=0}^{b[1]}{\sum }_{j=i+1}^{b[2}{\sum }_{k=j+1}^{b[3]}$

代码微积分计算部分过程原理：

首先，让我们模拟一遍，
假设 $ans={\int }_0^{b[1]}{\int }_i^{b[2]}{\int }_j^{b[3]}1dkdjdi$ _{最外层是0-b[1]哦}

对k层积分，
${\int }_j^{b[3]}1dk=[k{]}_j^{b[3]}=[c_0+c_{1}k{]}_j^{b[3]}$ _{这里的cx表示的是幂次数为x的项的系数，只是作为符号，值不一定相等，这里c0=0,c1=1}

得到j层的函数，
$[c_0+c_{1}k{]}_j^{b[3]}=(c_0+c_{1}b[3])-(c_0+c_{1}j)=c_{1}b[3]-c_{1}j=c_0+c_{1}j$

然后，对j层积分
${\int }_i^{b[2]}{c}_0+c_{1}jdj=[c_{1}j+c_2{j}^2{]}_i^{b[2]}$

得到i层的函数
$[c_{1}j+c_2{j}^2{]}_i^{b[2]}=(c_{1}b[2]+c_{2}b[2{]}^2)-(c_{1}i+c_2{j}^2)=c_0+c_{1}i+c_2{i}^2$

对i层积分得到最终答案
${\int }_0^{b[1]}{c}_0+c_{1}i+c_2{i}^{2}di=[c_{1}i+c_2{i}^2+c_3{i}^3{]}_0^{b[1]}=c_{1}b[1]+c_{2}b[1{]}^2+c_{3}b[1{]}^3$

通过以上式子观察得出，
微积分应该由内向外积分，
第i次积分，该层的 $c_0={\sum }_{k=1}^{i-1}$上一层的 $c_k\ast b[$本层积分的实际上限${]}^k$，求c₀是为了求出下一层的c₁
而下一层的 $c_k$ 都是由本层的$-c_{k-1}\ast \frac{1}{k}$ 转化而来（积分嘛）
最终所有计算结果统计在最外层中，

整理一下我们会用到哪些东西呢？
原上限，实际上限，各种次幂的系数，各个数的逆元

---

---

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e6+10;
const int mod=1e9+7;

ll qpow(ll a,ll b){
    ll res=1;
    while(b){
        if(b&1) res=res*a%mod;
        b>>=1;
        a=a*a%mod;
    }
    return res;
}

int n,m,k;
ll a[N],b[N],c[N];//原上限 实际上限 指数为i次的项的系数
ll inv[N];//各个数的逆元

int main(){
    ios::sync_with_stdio(0);

    inv[1]=1;
    for (int i = 2; i <= 2005; ++i) {
        inv[i]=qpow(i,mod-2)%mod;
    }

    int T;
    cin>>T;
    for (int cs = 1; cs <= T; ++cs) {
        memset(c, 0, sizeof(c));

        cin>>n;
        ll div=1;
        for (int i = 1; i <= n; ++i) {
            cin>>a[i];
            b[i]=a[i];
            div=div*a[i]%mod;
        }
        //处理实际上限
        for (int i = n; i ; --i) {
            if(b[i-1]>b[i])
                b[i-1]=b[i];
        }

        c[1]=1;c[0]=0;
        for (int i = n,p=2; i>1; --i,p++) {//p为本层函数内的最高次项的次数
            //处理本层的c[0]
            for (int j = 0; j < p; ++j) {
                c[0]+=c[j]*qpow(b[i],j)%mod;
                c[0]%=mod;
                c[j]=-c[j];
            }
            //再处理下一层的各个系数
            for (int j = p; j; --j) {
                c[j]=c[j-1]*inv[j]%mod;
            }
            c[0]=0;//记得清0 每一层的原函数的c[0]都是0
        }

        ll ans=0;
        for (int i = 0; i <= n; ++i) {
            ans+=c[i]*qpow(b[1],i)%mod;
            ans%=mod;
        }
        ans=(ans*qpow(div,mod-2)%mod+mod)%mod;
        printf("%lld\n", ans);
    }
    return 0;
}