伯努利数（详解 + 例题 ：P3711 仓鼠的数学题）

最新推荐文章于 2021-06-28 18:48:54 发布
最新推荐文章于 2021-06-28 18:48:54 发布
阅读量1.5k 收藏 3
点赞数 5
分类专栏： 多项式及生成函数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_45483201/article/details/114414524
版权

伯努利数

定义 S k ( n ) = ∑ i = 0 n − 1 i k S_k(n) = \sum\limits_{i = 0} ^{n - 1} i ^ k Sk​(n)=i=0∑n−1​ik。

从二项式出发

( 0 + 1 ) k + 1 = ∑ i = 0 k C k + 1 i 0 i + 0 k + 1 ⋮ ( n − 1 + 1 ) k + 1 = ∑ i = 0 k C k + 1 i ( n − 1 ) i + ( n − 1 ) k + 1 把 次 方 k + 1 的 移 项 ， 再 整 体 相 加 ， 得 n k + 1 = ∑ i = 0 k C k + 1 i S i ( n ) n k + 1 = ∑ i = 0 k − 1 C k + 1 i S i ( n ) + ( k + 1 ) S k ( n ) S k ( n ) = 1 k + 1 ( n k + 1 − ∑ i = 0 k − 1 C k + 1 i S i ( n ) ) (0 + 1) ^ {k + 1} = \sum_{i = 0} ^{k} C_{k + 1} ^{i} 0 ^ {i} + 0 ^{k + 1}\\ \vdots\\ (n - 1 + 1) ^{k + 1} = \sum_{i = 0} ^{k} C_{k + 1} ^{i} (n - 1) ^{i} + (n - 1) ^{k + 1}\\ 把次方k + 1的移项，再整体相加，得\\ n ^{k + 1} = \sum_{i = 0} ^{k} C_{k + 1} ^{i} S_i(n)\\ n ^{k + 1} = \sum_{i = 0} ^{k - 1} C_{k + 1} ^{i} S_i(n) + (k + 1) S_k(n)\\ S_k(n) = \frac{1}{k + 1}\left(n ^{k + 1} - \sum_{i = 0} ^{k - 1} C_{k + 1} ^{i} S_i(n) \right)\\ (0+1)k+1=i=0∑k​Ck+1i​0i+0k+1⋮(n−1+1)k+1=i=0∑k​Ck+1i​(n−1)i+(n−1)k+1把次方k+1的移项，再整体相加，得nk+1=i=0∑k​Ck+1i​Si​(n)nk+1=i=0∑k−1​Ck+1i​Si​(n)+(k+1)Sk​(n)Sk​(n)=k+11​(nk+1−i=0∑k−1​Ck+1i​Si​(n))
于是我们有了一个 O ( n 2 ) O(n ^ 2) O(n2)递推求 S i ( n ) S_i(n) Si​(n)得算法，

对二项式再进一步考虑：

