伯努利求幂和模版

最新推荐文章于 2021-12-26 20:51:22 发布

自ran而然

最新推荐文章于 2021-12-26 20:51:22 发布

阅读量338

点赞数

分类专栏：模版数论

本文链接：https://blog.csdn.net/S1766434333/article/details/81132911

版权

模版同时被 2 个专栏收录

39 篇文章 0 订阅

订阅专栏

数论

15 篇文章 0 订阅

订阅专栏

第二伯努利数 $\sum_{i=1}^n{i^k}=\frac{{(n+B^+)}^{k+1}-(B^+)^{k+1}}{k+1}$
第一伯努利数 $\sum_{i=1}^n{i^k}=\frac{1}{k+1}\sum_{i=0}^{k}{(-1)^i \binom{k+1}{i}}B_i^-n^{k+1-i}$
这里模版求的是 $B^-$
模版借鉴 acdreamers 的博客
代码如下：

const ll mod = 1e9+7;
const int N = 1010;

ll C[N][N];
ll B[N],Inv[N];
ll Tmp[N];
ll n;
ll pre[N];

void Init()
{
    //预处理组合数
    for(int i=0; i<N; i++)
    {
        C[i][0] = C[i][i] = 1;
        if(i == 0) continue;
        for(int j=1; j<i; j++)
            C[i][j] = (C[i-1][j] % mod + C[i-1][j-1] % mod) % mod;
    }
    //预处理逆元
    Inv[1] = 1;
    for(int i=2; i<N; i++)
        Inv[i] = (mod - mod / i) * Inv[mod % i] % mod;
    //预处理伯努利数
    B[0] = 1;
    for(int i=1; i<N; i++)
    {
        ll ans = 0;
        if(i == N - 1) break;
        for(int j=0; j<i; j++)
        {
            ans += C[i+1][j] * B[j];
            ans %= mod;
        }
        ans *= -Inv[i+1];
        ans = (ans % mod + mod) % mod;
        B[i] = ans;
    }
}

ll Work(int k)  //求k次幂的和
{
    ll ans = Inv[k+1];
    ll sum = 0;
    for(int i=0; i<=k; i++)
    {
        if(i&1)
            sum=(sum-C[k+1][i] * Tmp[k+1-i] % mod * B[i] % mod+mod)%mod;
        else
            sum += C[k+1][i] * Tmp[k+1-i] % mod * B[i] % mod;
        sum %= mod;
    }
    ans *= sum;
    ans %= mod;
    return ans;
}