多项式全家桶（烂尾不修了，不过板子常数是小的）

最新推荐文章于 2024-05-20 20:08:58 发布
最新推荐文章于 2024-05-20 20:08:58 发布
阅读量102 收藏
点赞数
文章标签： 算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_61179348/article/details/134000743
版权

一些约定：下面 f i ( x ) f^i(x) fi(x) 表示 i i i 阶导数。 f ( x ) i f(x)^i f(x)i 表示幂次。若不说明绝大部分除一个多项式时都是代表乘上它的逆。

多项式加减，求导积分

过于简单不讲。 ( x a ) ′ = a x a − 1 , ∫ x a d x = x a + 1 a + 1 + C (x^a)'=ax^{a-1},\displaystyle\int x^a {\rm d}x=\dfrac{x^{a+1}}{a+1}+C (xa)′=axa−1,∫xadx=a+1xa+1​+C。实际操作时 C = 0 C=0 C=0，正确性是好证的。（如果你实在闲可以每次证一下）

多项式乘法，任意模数多项式乘法

目前常数很小而且代码不长的 NTT：参考这里。

任意模数请自行学习，基本 useless。

普通 $\texttt{NTT code}$ 
#include<bits/stdc++.h>
#define fr(x) freopen(#x".in","r",stdin);freopen(#x".out","w",stdout);
using namespace std;
const int mod=998244353,N=4e6+5;
int n,m,a[N],b[N],w[N],mmax;
inline int rd()
{
    int x=0,zf=1;
    char ch=getchar();
    while(ch<'0'||ch>'9') (ch=='-')and(zf=-1),ch=getchar();
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*zf;
}
inline void wr(int x)
{
    if(x==0) return putchar('0'),putchar(' '),void();
    int num[35],len=0;
    while(x) num[++len]=x%10,x/=10;
    for(int i=len;i>=1;i--) putchar(num[i]+'0');
    putchar(' ');
}
inline int bger(int x){return x|=x>>1,x|=x>>2,x|=x>>4,x|=x>>8,x|=x>>16,x+1;}
inline int md(int x){return x>=mod?x-mod:x;}
inline int ksm(int x,int p){int s=1;for(;p;(p&1)&&(s=1ll*s*x%mod),x=1ll*x*x%mod,p>>=1);return s;}
inline void init(int mmax)
{
	for(int i=1,j,k;i<mmax;i<<=1)
		for(w[j=i]=1,k=ksm(3,(mod-1)/(i<<1)),j++;j<(i<<1);j++)
			w[j]=1ll*w[j-1]*k%mod;
}
inline void DNT(int *a,int mmax)
{
	for(int i,j,k=mmax>>1,L,*W,*x,*y,z;k;k>>=1)
		for(L=k<<1,i=0;i<mmax;i+=L)
			for(j=0,W=w+k,x=a+i,y=x+k;j<k;j++,W++,x++,y++)
				*y=1ll*(*x+mod-(z=*y))* *W%mod,*x=md(*x+z);
}
inline void IDNT(int *a,int mmax)
{
	for(int i,j,k=1,L,*W,*x,*y,z;k<mmax;k<<=1)
		for(L=k<<1,i=0;i<mmax;i+=L)
			for(j=0,W=w+k,x=a+i,y=x+k;j<k;j++,W++,x++,y++)
				z=1ll* *W* *y%mod,*y=md(*x+mod),*x=md(*x);
	reverse(a+1,a+mmax);
	for(int inv=ksm(mmax,mod-2),i=0;i<mmax;i++) a[i]=1ll*a[i]*inv%mod;
}
inline void NTT(int *a,int *b,int n,int m)
{
	mmax=bger(n+m);init(mmax);DNT(a,mmax);DNT(b,mmax);
	for(int i=0;i<mmax;i++) a[i]=1ll*a[i]*b[i]%mod;IDNT(a,mmax);
}
int main()
{
	n=rd();m=rd();
	for(int i=0;i<=n;i++) a[i]=rd();
	for(int i=0;i<=m;i++) b[i]=rd();
	NTT(a,b,n,m);
	for(int i=0;i<=n+m;i++) wr(a[i]);
    return 0;
}

多项式牛顿迭代

给定多项式 g ( x ) g(x) g(x)，已知 g ( f ( x ) ) ≡ 0 ( m o d x n ) g(f(x))\equiv 0\pmod{x^n} g(f(x))≡0(modxn)，要解出 f ( x ) f(x) f(x)。其中 g ( x ) g(x) g(x) 的形式较为简单（具体看下面求逆）。（注：这没有模板，只是一个前置知识，会应用到各个模板）

考虑倍增。

若已经求出 f 0 ( x ) f_0(x) f0​(x) 满足： g ( f 0 ( x ) ) ≡ 0 ( m o d x ⌈ n / 2 ⌉ ) g(f_0(x))\equiv 0\pmod{x^{\lceil n/2\rceil}} g(f0​(x))≡0(modx⌈n/2⌉)，下面考虑求出 f ( x ) f(x) f(x)。

把 g ( f ( x ) ) g(f(x)) g(f(x)) 在 f 0 ( x ) f_0(x) f0​(x) 出泰勒展开：

0 ≡ g ( f ( x ) ) = ∑ i ≥ 0 1 i ! g i ( f 0 ( x ) ) ( f ( x ) − f 0 ( x ) ) i ( m o d x n ) 0\equiv g(f(x))=\sum\limits_{i\ge 0} \dfrac{1}{i!}g^i(f_0(x))(f(x)-f_0(x))^i\pmod {x^n} 0≡g(f(x))=i≥0∑​i!1​gi(f0​(x))(f(x)−f0​(x))i(modxn)。

注意到 f ( x ) − f 0 ( x ) f(x)-f_0(x) f(x)−f0​(x) 的最低次项次数 ≥ ⌈ n / 2 ⌉ \ge \lceil n/2\rceil ≥⌈n/2⌉，于是 ∀ i ≥ 2 , ( f ( x ) − f 0 ( x ) ) i ≡ 0 ( m o d x n ) \forall i\ge 2,(f(x)-f_0(x))^i\equiv 0\pmod {x^n} ∀i≥2,(f(x)−f0​(x))i≡0(modxn)。

于是 g ′ ( f 0 ( x ) ) ( f ( x ) − f 0 ( x ) ) + g ( f 0 ( x ) ) ≡ 0 ( m o d x n ) ⇒ f ( x ) ≡ f 0 ( x ) − g ( f 0 ( x ) ) g ′ ( f 0 ( x ) ) ( m o d x n ) g'(f_0(x))(f(x)-f_0(x))+g(f_0(x))\equiv 0\pmod {x^n}\Rightarrow f(x)\equiv f_0(x)-\dfrac{g(f_0(x))}{g'(f_0(x))}\pmod {x^n} g′(f0​(x))(f(x)−f0​(x))+g(f0​(x))≡0(modxn)⇒f(x)≡f0​(x)−g′(f0​(x))g(f0​(x))​(modxn)。

你可能觉得这不是还要求逆吗？你为啥不先讲求逆？你先别急。

多项式求逆

给定一个多项式 f ( x ) f(x) f(x)，请求出一个多项式 g ( x ) g(x) g(x)，满足 f ( x ) × g ( x ) ≡ 1 ( m o d x n ) f(x)\times g(x)\equiv1\pmod{x^n} f(x)×g(x)≡1(modxn)。系数对 998244353 998244353 998244353 取模。

考虑牛顿迭代。

G ( X ) = f ( x ) − 1 X G(X)=f(x)-\dfrac{1}{X} G(X)=f(x)−X1​，则 G ( g ( x ) ) ≡ 0 ( m o d x n ) G(g(x))\equiv0 \pmod {x^n} G(g(x))≡0(modxn)。

假设求出了 g 0 ( x ) g_0(x) g0​(x) 使得 G ( g 0 ( x ) ) ≡ 0 ( m o d x ⌈ n / 2 ⌉ ) G(g_0(x))\equiv 0\pmod{x^{\lceil n/2\rceil}} G(g0​(x))≡0(modx⌈n/2⌉)，则：

