伯努利数与自然数幂和
今天我们讨论的问题是如何有效地求自然数的幂和。接下来以3个经典题目为例来讲解。
题目:http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=1864
分析:其实求自然数的幂和方法有很多种,先来看看普通的递推求法,由于
那么对于所有的累加得到
进一步得到
可以看出这是一个递推式,如果我们记
那么得到如下递归式
递归出口是
为了提高效率,在递归的时候需要记忆化。由于要用到大数,用Java实现。
代码:
import java.math.*;
import java.util.*;public class Main {public static final int N = 105;public static final BigInteger FLAG = (BigInteger.ZERO).subtract(BigInteger.ONE);public static BigInteger[][] C = new BigInteger[N][N];public static BigInteger[] ans = new BigInteger[N];public static void Init(){for(int i=0; i<N; i++){C[i][0] = C[i][i] = BigInteger.ONE;if(i == 0) continue;for(int j=1; j<i; j++)C[i][j] = C[i-1][j].add(C[i-1][j-1]);}}public static BigInteger Solve(BigInteger n, int k){if(ans[k].compareTo(FLAG) != 0){return ans[k];}if(k == 1){ans[k] = ((n.add(BigInteger.ONE)).multiply(n)).divide(BigInteger.valueOf(2));return ans[k];}BigInteger tmp = BigInteger.ONE;for(int i=0; i<k+1; i++){tmp = tmp.multiply(n.add(BigInteger.ONE));}tmp = tmp.subtract(n.add(BigInteger.ONE));BigInteger sum = BigInteger.ZERO;for(int i=1; i<k; i++){BigInteger t = C[k+1][i+1].multiply(Solve(n, k-i));sum = sum.add(t);}ans[k] = (tmp.subtract(sum)).divide(BigInteger.valueOf(k+1));return ans[k];}public static void main(String[] args){Init();Scanner cin = new Scanner(System.in);while(cin.hasNext()){BigInteger n = cin.nextBigInteger();int k = cin.nextInt();for(int i=0; i<N; i++){ans[i] = FLAG;}System.out.println(Solve(n, k));}}
}
题目:http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1228
分析:本题题意就是求自然数的幂和,但是它的case比较多。对于求幂和本身就需要的时间复杂度,如果继
续用上述方法来求自然数的幂和,5000个case会TLE,接下来介绍另一个求自然数幂和的方法,它是基于伯
努利数的,公式描述如下
可以看出只要我们预处理出每一项,就可以在线性时间内求得自然数的幂和。前面的倒数可以用递推法求逆元
预处理,组合数也可以预处理,也可以先预处理,现在关键是如何预处理伯努利数。
伯努利数满足条件,且有
那么继续得到
这就是伯努利数的递推式,逆元部分同样可以预处理。
代码:
#include <iostream>
#include <string.h>
#include <stdio.h>using namespace std;
typedef long long LL;
const LL MOD = 1000000007;
const int N = 2005;LL C[N][N];
LL B[N],Inv[N];
LL Tmp[N];
LL 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)
{LL ans = Inv[k+1];LL sum = 0;for(int i=1; i<=k+1; i++){sum += C[k+1][i] * Tmp[i] % MOD * B[k+1-i] % MOD;sum %= MOD;}ans *= sum;ans %= MOD;return ans;
}int main()
{int T;Init();scanf("%d", &T);while(T--){int k;scanf("%I64d %d", &n, &k);n %= MOD;Tmp[0] = 1;for(int i=1; i<N; i++)Tmp[i] = Tmp[i-1] * (n + 1) % MOD;printf("%I64d\n", Work(k));}return 0;
}
题目:http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1258
分析:本题与上题不同的是值比较大,达到50000,如果采用同样的方法,会TLE的。那么必定要进行优化。
参考Pick大神的方法,如下
链接:http://picks.logdown.com/posts/189620-the-inverse-element-of-polynomial
伯努利数与自然数幂和相关推荐
- 【学习笔记】自然数幂和
温馨提示: 本文文档大小约\(11KB\). 引入 自然数幂和是一个我们从小就耳熟能详的经典问题.定义\(S(n,k)=\sum^{n}_{i=0} i^k\), 显然\(S(n,k)\)为关于\(n ...
- 自然数幂和伯努利数(Bernoulli)
二项式定理求自然数幂和 由二项式定理展开得 \[ (n+1)^{k+1}-n^{k+1}=\binom {k+1}1n^k+\binom {k+1}2n^{k-1}+\cdots+\binom {k+ ...
- 自然数幂和取模问题进一步探究
在上次的文章 http://blog.csdn.net/acdreamers/article/details/38929067 中,学习了求自然数幂 和的有效方法,并且了解了求伯努利数的优美算法.今天 ...
- Easy Summation 假的自然数幂的和
easy summation 题意:就是求自然数幂的和.注意平常求自然数幂的和当然是用伯努利数,但是本题n很大,k很小,用伯努利数来做会内存超限. 用快速幂暴力求即可. #include<bit ...
- 自然数幂和 斯特林数
2018 UPD: 其实第二类斯特林数做自然数幂和更简单,这里简单写一下: 由一个基本式子出发 nk=∑i=0k{k i}[n]i n^k=\sum_{i=0}^k\big\{^k_{\ i}\big ...
- 专题——自然数幂和 详解
自然数幂和,有许多种解决方法,其中最典型的有下面几种. Description 求 ∑ i = 1 n i k \sum_{i=1}^n i^k i=1∑nik 由于答案可能过大,请将其对 1 0 ...
- 自然数幂求和方法1:扰动法(求两次)
自然数幂求和方法1:扰动法(求两次) 先来搞一搞等比数列 标号从1开始,\(a_n=a1*q^{n-1}\) \(S_n=\sum_{k=1}^n a_k\) \[\begin{aligned}S_n ...
- jzoj5224 [GDOI2018模拟7.12]C 杜教筛+自然数幂和
Description 求 ∑i=1n∑j=1ngcd(i,j)k且n≤1010,k≤5 ∑ i = 1 n ∑ j = 1 n gcd ( i , j ) k 且 n ≤ 10 10 , k ≤ 5 ...
- 单片机c语言如何实现幂方,自然数幂方和二项系数表示的系数公式.pdf
维普资讯 山西师范大学学报 (自然科学版) 第22卷第3期 JournalofShanxiNormalUniversity V01.22 No.3 2008年9月 NaturalScienceEdit ...
最新文章
- Oracle 计算两个日期间隔的天数、月数和年数
- script标签中的async和defer
- Linux open函数使用方法记录
- python问题:NameError: name 'reload' is not defined
- java arraylist 字符串数组_java中String,数组,ArrayList三者之间的转换
- JQuery实现步骤条效果(上一步,下一步)
- detectron2 学习笔记
- springboot+党员信息管理系统 毕业设计-附源码161528
- 最低成本的ARM调试解决方案——有关于Wiggler、H-Jtag、OpenOCD、GDB
- python上传文件到onedrive_python-onedrive使用教程【linux备份至onedrive】 | C/C++程序员之家...
- 微信公众号获取微信用户共享收货地址
- 此实现不是 Windows 平台 FIPS 验证的加密算法的一部
- 采药问题 c语言程序,采药-动态规划题解(C语言代码)
- 9.9实例操作:MVC设计模式应用遇到的问题和解决办法
- Linux磁盘管理练习题(习题+演示)
- 打破界限,获得新生——观电影《云图》有感
- 智能电灯代码_智能电灯开关与智能电灯泡:您应该购买哪个?
- 一元线性回归及Excel回归分析
- python编程用什么软件?
- JDK1.8 or JDK8傻傻分不清楚