( 0 + m ) k = ∑ i = 0 k C k i 0 i m k − i ⋮ ( n − 1 + m ) k = ∑ i = 0 k C k i ( n − 1 ) i m k − i 左 右 两 边 同 时 相 加 S k ( n + m ) − S k ( m ) = ∑ i = 0 k C k i S i ( n ) m k − i 左 右 两 边 对 n 求 导 S k ′ ( n + m ) = ∑ i = 0 k C k i S i ′ ( n ) m k − i n = 0 ， 代 入 S k ′ ( m ) = ∑ i = 0 k C k i S i ′ ( 0 ) m k − i 两 边 对 m 从 0 − > n 求 积 分 S k ( n ) − S k ( 0 ) = ∑ i = 0 k C k i S i ′ ( 0 ) n k − i + 1 k − i + 1 − ∑ i = 0 k C k i S i ′ ( 0 ) 0 k − i + 1 k − i + 1 S k ( 0 ) = 0 , 对 于 最 右 侧 i ≠ k + 1 , 所 以 0 k − i + 1 = 0 S k ( n ) = ∑ i = 0 k C k i S i ′ ( 0 ) n k − i + 1 k − i + 1 = 1 k + 1 ∑ i = 0 k C k + 1 i S i ′ ( 0 ) n k + i − 1 (0 + m) ^{k} = \sum_{i = 0} ^{k} C_{k} ^{i} 0 ^{i} m ^{k - i}\\ \vdots\\ (n - 1 + m) ^{k} = \sum_{i = 0} ^{k} C_{k} ^{i}(n - 1) ^{i} m ^{k - i}\\ 左右两边同时相加\\ S_k(n + m) - S_k(m) = \sum_{i = 0} ^{k} C_{k} ^{i} S_{i} (n) m ^{k - i}\\ 左右两边对n求导\\ S'_k(n + m) = \sum_{i = 0} ^{k} C_{k} ^{i} S'_i(n) m^{k - i}\\ \\ n = 0，代入\\ S'_k(m) = \sum_{i = 0} ^{k} C_{k} ^{i} S'_i(0) m^{k - i}\\ 两边对m从0->n求积分\\ S_k(n) - S_k(0) = \sum_{i = 0} ^{k} C_{k} ^{i} S'_i(0) \frac{n ^{k - i + 1}}{k - i + 1} - \sum_{i = 0} ^{k} C_{k} ^{i} S'_i(0) \frac{0 ^{k - i + 1}}{k - i + 1}\\ S_k(0) = 0, 对于最右侧i \neq k + 1,所以0 ^{k - i + 1} = 0 \\ S_k(n) = \sum_{i = 0} ^{k} C_{k} ^{i} S'_i(0) \frac{n ^{k - i + 1}}{k - i + 1} = \frac{1}{k + 1} \sum_{i = 0} ^{k} C_{k + 1} ^{i} S'_i(0) n ^{k + i - 1}\\ (0+m)k=i=0∑k​Cki​0imk−i⋮(n−1+m)k=i=0∑k​Cki​(n−1)imk−i左右两边同时相加Sk​(n+m)−Sk​(m)=i=0∑k​Cki​Si​(n)mk−i左右两边对n求导Sk′​(n+m)=i=0∑k​Cki​Si′​(n)mk−in=0，代入Sk′​(m)=i=0∑k​Cki​Si′​(0)mk−i两边对m从0−>n求积分Sk​(n)−Sk​(0)=i=0∑k​Cki​Si′​(0)k−i+1nk−i+1​−i=0∑k​Cki​Si′​(0)k−i+10k−i+1​Sk​(0)=0,对于最右侧i​=k+1,所以0k−i+1=0Sk​(n)=i=0∑k​Cki​Si′​(0)k−i+1nk−i+1​=k+11​i=0∑k​Ck+1i​Si′​(0)nk+i−1
其实，这里的 S i ( 0 ) S_i(0) Si​(0)就是伯努利数，为了方便书写， B i = S i ′ ( 0 ) B_i = S'_i(0) Bi​=Si′​(0)。
我 们 先 将 n = 1 代 入 上 式 S k ( 1 ) = 1 k + 1 ∑ i = 0 k C k + 1 i B i 同 样 的 我 们 把 右 边 的 最 高 项 给 提 出 来 S k ( 1 ) = 1 k + 1 ∑ i = 0 k − 1 C k + 1 i B i + B k 则 有 B k − S k ( 1 ) = − 1 k + 1 ∑ i = 0 k − 1 C k + 1 i B i 特 殊 地 ， 考 虑 k = 0 , 有 S 0 ( 1 ) = 1 , 所 以 B 0 = 1 k ≥ 1 , S k ( 1 ) = 0 , 所 以 B k = − 1 k + 1 ∑ i = 0 k − 1 C k + 1 i B i 我们先将n = 1代入上式\\ S_k(1) = \frac{1}{k + 1} \sum_{i = 0} ^{k} C_{k + 1} ^{i} B_i\\ 同样的我们把右边的最高项给提出来\\ S_k(1) = \frac{1}{k + 1} \sum_{i = 0} ^{k - 1} C_{k + 1} ^{i} B_i + B_k\\ 则有B_k - S_k(1) = - \frac{1}{k + 1} \sum_{i = 0} ^{k - 1} C_{k + 1} ^{i} B_i\\ 特殊地，考虑k = 0,有S_0(1) = 1, 所以B_0 = 1\\ k \geq 1, S_k(1) = 0, 所以B_k = - \frac{1}{k + 1} \sum_{i = 0} ^{k - 1} C_{k + 1} ^{i} B_i\\ 我们先将n=1代入上式Sk​(1)=k+11​i=0∑k​Ck+1i​Bi​同样的我们把右边的最高项给提出来Sk​(1)=k+11​i=0∑k−1​Ck+1i​Bi​+Bk​则有Bk​−Sk​(1)=−k+11​i=0∑k−1​Ck+1i​Bi​特殊地，考虑k=0,有S0​(1)=1,所以B0​=1k≥1,Sk​(1)=0,所以Bk​=−k+11​i=0∑k−1​Ck+1i​Bi​
由此我们可以 O ( k 2 ) O(k ^ 2) O(k2)递推出伯努利数出来了。

