\(k\leq 200000\)

考虑转化成枚举 \(k\) 的形式

我们错位相减!

\[A_k=\sum_{i=1}^N i^K\times R^i
\\
RA_k=\sum_{i=2}^{N+1} (i-1)^KR^i
\\
(R-1)A_k=N^kR^{N+1}+\sum_{i=1}^{N}[(i-1)^k-i^k]R^i
\]

二项式展开!

\[(R-1)A_k=N^kR^{N+1}+\sum_{i=1}^{N}[\sum_{j=0}^k(-1)^{k-j}i^{j}-i^k]R^i
\\
=N^kR^{N+1}+\sum_{i=1}^{N}\sum_{j=0}^{k-1}\binom{k}{j}(-1)^{k-j}i^{j}R^i
\\
=N^kR^{N+1}+\sum_{j=0}^{k-1}\binom{k}{j}(-1)^{k-j}\sum_{i=1}^{N}i^{j}R^i
\\
=N^kR^{N+1}+\sum_{j=0}^{k-1}\binom{k}{j}(-1)^{k-j}A_j
\]

这个时候我们递归计算可以得到 \(O(k^2)\) 的复杂度

考虑能否插值出一个 \(k\) 次多项式?

先暴力展开第零项

\[A_0=\sum_{i=1}^N \times R^i
\]

等比数列求和公式套一下

\[S_n=\frac{a_1(1-q^n)}{1-q}
\\
A_0=\frac{R(1-R^n)}{1-R}=\frac{R^{n+1}-R}{R-1}=R^{N + 1}\cdot (\frac{1}{R - 1}) - R\cdot (\frac{1}{R - 1})
\]

暴力展开第一项

\[A_1=\frac{NR^{N+1}-A_0}{R-1}
\\
A_1=\frac{NR^{N+1}-\frac{R^{n+1}-R}{R-1}}{R-1}
\\
A_1=\frac{NR^{N+1}}{R-1}-\frac{\frac{R^{n+1}-R}{R-1}}{R-1}
\\
A_1=\frac{NR^{N+1}}{R-1}-\frac{R^{n+1}-R}{(R-1)^2}
\\
A_1=R^{N + 1}(\frac{N}{R-1}-\frac{1}{(R-1)^2}) + R\times \frac{1}{(R-1)^2}
\]

我们的 \(A_k\) 可以插值,次数可以保证,我们直接把前 \(O(k)\) 项求出来插值,线性拉插一下就好了

