1822 序列求和V5

于 2021-08-19 11:42:24 发布
阅读量337 收藏
点赞数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/t56346/article/details/119798584
版权

若R≡0R≡0，那么显然答案就为0.
　　接下来，我们用符号fk(n)fk(n)代表我们所要求的值。首先我们要预处理fk(1),⋯,fk(k+1)fk(1),⋯,fk(k+1)的值，由于xkxk是一个完全积性函数，利用线性筛，这一预处理可以很容易在O(k)O(k)的时间内完成。
　　若R≡1R≡1，则原问题转化为求∑ni=1ik∑i=1nik的值，我们知道，这时fk(n)fk(n)是关于nn的k+1k+1次多项式。同时，我们也知道，n≥k+1n≥k+1时，每一个CinCni都是关于nn的ii次多项式，并且它们线性无关，而n≤kn≤k时，可以O(1)O(1)直接查询。对于n≥k+1n≥k+1的情况，我们用C0n,⋯,Ck+1nCn0,⋯,Cnk+1这k+2k+2个数来线性表出fk(n)fk(n)，设fk(n)=∑k+1i=0Cinaifk(n)=∑i=0k+1Cniai，根据二项式反演，有an=∑ni=0(−1)n−iCinfk(i)an=∑i=0n(−1)n−iCnifk(i)。
　　这里先写一个引理：∑ki=0(−1)iCin=(−1)kCkn−1∑i=0k(−1)iCni=(−1)kCn−1k，这个很好证我就不写了。
　　将anan代回原式得：
fk(n)=∑i=0k+1Cin∑j=0i(−1)i−jCjifk(j)=∑i=0k+1fk(i)∑j=ik+1(−1)j−iCjnCij=∑i=0k+1fk(i)Cin∑j=ik+1(−1)j−iCj−in−i=∑i=0k+1fk(i)Cin∑j=0k−i+1(−1)jCjn−i=∑i=0k+1(−1)k−i+1CinCk−i+1n−i−1=∑i=0k+1∏0≤j≤k+1,j≠in−ji!(k−i+1)!
fk(n)=∑i=0k+1Cni∑j=0i(−1)i−jCijfk(j)=∑i=0k+1fk(i)∑j=ik+1(−1)j−iCnjCji=∑i=0k+1fk(i)Cni∑j=ik+1(−1)j−iCn−ij−i=∑i=0k+1fk(i)Cni∑j=0k−i+1(−1)jCn−ij=∑i=0k+1(−1)k−i+1CniCn−i−1k−i+1=∑i=0k+1∏0≤j≤k+1,j≠in−ji!(k−i+1)!
　　这样，分子和分母的部分都可以预处理出来，复杂度O(k)O(k)。
　　若R≢1R≢1，有：
(R−1)fk(n)∴fk(n)=(R−1)∑i=1nikRi=∑i=1nikRi+1−∑i=1nikRi=∑i=1n+1(i−1)kRi−∑i=1nikRi=nkRn+1−∑i=1n∑j=0k−1(−1)k−j−1CjkijRi=nkRn+1−∑j=0k−1(−1)k−j+1Cjkfj(n)=nkRn+1−∑k−1j=0(−1)k−j+1Cjkfj(n)R−1
(R−1)fk(n)=(R−1)∑i=1nikRi=∑i=1nikRi+1−∑i=1nikRi=∑i=1n+1(i−1)kRi−∑i=1nikRi=nkRn+1−∑i=1n∑j=0k−1(−1)k−j−1CkjijRi=nkRn+1−∑j=0k−1(−1)k−j+1Ckjfj(n)∴fk(n)=nkRn+1−∑j=0k−1(−1)k−j+1Ckjfj(n)R−1
　　下面我们证明一个结论，对∀fi(n),0≤i≤k,∃Fi(n)∈Qi[x]∀fi(n),0≤i≤k,∃Fi(n)∈Qi[x]，有fi(n)=Rn+1Fi(n)−RFi(0)fi(n)=Rn+1Fi(n)−RFi(0)。
证明：当i当ifk(n)这样，=0时,结论已经成立(请自行验证)=k时，假设结论对0,⋯,k−1均成立=nkRn+1−∑k−1j=0(−1)k−j+1Cjkfj(n)R−1=nkRn+1−∑k−1j=0(−1)k−j+1Cjk[Rn+1Fj(n)−RFj(0)]R−1=Rn+1[nk−∑k−1j=0(−1)k−j+1CkjFj(n)]+R[∑k−1j=0(−1)k−j+1CjkFj(0)]R−1我们令Fk(n)=nk−∑j=0k−1(−1)k−j+1CkjFj(n)，容易验证它满足题目要求。
证明：当i=0时,结论已经成立(请自行验证)当i=k时，假设结论对0,⋯,k−1均成立fk(n)=nkRn+1−∑j=0k−1(−1)k−j+1Ckjfj(n)R−1=nkRn+1−∑j=0k−1(−1)k−j+1Ckj[Rn+1Fj(n)−RFj(0)]R−1=Rn+1[nk−∑j=0k−1(−1)k−j+1CjkFj(n)]+R[∑j=0k−1(−1)k−j+1CkjFj(0)]R−1这样，我们令Fk(n)=nk−∑j=0k−1(−1)k−j+1CjkFj(n)，容易验证它满足题目要求。