F a u l h a b e r Faulhaber Faulhaber公式

S k ( x ) = 1 k + 1 ∑ i = 0 k C k + 1 i B i x k − i + 1 S_k(x) = \frac{1}{k + 1} \sum\limits_{i = 0} ^{k}C_{k + 1} ^{i}B_i x ^{k - i + 1} Sk​(x)=k+11​i=0∑k​Ck+1i​Bi​xk−i+1

一个自然数幂次求和的公式，其实就是我们推导过程中的一步。

#include <bits/stdc++.h>

using namespace std;

const int N = 3e3 +10, mod = 1e9 + 7;

int C[N][N], B[N], inv[N];

void init() {
  for (int i = 0; i < N; i++) {
    C[i][0] = C[i][i] = 1;
    for (int j = 1; j < i; j++) {
      C[i][j] = (C[i - 1][j] + C[i - 1][j - 1]) % mod;
    }
  }
  inv[1] = 1;
  for (int i = 2; i < N; i++) {
    inv[i] = 1ll * (mod - mod / i) * inv[mod % i] % mod;
  }
  B[0] = 1;
  for (int i = 1; i < N - 1; i++) {
    int cur = 0;
    for (int j = 0; j < i; j++) {
      cur = (cur + 1ll * C[i + 1][j] * B[j] % mod) % mod;
    }
    cur = 1ll * cur * inv[i + 1] % mod;
    B[i] = (mod - cur) % mod;
  }
}

int main() {
  // freopen("in.txt", "r", stdin);
  // freopen("out.txt", "w", stdout);
  // ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
  init();
  int T, k;
  scanf("%d", &T);
  while (T--) {
    long long n, ans = 0, cur = 1;
    scanf("%lld %d", &n, &k);
    n++;
    n %= mod;
    for (int i = k; i >= 0; i--) {
      cur = cur * n % mod;
      ans = (ans + 1ll * C[k + 1][i] * B[i] % mod * cur % mod) % mod;
    }
    printf("%lld\n", ans * inv[k + 1] % mod);
  }
  return 0;
}

考虑生成函数

∑ i = 0 n C n + 1 i B i = [ n = 0 ] 对 两 边 加 上 B n + 1 ∑ i = 0 n + 1 C n + 1 i B i = [ n = 0 ] + B n + 1 ∑ i = 0 n C n i B i = [ n = 1 ] + B n ∑ i = 0 n 1 i ! ( n − i ) ! B i = [ n = 1 ] + B n n ! 设 指 数 型 生 成 函 数 B ( x ) = ∑ n ≥ 0 B n x n 有 B ( x ) e x = x + B ( x ) 从 而 有 B ( x ) = x e x − 1 B ( x ) = x ∑ n ≥ 0 1 n ! x n − 1 x ∑ n ≥ 1 1 n ! x n 1 ∑ n ≥ 0 1 ( n + 1 ) ! x n \sum_{i = 0} ^{n} C_{n + 1} ^{i} B_i = [n = 0]\\ 对两边加上B_{n + 1}\\ \sum_{i = 0} ^{n + 1} C_{n + 1} ^{i} B_i = [n = 0] + B_{n + 1}\\ \sum_{i = 0} ^{n} C_{n} ^{i} B_i = [n = 1] + B_{n}\\ \sum_{i = 0} ^{n} \frac{1}{i!(n - i)!} B_{i} = [n = 1] + \frac{B_{n}}{n!}\\ 设指数型生成函数B(x) = \sum_{n \geq 0} B_n x ^ n\\ 有B(x) e ^x = x + B(x)\\ 从而有B(x) = \frac{x}{e ^ x - 1}\\ B(x) = \frac{x}{\sum\limits_{n \geq 0} \frac{1}{n!} x ^ n - 1}\\ \frac{x}{\sum\limits_{n \geq 1} \frac{1}{n!} x ^ n}\\ \frac{1}{\sum\limits_{n \geq 0} \frac{1}{(n + 1)!}x ^ n}\\ i=0∑n​Cn+1i​Bi​=[n=0]对两边加上Bn+1​i=0∑n+1​Cn+1i​Bi​=[n=0]+Bn+1​i=0∑n​Cni​Bi​=[n=1]+Bn​i=0∑n​i!(n−i)!1​Bi​=[n=1]+n!Bn​​设指数型生成函数B(x)=n≥0∑​Bn​xn有B(x)ex=x+B(x)从而有B(x)=ex−1x​B(x)=n≥0∑​n!1​xn−1x​n≥1∑​n!1​xnx​n≥0∑​(n+1)!1​xn1​
只需多项式求逆即可解决。

