多项式

多项式

形如

\[F(x) = \sum_{i\ge 0} f_i x^i
\]

这里的多项式实际上指的是形式幂级数, 也就是 \(\mathrm R[[x]]\), 因为我们只关注它每一项的系数.

牛顿迭代法

给定定义域为多项式的函数 \(G(x)\), 求多项式 \(F(x)\), 使得

\[G(F(x)) \equiv 0 \pmod a
\]

考虑倍增求出. 如果现在求出了低 \({\lceil \frac n2 \rceil}\) 项的系数, 即

\[G(F_0(x)) \equiv 0 \pmod{x^{\lceil \frac n2 \rceil}}
\]

, 那么利用泰勒展开可以得到

\[F(x)\equiv F_0(x)-\frac{G(F_0(x))}{G'(F_0(x))}\pmod{x^n} \tag{1}
\]

求逆

求 \(F(x)\), 满足 \(F(x)A(x) \equiv 1 \pmod{x^n}\).

设 \(A(x)\) 为原函数, 那么 \(G(F(x)) = A(x) - \frac{1}{F(x)}\).

带入 \((1)\),

\[\begin{aligned}
F(x) & \equiv F_0(x)-\frac{A(x) - F(x)^{-1}}{F(x)^{-2}} & \pmod{x^n} \\
& \equiv F_0(x)-A(x)F_0^2(x)+F_0(x) \\
& \equiv F_0(x)(2-(A(x)F_0(x)) \\
\end{aligned}
\]

多项式有逆的充要条件是常数项 \(a_0 \not = 0\), 递归终点 \(f_0 = a_0^{-1}\).

还可以设 \(G(F(x)) = F(x)A(x) - 1\), 推导略有不同:

带入 \((1)\), 即

\[F(x) \equiv F_0(x)-\frac{G(F_0(x))}{G'(F_0(x))}\equiv F_0(x)-\frac{A(x)F_0(x)-1}{A(x)} \pmod{x^n}
\]

由于\(A(x)F_0(x) - 1\) 的低 \(\frac n2\) 位均为 \(0\), 那么在模 \(x^n\) 意义下, 它乘上的数只需考虑低 \(\frac n2\) 位即可.

而 \(A(x)F_0(x)\) 的低 \(\frac n2\) 位中, 只有常数项为 \(1\), 其他项为 \(0\). 那么

\[\frac {A(x)F_0(x) - 1} {A(x)F_0(x)} \equiv {A(x)F_0(x) - 1} \pmod{x^{n}}
\]

上面的式子就可以化为

\[\begin{aligned}
F(x) & \equiv F_0(x)-\frac{A(x)F_0(x)-1}{A(x)} & \pmod{x^n} \\
& \equiv F_0(x)-\frac{(A(x)F_0(x)-1)F_0(x)}{A(x)F_0(x)} \\
& \equiv F_0(x)-{(A(x)F_0(x)-1)F_0(x)} \\
& \equiv F_0(x)(2-(A(x)F_0(x)) \\
\end{aligned}
\]

倍增即可求解.

求对数

求 \(F(x)\), 满足 \(F(x) \equiv \ln A(x) \pmod{x^n}\).

\[\begin{aligned}
F(x)
& \equiv \ln A(x) & \pmod{x^n} \\
& \equiv \int \frac{A'(x)}{A(x)} \, \mathrm d x
\end{aligned}
\]

多项式可以求对数的充要条件是常数项 \(a_0 \not = 0\), 经过积分得到常数项为 \(0\).

容易发现在求 \(\ln\) 过程中 \(A(x)\) 的常数项被消成了 \(1\) (分子分母约分).

这里的求对数其实是假设 \(a_0 = 1\), 因为在模意义下, \(\ln x\) 在 \(x \not = 1\) 处无定义.

求指数

求 \(F(x)\), 满足 \(F(x) \equiv e^{A(x)} \pmod{x^n}\).

\[G(x) = A(x) - \ln F(x)
\]

那么, 通过牛顿迭代

\[F(x)\equiv F_0(x)(1-\ln F_0(x)+A(x))\pmod{x^n}
\]

多项式可以求指数的充要条件是常数项 \(a_0 = 0\), 递归终点 \(f_0 = 1\).

快速幂

求 \(F(x)\), 满足 \(F(x) \equiv A^k(x) \pmod{x^n}\). \(k\) 可为 (模意义下有意义的) 任意有理数.

\[\begin{aligned}
F(x)
& \equiv A^k(x) & \pmod{x^n} \\
& \equiv e^{k \ln A(x)}
\end{aligned}
\]

由于 \(\ln\) 过程的限制, 这个过程其实假设了 \(a_0 = 1\).

如果 \(a_0\) 为 \(> 1\) 正整数, 那么设多项式 \(\exp\) 的递归终点为 \(f_0 = a_0^k\) (模意义下; 如果 \(k\) 为分数或负数需要求高次剩余/逆元) 即可;

如果 \(a_0 = 0\), 需要化为 \(A'(x) = \frac{A(x)}{x^l}\), 其中 \(A'(x)\) 常数项不为 \(0\). 答案则为 \((A'(x))^k \cdot x^{lk}\).

多点求值

过于难写... 推荐秦九韶算法, 时间复杂度 \(O(n^2)\).

快速插值

给定 \((x_i, y_i) , \forall i \in \{1,2,\dotsc n\}\), 求一个 \(n-1\) 次多项式 \(F(x)\), 满足 \(F(x_i) \equiv y_i \pmod p\).

同上... 不想写快速插值...

拉格朗日插值

\[F(x)=\sum_{i=1}^ny_i \cdot \frac{\prod_{1\le j\le n,j!=i}(x-x_j)}{\prod_{1\le j\le n,j!=i}(x_i-x_j)}
\]

求单点值的时间复杂度 \(O(n^2)\).

拉格朗日反演

设\(F(x)\), \(G(x)\)是形式幂级数, 满足 \(F(G(x))=x\), 称 \(F(x)\) 为 \(G(x)\) 的复合逆.

根据形式幂级数的性质, 可以得到 \(F(x)\) 常数项为 \(0\), 并且 \(G(F(x)) = x\).

可以证明

\[[x^n]G(x)=\frac{1}{n}[x^{-1}]\frac{1}{F(x)^n}
\]

如果有 \(F(x) = \frac{x}{H(x)}\), 还有

\[[x^n]G(x)=\frac{1}{n}[x^{n-1}](\frac{x}{F(x)})^n
\]

这个过程实际上是求 \(\frac{F(x)}{x}\) 的逆的 \(n\) 次幂. 多项式求逆, 快速幂即可.

另外, 给定另一个多项式 \(A(x)\), 还可以求多项式的复合:

\[[x^n]A(G(x))=\frac{1}{n}[x^{n-1}](A'(x)(\frac{x}{F(x)})^n)
\]

实现的细节

  • 长度为 n, 次数为 n-1
  • 数组长度为超过 2*n 的2的幂
  • 必要时要清空
  • 不要忘记取模!

代码

多项式全家桶:

const int nsz=4e5+50;
const ll nmod=998244353,g=3,ginv=332748118;
//998244353 g=3
//1004535809 g=3 int n;
ll a[nsz],b[nsz],ans[nsz]; ll qp(ll a,ll b=nmod-2){
ll res=1;
for(;b;b>>=1,a=a*a%nmod)if(b&1)res=res*a%nmod;
return res;
}
namespace npoly{
//l means length of array
ll len,l2,rev[nsz]; il void cp(ll *a,ll *b,int l0,int l1){
memcpy(a,b,l0*sizeof(ll));
memset(a+l0,0,(l1-l0)*sizeof(ll));
} il void fftinit(int l0){
l2=0,len=1;
while(len<l0)++l2,len<<=1;
rep(i,0,len-1)rev[i]=(rev[i>>1]>>1)|((i&1)<<(l2-1));//
}
void dft(ll *a,int l,int fl){
rep(i,0,l-1)if(i<rev[i])swap(a[i],a[rev[i]]);
for(int i=1;i<l;i<<=1){
ll wn=qp(fl==1?g:ginv,(nmod-1)/(i<<1));
for(int j=0,p=i<<1;j<l;j+=p){
ll w=1;
for(int k=0;k<i;++k,w=w*wn%nmod){
ll x=a[j+k],y=w*a[j+k+i]%nmod;
a[j+k]=(x+y)%nmod,a[j+k+i]=(x-y)%nmod;
}
}
}
if(fl==-1){
ll v=qp(l);
rep(i,0,l-1)a[i]=a[i]*v%nmod;
}
}
void mul(ll *a,int l1,ll *b,int l2,ll *c,int l3=-1){
if(l3==-1)l3=l1+l2-1;
static ll c1[nsz],c2[nsz];
fftinit(l1+l2-1);
cp(c1,a,l1,len),cp(c2,b,l2,len);
dft(c1,len,1),dft(c2,len,1);
rep(i,0,len-1)c1[i]=c1[i]*c2[i]%nmod;
dft(c1,len,-1);
cp(c,c1,l3,l3);
}
void inv(ll *a,int l,ll *b){
if(l==1){b[0]=qp(a[0]);return;}
int l0=(l+1)>>1;
inv(a,l0,b); static ll c1[nsz],c2[nsz];
fftinit(l<<1);//需要两倍长度dft保证乘法正确
cp(c1,a,l,len),cp(c2,b,l0,len);
dft(c1,len,1),dft(c2,len,1);
rep(i,0,len-1)c1[i]=c2[i]*(2-c1[i]*c2[i]%nmod)%nmod;
dft(c1,len,-1);
cp(b,c1,l,l);
} void deriv(ll *a,int l,ll *b){ //b could be eq to a
rep(i,0,l-2)b[i]=a[i+1]*(i+1)%nmod;
b[l-1]=0;
}
void integ(ll *a,int l,ll *b){
repdo(i,l-2,0)b[i+1]=a[i]*qp(i+1)%nmod;
b[0]=0;
}
//a[0] should be 1
void ln(ll *a,int l,ll *b){
static ll c1[nsz],c2[nsz];
inv(a,l,c1),deriv(a,l,c2);
mul(c1,l,c2,l,b,l);
integ(b,l,b);
} //a[0] should be 0
ll expzero=1;
void exp(ll *a,int l,ll *b){
if(l==1){b[0]=expzero;return;}
int l0=(l+1)>>1;
exp(a,l0,b); static ll c1[nsz],c2[nsz];
fftinit(l<<1);
rep(i,l0,l-1)b[i]=0;
ln(b,l,c1);
rep(i,0,l-1)c1[i]=(a[i]-c1[i]+(i==0))%nmod;
rep(i,l,len-1)c1[i]=0;
cp(c2,b,l0,len);
dft(c1,len,1),dft(c2,len,1);
rep(i,0,len-1)c1[i]=c1[i]*c2[i]%nmod;
dft(c1,len,-1);
cp(b,c1,l,l);
} void pow(ll *a,int l,ll p,ll *b){//suppose a[0]!=0;
static ll c[nsz];
ln(a,l,c);
rep(i,0,l-1)c[i]=c[i]*p%nmod;
// expzero=qp(a[0],p);
exp(c,l,b);
}
void sqrt(ll *a,int l,ll *b){pow(a,l,qp(2),b);} void dncfft(ll *a,int l0,ll *b){
static ll c1[nsz];
rep(i,0,l0-1)c1[i]=-a[i];
c1[0]=(c1[0]+1)%nmod;
inv(c1,l0,b);
} //tests
void testfft(){
int n,m;
cin>>n>>m,++n,++m;
rep(i,0,n-1)cin>>a[i];
rep(i,0,m-1)cin>>b[i];
mul(a,n,b,m,ans);
rep(i,0,n+m-2)cout<<(ans[i]+nmod)%nmod<<' ';
cout<<'\n';
}
//3 1 2 1
//1 -2 3
void testinv(){
int n;
cin>>n;
rep(i,0,n-1)cin>>a[i];
inv(a,n,ans);
rep(i,0,n-1)cout<<(ans[i]+nmod)%nmod<<' ';
cout<<'\n';
}
void testdnc(){
int n;
cin>>n;
rep(i,1,n-1)cin>>a[i];
dncfft(a,n,ans);
rep(i,0,n-1)cout<<(ans[i]+nmod)%nmod<<' ';
cout<<'\n';
}
//3 1 2 1
//0 2 -1
void testln(){
int n;
cin>>n;
rep(i,0,n-1)cin>>a[i];
ln(a,n,ans);
rep(i,0,n-1)cout<<(ans[i]+nmod)%nmod<<' ';
cout<<'\n';
}
//3 0 1 2
//1 1 499122179(5/2)
void testexp(){
int n;
cin>>n;
rep(i,0,n-1)cin>>a[i];
exp(a,n,ans);
rep(i,0,n-1)cout<<(ans[i]+nmod)%nmod<<' ';
cout<<'\n';
}
void testsq(){
int n;
cin>>n;
rep(i,0,n-1)cin>>a[i];
sqrt(a,n,ans);
rep(i,0,n-1)cout<<(ans[i]+nmod)%nmod<<' ';
cout<<'\n';
}
void testpow(){
int n,k;
cin>>n>>k;
rep(i,0,n-1)cin>>a[i];
pow(a,n,k,ans);
rep(i,0,n-1)cout<<(ans[i]+nmod)%nmod<<' ';
cout<<'\n';
}
} int main(){
ios::sync_with_stdio(0),cin.tie(0);
npoly::testexp();
return 0;
}

拉格朗日插值:

pair<ll,ll> li[nsz];
ll lag(ll x){
sort(li+1,li+n+1);
n=unique(li+1,li+n+1)-li-1;
ll res=0;
rep(i,1,n){
ll tmp=1;
rep(j,1,n)if(i!=j)tmp=tmp*(li[i].first-li[j].first)%nmod;
tmp=inv(tmp)*li[i].second%nmod;
rep(j,1,n)if(i!=j)tmp=tmp*(x-li[j].first)%nmod;
res=(res+tmp)%nmod;
}
return res;
}

[模板] 多项式: 乘法/求逆/分治fft/微积分/ln/exp/幂的更多相关文章

  1. 多项式求逆/分治FFT 学习笔记

    一.多项式求逆 给定一个多项式 \(F(x)\),请求出一个多项式 \(G(x)\), 满足 \(F(x) * G(x) \equiv 1 ( \mathrm{mod\:} x^n )\).系数对 \ ...

  2. bzoj 3456 城市规划 多项式求逆+分治FFT

    城市规划 Time Limit: 40 Sec  Memory Limit: 256 MBSubmit: 1091  Solved: 629[Submit][Status][Discuss] Desc ...

  3. 【Uoj34】多项式乘法(NTT,FFT)

    [Uoj34]多项式乘法(NTT,FFT) 题面 uoj 题解 首先多项式乘法用\(FFT\)是一个很久很久以前就写过的东西 直接贴一下代码吧.. #include<iostream> # ...

  4. CF848E Days of Floral Colours——DP+多项式求逆/分治NTT

    官方题解:http://codeforces.com/blog/entry/54233 就是由简入繁 1.序列处理,只考虑一个半圆 2.环形处理(其实这个就是多了旋转同构) 然后基于分割线邻居的跨越与 ...

  5. 洛谷.3803.[模板]多项式乘法(FFT)

    题目链接:洛谷.LOJ. FFT相关:快速傅里叶变换(FFT)详解.FFT总结.从多项式乘法到快速傅里叶变换. 5.4 又看了一遍,这个也不错. 2019.3.7 叕看了一遍,推荐这个. #inclu ...

  6. P3803 [模板] 多项式乘法 (FFT)

    Rt 注意len要为2的幂 #include <bits/stdc++.h> using namespace std; const double PI = acos(-1.0); inli ...

  7. 洛谷.3803.[模板]多项式乘法(NTT)

    题目链接:洛谷.LOJ. 为什么和那些差那么多啊.. 在这里记一下原根 Definition 阶 若\(a,p\)互质,且\(p>1\),我们称使\(a^n\equiv 1\ (mod\ p)\ ...

  8. hdu 5730 Shell Necklace [分治fft | 多项式求逆]

    hdu 5730 Shell Necklace 题意:求递推式\(f_n = \sum_{i=1}^n a_i f_{n-i}\),模313 多么优秀的模板题 可以用分治fft,也可以多项式求逆 分治 ...

  9. FFT模板(多项式乘法)

    FFT模板(多项式乘法) 标签: FFT 扯淡 一晚上都用来捣鼓这个东西了...... 这里贴一位神犇的博客,我认为讲的比较清楚了.(刚好适合我这种复数都没学的) http://blog.csdn.n ...

随机推荐

  1. Java 重建二叉树 根据前序中序重建二叉树

    题目:输入某二叉树的前序遍历和中序遍历的结果,请重建出该二叉树.假设输入的前序遍历和中序遍历的结果中都不含重复的数字.例如输入前序遍历序列{1,2,4,7,3,5,6,8}和中序遍历序列{4,7,2, ...

  2. Xamarin for Visual Studio下载后的文件路径

    Xamarin for Visual Studio的下载很纠结,在官网上不知道如何下载?现在找到一个办法:可以先在网上找一个低版本的之后安装,然后利用VS更新.利用VS更新这里也遇到了问题,下载成功之 ...

  3. 章节九、4-ChromDriver介绍

    一.首先下载Chrom浏览器驱动,将驱动解压到存放火狐浏览器驱动文件路径中(请观看前面的章节) 1.进入该网址下载匹配本地浏览器版本的驱动 http://chromedriver.storage.go ...

  4. PJProject(2.6) 工程介绍

    pjlib pjlib\build\pjlib.vcproj pjlib_test pjlib\build\pjlib_test.vcproj pjsip_core pjsip\build\pjsip ...

  5. C#基础学习第一天

    ..net与C# .NET是一个框架.一种平台.一种技术 C#是一种编程语言,可以开发基于.NET平台的应用 .NET能干什么 Winform ASP.NET Wwb wphone Unity3D游戏 ...

  6. Linux分页机制之分页机制的演变--Linux内存管理(七)

    1 页式管理 1.1 分段机制存在的问题 分段,是指将程序所需要的内存空间大小的虚拟空间,通过映射机制映射到某个物理地址空间(映射的操作由硬件完成).分段映射机制解决了之前操作系统存在的两个问题: 地 ...

  7. c/c++ 头文件的血案

    头文件的血案 不小心在一个头文件里,加了函数的定义,结果导致编译时,提示这个函数被重复定义:( Quote.h #ifndef __QUOTE_H__ #define __QUOTE_H__ #inc ...

  8. mysql解决select * from 表名 (where + 约束条件为空)

    mysql解决select * from 表名 (where + 约束条件为空),示例如下: SELECT * from tableName WHERE name is NULL; 从 tableNa ...

  9. SQLServer之修改FOREIGN KEY约束

    使用SSMS数据库管理工具修改FOREIGN KEY约束 1.连接数据库,选择数据表->右键点击->选择设计(或者展开键,选择要修改的外键,右键点击,选择修改,后面修改步骤相同). 2.在 ...

  10. 正益移动推出新产品正益工作 实现PaaS+SaaS新组合

    近期,正益移动不仅将其AppCan 移动平台云化,还通过发布全新 SaaS 产品 -- 正益工作,这款集合了企业信息聚合.应用聚合.社交聚合为一体的企业移动综合门户,与 AppCan 平台一起实现了P ...