　　我们再来证明一个引理：对 ∀0≤m≤k,∑ki=0(−1)iCikCmi=0∀0≤m≤k,∑i=0k(−1)iCkiCim=0
证明：∑i=0k(−1)iCikCmi=∑i=0k(−1)iCmkCi−mk−m=Cmk∑i=0k−m(−1)i+mCik−m=Cmk(−1)m(1−1)(k−m)=0
证明：∑i=0k(−1)iCkiCim=∑i=0k(−1)iCkmCk−mi−m=Ckm∑i=0k−m(−1)i+mCk−mi=Ckm(−1)m(1−1)(k−m)=0

　　在上一部分的讨论中，我们已经知道，作为一个 kk次多项式，Fk(n)Fk(n)可以表示为 C0n,⋯,CknCn0,⋯,Cnk的线性组合。设 Fk(n)=∑ki=0CinaiFk(n)=∑i=0kCniai，则
∑i=0k(−1)iCikFk(n)=∑i=0k(−1)iCik∑j=0kCjiaj=∑j=0k∑i=0k(−1)iCikCjiaj=∑j=0k0=0
∑i=0k(−1)iCkiFk(n)=∑i=0k(−1)iCki∑j=0kCijaj=∑j=0k∑i=0k(−1)iCkiCijaj=∑j=0k0=0

　　我们又有：
(n+1)kRn+1∴Fk(n+1)=fk(n+1)−fk(n)=Rn+2Fk(n+1)−Rn+1+RFk(0)=(n+1)k+Fk(n)R
(n+1)kRn+1=fk(n+1)−fk(n)=Rn+2Fk(n+1)−Rn+1+RFk(0)∴Fk(n+1)=(n+1)k+Fk(n)R

　　通过上面两个等式，我们就可以通过解一个一元一次方程来得到 Fk(0),⋯,Fk(k)Fk(0),⋯,Fk(k)的值，之后我们用拉格朗日插值公式就可以 O(k)O(k)地得到 Fk(n)Fk(n)的值，再代回原式，这道题就解决了。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <ctime>
#include <vector>
#include <set>
//#define ivorysi
#define eps 1e-8
#define mo 974711
#define pb push_back
#define mp make_pair
#define pii pair<int,int>
#define fi first
#define se second
#define MAXN 200005
#define space putchar(' ')
#define enter putchar('\n')
using namespace std;
typedef long long int64;
typedef unsigned int u32;
typedef unsigned long long u64;
typedef double db;
template<class T>
void read(T &res) {
    res = 0;char c = getchar();T f = 1;
    while(c < '0' || c > '9') {
    if(c == '-') f = -1;
    c = getchar();
    }
    while(c >= '0' && c <= '9') {
    res = res * 10 + c - '0';
    c = getchar();
    }
    res *= f;
}
template<class T>
void out(T x) {
    if(x < 0) putchar('-');
    if(x >= 10) {
    out(x / 10);
    }
    putchar('0' + x % 10);
}
const int64 MOD = 985661441;
int64 fpow(int64 x,int64 c) {
    int64 res = 1,t = x;
    while(c) {
    if(c & 1) res = res * t % MOD;
    t = t * t % MOD;
    c >>= 1;
    }
    return res;
}