#include <bits/stdc++.h>

using namespace std;

const int N = 1e6 + 10, mod = 998244353;

int r[N], b[N], t[N];

int quick_pow(int a, int n) {
  int ans = 1;
  while (n) {
    if (n & 1) {
      ans = 1ll * a * ans % mod;
    }
    a = 1ll * a * a % mod;
    n >>= 1;
  }
  return ans;
}

void get_r(int lim) {
  for (int i = 0; i < lim; i++) {
    r[i] = (i & 1) * (lim >> 1) + (r[i >> 1] >> 1);
  }
}

void NTT(int *f, int lim, int rev) {
  for (int i = 0; i < lim; i++) {
    if (i < r[i]) {
      swap(f[i], f[r[i]]);
    }
  }
  for (int mid = 1; mid < lim; mid <<= 1) {
    int wn = quick_pow(3, (mod - 1) / (mid << 1));
    for (int len = mid << 1, cur = 0; cur < lim; cur += len) {
      int w = 1;
      for (int k = 0; k < mid; k++, w = 1ll * w * wn % mod) {
        int x = f[cur + k], y = 1ll * w * f[cur + mid + k] % mod;
        f[cur + k] = (x + y) % mod, f[cur + mid + k] = (x - y + mod) % mod;
      }
    }
  }
  if (rev == -1) {
    int inv = quick_pow(lim, mod - 2);
    reverse(f + 1, f + lim);
    for (int i = 0; i < lim; i++) {
      f[i] = 1ll * f[i] * inv % mod;
    }
  }
}

void polyinv(int *f, int *g, int n) {
  if (n == 1) {
    g[0] = quick_pow(f[0], mod - 2);
    return ;
  }
  polyinv(f, g, n + 1 >> 1);
  for (int i = 0; i < n; i++) {
    t[i] = f[i];
  }
  int lim = 1;
  while (lim < 2 * n) {
    lim <<= 1;
  }
  get_r(lim);
  NTT(t, lim, 1);
  NTT(g, lim, 1);
  for (int i = 0; i < lim; i++) {
    int cur = (2 - 1ll * g[i] * t[i] % mod + mod) % mod;
    g[i] = 1ll * g[i] * cur % mod;
    t[i] = 0;
  }
  NTT(g, lim, -1);
  for (int i = n; i < lim; i++) {
    g[i] = 0;
  }
}

int fac[N], ifac[N], B[N], n;

void init() {
  fac[0] = 1;
  for (int i = 1; i < N; i++) {
    fac[i] = 1ll * fac[i - 1] * i % mod;
  }
  ifac[N - 1] = quick_pow(fac[N - 1], mod - 2);
  for (int i = N - 2; i >= 0; i--) {
    ifac[i] = 1ll * ifac[i + 1] * (i + 1) % mod;
  }
}

int main() {
  // freopen("in.txt", "r", stdin);
  // freopen("out.txt", "w", stdout);
  // ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
  init();
  n = 100000;
  for (int i = 0; i <= n + 1; i++) {
    B[i] = ifac[i + 1];
  }
  polyinv(B, b, n + 2);
  for (int i = 0; i <= n + 1; i++) {
    B[i] = 1ll * b[i] * fac[i] % mod;
  }
  for (int i = 0; i <= n + 1; i++) {
    printf("%d%c", B[i], i == n + 1 ? '\n' : ' ');
  }
  return 0;
}

P3711 仓鼠的数学题

