UVA 766 伯努利公式

UVA 766 伯努利公式

题意

$$S_k(n)=\sum _{i=1}^n{i}^k$$

找到最小的M，以及 $a_k,...a_1$ 使得上式变成下式
$$S_k(n)=\frac{1}{M}(a_{k+1}\ast n^{k+1}+a_k{n}^k+...+a_1\ast n+a_0)$$

题解

主要思路

该题就是套公式，公式理解都来自于这两个链接。

wiki上关于伯努利数的介绍

知乎上对伯努利数的推导

定义等幂和如下， 其中$m,n\ge$ 0:
$$S_k(n)=\sum _{i=1}^n{i}^k=1^k+2^k+...n^k$$
这数列和的公式必定是变量为n，次数为m+1次的多项式，称为伯努利多项式。

伯努利多项式的系数与 伯努利数$B_n$ 有密切关系如下：
$$S_k(n)=\frac{1}{k+1}\sum _{i=0}^k{C}_{k+1}^i{B}_i{n}^{k+1-i}$$
而伯努利$B_i$ 可以通过递推得到
$$B_m=-\frac{1}{m+1}\sum _{i=0}^{m-1}{C}_{m+1}^i{B}_i$$
代入 $B_0=1$ ，就可以得到一系列伯努利数了

算法过程

1. 代入公式，把系数用分数求出来
2. 再把这些分数形式的系数乘上他们分母的公倍数，就能让系数都为整数，并且$M$ 最小。

代码

#include <stdio.h>
#include <string.h>

long long gcd(long long a, long long b) {
    if (!b) return a;
    return gcd(b, a % b);
}

long long lcm(long long a, long long b) {
    a = a / gcd(a, b) * b;
    if (a < 0) a = -a;
    return a;
}

struct Fraction {
    long long a, b;
    Fraction() {a = 0; b = 1;}

    Fraction(long long x) {
        a = x; b = 1;
    }

    Fraction(long long x, long long y) {
        a = x; b = y;
    }

    void deal() {
        if (b < 0) {b = -b; a = -a;}
        long long k = gcd(a, b);
        if (k < 0) k = -k;
        a /= k; b /= k;
    }

    Fraction operator+(Fraction p) {
        Fraction ans;
        ans.b = lcm(b, p.b);
        ans.a = ans.b / b * a + ans.b / p.b * p.a;
        ans.deal();
        return ans;
    }

    Fraction operator-(Fraction p) {
        Fraction ans;
        ans.b = lcm(b, p.b);
        ans.a = ans.b / b * a - ans.b / p.b * p.a;
        ans.deal();
        return ans;
    }

    Fraction operator*(Fraction p) {
        Fraction ans;
        ans.a = a * p.a;
        ans.b = b * p.b;
        ans.deal();
        return ans;
    }

    Fraction operator/(Fraction p) {
        Fraction ans;
        ans.a = a * p.b;
        ans.b = b * p.a;
        ans.deal();
        return ans;
    }

    void operator=(int x) {
        a = x;
        b = 1;
    }

    void print() {
        printf("%lld/%lld\n", a, b);
    }
};

const int N = 25;

Fraction B[N], C[N][N], a[N];
long long L;

void init() {
    for (int i = 0; i < N; i++) {
        C[i][0] = 1;
        C[i][i] = 1;
        for (int j = 1; j < i; j++) {
            C[i][j] = C[i - 1][j - 1] + C[i - 1][j];
        }
    }

    B[0] = 1;
    for (int i = 1; i <= 20; i++) {
        B[i] = 0;
        for (int j = 0; j < i; j++) B[i] = B[i] - C[i + 1][j] * B[j];
        B[i] = B[i] / (i+1);
    }
}

int main() {
//    freopen("in.text","r",stdin);

    int cas,k;
    init();
    scanf("%d", &cas);
    while (cas--) {
        L = 1;
        scanf("%d", &k);
        for (int i = 0; i <= k; i++){
            a[i] = C[k+ 1][i] * B[i] * Fraction(1, k+ 1);
            L = lcm(L, a[i].b);
        }
        printf("%lld ", L);
        a[1] = a[1] + 1;
        for (int i = 0; i <=k; i++)printf("%lld ", L / a[i].b * a[i].a);
        printf("0\n");
        if (cas) printf("\n");
    }
    return 0;
}