int T,K;
int prime[MAXN],tot;
int64 N,R,F[MAXN];
int64 a[MAXN],b[MAXN];
int64 MK[MAXN];
int64 fac[MAXN],invfac[MAXN],Le[MAXN],Ri[MAXN];
bool nonprime[MAXN];
int64 C(int n,int m) {
    if(n < m) return 0;
    return fac[n] * invfac[n - m] % MOD * invfac[m] % MOD;
}
void Solve() {
    read(K);read(N);read(R);
    R %= MOD;
    if(R == 0) {puts("0");return;}   
    MK[1] = 1;
    memset(nonprime,0,sizeof(nonprime));
    tot = 0;
    for(int i = 2 ; i <= K + 2 ; ++i) {
    if(!nonprime[i]) {
        prime[++tot] = i;
        MK[i] = fpow(i,K); 
    }
    for(int j = 1 ; j <= tot ; ++j) {
        if(prime[j] > (K + 2) / i) break;
        nonprime[i * prime[j]] = 1;
        MK[i * prime[j]] = MK[i] * MK[prime[j]] % MOD;
        if(i % prime[j] == 0) break;
    }
    }
    if(R == 1) {
    F[0] = 0;
    for(int i = 1 ; i <= K + 2 ; ++i) F[i] = (F[i - 1] + MK[i]) % MOD;
    if(N <= K + 2) {out(F[N]);enter;return;}
    int64 t = 1,ans = 0;
    Le[0] = 1;
    N %= MOD;
    for(int i = 1 ; i <= K + 2 ; ++i) {
        Le[i] = Le[i - 1] * (N + MOD - i) % MOD;
    }
    Ri[K + 3] = 1;
    for(int i = K + 2 ; i >= 1 ; --i) {
        Ri[i] = Ri[i + 1] * (N + MOD - i) % MOD;
    }
    for(int i = K + 2 ; i >= 1 ; --i) {
        ans += t * invfac[i - 1] % MOD * invfac[K + 2 - i] % MOD * Le[i - 1] % MOD * Ri[i + 1] % MOD * F[i] % MOD;
        t = t * (MOD - 1) % MOD;
    }
    ans %= MOD;
    out(ans);enter;return;
    }
    a[0] = 1,b[0] = 0;
    int64 InvR = fpow(R,MOD - 2);
    for(int i = 1 ; i <= K + 1; ++i) {
    a[i] = 1LL * a[i - 1] * InvR % MOD;
    b[i] = (1LL * b[i - 1] * InvR + MK[i]) % MOD;
    }
    int64 suma = 0,sumb = 0,t = 1;
    if(K & 1) t = MOD - 1;
    for(int i = 0 ; i <= K + 1; ++i) {
    int64 h = t * C(K + 1,i) % MOD;
    suma = (suma + a[i] * h) % MOD;
    sumb = (sumb + b[i] * h) % MOD;
    t = t * (MOD - 1) % MOD;
    }
    F[0] = (MOD - sumb) * fpow(suma,MOD - 2) % MOD;
    for(int i = 1 ; i <= K + 1; ++i) F[i] = (a[i] * F[0] + b[i]) % MOD;
    int64 ans = 0;
    t = 1;
    if(N <= K) {
    ans = (fpow(R,N) * F[N] - F[0] + MOD) % MOD;
    out(ans);enter;
    return ;
    }
    int64 T = N % (MOD - 1);
    N %= MOD;
    Le[0] = N;
    for(int i = 1 ; i <= K ; ++i) {
    Le[i] = Le[i - 1] * (N + MOD - i) % MOD;
    }
    Ri[K + 1] = 1;
    for(int i = K ; i >= 0 ; --i) {
    Ri[i] = Ri[i + 1] * (N + MOD - i) % MOD;
    }
    if(K & 1) t = MOD - 1;
    for(int i = 0 ; i <= K; ++i) {
    int64 h;
    if(i == 0) h = t * invfac[i] % MOD * invfac[K - i] % MOD * Ri[i + 1] % MOD * F[i] % MOD;
    else h = t * invfac[i] % MOD * invfac[K - i] % MOD * Le[i - 1] % MOD * Ri[i + 1] % MOD * F[i] % MOD;
    t = t * (MOD - 1) % MOD;
    ans += h;
    }
    ans %= MOD;
    ans = (ans * fpow(R,T) - F[0] + MOD) % MOD;
    out(ans);enter;
}
int main() {
#ifdef ivorysi
    freopen("f1.in","r",stdin);
#endif
    fac[0] = 1;
    for(int i = 1 ; i <= 200001 ; ++i) fac[i] = fac[i - 1] * i % MOD;
    invfac[200001] = fpow(fac[200001],MOD - 2);
    for(int i = 200000 ; i >= 0 ; --i) invfac[i] = invfac[i + 1] * (i + 1) % MOD;
    read(T);
    while(T--) {
    Solve();
    }
    return 0;
}