F ( x ) = ∑ k = 0 n S k ( x ) a k 原 本 定 义 的 S k ( x ) = ∑ i = 0 x i k 根 据 伯 努 利 数 的 定 义 S k ′ ( x ) = ∑ i = 0 x − 1 i k 则 我 们 求 F ( x ) = ∑ k = 0 n S k ′ ( x ) a k , 答 案 即 为 F ( x + 1 ) 考 虑 先 求 F ( x ) = ∑ k = 0 n a k 1 k + 1 ∑ i = 0 k C k + 1 i B i x k − i + 1 ∑ k = 0 n a k k + 1 ∑ i = 0 k ( k + 1 ) ! i ! ( k + 1 − i ) ! B i x k − i + 1 ∑ k = 0 n a k k ! ∑ i = 0 k B i i ! x k − i + 1 ( k − i + 1 ) ! ∑ i = 1 n + 1 x i i ! ∑ k = i − 1 n ( a k k ! ) ( B k + 1 − i ( k + 1 − i ) ! ) 设 f ( n ) = a k k ! , g ( n ) = B n n ! , 另 g ( n ) = g ( n − i ) 即 翻 转 一 下 ∑ i = 1 n + 1 x i i ! ∑ k = i − 1 n f ( k ) g ( n + i − k ) 容 易 发 现 后 面 是 一 个 卷 积 ， 所 以 可 以 优 化 F(x) = \sum_{k = 0} ^{n} S_k(x) a_k\\ 原本定义的S_k(x) = \sum_{i = 0} ^{x} i ^ k\\ 根据伯努利数的定义S'_k(x) = \sum_{i = 0} ^{x - 1} i ^ k\\ 则我们求F(x) = \sum_{k = 0} ^{n} S'_k(x) a_k, 答案即为F(x + 1)\\ 考虑先求F(x) = \sum_{k = 0} ^{n} a_k \frac{1}{k + 1} \sum\limits_{i = 0} ^{k}C_{k + 1} ^{i}B_i x ^{k - i + 1}\\ \sum_{k = 0} ^{n} \frac{a_k}{k + 1} \sum_{i = 0} ^{k} \frac{(k + 1)!}{i!(k + 1 - i)!} B_i x ^{k - i + 1}\\ \sum_{k = 0} ^{n} a_k k! \sum_{i = 0} ^{k} \frac{B_i}{i!}\frac{x ^{k - i + 1}}{(k - i + 1)!}\\ \sum_{i = 1} ^{n + 1} \frac{x ^{i}}{i!} \sum_{k = i - 1} ^{n} \left(a_k k!\right) \left( \frac{B_{k + 1 - i}}{(k + 1 - i)!} \right)\\ 设f(n) = a_k k!, g(n) = \frac{B_{n}}{n!}, 另g(n) = g(n - i)即翻转一下\\ \sum_{i = 1} ^{n + 1} \frac{x ^ i}{i!} \sum_{k = i - 1} ^{n} f(k) g(n + i - k)\\ 容易发现后面是一个卷积，所以可以优化\\ F(x)=k=0∑n​Sk​(x)ak​原本定义的Sk​(x)=i=0∑x​ik根据伯努利数的定义Sk′​(x)=i=0∑x−1​ik则我们求F(x)=k=0∑n​Sk′​(x)ak​,答案即为F(x+1)考虑先求F(x)=k=0∑n​ak​k+11​i=0∑k​Ck+1i​Bi​xk−i+1k=0∑n​k+1ak​​i=0∑k​i!(k+1−i)!(k+1)!​Bi​xk−i+1k=0∑n​ak​k!i=0∑k​i!Bi​​(k−i+1)!xk−i+1​i=1∑n+1​i!xi​k=i−1∑n​(ak​k!)((k+1−i)!Bk+1−i​​)设f(n)=ak​k!,g(n)=n!Bn​​,另g(n)=g(n−i)即翻转一下i=1∑n+1​i!xi​k=i−1∑n​f(k)g(n+i−k)容易发现后面是一个卷积，所以可以优化
下面的 n = n + 1 n = n + 1 n=n+1，由上易知 F ( x ) F(x) F(x)是一个 n + 1 n + 1 n+1次多项式。
设 F ( x ) = ∑ i = 0 n a i x i F ( x + 1 ) = ∑ i = 0 n a i ( x + 1 ) i ∑ i = 0 n a i ∑ j = 0 i C j i x j ∑ i = 0 n a i ∑ j = 0 i i ! j ! ( i − j ) ! x j ∑ i = 0 n x i i ! ∑ j = i n a j j ! 1 ( j − i ) ! 另 f ( n ) = a n n ! , g ( n ) = 1 n ! , 再 翻 转 g ( n ) ∑ i = 0 n x i i ! ∑ j = i n f ( j ) g ( n + i − j ) 后 面 也 同 样 是 一 个 卷 积 ， 所 以 可 以 优 化 设F(x) = \sum_{i = 0} ^{n} a_i x ^ i\\ F(x + 1) = \sum_{i = 0} ^{n} a_i (x + 1) ^ i\\ \sum_{i = 0} ^{n} a_i \sum_{j = 0} ^{i} C_{j} ^{i} x ^ j\\ \sum_{i = 0} ^{n} a_i \sum_{j = 0} ^{i} \frac{i!}{j!(i - j)!} x ^ j\\ \sum_{i = 0} ^{n} \frac{x ^ i}{i!} \sum_{j = i} ^{n} a_j j! \frac{1}{(j - i)!}\\ 另f(n) = a_n n!, g(n) = \frac{1}{n!},再翻转g(n)\\ \sum_{i = 0} ^{n} \frac{x ^{i}}{i!} \sum_{j = i} ^{n} f(j) g(n + i - j)\\ 后面也同样是一个卷积，所以可以优化\\ 设F(x)=i=0∑n​ai​xiF(x+1)=i=0∑n​ai​(x+1)ii=0∑n​ai​j=0∑i​Cji​xji=0∑n​ai​j=0∑i​j!(i−j)!i!​xji=0∑n​i!xi​j=i∑n​aj​j!(j−i)!1​另f(n)=an​n!,g(n)=n!1​,再翻转g(n)i=0∑n​i!xi​j=i∑n​f(j)g(n+i−j)后面也同样是一个卷积，所以可以优化