g ( x ) ≡ g 0 ( x ) − G ( g 0 ( x ) ) G ′ ( g 0 ( x ) ) = g 0 ( x ) − f ( x ) − 1 g 0 ( x ) 1 g 0 ( x ) 2 ≡ 2 g 0 ( x ) − f ( x ) g 0 ( x ) 2 ( m o d x n ) g(x)\equiv g_0(x)-\dfrac{G(g_0(x))}{G'(g_0(x))}=g_0(x)-\dfrac{f(x)-\dfrac{1}{g_0(x)}}{\dfrac{1}{g_0(x)^2}}\equiv 2g_0(x)-f(x)g_0(x)^2\pmod {x^n} g(x)≡g0​(x)−G′(g0​(x))G(g0​(x))​=g0​(x)−g0​(x)21​f(x)−g0​(x)1​​≡2g0​(x)−f(x)g0​(x)2(modxn)。

用两次多项式乘法即可，但是可以通过点值进行特殊变换做到一次多项式乘法。

即一般的多项式乘法设两个对应的点值为 x , y x,y x,y，则变换 f ( x , y ) = x y f(x,y)=xy f(x,y)=xy，这里 f ( x , y ) = 2 y − x y 2 f(x,y)=2y-xy^2 f(x,y)=2y−xy2 即可。

复杂度 T ( n ) = T ( n / 2 ) + O ( n log ⁡ n ) = O ( n log ⁡ n ) T(n)=T(n/2)+O(n\log n)=O(n\log n) T(n)=T(n/2)+O(nlogn)=O(nlogn)。

$\texttt{code}$ 
#include<bits/stdc++.h>
#define LL long long
using namespace std;
const int mod=998244353,N=4e5+5;
int n,m,a[N],b[N],w[N],mmax;
inline int rd()
{
    int x=0,zf=1;
    char ch=getchar();
    while(ch<'0'||ch>'9') (ch=='-')and(zf=-1),ch=getchar();
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*zf;
}
inline void wr(int x)
{
    if(x==0) return putchar('0'),putchar(' '),void();
    int num[35],len=0;
    while(x) num[++len]=x%10,x/=10;
    for(int i=len;i>=1;i--) putchar(num[i]+'0');
    putchar(' ');
}
inline int bger(int x){return x|=x>>1,x|=x>>2,x|=x>>4,x|=x>>8,x|=x>>16,x+1;}
inline int md(int x){return x>=mod?x-mod:x;}
inline int ksm(int x,int p){int s=1;for(;p;(p&1)&&(s=1ll*s*x%mod),x=1ll*x*x%mod,p>>=1);return s;}
inline void init(int mmax)
{
	for(int i=1,j,k;i<mmax;i<<=1)
		for(w[j=i]=1,k=ksm(3,(mod-1)/(i<<1)),j++;j<(i<<1);j++)
			w[j]=1ll*w[j-1]*k%mod;
}
inline void DNT(int *a,int mmax)
{
	for(int i,j,k=mmax>>1,L,*W,*x,*y,z;k;k>>=1)
		for(L=k<<1,i=0;i<mmax;i+=L)
			for(j=0,W=w+k,x=a+i,y=x+k;j<k;j++,W++,x++,y++)
				*y=1ll*(*x+mod-(z=*y))* *W%mod,*x=md(*x+z);
}
inline void IDNT(int *a,int mmax)
{
	for(int i,j,k=1,L,*W,*x,*y,z;k<mmax;k<<=1)
		for(L=k<<1,i=0;i<mmax;i+=L)
			for(j=0,W=w+k,x=a+i,y=x+k;j<k;j++,W++,x++,y++)
				z=1ll* *W* *y%mod,*y=md(*x+mod-z),*x=md(*x+z);
	reverse(a+1,a+mmax);
	for(int inv=ksm(mmax,mod-2),i=0;i<mmax;i++) a[i]=1ll*a[i]*inv%mod;
}
void INV(int num,int *a,int *b)
{
	if(num==1) return b[0]=ksm(a[0],mod-2),void();
	INV((num+1)>>1,a,b);
	int mmax=bger(num<<1);init(mmax);
	static int c[N];
	for(int i=0;i<num;i++) c[i]=a[i];for(int i=num;i<mmax;i++) c[i]=0;
	DNT(c,mmax);DNT(b,mmax);
	for(int i=0;i<mmax;i++) b[i]=1ll*(2-1ll*c[i]*b[i]%mod+mod)%mod*b[i]%mod;
	IDNT(b,mmax);
	for(int i=num;i<mmax;i++) b[i]=0;
}
signed main()
{
	n=rd();
	for(int i=0;i<n;i++) a[i]=rd();
	INV(n,a,b);
	for(int i=0;i<n;i++) wr(b[i]);
	return 0;
}

分治 FFT

给定序列 g 1 , 2 … n − 1 g_{1,2\dots n - 1} g1,2…n−1​，求序列 f 0 , 1 … n − 1 f_{0,1\dots n - 1} f0,1…n−1​。

其中 f i = ∑ j = 1 i f i − j g j f_i=\sum\limits_{j=1}^if_{i-j}g_j fi​=j=1∑i​fi−j​gj​，边界为 f 0 = 1 f_0=1 f0​=1。答案对 998244353 998244353 998244353 取模。

设 F ( x ) = ∑ i ≥ 0 f i x i , G ( x ) = ∑ i ≥ 0 g i x i F(x)=\sum\limits_{i\ge 0}f_ix^i,G(x)=\sum\limits_{i\ge 0}g_ix^i F(x)=i≥0∑​fi​xi,G(x)=i≥0∑​gi​xi，其中令 g 0 = 0 g_0=0 g0​=0。

则 F ( x ) × G ( x ) = ∑ i ≥ 0 x i ∑ j = 0 i f i − j g j = ∑ i ≥ 0 x i ∑ j = 0 i f i − j g j = ∑ i ≥ 1 x i f i = F ( x ) − f 0 = F ( x ) − 1 F(x)\times G(x)=\sum\limits_{i\ge 0}x^i\sum\limits_{j=0}^i f_{i-j}g_j=\sum\limits_{i\ge 0}x^i\sum\limits_{j=0}^i f_{i-j}g_j=\sum\limits_{i\ge 1}x^if_i=F(x)-f_0=F(x)-1 F(x)×G(x)=i≥0∑​xij=0∑i​fi−j​gj​=i≥0∑​xij=0∑i​fi−j​gj​=i≥1∑​xifi​=F(x)−f0​=F(x)−1。

于是 G ( x ) = 1 1 − F ( x ) G(x)=\dfrac{1}{1-F(x)} G(x)=1−F(x)1​。多项式求逆即可，复杂度 O ( n log ⁡ n ) O(n\log n) O(nlogn)。

$\texttt{code}$ 
#include<bits/stdc++.h>
#define fr(x) freopen(#x".in","r",stdin);freopen(#x".out","w",stdout);
using namespace std;
const int mod=998244353,N=4e5+5;
int n,m,a[N],b[N],w[N],mmax;
inline int rd()
{
    int x=0,zf=1;
    char ch=getchar();
    while(ch<'0'||ch>'9') (ch=='-')and(zf=-1),ch=getchar();
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*zf;
}
inline void wr(int x)
{
    if(x==0) return putchar('0'),putchar(' '),void();
    int num[35],len=0;
    while(x) num[++len]=x%10,x/=10;
    for(int i=len;i>=1;i--) putchar(num[i]+'0');
    putchar(' ');
}
inline int bger(int x){return x|=x>>1,x|=x>>2,x|=x>>4,x|=x>>8,x|=x>>16,x+1;}
inline int md(int x){return x>=mod?x-mod:x;}
inline int ksm(int x,int p){int s=1;for(;p;(p&1)&&(s=1ll*s*x%mod),x=1ll*x*x%mod,p>>=1);return s;}
inline void init(int mmax)
{
	for(int i=1,j,k;i<mmax;i<<=1)
		for(w[j=i]=1,k=ksm(3,(mod-1)/(i<<1)),j++;j<(i<<1);j++)
			w[j]=1ll*w[j-1]*k%mod;
}
inline void DNT(int *a,int mmax)
{
	for(int i,j,k=mmax>>1,L,*W,*x,*y,z;k;k>>=1)
		for(L=k<<1,i=0;i<mmax;i+=L)
			for(j=0,W=w+k,x=a+i,y=x+k;j<k;j++,W++,x++,y++)
				*y=1ll*(*x+mod-(z=*y))* *W%mod,*x=md(*x+z);
}
inline void IDNT(int *a,int mmax)
{
	for(int i,j,k=1,L,*W,*x,*y,z;k<mmax;k<<=1)
		for(L=k<<1,i=0;i<mmax;i+=L)
			for(j=0,W=w+k,x=a+i,y=x+k;j<k;j++,W++,x++,y++)
				z=1ll* *W* *y%mod,*y=md(*x+mod-z),*x=md(*x+z);
	reverse(a+1,a+mmax);
	for(int inv=ksm(mmax,mod-2),i=0;i<mmax;i++) a[i]=1ll*a[i]*inv%mod;
}
void dfs(int num,int *a,int *b)
{
	if(num==1) return b[0]=ksm(a[0],mod-2),void();
	dfs((num+1)>>1,a,b);
	int mmax=bger(num<<1);init(mmax);
	static int c[N];
	for(int i=0;i<num;i++) c[i]=a[i];for(int i=num;i<mmax;i++) c[i]=0;
	DNT(c,mmax);DNT(b,mmax);
	for(int i=0;i<mmax;i++) b[i]=1ll*(2-1ll*c[i]*b[i]%mod+mod)%mod*b[i]%mod;
	IDNT(b,mmax);
	for(int i=num;i<mmax;i++) b[i]=0;
}
signed main()
{
	n=rd();
	for(int i=1;i<n;i++) a[i]=(mod-rd())%mod;a[0]=1;
	dfs(n,a,b);
	for(int i=0;i<n;i++) wr(b[i]);
	return 0;
}

多项式除法

给定一个 n n n 次多项式 F ( x ) F(x) F(x) 和一个 m m m 次多项式 G ( x ) G(x) G(x) ，请求出多项式 Q ( x ) Q(x) Q(x), R ( x ) R(x) R(x)，满足以下条件：

  • Q ( x ) Q(x) Q(x) 次数为 n − m n-m n−m， R ( x ) R(x) R(x) 次数小于 m m m
  • F ( x ) = Q ( x ) × G ( x ) + R ( x ) F(x) = Q(x) \times G(x) + R(x) F(x)=Q(x)×G(x)+R(x)

所有的运算在模 998244353 998244353 998244353 意义下进行。保证 m < n m<n m<n。

这种 trick 基本不会用第二次了。注意到我们只需求出 Q ( x ) Q(x) Q(x) 即可一次多项式乘法解决。

设 F R ( x ) F_R(x) FR​(x) 表示 F ( x ) F(x) F(x) 翻转系数形成的多项式，即 F R ( x ) = x deg ⁡ f F ( 1 x ) F_R(x)=x^{\deg f}F(\frac{1}{x}) FR​(x)=xdegfF(x1​)。而 deg ⁡ F R ( x ) = F ( x ) \deg F_R(x)=F(x) degFR​(x)=F(x)。

F ( x ) = Q ( x ) × G ( x ) + R ( x ) ⇒ x n F ( 1 x ) = ( x m G ( 1 x ) ( x n − m Q ( 1 x ) ) + R ( x ) ⇒ F R ( x ) = G R ( x ) Q R ( x ) + x n − m + 1 R R ( x ) F(x) = Q(x) \times G(x) + R(x)\Rightarrow x^nF(\frac{1}{x})=(x^mG(\frac{1}{x})(x^{n-m}Q(\frac{1}{x}))+R(x)\Rightarrow F_R(x)=G_R(x)Q_R(x)+x^{n-m+1}R_R(x) F(x)=Q(x)×G(x)+R(x)⇒xnF(x1​)=(xmG(x1​)(xn−mQ(x1​))+R(x)⇒FR​(x)=GR​(x)QR​(x)+xn−m+1RR​(x)。

我们消去 R R ( x ) R_R(x) RR​(x) 的影响： F R ( x ) ≡ G R ( x ) Q R ( x ) ( m o d x n − m + 1 ) F_R(x)\equiv G_R(x)Q_R(x)\pmod {x^{n-m+1}} FR​(x)≡GR​(x)QR​(x)(modxn−m+1)。

于是多项式求逆即可，注意细节。复杂度 O ( n log ⁡ n ) O(n\log n) O(nlogn)。

$\texttt{code}$ 
#include<bits/stdc++.h>
#define fr(x) freopen(#x".in","r",stdin);freopen(#x".out","w",stdout);
#define LL long long
using namespace std;
const int mod=998244353,N=4e5+5;
int n,m,a[N],b[N],c[N],d[N],e[N],w[N],mmax;
inline int rd()
{
    int x=0,zf=1;
    char ch=getchar();
    while(ch<'0'||ch>'9') (ch=='-')and(zf=-1),ch=getchar();
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*zf;
}
inline void wr(int x)
{
    if(x==0) return putchar('0'),putchar(' '),void();
    int num[35],len=0;
    while(x) num[++len]=x%10,x/=10;
    for(int i=len;i>=1;i--) putchar(num[i]+'0');
    putchar(' ');
}
inline int bger(int x){return x|=x>>1,x|=x>>2,x|=x>>4,x|=x>>8,x|=x>>16,x+1;}
inline int md(int x){return x>=mod?x-mod:x;}
inline int ksm(int x,int p){int s=1;for(;p;(p&1)&&(s=1ll*s*x%mod),x=1ll*x*x%mod,p>>=1);return s;}
inline void init(int mmax)
{
	for(int i=1,j,k;i<mmax;i<<=1)
		for(w[j=i]=1,k=ksm(3,(mod-1)/(i<<1)),j++;j<(i<<1);j++)
			w[j]=1ll*w[j-1]*k%mod;
}
inline void DNT(int *a,int mmax)
{
	for(int i,j,k=mmax>>1,L,*W,*x,*y,z;k;k>>=1)
		for(L=k<<1,i=0;i<mmax;i+=L)
			for(j=0,W=w+k,x=a+i,y=x+k;j<k;j++,W++,x++,y++)
				*y=1ll*(*x+mod-(z=*y))* *W%mod,*x=md(*x+z);
}
inline void IDNT(int *a,int mmax)
{
	for(int i,j,k=1,L,*W,*x,*y,z;k<mmax;k<<=1)
		for(L=k<<1,i=0;i<mmax;i+=L)
			for(j=0,W=w+k,x=a+i,y=x+k;j<k;j++,W++,x++,y++)
				z=1ll* *W* *y%mod,*y=md(*x+mod-z),*x=md(*x+z);
	reverse(a+1,a+mmax);
	for(int inv=ksm(mmax,mod-2),i=0;i<mmax;i++) a[i]=1ll*a[i]*inv%mod;
}
inline void NTT(int *a,int *b,int n,int m)
{
	mmax=bger(n+m);init(mmax);
	DNT(a,mmax);DNT(b,mmax);
	for(int i=0;i<mmax;i++) a[i]=1ll*a[i]*b[i]%mod;
	IDNT(a,mmax);
}
void INV(int num,int *a,int *b)
{
	if(num==1) return b[0]=ksm(a[0],mod-2),void();
	INV((num+1)>>1,a,b);
	int mmax=bger(num<<1);init(mmax);
	static int c[N];
	for(int i=0;i<num;i++) c[i]=a[i];for(int i=num;i<mmax;i++) c[i]=0;
	DNT(c,mmax);DNT(b,mmax);
	for(int i=0;i<mmax;i++) b[i]=1ll*(2-1ll*c[i]*b[i]%mod+mod)%mod*b[i]%mod;
	IDNT(b,mmax);
	for(int i=num;i<mmax;i++) b[i]=0;
}
inline void div(int *a,int *b,int n,int m,int *d)//余式是a 
{
	reverse(a,a+1+n);reverse(b,b+1+m);INV(n-m+1,b,c);
	int mmax=bger(2*n-m+1);static int e[N];init(mmax);
	for(int i=0;i<=n;i++) d[i]=a[i];
	DNT(d,mmax);DNT(c,mmax);
	for(int i=0;i<mmax;i++) d[i]=1ll*d[i]*c[i]%mod;
	IDNT(d,mmax);reverse(d,d+1+n-m);for(int i=n-m+1;i<=n;i++) d[i]=0;
	for(int i=0;i<n-m+1;i++) e[i]=d[i];
	mmax=bger(n<<1);init(mmax);
	reverse(a,a+1+n);reverse(b,b+1+m);
	DNT(b,mmax);DNT(e,mmax);
	for(int i=0;i<mmax;i++) b[i]=1ll*b[i]*e[i]%mod;
	IDNT(b,mmax);
	for(int i=0;i<m;i++) a[i]=(mod+a[i]-b[i])%mod;
}
signed main()
{
	n=rd();m=rd();
	for(int i=0;i<=n;i++) a[i]=rd();
	for(int i=0;i<=m;i++) b[i]=rd();
	div(a,b,n,m,d);
	for(int i=0;i<n-m+1;i++) wr(d[i]);puts("");
	for(int i=0;i<m;i++) wr(a[i]);
	return 0;
}

多项式 ln/exp

多项式 ln

给出 n − 1 n-1 n−1 次多项式 A ( x ) A(x) A(x)，求一个   m o d     x n \bmod{\:x^n} modxn 下的多项式 B ( x ) B(x) B(x)，满足 B ( x ) ≡ ln ⁡ A ( x ) B(x) \equiv \ln A(x) B(x)≡lnA(x)

在 mod  998244353 \text{mod } 998244353 mod 998244353 下进行，且 a i ∈ [ 0 , 998244353 ) ∩ Z a_i \in [0, 998244353) \cap \mathbb{Z} ai​∈[0,998244353)∩Z

注意到 ( ln ⁡ f ( x ) ) ′ = f ′ ( x ) f ( x ) (\ln f(x))'=\dfrac{f'(x)}{f(x)} (lnf(x))′=f(x)f′(x)​，于是 ln ⁡ ( A ( x ) ) = ∫ A ′ ( x ) A ( x ) d x \ln(A(x))=\displaystyle\int \dfrac{A'(x)}{A(x)} {\rm d}x ln(A(x))=∫A(x)A′(x)​dx。

只需一次求导，一次求逆，一次积分即可。复杂度 O ( n log ⁡ n ) O(n\log n) O(nlogn)。

$\texttt{code}$ 
//洛谷 P4725
//https://www.luogu.com.cn/problem/P4725
#include<bits/stdc++.h>
#define fr(x) freopen(#x".in","r",stdin);freopen(#x".out","w",stdout);
using namespace std;
const int mod=998244353,N=4e5+5;
int n,a[N],w[N],mmax;
inline int rd()
{
    int x=0,zf=1;
    char ch=getchar();
    while(ch<'0'||ch>'9') (ch=='-')and(zf=-1),ch=getchar();
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*zf;
}
inline void wr(int x)
{
    if(x==0) return putchar('0'),putchar(' '),void();
    int num[35],len=0;
    while(x) num[++len]=x%10,x/=10;
    for(int i=len;i>=1;i--) putchar(num[i]+'0');
    putchar(' ');
}
inline int bger(int x){return x|=x>>1,x|=x>>2,x|=x>>4,x|=x>>8,x|=x>>16,x+1;}
inline int md(int x){return x>=mod?x-mod:x;}
inline int ksm(int x,int p){int s=1;for(;p;(p&1)&&(s=1ll*s*x%mod),x=1ll*x*x%mod,p>>=1);return s;}
inline void dao(int *a,int n){for(int i=1;i<n;i++) a[i-1]=1ll*i*a[i]%mod;a[n-1]=0;}
inline void ji(int *a,int n){for(int i=n-1;i>=1;i--) a[i]=1ll*ksm(i,mod-2)*a[i-1]%mod;a[0]=0;}
inline void init(int mmax)
{
	for(int i=1,j,k;i<mmax;i<<=1)
		for(w[j=i]=1,k=ksm(3,(mod-1)/(i<<1)),j++;j<(i<<1);j++)
			w[j]=1ll*w[j-1]*k%mod;
}
inline void DNT(int *a,int mmax)
{
	for(int i,j,k=mmax>>1,L,*W,*x,*y,z;k;k>>=1)
		for(L=k<<1,i=0;i<mmax;i+=L)
			for(j=0,W=w+k,x=a+i,y=x+k;j<k;j++,W++,x++,y++)
				*y=1ll*(*x+mod-(z=*y))* *W%mod,*x=md(*x+z);
}
inline void IDNT(int *a,int mmax)
{
	for(int i,j,k=1,L,*W,*x,*y,z;k<mmax;k<<=1)
		for(L=k<<1,i=0;i<mmax;i+=L)
			for(j=0,W=w+k,x=a+i,y=x+k;j<k;j++,W++,x++,y++)
				z=1ll* *W* *y%mod,*y=md(*x+mod-z),*x=md(*x+z);
	reverse(a+1,a+mmax);
	for(int inv=ksm(mmax,mod-2),i=0;i<mmax;i++) a[i]=1ll*a[i]*inv%mod;
}
inline void NTT(int *a,int *b,int n,int m)
{
	mmax=bger(n+m);init(mmax);
	DNT(a,mmax);DNT(b,mmax);
	for(int i=0;i<mmax;i++) a[i]=1ll*a[i]*b[i]%mod;
	IDNT(a,mmax);
}
void INV(int num,int *a,int *b)
{
	if(num==1) return b[0]=ksm(a[0],mod-2),void();
	INV((num+1)>>1,a,b);
	int mmax=bger(num<<1);init(mmax);
	static int c[N];
	for(int i=0;i<num;i++) c[i]=a[i];for(int i=num;i<mmax;i++) c[i]=0;
	DNT(c,mmax);DNT(b,mmax);
	for(int i=0;i<mmax;i++) b[i]=1ll*(2-1ll*c[i]*b[i]%mod+mod)%mod*b[i]%mod;
	IDNT(b,mmax);
	for(int i=num;i<mmax;i++) b[i]=0;
}
inline void Ln(int *a,int n){static int b[N];for(int i=0;i<bger(n<<1);i++) b[i]=0;INV(n,a,b);dao(a,n);NTT(a,b,n,n);ji(a,n);for(int i=n;i<bger(n<<1);i++) a[i]=0;}
int main()
{
	n=rd();for(int i=0;i<n;i++) a[i]=rd();
	Ln(a,n);for(int i=0;i<n;i++) wr(a[i]);
    return 0;
}

多项式 exp

给出 n − 1 n-1 n−1 次多项式 A ( x ) A(x) A(x)，求一个   m o d     x n \bmod{\:x^n} modxn 下的多项式 B ( x ) B(x) B(x)，满足 B ( x ) ≡ e A ( x ) B(x) \equiv \text e^{A(x)} B(x)≡eA(x)。系数对 998244353 998244353 998244353 取模。

考虑牛顿迭代。

构造 g ( X ) = ln ⁡ ( X ) − A ( x ) g(X)=\ln(X)-A(x) g(X)=ln(X)−A(x)，则 g ( B ( x ) ) ≡ 0 ( m o d x n ) g(B(x))\equiv 0\pmod {x^n} g(B(x))≡0(modxn)。

设 B 0 ( x ) B_0(x) B0​(x) 满足 g ( B 0 ( x ) ) ≡ 0 ( m o d x ⌈ x / 2 ⌉ ) g(B_0(x))\equiv 0\pmod {x^{\lceil x/2\rceil}} g(B0​(x))≡0(modx⌈x/2⌉)

则 B ( x ) ≡ B 0 ( x ) − g ( B 0 ( x ) ) g ′ ( B 0 ( x ) ) = B 0 ( x ) − ln ⁡ ( B 0 ( x ) ) − A ( x ) 1 B 0 ( x ) = B 0 ( x ) ( 1 − ln ⁡ ( B 0 ( x ) ) + A ( x ) ) ( m o d x n ) B(x)\equiv B_0(x)-\dfrac{g(B_0(x))}{g'(B_0(x))}=B_0(x)-\dfrac{\ln(B_0(x))-A(x)}{\dfrac{1}{B_0(x)}}=B_0(x)(1-\ln(B_0(x))+ A(x))\pmod {x^n} B(x)≡B0​(x)−g′(B0​(x))g(B0​(x))​=B0​(x)−B0​(x)1​ln(B0​(x))−A(x)​=B0​(x)(1−ln(B0​(x))+A(x))(modxn)

这其中需要一次多项式乘法，一次多项式 ln ⁡ \ln ln，于是复杂度 T ( n ) = T ( n / 2 ) + O ( n log ⁡ n ) = O ( n log ⁡ n ) T(n)=T(n/2)+O(n\log n)=O(n\log n) T(n)=T(n/2)+O(nlogn)=O(nlogn)。

$\texttt{code}$ 
//落谷 P4726
//https://www.luogu.com.cn/problem/P4726
#include<bits/stdc++.h>
#define fr(x) freopen(#x".in","r",stdin);freopen(#x".out","w",stdout);
using namespace std;
const int mod=998244353,N=4e5+5;
int n,a[N],b[N],w[N],mmax;
inline int rd()
{
    int x=0,zf=1;
    char ch=getchar();
    while(ch<'0'||ch>'9') (ch=='-')and(zf=-1),ch=getchar();
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*zf;
}
inline void wr(int x)
{
    if(x==0) return putchar('0'),putchar(' '),void();
    int num[35],len=0;
    while(x) num[++len]=x%10,x/=10;
    for(int i=len;i>=1;i--) putchar(num[i]+'0');
    putchar(' ');
}
inline int bger(int x){return x|=x>>1,x|=x>>2,x|=x>>4,x|=x>>8,x|=x>>16,x+1;}
inline int md(int x){return x>=mod?x-mod:x;}
inline int ksm(int x,int p){int s=1;for(;p;(p&1)&&(s=1ll*s*x%mod),x=1ll*x*x%mod,p>>=1);return s;}
inline void dao(int *a,int n){for(int i=1;i<n;i++) a[i-1]=1ll*i*a[i]%mod;a[n-1]=0;}
inline void ji(int *a,int n){for(int i=n-1;i>=1;i--) a[i]=1ll*ksm(i,mod-2)*a[i-1]%mod;a[0]=0;}
inline void init(int mmax)
{
	for(int i=1,j,k;i<mmax;i<<=1)
		for(w[j=i]=1,k=ksm(3,(mod-1)/(i<<1)),j++;j<(i<<1);j++)
			w[j]=1ll*w[j-1]*k%mod;
}
inline void DNT(int *a,int mmax)
{
	for(int i,j,k=mmax>>1,L,*W,*x,*y,z;k;k>>=1)
		for(L=k<<1,i=0;i<mmax;i+=L)
			for(j=0,W=w+k,x=a+i,y=x+k;j<k;j++,W++,x++,y++)
				*y=1ll*(*x+mod-(z=*y))* *W%mod,*x=md(*x+z);
}
inline void IDNT(int *a,int mmax)
{
	for(int i,j,k=1,L,*W,*x,*y,z;k<mmax;k<<=1)
		for(L=k<<1,i=0;i<mmax;i+=L)
			for(j=0,W=w+k,x=a+i,y=x+k;j<k;j++,W++,x++,y++)
				z=1ll* *W* *y%mod,*y=md(*x+mod-z),*x=md(*x+z);
	reverse(a+1,a+mmax);
	for(int inv=ksm(mmax,mod-2),i=0;i<mmax;i++) a[i]=1ll*a[i]*inv%mod;
}
inline void NTT(int *a,int *b,int n,int m)
{
	mmax=bger(n+m);init(mmax);
	DNT(a,mmax);DNT(b,mmax);
	for(int i=0;i<mmax;i++) a[i]=1ll*a[i]*b[i]%mod;
	IDNT(a,mmax);
}
void INV(int num,int *a,int *b)
{
	if(num==1) return b[0]=ksm(a[0],mod-2),void();
	INV((num+1)>>1,a,b);
	int mmax=bger(num<<1);init(mmax);
	static int c[N];
	for(int i=0;i<num;i++) c[i]=a[i];for(int i=num;i<mmax;i++) c[i]=0;
	DNT(c,mmax);DNT(b,mmax);
	for(int i=0;i<mmax;i++) b[i]=1ll*(2-1ll*c[i]*b[i]%mod+mod)%mod*b[i]%mod;
	IDNT(b,mmax);
	for(int i=num;i<mmax;i++) b[i]=0;
}
inline void Ln(int *a,int n){static int b[N];for(int i=0;i<bger(n<<1);i++) b[i]=0;INV(n,a,b);dao(a,n);NTT(a,b,n,n);ji(a,n);for(int i=n;i<bger(n<<1);i++) a[i]=0;}
inline void Exp(int *a,int *b,int n)
{
	if(n==1) return b[0]=1,void();
	Exp(a,b,(n+1)>>1);static int c[N];for(int i=0;i<bger(n<<1);i++) c[i]=0;
	for(int i=0;i<n;i++) c[i]=b[i];Ln(c,n);
	for(int i=0;i<n;i++) c[i]=md(mod-c[i]+a[i]);c[0]=md(c[0]+1);
	NTT(b,c,n,n);for(int i=n;i<bger(n<<1);i++) b[i]=0;
}
int main()
{
	n=rd();for(int i=0;i<n;i++) a[i]=rd();
	Exp(a,b,n);for(int i=0;i<n;i++) wr(b[i]);
    return 0;
}

多项式快速幂/加强版

给定一个 n − 1 n-1 n−1 次多项式 A ( x ) A(x) A(x)，求一个在   m o d     x n \bmod\ x^n mod xn 意义下的多项式 B ( x ) B(x) B(x)，使得 B ( x ) ≡ ( A ( x ) ) k   (   m o d     x n ) B(x) \equiv (A(x))^k \ (\bmod\ x^n) B(x)≡(A(x))k (mod xn)。

多项式的系数在   m o d     998244353 \bmod\ 998244353 mod 998244353 的意义下进行运算。

普通版保证 a 0 = 1 a_0=1 a0​=1，加强版则不保证。

注意到 ( A ( x ) ) k = exp ⁡ ( k ln ⁡ ( A ( x ) ) ) (A(x))^k=\exp(k\ln(A(x))) (A(x))k=exp(kln(A(x)))，于是一次 exp 一次 ln 即可，复杂度 O ( n log ⁡ n ) O(n\log n) O(nlogn)。

$\texttt{code}$ 
//洛谷 P4726
//https://www.luogu.com.cn/problem/P4726
#include<bits/stdc++.h>
#define fr(x) freopen(#x".in","r",stdin);freopen(#x".out","w",stdout);
using namespace std;
const int mod=998244353,N=4e6+5;
int n,m,a[N],b[N],w[N],mmax;
inline int rd()
{
    int x=0,zf=1;
    char ch=getchar();
    while(ch<'0'||ch>'9') (ch=='-')and(zf=-1),ch=getchar();
    while(ch>='0'&&ch<='9') x=(10ll*x%mod+ch-'0')%mod,ch=getchar();
    return x*zf;
}
inline void wr(int x)
{
    if(x==0) return putchar('0'),putchar(' '),void();
    int num[35],len=0;
    while(x) num[++len]=x%10,x/=10;
    for(int i=len;i>=1;i--) putchar(num[i]+'0');
    putchar(' ');
}
inline int bger(int x){return x|=x>>1,x|=x>>2,x|=x>>4,x|=x>>8,x|=x>>16,x+1;}
inline int md(int x){return x>=mod?x-mod:x;}
inline int ksm(int x,int p){int s=1;for(;p;(p&1)&&(s=1ll*s*x%mod),x=1ll*x*x%mod,p>>=1);return s;}
inline void dao(int *a,int n){for(int i=1;i<n;i++) a[i-1]=1ll*i*a[i]%mod;a[n-1]=0;}
inline void ji(int *a,int n){for(int i=n-1;i>=1;i--) a[i]=1ll*ksm(i,mod-2)*a[i-1]%mod;a[0]=0;}
inline void init(int mmax)
{
	for(int i=1,j,k;i<mmax;i<<=1)
		for(w[j=i]=1,k=ksm(3,(mod-1)/(i<<1)),j++;j<(i<<1);j++)
			w[j]=1ll*w[j-1]*k%mod;
}
inline void DNT(int *a,int mmax)
{
	for(int i,j,k=mmax>>1,L,*W,*x,*y,z;k;k>>=1)
		for(L=k<<1,i=0;i<mmax;i+=L)
			for(j=0,W=w+k,x=a+i,y=x+k;j<k;j++,W++,x++,y++)
				*y=1ll*(*x+mod-(z=*y))* *W%mod,*x=md(*x+z);
}
inline void IDNT(int *a,int mmax)
{
	for(int i,j,k=1,L,*W,*x,*y,z;k<mmax;k<<=1)
		for(L=k<<1,i=0;i<mmax;i+=L)
			for(j=0,W=w+k,x=a+i,y=x+k;j<k;j++,W++,x++,y++)
				z=1ll* *W* *y%mod,*y=md(*x+mod-z),*x=md(*x+z);
	reverse(a+1,a+mmax);
	for(int inv=ksm(mmax,mod-2),i=0;i<mmax;i++) a[i]=1ll*a[i]*inv%mod;
}
inline void NTT(int *a,int *b,int n,int m)
{
	mmax=bger(n+m);init(mmax);
	DNT(a,mmax);DNT(b,mmax);
	for(int i=0;i<mmax;i++) a[i]=1ll*a[i]*b[i]%mod;
	IDNT(a,mmax);
}
void INV(int num,int *a,int *b)
{
	if(num==1) return b[0]=ksm(a[0],mod-2),void();
	INV((num+1)>>1,a,b);
	int mmax=bger(num<<1);init(mmax);
	static int c[N];
	for(int i=0;i<num;i++) c[i]=a[i];for(int i=num;i<mmax;i++) c[i]=0;
	DNT(c,mmax);DNT(b,mmax);
	for(int i=0;i<mmax;i++) b[i]=1ll*(2-1ll*c[i]*b[i]%mod+mod)%mod*b[i]%mod;
	IDNT(b,mmax);
	for(int i=num;i<mmax;i++) b[i]=0;
}
inline void Ln(int *a,int n){static int b[N];for(int i=0;i<bger(n<<1);i++) b[i]=0;INV(n,a,b);dao(a,n);NTT(a,b,n,n);ji(a,n);for(int i=n;i<bger(n<<1);i++) a[i]=0;}
inline void Exp(int *a,int *b,int n)
{
	if(n==1) return b[0]=1,void();
	Exp(a,b,(n+1)>>1);static int c[N];for(int i=0;i<bger(n<<1);i++) c[i]=0;
	for(int i=0;i<n;i++) c[i]=b[i];Ln(c,n);
	for(int i=0;i<n;i++) c[i]=md(mod-c[i]+a[i]);c[0]=md(c[0]+1);
	NTT(b,c,n,n);for(int i=n;i<bger(n<<1);i++) b[i]=0;
}
int main()
{
	n=rd();m=rd();for(int i=0;i<n;i++) a[i]=rd();
	Ln(a,n);for(int i=0;i<n;i++) a[i]=1ll*a[i]*m%mod;
	Exp(a,b,n);for(int i=0;i<n;i++) wr(b[i]);
    return 0;
}

加强版只需注意到：

  • 由于费马小定理，幂次可以对 998244352 998244352 998244352 取模。

  • A ( x ) k = c k ( A ( x ) c ) k A(x)^k=c^k\left(\dfrac{A(x)}{c}\right)^k A(x)k=ck(cA(x)​)k，而后求后面的幂次。

  • 设 A ( x ) = ∑ i = 0 n − 1 a i x i A(x)=\sum\limits_{i=0}^{n-1} a_ix^i A(x)=i=0∑n−1​ai​xi，第一个非零次为 x t x^t xt，则 A ( x ) k = x t k ( ∑ i = t n − 1 a i x i − t ) k A(x)^k=x^{tk}\left(\sum\limits_{i=t}^{n-1}a_ix^{i-t}\right)^k A(x)k=xtk(i=t∑n−1​ai​xi−t)k，而后求后面的幂次。

于是就能转化为 a 0 = 1 a_0=1 a0​=1 的情况啦。复杂度不变。当然加强版基本 useless 啦！

$\texttt{code}$ 
//落谷 P4726
//https://www.luogu.com.cn/problem/P4726
#include<bits/stdc++.h>
#define fr(x) freopen(#x".in","r",stdin);freopen(#x".out","w",stdout);
using namespace std;
const int mod=998244353,N=4e5+5;
int n,m,a[N],b[N],w[N],mmax;
char str[N];
inline int rd()
{
    int x=0,zf=1;
    char ch=getchar();
    while(ch<'0'||ch>'9') (ch=='-')and(zf=-1),ch=getchar();
    while(ch>='0'&&ch<='9') x=(10*x+ch-'0'),ch=getchar();
    return x*zf;
}
inline void wr(int x)
{
    if(x==0) return putchar('0'),putchar(' '),void();
    int num[35],len=0;
    while(x) num[++len]=x%10,x/=10;
    for(int i=len;i>=1;i--) putchar(num[i]+'0');
    putchar(' ');
}
inline int bger(int x){return x|=x>>1,x|=x>>2,x|=x>>4,x|=x>>8,x|=x>>16,x+1;}
inline int md(int x){return x>=mod?x-mod:x;}
inline int ksm(int x,int p){int s=1;for(;p;(p&1)&&(s=1ll*s*x%mod),x=1ll*x*x%mod,p>>=1);return s;}
inline void dao(int *a,int n){for(int i=1;i<n;i++) a[i-1]=1ll*i*a[i]%mod;a[n-1]=0;}
inline void ji(int *a,int n){for(int i=n-1;i>=1;i--) a[i]=1ll*ksm(i,mod-2)*a[i-1]%mod;a[0]=0;}
inline void init(int mmax)
{
	for(int i=1,j,k;i<mmax;i<<=1)
		for(w[j=i]=1,k=ksm(3,(mod-1)/(i<<1)),j++;j<(i<<1);j++)
			w[j]=1ll*w[j-1]*k%mod;
}
inline void DNT(int *a,int mmax)
{
	for(int i,j,k=mmax>>1,L,*W,*x,*y,z;k;k>>=1)
		for(L=k<<1,i=0;i<mmax;i+=L)
			for(j=0,W=w+k,x=a+i,y=x+k;j<k;j++,W++,x++,y++)
				*y=1ll*(*x+mod-(z=*y))* *W%mod,*x=md(*x+z);
}
inline void IDNT(int *a,int mmax)
{
	for(int i,j,k=1,L,*W,*x,*y,z;k<mmax;k<<=1)
		for(L=k<<1,i=0;i<mmax;i+=L)
			for(j=0,W=w+k,x=a+i,y=x+k;j<k;j++,W++,x++,y++)
				z=1ll* *W* *y%mod,*y=md(*x+mod-z),*x=md(*x+z);
	reverse(a+1,a+mmax);
	for(int inv=ksm(mmax,mod-2),i=0;i<mmax;i++) a[i]=1ll*a[i]*inv%mod;
}
inline void NTT(int *a,int *b,int n,int m)
{
	mmax=bger(n+m);init(mmax);
	DNT(a,mmax);DNT(b,mmax);
	for(int i=0;i<mmax;i++) a[i]=1ll*a[i]*b[i]%mod;
	IDNT(a,mmax);
}
void INV(int num,int *a,int *b)
{
	if(num==1) return b[0]=ksm(a[0],mod-2),void();
	INV((num+1)>>1,a,b);
	int mmax=bger(num<<1);init(mmax);
	static int c[N];
	for(int i=0;i<num;i++) c[i]=a[i];for(int i=num;i<mmax;i++) c[i]=0;
	DNT(c,mmax);DNT(b,mmax);
	for(int i=0;i<mmax;i++) b[i]=1ll*(2-1ll*c[i]*b[i]%mod+mod)%mod*b[i]%mod;
	IDNT(b,mmax);
	for(int i=num;i<mmax;i++) b[i]=0;
}
inline void Ln(int *a,int n){static int b[N];for(int i=0;i<bger(n<<1);i++) b[i]=0;INV(n,a,b);dao(a,n);NTT(a,b,n,n);ji(a,n);for(int i=n;i<bger(n<<1);i++) a[i]=0;}
inline void Exp(int *a,int *b,int n)
{
	if(n==1) return b[0]=1,void();
	Exp(a,b,(n+1)>>1);static int c[N];for(int i=0;i<bger(n<<1);i++) c[i]=0;
	for(int i=0;i<n;i++) c[i]=b[i];Ln(c,n);
	for(int i=0;i<n;i++) c[i]=md(mod-c[i]+a[i]);c[0]=md(c[0]+1);
	NTT(b,c,n,n);for(int i=n;i<bger(n<<1);i++) b[i]=0;
}
inline void ksm(int *a,char *str,int n)
{
	int t=n,t1,t2,t3,L=strlen(str),X=0,XX=0;
	static int A[N],b[N];for(int i=0;i<n;i++) A[i]=b[i]=0;
	for(int i=0;i<L;i++) m=(10ll*m%mod+str[i]-'0')%mod;
	for(int i=0;i<L;i++) X=(10ll*X%(mod-1)+str[i]-'0')%(mod-1);
	for(int i=0;i<min(L,7);i++) XX=(10*XX+str[i]-'0');
	if(a[0]==0&&XX>=n){for(int i=0;i<n;i++) a[i]=0;return;}
	for(int i=0;i<n;i++) if(a[i]!=0){t=i;break;}
	for(int i=t;i<n;i++) A[i-t]=a[i];t1=ksm(A[0],mod-2),t2=ksm(A[0],X);
	for(int i=0;i<n-t;i++) A[i]=1ll*A[i]*t1%mod;
	Ln(A,n-t);for(int i=0;i<n-t;i++) A[i]=1ll*A[i]*m%mod;t3=min(1ll*t*m,1ll*n);
	Exp(A,b,n-t);for(int i=0;i<t3;i++) A[i]=0;for(int i=t3;i<n;i++) A[i]=1ll*b[i-t3]*t2%mod;
	for(int i=0;i<n;i++) a[i]=A[i];
}
int main()
{
	n=rd();scanf("%s",str);
	for(int i=0;i<n;i++) a[i]=rd();
	ksm(a,str,n);for(int i=0;i<n;i++) wr(a[i]);
    return 0;
}

多项式开根/加强版

给定一个 n − 1 n-1 n−1 次多项式 A ( x ) A(x) A(x) ，求一个在   m o d     x n \bmod\ {x^n} mod xn 意义下的多项式 B ( x ) B(x) B(x) ，使得 B 2 ( x ) ≡ A ( x ) ( m o d x n ) B^2(x)\equiv A(x)\pmod {x^n} B2(x)≡A(x)(modxn)。

多项式的系数在   m o d     998244353 \bmod\ 998244353 mod 998244353 的意义下进行运算。

普通版保证 a 0 = 1 a_0=1 a0​=1，加强版仅保证 a 0 a_0 a0​ 是   m o d     998244353 \bmod\ 998244353 mod 998244353 的二次剩余。

直接讲加强版，设我们求出 a 0 a_0 a0​ 的二次剩余为 r r r。

考虑牛顿迭代。 g ( X ) = X 2 − A ( x ) g(X)=X^2-A(x) g(X)=X2−A(x)，则 g ( B ( x ) ) ≡ 0 ( m o d x n ) g(B(x))\equiv 0\pmod {x^n} g(B(x))≡0(modxn)。

设 B 0 ( x ) B_0(x) B0​(x) 满足 g ( B 0 ( x ) ) ≡ 0 ( m o d x ⌈ x / 2 ⌉ ) g(B_0(x))\equiv 0\pmod {x^{\lceil x/2\rceil}} g(B0​(x))≡0(modx⌈x/2⌉)

则 B ( x ) ≡ B 0 ( x ) − g ( B 0 ( x ) ) g ′ ( B 0 ( x ) ) = B 0 ( x ) − B 0 ( x ) 2 − A ( x ) 2 B 0 ( x ) = B 0 ( x ) 2 + A ( x ) 2 B 0 ( x ) ( m o d x n ) B(x)\equiv B_0(x)-\dfrac{g(B_0(x))}{g'(B_0(x))}=B_0(x)-\dfrac{B_0(x)^2-A(x)}{2B_0(x)}=\dfrac{B_0(x)}{2}+\dfrac{A(x)}{2B_0(x)}\pmod {x^n} B(x)≡B0​(x)−g′(B0​(x))g(B0​(x))​=B0​(x)−2B0​(x)B0​(x)2−A(x)​=2B0​(x)​+2B0​(x)A(x)​(modxn)。

于是这里一次多项式求逆即可，边界时 n = 1 n=1 n=1， B ( x ) = r B(x)=r B(x)=r。

复杂度 T ( n ) = T ( n / 2 ) + O ( n log ⁡ n ) = O ( n log ⁡ n ) T(n)=T(n/2)+O(n\log n)=O(n\log n) T(n)=T(n/2)+O(nlogn)=O(nlogn)。

加强版 $\texttt{code}$ 
//洛谷 P5277
//https://www.luogu.com.cn/problem/P5277
#include<bits/stdc++.h>
#define fr(x) freopen(#x".in","r",stdin);freopen(#x".out","w",stdout);
#define LL long long
using namespace std;
const int mod=998244353,N=4e6+5;
int n,m,a[N],b[N],c,w[N],jc[N],inv[N],mmax;
inline int rd()
{
    int x=0,zf=1;char ch=getchar();
    while(ch<'0'||ch>'9') (ch=='-')and(zf=-1),ch=getchar();
    while(ch>='0'&&ch<='9') x=((x<<3)+(x<<1)+ch-'0')%mod,ch=getchar();
    return x*zf;
}
inline void wr(int x)
{
    if(x==0) return putchar('0'),putchar(' '),void();
    int num[35],len=0;
    while(x) num[++len]=x%10,x/=10;
    for(int i=len;i>=1;i--) putchar(num[i]+'0');
    putchar(' ');
}
inline int bger(int x){return x|=x>>1,x|=x>>2,x|=x>>4,x|=x>>8,x|=x>>16,x+1;}
inline int md(int x){return x>=mod?x-mod:x;}
inline int yy(int x){if(x&1) return (x+mod)>>1;return x>>1;}
inline int ksm(int x,int p){int s=1;for(;p;(p&1)&&(s=1ll*s*x%mod),x=1ll*x*x%mod,p>>=1);return s;}
inline void dao(int *a,int n){for(int i=1;i<n;i++) a[i-1]=1ll*i*a[i]%mod;a[n-1]=0;}
inline void ji(int *a,int n){for(int i=n-1;i>=1;i--) a[i]=1ll*ksm(i,mod-2)*a[i-1]%mod;a[0]=0;}
namespace er
{
	int n,w;
	struct comp
	{
		int x,y;
		inline comp operator *(comp X){return {(1ll*x*X.x+1ll*y*X.y%mod*w)%mod,(1ll*x*X.y+1ll*y*X.x)%mod};}
	};
	inline comp ksm1(comp x,int p){comp s={1,0};for(;p;(p&1)&&(s=s*x,1),x=x*x,p>>=1);return s;}
	inline int Find()
	{
		int x;
		while(1)
		{
			x=rand()%mod;w=((LL)x*x%mod-n+mod)%mod;
			if(ksm(w,(mod-1)>>1)==mod-1) break;
		}
		return ksm1({x,1},(mod+1)>>1).x;
	}	
}using er::Find;
inline void init(int mmax)
{
	for(int i=1,j,k;i<mmax;i<<=1)
		for(w[j=i]=1,k=ksm(3,(mod-1)/(i<<1)),j++;j<(i<<1);j++)
			w[j]=1ll*w[j-1]*k%mod;
}
inline void DNT(int *a,int mmax)
{
	for(int i,j,k=mmax>>1,L,*W,*x,*y,z;k;k>>=1)
		for(L=k<<1,i=0;i<mmax;i+=L)
			for(j=0,W=w+k,x=a+i,y=x+k;j<k;j++,W++,x++,y++)
				*y=1ll*(*x+mod-(z=*y))* *W%mod,*x=md(*x+z);
}
inline void IDNT(int *a,int mmax)
{
	for(int i,j,k=1,L,*W,*x,*y,z;k<mmax;k<<=1)
		for(L=k<<1,i=0;i<mmax;i+=L)
			for(j=0,W=w+k,x=a+i,y=x+k;j<k;j++,W++,x++,y++)
				z=1ll* *W* *y%mod,*y=md(*x+mod-z),*x=md(*x+z);
	reverse(a+1,a+mmax);
	for(int inv=ksm(mmax,mod-2),i=0;i<mmax;i++) a[i]=1ll*a[i]*inv%mod;
}
inline void NTT(int *a,int *b,int n,int m)
{
	mmax=bger(n+m);init(mmax);
	DNT(a,mmax);DNT(b,mmax);
	for(int i=0;i<mmax;i++) a[i]=1ll*a[i]*b[i]%mod;
	IDNT(a,mmax);
}
void INV(int num,int *a,int *b)
{
	if(num==1) return b[0]=ksm(a[0],mod-2),void();INV((num+1)>>1,a,b);
	int mmax=bger(num<<1);init(mmax);static int c[N];
	for(int i=0;i<num;i++) c[i]=a[i];for(int i=num;i<mmax;i++) c[i]=0;DNT(c,mmax);DNT(b,mmax);
	for(int i=0;i<mmax;i++) b[i]=1ll*(2-1ll*c[i]*b[i]%mod+mod)%mod*b[i]%mod;IDNT(b,mmax);
	for(int i=num;i<mmax;i++) b[i]=0;
}
void SQ(int n,int *a,int *b)
{
	if(n==1) return b[0]=c,void();SQ((n+1)>>1,a,b);static int c[N],d[N];
	for(int i=0;i<n;i++) c[i]=b[i];INV(n,c,d);for(int i=n;i<mmax;i++) c[i]=d[i]=0;
	for(int i=0;i<n;i++) c[i]=a[i];NTT(c,d,n-1,n-1);
	for(int i=0;i<n;i++) b[i]=yy(md(c[i]+b[i]));
	for(int i=0;i<mmax;i++) c[i]=d[i]=0;
}
int main()
{
	srand(time(NULL));n=rd();for(int i=0;i<n;i++) a[i]=rd();
	er::n=a[0];c=Find();SQ(n,a,b);
	for(int i=0;i<n;i++) wr((b[0]<=mod-b[0])?b[i]:mod-b[i]);
    return 0;
}

挑战多项式

按照惯例需要给出挑战多项式的代码作为小全家桶。代码。

常系数齐次线性递推

求一个满足 k k k 阶齐次线性递推数列 a i {a_i} ai​ 的第 n n n 项，即给定 f 1 , 2 … k , a 0 , 1 … k − 1 f_{1,2\dots k},a_{0,1\dots k-1} f1,2…k​,a0,1…k−1​，以及递推式： a n = ∑ i = 1 k f i × a n − i a_n=\sum\limits_{i=1}^{k}f_i \times a_{n-i} an​=i=1∑k​fi​×an−i​，求 a n a_n an​

k ≤ 3.2 × 1 0 4 , n ≤ 1 0 9 k\le 3.2\times10^4,n\le 10^9 k≤3.2×104,n≤109.

前置知识-- Bostan-mori 算法

我们有一个生成函数 H ( x ) = F ( x ) G ( x ) H(x)=\dfrac{F(x)}{G(x)} H(x)=G(x)F(x)​ 。其中 deg ⁡ ( F ) , deg ⁡ ( G ) ≤ k ≤ 32000 \deg(F),\deg(G)\le k\le32000 deg(F),deg(G)≤k≤32000。要求其第 n ≤ 1 0 9 n\le10^9 n≤109 项。即求 [ x n ] F ( x ) G ( x ) [x^n]\dfrac{F(x)}{G(x)} [xn]G(x)F(x)​
。

推柿子： [ x n ] F ( x ) G ( x ) = [ x n ] F ( x ) G ( − x ) G ( x ) G ( − x ) = [ x n ] F 0 ( x 2 ) + x F 1 ( x 2 ) G 0 ( x 2 ) = { [ x n / 2 ] F 0 ( x ) G 0 ( x ) ( 2 ∣ n ) [ x n / 2 ] F 1 ( x ) G 0 ( x ) ( 2 ∤ n ) [x^n]\dfrac{F(x)}{G(x)}=[x^n]\dfrac{F(x)G(-x)}{G(x)G(-x)}=[x^n]\dfrac{F_0(x^2)+xF_1(x^2)}{G_0(x^2)}=\begin{cases}[x^{n/2}]\dfrac{F_0(x)}{G_0(x)}(2\mid n)\\\\ [x^{n/2}]\dfrac{F_1(x)}{G_0(x)}(2\nmid n)\end{cases} [xn]G(x)F(x)​=[xn]G(x)G(−x)F(x)G(−x)​=[xn]G0​(x2)F0​(x2)+xF1​(x2)​=⎩ ⎨ ⎧​[xn/2]G0​(x)F0​(x)​(2∣n)[xn/2]G0​(x)F1​(x)​(2∤n)​

又 [ x 0 ] F ( x ) G ( x ) = [ x 0 ] F ( x ) [ x 0 ] G ( x ) [x^0]\dfrac{F(x)}{G(x)}=\dfrac{[x^0]F(x)}{[x^0]G(x)} [x0]G(x)F(x)​=[x0]G(x)[x0]F(x)​。于是就可以求啦！

有一个细节：就是 G ( − x ) G(-x) G(−x) 要和两个多项式分别做卷积，于是要开两个数组（下面的 c , d c,d c,d 数组）分别和两个做卷积，不能只用一个数组分别和两个做卷积。

$\texttt{function}$ 
inline int genf(int *a,int *b,int n,int m,int k)
{
	static int c[N],d[N];
	for(;k;k>>=1)
	{
		memset(c,0,sizeof(c));for(int i=0;i<=m;i++) c[i]=b[i];
		for(int i=1;i<=m;i+=2) c[i]=(mod-c[i])%mod;
		memset(d,0,sizeof(d));for(int i=0;i<=m;i++) d[i]=c[i];
		NTT(a,c,n,m);NTT(b,d,m,m);int j=0;
		for(int i=(k&1);i<=n+m;i+=2,j++) a[i/2]=a[i];n=j-1;j=0;
		for(int i=0;i<=2*m;i+=2,j++) b[i/2]=b[i];m=j-1;
		for(int i=n+1;i<=N-5;i++) a[i]=0;
		for(int i=m+1;i<=N-5;i++) b[i]=0;
	}
	return 1ll*a[0]*ksm(b[0],mod-2)%mod;
}

计算递推部分

其实内核不难发现。只是名字高端一点。

已知： a 0 , 1 , . . . , k − 1 , f 1 , 2 , . . . , k , a n = ∑ i = 1 k f i × a n − i a_{0,1,...,k-1},f_{1,2,...,k},a_n=\sum\limits_{i=1}^{k}f_i \times a_{n-i} a0,1,...,k−1​,f1,2,...,k​,an​=i=1∑k​fi​×an−i​。要求 a a a 的生成函数。

设 a a a 表示 ∑ i ≥ 0 a i x i \sum\limits_{i\ge 0} a_ix^i i≥0∑​ai​xi，由递推式可得： a = a × ∑ i = 1 k f i x i + t ( x ) ( deg ⁡ ( t ) ≤ k − 1 ) a=a\times\sum\limits_{i=1}^k f_ix^i +t(x)(\deg(t)\le k-1) a=a×i=1∑k​fi​xi+t(x)(deg(t)≤k−1)。其中 a a a 是原来的数列， t ( x ) = ∑ i = 0 k − 1 a i x i − ∑ i = 0 k − 2 a i ∑ j = 1 k − 1 − i f j x i + j t(x)=\sum\limits_{i=0}^{k-1}a_ix^i-\sum\limits_{i=0}^{k-2}a_i\sum\limits_{j=1}^{k-1-i}f_jx^{i+j} t(x)=i=0∑k−1​ai​xi−i=0∑k−2​ai​j=1∑k−1−i​fj​xi+j。是一个幂次小于 k − 1 k-1 k−1 的一个多项式（很难求，于是我们不需要知道它是多少）。

于是 a = t ( x ) T ( x ) ( T ( x ) = 1 − ∑ i = 1 k f i x i ) a=\dfrac{t(x)}{T(x)}(T(x)=1-\sum\limits_{i=1}^k f_ix^i) a=T(x)t(x)​(T(x)=1−i=1∑k​fi​xi)。

由于我们知道 a a a 的前 k − 1 k-1 k−1 项 , T ( x ) T(x) T(x) 又有足足 k k k 次，于是 t ( x ) t(x) t(x) 的前 k − 1 k-1 k−1 项可以通过 A × T A\times T A×T 做一个卷积（其中 A A A 为 a a a 的前 k k k 项构成的多项式）模 x k x^k xk 唯一确定。

最后套用 Bostan-mori 算法计算这个生成函数的第 n n n 项即可，复杂度 O ( k log ⁡ k log ⁡ n ) O(k\log k\log n) O(klogklogn)。

$\texttt{code}$ 
#include<bits/stdc++.h>
#define LL long long
#define fr(x) freopen(#x".in","r",stdin);freopen(#x".out","w",stdout);
using namespace std;
const int mod=998244353,N=4e5+5;
int n,k,a[N],b[N],f[N],w[N],mmax;
inline int rd()
{
    int x=0,zf=1;
    char ch=getchar();
    while(ch<'0'||ch>'9') (ch=='-')and(zf=-1),ch=getchar();
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*zf;
}
inline void wr(int x)
{
    if(x==0) return putchar('0'),putchar(' '),void();
    int num[35],len=0;
    while(x) num[++len]=x%10,x/=10;
    for(int i=len;i>=1;i--) putchar(num[i]+'0');
    putchar(' ');
}
inline int bger(int x){return x|=x>>1,x|=x>>2,x|=x>>4,x|=x>>8,x|=x>>16,x+1;}
inline int md(int x){return x>=mod?x-mod:x;}
inline int ksm(int x,int p){int s=1;for(;p;(p&1)&&(s=1ll*s*x%mod),x=1ll*x*x%mod,p>>=1);return s;}
inline void init(int mmax)
{
	for(int i=1,j,k;i<mmax;i<<=1)
		for(w[j=i]=1,k=ksm(3,(mod-1)/(i<<1)),j++;j<(i<<1);j++)
			w[j]=1ll*w[j-1]*k%mod;
}
inline void DNT(int *a,int mmax)
{
	for(int i,j,k=mmax>>1,L,*W,*x,*y,z;k;k>>=1)
		for(L=k<<1,i=0;i<mmax;i+=L)
			for(j=0,W=w+k,x=a+i,y=x+k;j<k;j++,W++,x++,y++)
				*y=1ll*(*x+mod-(z=*y))* *W%mod,*x=md(*x+z);
}
inline void IDNT(int *a,int mmax)
{
	for(int i,j,k=1,L,*W,*x,*y,z;k<mmax;k<<=1)
		for(L=k<<1,i=0;i<mmax;i+=L)
			for(j=0,W=w+k,x=a+i,y=x+k;j<k;j++,W++,x++,y++)
				z=1ll* *W* *y%mod,*y=md(*x+mod-z),*x=md(*x+z);
	reverse(a+1,a+mmax);
	for(int inv=ksm(mmax,mod-2),i=0;i<mmax;i++) a[i]=1ll*a[i]*inv%mod;
}
inline void NTT(int *a,int *b,int n,int m)
{
	mmax=bger(n+m);init(mmax);
	DNT(a,mmax);DNT(b,mmax);
	for(int i=0;i<mmax;i++) a[i]=1ll*a[i]*b[i]%mod;
	IDNT(a,mmax);
}
inline int genf(int *a,int *b,int n,int m,int k)
{
	static int c[N],d[N];
	for(;k;k>>=1)
	{
		memset(c,0,sizeof(c));for(int i=0;i<=m;i++) c[i]=b[i];
		for(int i=1;i<=m;i+=2) c[i]=md(mod-c[i]);
		memset(d,0,sizeof(d));for(int i=0;i<=m;i++) d[i]=c[i];
		NTT(a,c,n,m);NTT(b,d,m,m);int j=0;
		for(int i=(k&1);i<=n+m;i+=2,j++) a[i/2]=a[i];n=j-1;j=0;
		for(int i=0;i<=2*m;i+=2,j++) b[i/2]=b[i];m=j-1;
		for(int i=n+1;i<=N-5;i++) a[i]=0;
		for(int i=m+1;i<=N-5;i++) b[i]=0;
	}
	return 1ll*a[0]*ksm(b[0],mod-2)%mod;
}
inline int rec(int *f,int *a,int n,int k)
{
	for(int i=1;i<=k;i++) f[i]=md(mod-f[i]);f[0]=1;
	static int c[N];memset(c,0,sizeof(c));
	for(int i=0;i<=k;i++) c[i]=f[i];
	NTT(c,a,k,k);for(int i=k;i<=N-5;i++) c[i]=0;
	return genf(c,f,k,k,n);
}
int main()
{
	n=rd();k=rd();
	for(int i=1;i<=k;i++) f[i]=md(rd()%mod+mod);
	for(int i=0;i<k;i++) a[i]=md(rd()%mod+mod);
	wr(rec(f,a,n,k));
	return 0;
}
masterhuanghuang
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
Python 最小二乘法拟合多项式
weixin_43896283的博客
06-15 2298
Python 最小二乘法拟合多项式
Open3D点云数据处理（二十一）：最小二乘多项式拟合
孙悟空
08-23 3346
本文介绍了最小二乘多项式拟合的原理与方法，同时介绍了多项式拟合的评估指标，包括残差平方和、均方误差、均方根误差、平均绝对误差和决定系数。文章介绍了使用numpy库中的np.polyfit()函数实现多项式拟合，并介绍了如何使用np.poly1d()函数将多项式系数转换为拟合方程。最后还提醒了在进行多项式拟合时需要注意的事项，如数据的准备、多项式阶数的选择、评估拟合效果等。
【学习笔记】多项式全家桶（包含全套证明）
繁凡さん的博客
12-29 1万+
目录0x00.多项式基本概述0x01.多项式的基础概念0x02. 多项式的度0x03. 多项式的乘法0x04. 多项式的逆元0x10.拉格朗日插值0x20.多项式乘法0x21.快速傅里叶变换（FFT）0x22.快速数论变换（NTT）0x23.快速沃尔什变换（FWT）0x24.快速莫比乌斯变换（FMT）0x30.多项式求逆0x40.多项式开方0x50.多项式除法 / 取模0x60.多项式牛顿迭代0x70.多项式多点求值 / 快速插值0x80.多项式三角函数0x90.多项式反...
多项式全家桶 基础知识+各种操作+例题（细节证明）
DD(XYX)的硫氰化氘博客
02-19 6054
只讲了一些基本操作，过于深奥的多项式三角函数、快速插值可以出门看 OneInDark 的博客
多项式全家桶（二）：多项式的逆，导，积
周道-Althen的博客
04-03 2万+
      &ThinSpace;\ \ \ \ \ \ \,      这一部分麻烦一点了，但是很重要，几乎所有多项式都得用的，不过好在也不是很复杂的。 1. 多项式求逆       &ThinSpa...
进军多项式（一）：多项式全家桶模板
jwg2732的专栏
02-09 4037
多项式乘法(NTT) 默认大家会了FFT…… 我们知道FFT中因为要用三角函数和复数，所以对精度要求高 并且常数比较大 那么有没有可以代替单位根的呢？ 我们思考一下在FFT中应用的单位根的性质 1:wn0=1w_n^0=1wn0​=1 保证了有初始值 2:wni+k=wniwnjw_n^{i+k}=w_n^iw_n^jwni+k​=wni​wnj​,也就是wnk=(wn1)kw_n^k=(w_n^1)^kwnk​=(wn1​)k 因此只要求出一个就可以递推 3:wnk+n/2=−wnkw_n^{k+n/2}
最小二乘法—多项式拟合非线性函数
纸上得来终觉浅，绝知此事要躬行
07-10 1万+
原文地址：https://www.jianshu.com/p/af0a4f71c05a 目录 一、函数的近似表示—高次多项式 二、误差函数—最小二乘法 三、引出案例函数曲线 四、目标函数 五、优化目标函数 六、优化目标函数—梯度下降法 七、优化目标函数—求解线性方程组 八、python编程实现拟合曲线函数 九、结果分析 一、函数的近似表示—高次多项式 为了研究一些复杂...
繁凡的ACM算法全家桶（全新的模板整合计划）
繁凡さん的博客
01-21 2万+
整理的算法模板合集： ACM模板 点我看算法全家桶系列！！！ 实际上是一个全新的模板整合计划 目录一、多项式全家桶 算法全家桶！！！ 在搞完多项式全家桶之后就想着要不要把我学过的数论，组合数学，图论，数据结构，线性代数什么的都搞一个全家桶，实际上就是一个模板hhh，只不过是整合到一块的，打算把我之前搞的ACM模板全部再整理一下，完善一下，正好也全部复习一遍，把之前冗余的模板全部精炼一下。 大概需要10天的时间，明天先从数论全家桶开始 现在只有一个多项式全家桶撑场面hhh 一、多项式全家桶 点我看多项式
从寻找谷神星的过程，谈最小二乘法实现多项式拟合
Python作业辅导员 - 天元浪子
09-24 5760
科学史上众星云集，璨若星河。这些牛人基本上都是天才，但也不乏无名之辈凭借匪夷所思、骇世惊俗的猜想而立足于巨星之列。比如，门捷列夫，整了一张留空的元素周期表，引得全世界的化学家去做填空题。还有一位德国的中学老师，名唤约翰·提丢斯（Johann Daniel Titius）的，在1766年写下了这么一个数列：
最小二乘法多项式曲线拟合及其python实现
热门推荐
wolfcsharp的博客
04-20 3万+
最小二乘法多项式曲线拟合及其python实现多项式曲线拟合问题描述最小二乘法针对overfitting，加入正则项python实现运行结果 多项式曲线拟合问题描述 问题描述：给定一些数据点，用一个多项式尽可能好的拟合出这些点排布的轨迹，并给出解析解 判断拟合的好坏常用的误差衡量方法是均方根误差，要求均方根误差先要求平方和误差： 然后计算均方根误差： 多项式拟合问题本质是一个优化问题，目标函数是使R...
C++实现最小二乘法一元回归和多项式拟合
03-23
在进行曲线拟合时用的最多的是最小二乘法，其中以一元函数（线性）和多元函数（多项式）居多，下面这个类专门用于进行多项式拟合，可以根据用户输入的阶次进行多项式拟合，算法来自于网上，和GSL的拟合算法对比过，...
最小二乘法拟合多项式_多项式拟合_最小二乘法_最小二乘法拟合多项式_
10-03
最小二乘法拟合多项式：给定数据点（xi
基于最小二乘多项式的焊缝图像拟合研究
01-12
针对不同类型焊缝线提取拟合精度低的问题，提出基于最小二乘多项式的拟合方法。利用实验跟踪系统，将采集后的焊缝图像进行中值滤波和二值化，在二值化图像的基础上提取焊缝的位置坐标。建立拟合焊缝的多项式函数模型...
C++最小二乘法求多项式拟合曲线.doc
11-15
C++最小二乘法求多项式拟合曲线.doc
zuixiaoercheng.rar_多项式拟合_最小二乘法
09-21
利用给定的点，用最小二乘法求拟合多项式。
【算法】二分算法——x的平方根
2302_79031646的博客
05-17 789
用二分边界算法来求解——x的平方根，从暴力求解到二分的优化细说
【算法】【二叉树，DFS，哈希集合，分类讨论】力扣1110. 删点成林
未来就在脚下。
05-20 496
本题通过DFS遍历二叉树，结合分类讨论的方法，逐步删除指定节点并生成新的森林。该算法有效地处理了节点删除后树结构的调整问题，并通过哈希集合优化删除判断的时间复杂度，最终实现了高效的解决方案。
力扣62 不同路径 Java版本
最新发布
m0_47066863的博客
05-20 186
机器人每次只能向下或者向右移动一步。机器人试图达到网格的右下角（在下图中标记为 “Finish” ）。一个机器人位于一个 m x n 网格的左上角 （起始点在下图中标记为 “Start” ）。从左上角开始，总共有 3 条路径可以到达右下角。题目数据保证答案小于等于 2 * 109。输入：m = 3, n = 7。输入：m = 3, n = 2。输入：m = 7, n = 3。输入：m = 3, n = 3。问总共有多少条不同的路径？
【二叉树算法题记录】从中序与后序遍历序列构造二叉树，从前序与中序遍历序列构造二叉树
weixin_43627638的博客
05-15 379
这样我们就能确定整棵树的根节点在。
多项式拟合 最小二乘拟合 分别是什么意思
08-19
多项式拟合是一种拟合数据的方法，它使用多项式函数来逼近给定数据集中的点。多项式函数由若干个项组成，每个项包含一个系数和一个变量的幂次。 最小二乘拟合是一种统计方法，用于拟合一个数学模型与观测数据之间的...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值