洛谷 P3711 - 仓鼠的数学题(多项式)
提供一种不太一样的做法。
假设要求的多项式为 \(f(x)\)。我们考察 \(f(x)-f(x-1)\),不难发现其等于 \(\sum\limits_{i=0}^na_ix^i\)
考虑设 \(f(x)=\sum\limits_{i=0}^{n+1}b_ix^i\),那么直接代入 \(x-1\) 并化简可以得到:
f(x-1)&=\sum\limits_{i=0}^{n+1}b_i(x-1)^i\\
&=\sum\limits_{i=0}^{n+1}b_i\sum\limits_{j=0}^i\dbinom{i}{j}(-1)^{i-j}x^j\\
&=\sum\limits_{i=0}^{n+1}(\sum\limits_{j=i}^{n+1}\dbinom{j}{i}(-1)^{j-i}b_j)x^i
\end{aligned}
\]
故
f(x)-f(x-1)&=\sum\limits_{i=0}^{n+1}(\sum\limits_{j=i+1}^{n+1}\dbinom{j}{i}(-1)^{j-i}b_j)x^i
\end{aligned}
\]
根据 \(f(x)-f(x-1)=\sum\limits_{i=0}^na_ix^i\) 可得
\]
按照套路拆组合数:
\]
设
\]
\]
\]
那么上式可以写作:
\]
喜闻乐见的减法卷积,按照套路设
\]
那么有 \(h'\) 为 \(f'\) 与 \(g\),按照常理是一遍求逆就可以解决的事,但是非常悲催的是 \(g_0=0\),因此无法直接求逆,不过注意到对于 \(f'\) 和 \(h'\) 同样有它们的第 \(0\) 项为 \(0\),因此可以将数组全部向前平移一位然后求逆。还有就是原数组是 \(n+1\) 次多项式,平移以后变成了 \(n\) 次多项式,因此我们可以大致确定的 \(f_i\) 只有 \(n+1\) 位,怎么办呢?首先注意到 \(f'_{n+1}=b_0\),也就是待求多项式的常数项,因此我们可以代入 \(x=0\) 求得 \(f'_{n+1}=a_{0}\)。还有就是由于你对 \(f',g\) 卷积是在 \(\bmod x^{n+1}\) 意义下进行,并且 \(g_0=0\),因此理论上来说 \(f'_n\) 是不能通过 \(f'\) 与 \(g\) 卷积为 \(h'\) 这一条件确定的,不过注意到多项式求逆的过程中如果我们固定住 \(f'_0\),那么所有数都可表示为 \(f'_i=v_i·\dfrac{1}{f'_0}\) 的形式,其中 \(v_i\) 为常数,也就是说我们求出来的 \(f'\) 与真正的 \(f'\) 是存在比例关系的,因此我们考虑代入 \(x=1\) 求出待求多项式的系数和,这样即可求出真正的 \(f'\) 相较于我们求出的 \(f'\) 缩放了多少倍,也就可以求出真正的 \(f'\)。
时间复杂度 \(n\log n\)。
const int pr=3;
const int ipr=332748118;
const int MAXN=255555;
const int MAXP=524288;
int n,fac[MAXN+5],ifac[MAXN+5];
void init_fac(int n){
for(int i=(fac[0]=ifac[0]=ifac[1]=1)+1;i<=n;i++) ifac[i]=1ll*ifac[MOD%i]*(MOD-MOD/i)%MOD;
for(int i=1;i<=n;i++) fac[i]=1ll*fac[i-1]*i%MOD,ifac[i]=1ll*ifac[i-1]*ifac[i]%MOD;
}
int qpow(int x,int e){
int ret=1;
for(;e;e>>=1,x=1ll*x*x%MOD) if(e&1) ret=1ll*ret*x%MOD;
return ret;
}
int rev[MAXP+5];
void NTT(vector<int> &a,int len,int type){
int lg=31-__builtin_clz(len);
for(int i=0;i<len;i++) rev[i]=(rev[i>>1]>>1)|((i&1)<<lg-1);
for(int i=0;i<len;i++) if(i<rev[i]) swap(a[i],a[rev[i]]);
for(int i=2;i<=len;i<<=1){
int W=qpow((type<0)?ipr:pr,(MOD-1)/i);
for(int j=0;j<len;j+=i){
for(int k=0,w=1;k<(i>>1);k++,w=1ll*w*W%MOD){
int X=a[j+k],Y=1ll*w*a[(i>>1)+j+k]%MOD;
a[j+k]=(X+Y)%MOD;a[(i>>1)+j+k]=(X-Y+MOD)%MOD;
}
}
} if(!~type){
int ivn=qpow(len,MOD-2);
for(int i=0;i<len;i++) a[i]=1ll*a[i]*ivn%MOD;
}
}
vector<int> conv(vector<int> a,vector<int> b){
int LEN=1;while(LEN<a.size()+b.size()) LEN<<=1;
a.resize(LEN,0);b.resize(LEN,0);NTT(a,LEN,1);NTT(b,LEN,1);
for(int i=0;i<LEN;i++) a[i]=1ll*a[i]*b[i]%MOD;
NTT(a,LEN,-1);return a;
}
vector<int> getinv(vector<int> a,int len){
vector<int> b(len,0);b[0]=qpow(a[0],MOD-2);
for(int i=2;i<=len;i<<=1){
vector<int> c(b.begin(),b.begin()+(i>>1));
vector<int> d(a.begin(),a.begin()+i);
c=conv(conv(c,c),d);
for(int j=0;j<i;j++) b[j]=(2ll*b[j]-c[j]+MOD)%MOD;
} return b;
}
vector<int> A,B;
int res[MAXN+5];
int main(){
scanf("%d",&n);init_fac(MAXN);A.resize(n+2);int sum=0;
for(int i=n;~i;i--) scanf("%d",&A[i]),sum=(sum+A[i])%MOD,A[i]=1ll*A[i]*fac[n-i]%MOD;
B.resize(n+2);for(int i=1;i<=n+1;i++) B[i-1]=(i&1)?(MOD-ifac[i]):ifac[i];
int LEN=1;while(LEN<=n+2) LEN<<=1;
A.resize(LEN,0);B.resize(LEN,0);
vector<int> C=conv(A,getinv(B,LEN));
// printf("A: ");for(int i=0;i<n+1;i++) printf("%d%c",A[i]," \n"[i==n]);
// printf("B: ");for(int i=0;i<n+1;i++) printf("%d%c",B[i]," \n"[i==n]);
for(int i=1;i<=n+1;i++) res[i]=1ll*C[n+1-i]*ifac[i]%MOD;
res[0]=A[n];int iv=qpow(res[n+1],MOD-2);
for(int i=1;i<=n+1;i++) res[i]=1ll*res[i]*iv%MOD;
int ssum=0;for(int i=1;i<=n+1;i++) ssum=(ssum+res[i])%MOD;
int mul=1ll*sum*qpow(ssum,MOD-2)%MOD;
for(int i=1;i<=n+1;i++) res[i]=1ll*res[i]*mul%MOD;
for(int i=0;i<=n+1;i++) printf("%d%c",res[i]," \n"[i==n+1]);
return 0;
}
洛谷 P3711 - 仓鼠的数学题(多项式)的更多相关文章
- 洛谷 P3711 仓鼠的数学题 [伯努利数 fft]
P3711 仓鼠的数学题 题意: \[ S_m(x) = \sum_{k=0}^x k^m, 0^0=1\quad 求 \sum_{m=0}^n S_m(x)a_m \] 的答案多项式\(\sum_{ ...
- 洛谷P3711 仓鼠的数学题(伯努利数+多项式求逆)
题面 传送门 题解 如果您不知道伯努利数是什么可以去看看这篇文章 首先我们把自然数幂和化成伯努利数的形式 \[\sum_{i=1}^{n-1}i^k={1\over k+1}\sum_{i=0}^k{ ...
- 洛谷 P3711 仓鼠的数学题【伯努利数+多项式科技】
有个东西叫伯努利数--一开始直接·用第一类斯特林推到自闭 式子来源:https://www.luogu.org/blog/ShadowassIIXVIIIIV/solution-p3711 https ...
- FFT/NTT总结+洛谷P3803 【模板】多项式乘法(FFT)(FFT/NTT)
前言 众所周知,这两个东西都是用来算多项式乘法的. 对于这种常人思维难以理解的东西,就少些理解,多背板子吧! 因此只总结一下思路和代码,什么概念和推式子就靠巨佬们吧 推荐自为风月马前卒巨佬的概念和定理 ...
- 洛谷P3412 仓鼠找$Sugar\ II$题解(期望+统计论?)
洛谷P3412 仓鼠找\(Sugar\ II\)题解(期望+统计论?) 标签:题解 阅读体验:https://zybuluo.com/Junlier/note/1327573 原题链接:洛谷P3412 ...
- 洛谷P4238【模板】多项式求逆
洛谷P4238 多项式求逆:http://blog.miskcoo.com/2015/05/polynomial-inverse 注意:直接在点值表达下做$B(x) \equiv 2B'(x) - A ...
- 洛谷p3398仓鼠找suger题解
我现在爱死树链剖分了 题目 具体分析的话在洛谷blog里 这里只是想放一下改完之后的代码 多了一个son数组少了一个for 少了找size最大的儿子的for #include <cstdio&g ...
- 伯努利数学习笔记&&Luogu P3711 仓鼠的数学题
新科技 Luogu P3711 题意 设$ S_{k,n}$表示$ \displaystyle\sum_{i=0}^n i^k$ 求多项式$\displaystyle\sum_{k=0}^n S_{k ...
- 洛谷P5050 【模板】多项式多点求值
传送门 人傻常数大.jpg 因为求逆的时候没清零结果调了几个小时-- 前置芝士 多项式除法,多项式求逆 什么?你不会?左转你谷模板区,包教包会 题解 首先我们要知道一个结论\[f(x_0)\equiv ...
随机推荐
- python的函数参数传递方式
python的一切数据类型都是对象.但是python的对象分为不可变对象和可变对象.python的变量是引用,对python变量的赋值是引用去绑定该对象. 可变对象的数据发生改变,例如列表和字典,引用 ...
- python join的用法
joinn其实就相当于用某个字符串来拼接列表或者元组中的元素 当然也可以将字符串以某一个str拼接起来 得出的结果自然也是字符串 ex1: results: 实例用处: 当我们从某个文件中读出内容时, ...
- Scrum Meeting 0602
零.说明 日期:2021-6-2 任务:简要汇报两日内已完成任务,计划后两日完成任务 一.进度情况 组员 负责 两日内已完成的任务 后两日计划完成的任务 困难 qsy PM&前端 完成后端管理 ...
- Noip模拟39 2021.8.14
T1 打地鼠 都切掉了的简单题 1 #include<bits/stdc++.h> 2 #define int long long 3 using namespace std; 4 con ...
- Vulnstack内网靶场3
Vulnstack内网靶场3 (qiyuanxuetang.net) 环境配置 打开虚拟机镜像为挂起状态,第一时间进行快照,部分服务未做自启,重启后无法自动运行. 挂起状态,账号已默认登陆,cento ...
- 数组中的逆序对 牛客网 剑指Offer
数组中的逆序对 牛客网 剑指Offer 题目描述 在数组中的两个数字,如果前面一个数字大于后面的数字,则这两个数字组成一个逆序对.输入一个数组,求出这个数组中的逆序对的总数P.并将P对10000000 ...
- Python NameError: name 'unicode' is not defined
Python2 的unicode 函数在 Python3 中被命名为 str.在 Python3 中使用 ·str 来代替 Python2 中的 unicode.
- hdu 1160 FatMouse's Speed(最长不下降子序列+输出路径)
题意: FatMouse believes that the fatter a mouse is, the faster it runs. To disprove this, you want to ...
- 在纯JaveScript中实现报表导出:从“PDF”到“JPG”
我们在前端报表中完成了各种工作数据的输入或内容处理之后,需要做什么? 数据的导出! 这些数据的常用导出格式有:PDF.Excel.HTML和图片几大类型. 但总有一些实际应用场景,需要的不仅仅是将现有 ...
- SpringCloud 2020.0.4 系列之 Gateway入门
1. 概述 老话说的好:做人要有幽默感,懂得幽默的人才会活的更开心. 言归正传,今天我们来聊聊 SpringCloud 的网关组件 Gateway,之前我们去访问 SpringCloud 不同服务的接 ...