#include <bits/stdc++.h>

using namespace std;

const int N = 1e6 + 10, mod = 998244353;

int r[N], t[N];

int quick_pow(int a, int n) {
  int ans = 1;
  while (n) {
    if (n & 1) {
      ans = 1ll * a * ans % mod;
    }
    a = 1ll * a * a % mod;
    n >>= 1;
  }
  return ans;
}

void get_r(int lim) {
  for (int i = 0; i < lim; i++) {
    r[i] = (i & 1) * (lim >> 1) + (r[i >> 1] >> 1);
  }
}

void NTT(int *f, int lim, int rev) {
  for (int i = 0; i < lim; i++) {
    if (i < r[i]) {
      swap(f[i], f[r[i]]);
    }
  }
  for (int mid = 1; mid < lim; mid <<= 1) {
    int wn = quick_pow(3, (mod - 1) / (mid << 1));
    for (int len = mid << 1, cur = 0; cur < lim; cur += len) {
      int w = 1;
      for (int k = 0; k < mid; k++, w = 1ll * w * wn % mod) {
        int x = f[cur + k], y = 1ll * w * f[cur + mid + k] % mod;
        f[cur + k] = (x + y) % mod, f[cur + mid + k] = (x - y + mod) % mod;
      }
    }
  }
  if (rev == -1) {
    int inv = quick_pow(lim, mod - 2);
    reverse(f + 1, f + lim);
    for (int i = 0; i < lim; i++) {
      f[i] = 1ll * f[i] * inv % mod;
    }
  }
}

void polyinv(int *f, int *g, int n) {
  if (n == 1) {
    g[0] = quick_pow(f[0], mod - 2);
    return ;
  }
  polyinv(f, g, n + 1 >> 1);
  for (int i = 0; i < n; i++) {
    t[i] = f[i];
  }
  int lim = 1;
  while (lim < 2 * n) {
    lim <<= 1;
  }
  get_r(lim);
  NTT(t, lim, 1);
  NTT(g, lim, 1);
  for (int i = 0; i < lim; i++) {
    int cur = (2 - 1ll * g[i] * t[i] % mod + mod) % mod;
    g[i] = 1ll * g[i] * cur % mod;
    t[i] = 0;
  }
  NTT(g, lim, -1);
  for (int i = n; i < lim; i++) {
    g[i] = 0;
  }
}

int fac[N], ifac[N], a[N], f[N], g[N], B[N], n;

void init() {
  fac[0] = 1;
  for (int i = 1; i < N; i++) {
    fac[i] = 1ll * fac[i - 1] * i % mod;
  }
  ifac[N - 1] = quick_pow(fac[N - 1], mod - 2);
  for (int i = N - 2; i >= 0; i--) {
    ifac[i] = 1ll * ifac[i + 1] * (i + 1) % mod;
  }
}