林子盛吧(✪ω✪)
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
Python实现分数序列求和
09-17
在本文中，我们将探讨如何使用Python来实现分数序列的求和。分数序列通常是指一系列分数，其中每个分数都是由两个正整数分子和分母组成的。在这个特定的问题中，我们面临的是一个特定的分数序列，它遵循斐波那契数列...
【51nod1822】序列求和 V5（拉格朗日插值）（组合数学推导）
zxyoi_dreamer的博客（不定期诈尸）
04-05 938
传送门 解析： 如果r≡0r\equiv 0r≡0显然答案就是0，如果r≡1r\equiv 1r≡1，则我们要求的就是自然数的幂和，可以看一下我的博客：自然数求幂和的几种方法。 现在考虑r̸≡1r\not \equiv 1r̸​≡1的情况。 设我们要求的是sk(n)=∑i=1nikris_k(n)=\sum\limits_{i=1}^{n}i^kr^isk​(n)=i=1∑n​ikri 利用序列...
nk 夹缝中求和 【二分】
ln2037的博客
12-13 260
链接：https://ac.nowcoder.com/acm/problem/213483 来源：牛客网 题目描述 给定一个数组a，以及两个正整数 和 ，求取两个数 和 aja_jaj​ ，满足的取法有多少种？ 注：只要两个取法有一个角标不同，则视为两种不同的取法。 输入描述: 第一行有三个正整数，，，分别代表数组长度，和的下界和上界。 第二行有个正整数 输出描述: 一个整数，代表最后的取法数量。 示例1 输入 复制 5 2 4 1 4 2 2 1 输出 复制 6 说明 共有以下6种取法： i=1，j=3
序列求和（gcd与逆元）
DoIdo
05-24 344
注释：Cii+Ci+1i+...+Cni=Cn+1i+1\small C_i^i+C_{i+1}^i+...+C_n^i=C_{n+1}^{i+1}Cii​+Ci+1i​+...+Cni​=Cn+1i+1​ 题面 题意：见题面。 解决思路：结论：12+22+32+...+n2=n(n+1)(2n+1)6\small 1^2+2^2+3^2+...+n^2=\frac{n(n+1)(2n+1)}{6}12+22+32+...+n2=6n(n+1)(2n+1)​ 证明 ∵i2=i2+i−i=i(i−1)+i..
常用组合数计算公式
shadandeajian的博客
08-26 6万+
Cmn=Cm−1n−1+Cmn−1Cnm=Cn−1m−1+Cn−1mC_n^m=C_{n-1}^{m-1}+C_{n-1}^{m} m∗Cmn=n∗Cm−1n−1m∗Cnm=n∗Cn−1m−1m*C_{n}^{m}=n*C_{n-1}^{m-1} C0n+C1n+C2n+⋯+Cnn=2nCn0+Cn1+Cn2+⋯+Cnn=2nC_{n}^{0}+C_{n}^{1}+C_{n}^{2}+\dots...
【numpy求和】numpy.sum()求和
热门推荐
TomRen
11-17 13万+
numpy.sum(a, axis=None, dtype=None, out=None, keepdims=, initial=)[source])用于计算array元素的和. python中常用的numpy进行数学计算，其中array的求和运算分为两种，一种是调用numpy.array自身的sum()方法,另一种是利用numpy的内建函数numpy.sum()使用。(tips:python...
数列求和各种方法总结归纳.ppt
08-17
数列求和
斐波那契数列求和_whyadm_斐波那契求和c_数列求和_poorbv2_
10-03
斐波那契数列是一个非常经典的数学序列，它在计算机科学和算法设计中有着广泛的应用。这个数列的定义是这样的：第一项F0为0，第二项F1为1，之后每一项Fi（i >= 2）都是前两项的和，即Fi = Fi-1 + Fi-2。数列的初始...
c语言数列求和c语言数列求和
11-18
等差数列是指序列中的每一项与前一项之间的差是恒定的。比如1, 3, 5, 7...，差值为2。C语言中求等差数列的公式可以表示为：`S = n/2 * (a1 + an)`，其中`S`是求和结果，`n`是项数，`a1`是首项，`an`是末项。可以...
数列求和-加强版 （20 分）
10-21
7-38 数列求和-加强版 （20 分） 给定某数字A（1≤A≤9）以及非负整数N（0≤N≤100000），求数列之和S=A+AA+AAA+⋯+AA⋯A（N个A）。例如A=1, N=3时，S=1+11+111=123。 输入格式： 输入数字A与非负整数N。 输出格式...
hdu5226 Tom and matrix 公式，Lucas
Every day
10-09 614
题意：给定x1,y1,x2,y2,求和C(i,j)，(x1 分析：因为\sum_{i=a}^{b}C_{i}^{k}=C_{b+1}^{k+1}-C_{a}^{k+1}∑​i=a​b​​C​i​k​​=C​b+1​k+1​​−C​a​k+1​​ 所以求同一列的数的和可以变成求两个组合数的差。由于p可能很小，当除数为p的倍数时就为0了，直接乘逆元会出问题，利用Lucas。 #include #i
自然数幂求和——第二类Strling数
A1847225889的博客
01-28 518
这个问题似乎有很多种求法，但感觉上第二类Strling数的做法是最方便的。 问题 求下面这个式子： ∑i=0nik\sum_{i=0}^n i^ki=0∑n​ik nnn的范围可以很大。 第二类Strling数 第二类Strling数记作S(n,m)S(n,m)S(n,m)、SnmS_n^mSnm​。 定义：将nnn个相同的球放在mmm个不同的箱子里的方案数（其中的每一个箱子至少有一个球）。 ...
二项式反演
jkchen's Haven
09-18 1799
反演： 设求解函数为G，当G直接求比较复杂时，可以用与G相关的函数F来转化，且一般F都比较好求，大多时候为一个未知量已知的且可以O(1)求出的公式 比如，G(d)=∑i=1n∑j=1m[gcd(i,j)==d]G(d)=\sum_{i=1}^n\sum_{j=1}^m[gcd(i,j)==d]G(d)=∑i=1n​∑j=1m​[gcd(i,j)==d]，就可以用莫比乌斯反演快速推出：F(d)=∑d...
几个重要的排列组合定理公式
不忘初心
02-01 4万+
1.排列的几个定理公式 .排列，一般地，从n个不同元素中取出m（m≤n）个元素，按照一定的顺序排成一列，叫做从n个元素中取出m个元素的一个排列(Arrangement)。特别地，当m=n时，这个排列被称作全排列(Permutation)。 .n个元素的循环r-排列的个数为上式除以r 循环，顾名思义，就是围成一圈，规定一个方向(顺或逆)，转一圈形成的为一种。如：若不围一圈，1234
组合数计算（从1000到1e9的组合数各类求法）
AEI的博客
09-29 2305
组合数涉及的问题很多，求解过程中涉及的定理判定也有很多，期间断断续续参考了很多大佬们的博客，最终进行了总结。在此我使用用问题引入算法的方法来一步一步了解组合数。 简单的组合数问题（运用快速幂及递归实现杨辉三角） 组合数求解 像关于组合数求和问题如∑ Cni（0<=i<=n） 可利用快速幂，此时∑ Cni = 2^n 此为二项式定理推导出的结果 练习题可见洛谷P3414 若求较小的组合数如C(n,m)（m<=n<=1000 保守估计1000，实际可能可以更大些，但不会超数量级
证明二项式系数平方和等于组合数C(2n,n)
qq_22125259的博客
07-02 7407
组合恒等式2 五个基本的组合恒等式 更复杂的技巧与例题
一个不愿透露姓名的博客
08-14 2248
组合恒等式2 五个基本的组合恒等式 复杂的例题使用两个数列的递推关系证明恒等式证明相反数与自身相等以说明恒等式为0组合数的倒数的裂项数学归纳法 这一讲继续讨论用五个基本的组合恒等式证明组合恒等式的方法。先复习一下基本恒等式： 等式一 组合数的对称性 Cnk=Cnn−kC_n^k = C_n^{n-k}Cnk​=Cnn−k​ 等式二 杨辉三角 Cnk=Cn−1k+Cn−1k−1C_n^k = C_{n-1}^k + C_{n-1}^{k-1}Cnk​=Cn−1k​+Cn−1k−1​ 等式三 多项式系数的
组合数有关的公式及常用求和
bigtiao097的博客
08-16 5万+
组合数
matlab数列求和
最新发布
06-18
在MATLAB中，求和数列是一个常见的数学运算，你可以使用内置函数`sum()`来实现。这个函数可以对向量、矩阵甚至复数数组求和。基本语法如下： ```matlab sum(vector) ``` 其中`vector`是你想要求和的数列。如果你想...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值