[51nod 1822]序列求和的更多相关文章

  1. 51nod 1258 序列求和 V4

    http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1258 1258 序列求和 V4  基准时间限制:8 秒 空间限制:131 ...

  2. 51NOD 1258 序列求和 V4 [任意模数fft 多项式求逆元 伯努利数]

    1258 序列求和 V4 题意:求\(S_m(n) = \sum_{i=1}^n i^m \mod 10^9+7\),多组数据,\(T \le 500, n \le 10^{18}, k \le 50 ...

  3. 51nod 1228 序列求和(伯努利数)

    1228 序列求和  题目来源: HackerRank 基准时间限制:3 秒 空间限制:131072 KB 分值: 160 难度:6级算法题  收藏  关注 T(n) = n^k,S(n) = T(1 ...

  4. 【51nod】1822 序列求和 V5

    题解 我是zz吧 nonprime[i * prime[j]] = 0 = = 还以为是要卡常,卡了半天就是过不掉 我们来说这道题-- 首先,我们考虑一个\(K^2\)做法 \(f_{k}(N) = ...

  5. 51nod 1228 序列求和 ( 1^k+2^k+3^k+...+n^k )

    C为组合数,B为伯努利数 具体推到过程略 参考博客:http://blog.csdn.net/acdreamers/article/details/38929067# (我的式子和博客中的不一样,不过 ...

  6. 51Nod - 1228 序列求和 (自然数幂和+伯努利数)

    https://vjudge.net/problem/51Nod-1228 Description T(n) = n^k,S(n) = T(1) + T(2) + ...... T(n).给出n和k, ...

  7. 51Nod 1228 序列求和

    T(n) = n^k,S(n) = T(1) + T(2) + ...... T(n).给出n和k,求S(n).   例如k = 2,n = 5,S(n) = 1^2 + 2^2 + 3^2 + 4^ ...

  8. [51nod]1229 序列求和 V2(数学+拉格朗日差值)

    题面 传送门 题解 这种颓柿子的题我可能死活做不出来-- 首先\(r=0\)--算了不说了,\(r=1\)就是个裸的自然数幂次和直接爱怎么搞怎么搞了,所以以下都假设\(r>1\) 设 \[s_p ...

  9. 【51Nod1258】序列求和V4(FFT)

    [51Nod1258]序列求和V4(FFT) 题面 51Nod 多组数据,求: \[Ans=\sum_{i=1}^ni^k,n\le 10^{18},k\le50000\] 题解 预处理伯努利数,时间 ...

随机推荐

  1. 142_Power BI之同比预测

    博客:www.jiaopengzi.com 焦棚子的文章目录 请点击下载附件 一.背景 最近刚好在做一个简单同比预测的模型,预测方法很简单,就是累计同比预测,把MTD展示出来. [video widt ...

  2. vue大型电商项目尚品汇(前台篇)day05终结篇

    前台部分到此结束,一路走来还挺怀念,今天主要是对整个项目的完成做一个最后的收尾工作,对于功能上的需求没有什么了,主要就是项目上线的一些注意事项. 一.个人中心二级路由 当我们点击查看订单应该跳转到个人 ...

  3. 《C Primer Plus》第六版笔记--7~10章

    目录 第七章 C控制语句:分支和跳转 第八章 字符输入/输出和输入验证 第九章 函数 第十章 数组和指针 第七章 C控制语句:分支和跳转 if else 用法 if (expression) //ex ...

  4. 数据分析工具Metabase--Metabase安装(最详细的安装教程)

    Meatabase介绍 Metabase 是一款开源的BI工具.主要可以实现在线的可视化分析,单独生成分析图标,定时刷新数据集,权限管理,报告分享等一系列功能. Metabase支持多种市面上主流的数 ...

  5. 入坑KeePass(四)KeePass通过坚果云WebDav同步方法

    参考博客:什么值得买平台-Keepass+坚果云:多平台密码无缝同步 .1.在坚果云上创建一个同步专用的文件夹(建议勾上"默认不同步到电脑"选项) 创建之后可以把刚才创建的数据库文 ...

  6. Java 基础常见知识点&面试题总结(下),2022 最新版!

    你好,我是 Guide.秋招即将到来,我对 JavaGuide 的内容进行了重构完善,同步一下最新更新,希望能够帮助你. 前两篇: Java 基础常见知识点&面试题总结(上),2022 最新版 ...

  7. SAP BDC 用户输入日期转系统日期格式: CONVERT_DATE_TO_EXTERNAL

    BDC中,日期输入格式不正确:可调用FM  CONVERT_DATE_TO_EXTERNAL DATA:l_bdcfield LIKE bdcdata-fval."BDC field val ...

  8. django--ORM表的多对一关系

    多对一关系是什么 Django使用django.db.models.ForeignKey定义多对一关系. ForeignKey需要一个位置参数:与该模型关联的类  class Info(models. ...

  9. Linux文本三剑客-sed

    sed工作原理: sed: Stream Editor.流编辑器 --- 属于行编辑工具 sed和vim一样都是文本编辑工具. 行编辑工具:一行一行处理文件内容 全屏编辑工具:一次性将文件内容加载到内 ...

  10. openGauss内核:SQL解析过程分析

    摘要:在传统数据库中SQL引擎一般指对用户输入的SQL语句进行解析.优化的软件模块.SQL的解析过程主要分为:词法.语法和语义分析. 本文分享自华为云社区< openGauss内核分析(三):S ...