int main() {
  // freopen("in.txt", "r", stdin);
  // freopen("out.txt", "w", stdout);
  // ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
  init();
  scanf("%d", &n);
  for (int i = 0; i <= n; i++) {
    scanf("%d", &a[i]);
  }
  for (int i = 0; i <= n + 1; i++) {
    B[i] = ifac[i + 1];
  }
  polyinv(B, g, n + 2);
  for (int i = 0; i <= n; i++) {
    f[i] = 1ll * a[i] * fac[i] % mod;
  }
  reverse(g, g + n + 2);
  int lim = 1;
  while (lim < 2 * n + 10) {
    lim <<= 1;
  }
  get_r(lim);
  NTT(f, lim, 1), NTT(g, lim, 1);
  for (int i = 0; i < lim; i++) {
    f[i] = 1ll * f[i] * g[i] % mod;
  }
  NTT(f, lim, -1);
  for (int i = 1; i <= n + 1; i++) {
    a[i] = 1ll * ifac[i] * f[n + i] % mod;
  }
  a[0] = 0;
  n++;
  for (int i = 0; i < lim; i++) {
    f[i] = g[i] = 0;
  }
  for (int i = 0; i <= n; i++) {
    f[i] = 1ll * a[i] * fac[i] % mod;
    g[i] = ifac[n - i];
  }
  NTT(f, lim, 1), NTT(g, lim, 1);
  for (int i = 0; i < lim; i++) {
    f[i] = 1ll * f[i] * g[i] % mod;
  }
  NTT(f, lim, -1);
  for (int i = 0; i <= n; i++) {
    printf("%lld%c", 1ll * ifac[i] * f[n + i] % mod, i == n ? '\n' : ' ');
  }
  return 0;
}
_lifehappy_
关注
  • 5
    点赞
  • 踩
  • 3
    收藏
    觉得还不错? 一键收藏
  • 3
    评论
专栏目录
python 伯努利分布详解
01-20
伯努利分布 是一种离散分布,有两种可能的结果。1表示成功，出现的概率为p(其中0<p<1)。0表示失败，出现的概率为q=1-p。这种分布在人工智能里很有用，比如你问机器今天某飞机是否起飞了，它的回复就是Yes或No，...
伯努利数与自然数幂和
weixin_45771864的博客
03-21 511
伯努利数与自然数幂和
伯努利数详解
Deep_Kevin的娱乐空间
06-28 981
正题 ​ 我们来考虑 ∑i=0nik\sum_{i=0}^n i^k∑i=0n​ik 如何用与 kkk 相关的时间快速计算。 我们记 S(n,k)=∑i=0nikS(n,k)=\sum_{i=0}^n i^kS(n,k)=∑i=0n​ik 。 我们构造其关于 kkk 的 EGFEGFEGF ，则有： Gn(x)=∑k=0∑i=0nikxkk!=∑i=0n∑kikxkk!=∑i=0nexi \\G_n(x)=\sum_{k=0}\sum_{i=0}^n i^k\frac{x^k}{k!} \\=\sum_{i
伯努利数（Bernoulli number）
weixin_30347009的博客
03-30 368
设B0＝1，当k＞0时，定义 这些Bi(i＝0, 1,…, k)被称为伯努利数。按定义，自然得出：B1＝－，B2＝，B3＝0，B4＝－，B5＝0，B6＝，B7＝0，B8＝－，…。伯努利数是瑞士数学家雅各布·伯努利引入的数，出自于他的著作《猜度术》(1713)。除了B1外，当k为奇数时，Bk＝0；当k为偶数时，B2, B6, B10,…是正分数；B4, B8, B12,…是负分数。雅各布·...
伯努利数
alpha202的博客
02-13 714
参考文献资料： 小蒟蒻yyb的博客 百度百科-伯努利数 伯努利数(Bernoulli)——学习笔记 康复计划#3 简单常用的几种计算自然数幂和的方法 Bernoulli number has been found for about 200 years.Its name is named after Bernoulli, its discoverer. Recursive definition： ...
[数论]伯努利数
wulalalawu的博客
07-04 381
线性求逆元快速版 #include<bits/stdc++.h> #define MOD mod #define MAX maxn using namespace std; typedef pair<int, int> pii; typedef long long ll; const double eps = 1e-6; const int maxn = 2e3 + 10, inf = 0x3f3f3f3f, mod = 1e9 + 7; ll quickpow(ll x, ll k
论文研究 - 幂加法，伯努利数，伯努利多项式被重新考虑
05-22
利用翻译算子，我们得到了作用于单项式上的微分算子形式的算术级数和阶伯努利多项式的幂和。... 为了完成这项工作，我们提出了一些表格，用于以最简单的方式计算伯努利数，伯努利多项式，幂和以及幂和的Faulhaber公式。
伯努利数求k次方幂
03-01
这是一篇有关用伯努利数求k次方幂的文章，不管是伯努利数还是最终涉及的求k次方幂的方法统统介绍得非常详细，本来是在BDWK里有的，但是看着一堆点点真心很不舒服，于是下载后放在这里……
伯努利过程的熵。：这个简短的 MATLAB 脚本显示了伯努利过程的熵。-matlab开发
05-30
伯努利过程是一个简单的随机实验，只有两个结果。 抛硬币是一个典型的伯努利过程。 此 MATLAB 代码生成并绘制伯努利过程的熵。
伯努利数计算
weixin_34248118的博客
11-09 1733
采用递推公式 上面(M,K)表示组合数他跟C表示法上下标数是相反的，上面表示从m个数里取k个的组合   以前老教材的数学书上教的是P，但是后来统一改成A了。现在组合用的是C。  ...
伯努利数公式
weixin_30855099的博客
08-04 1302
伯努利数公式： 伯努利数满足条件，且有 那么继续得到 这就是伯努利数的递推式，逆元部分同样可以预处理。 转载于:https://www.cnblogs.com/Profish/p/11299125.html...
【C++心路历程23】求幂的和 分治算法
ctf109的博客
02-09 1670
【问题描述】 　　题目很简单：请你计算(a^1+a^2+…+a^n) mod 1234567 的结果，其中（0 < a,n < 2^31 ）。【输入格式】 第一行T，表示数据组数，接下来的T行，每行包含a和n，表示一组数据。【输出格式】 对于每组数据，输出对应的答案。【输入样例】 5 1 7 3 10 5 8 9 20 17 100【输出样例】 7 88572 488280
UyHiP 往期趣题整理（6）2017 / 05~07
misaka_worst的博客
04-07 379
因为往后只剩 6 个月的题，有些题有好几小问，所以分两篇。2017-05这题是 2007-12 和 2010-06 题目的延伸，这两道题就等我哪年整理到那儿再说……（1）设 n&gt;1 是整数，请问对怎样的 n，存在 n 元交换群 (S, +) 以及 S→S 的一一映射 A, B, C，使得 A(x)+B(x)=C(x) 恒成立？（2）设 n&gt;1 是整数，请问对怎样的 n，存在 n 元群 ...
预处理伯努利数(多项式求逆)
yamiedie_
09-22 2903
伯努利数定义：  递推式： 将e^t展开得到： 然后就可以利用多项式逆元求出伯努利数了。。 多项式逆元就是给定一个n次多项式A(x),求另一个多项式B(x)，满足 其中mod x^n表示多项式除以x^n剩余的多项式，这里可以看成做卷积之后去掉大于等于n次的项。 假设已经求得B(x)满足： 要求另一个B'(x)满足： 详情见： Inverse Element of Pol
贝尔数&&斯特灵数&&调和数&&伯努利数
acmer的随记
05-06 995
整数序列
设{Xn: n∈N}是伯努利过程, 则以下随机过程中不具有马尔可夫性的是 :A,{Xn: n∈N};B,{(Xn,Xn+1): n∈N};C,{X0+```+Xn,n∈N};D,{Xn+Xn+1: n∈N}
最新发布
05-12
只有选项 C 不具有马尔可夫性。 一个随机过程具有马尔可夫性，当且仅当对于任意时刻 $n$，该时刻的状态 $X_n$ 在给定过去的状态 $X_0,X_1,\ldots,X_{n-1}$ 下，与未来的状态 $X_{n+1},X_{n+2},\ldots$ 是独立的。也就是说，给定过去的状态，未来的状态只与当前状态有关，与过去的状态无关。 对于选项 A、B 和 D，它们的状态集合中都包括了 $X_n$ 的值，因此它们具有马尔可夫性。 但是对于选项 C，$X_0+\cdots+X_n$ 的值不仅取决于当前状态 $X_n$，还取决于过去的状态 $X_{n-1},X_{n-2},\ldots,X_0$。因此，选项 C 不具有马尔可夫性。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 3
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值