多项式求ln,求exp,开方,快速幂 学习总结

最新推荐文章于 2021-07-23 14:06:04 发布
最新推荐文章于 2021-07-23 14:06:04 发布
阅读量384 收藏
点赞数
原文链接：http://www.cnblogs.com/joyouth/p/5604618.html
版权

按理说Po姐姐三月份来讲课的时候我就应该学了

但是当时觉得比较难加上自己比较懒，所以就QAQ了

现在不得不重新弄一遍了

 

首先说多项式求ln

设G(x)=lnF(x)

我们两边求导可以得到G'(x)=F‘(x)/F(x)

则G(x)就是F’(x)/F(x)的积分

我们知道多项式求导和积分是O(n)的,多项式求逆是O(nlogn)的

所以总时间复杂度O(nlogn)

多项式求ln一般解决的问题是这样的

设多项式f表示一些奇怪的东西，由一些奇怪的东西有序组成的方案为

f^1+f^2+f^3…… 化简之后得到多项式求逆的式子

而由这些奇怪的东西无序组成的方案为

f^1/1!+f^2/2！……

由泰勒展开我们知道化简后为e^f

则若我们知道e^f,求f就需要多项式求ln了

这样推导非常的优美

譬如说BZOJ 3456 城市规划

设答案为多项式f，设n个点的无向图个数为多项式g

不难发现无向图是由若干个无序联通块组成的

我们得到式子g=e^f

而多项式g是非常好构造的，f=ln(g)

多项式求ln即可

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define G 3
using namespace std;
 
typedef long long LL;
const int maxn=300010;
const int mod=1004535809;
int n,N,len,rev[maxn];
int jc[maxn],inv[maxn];
int g[maxn],f[maxn],w[maxn];
 
int pow_mod(int v,int p){
    int tmp=1;
    while(p){
        if(p&1)tmp=1LL*tmp*v%mod;
        v=1LL*v*v%mod;p>>=1;
    }return tmp;
}
void FFT(int *A,int n,int flag){
    for(int i=0;i<n;++i)if(i<rev[i])swap(A[i],A[rev[i]]);
    for(int k=2;k<=n;k<<=1){
        int mk=(k>>1);
        int wn=pow_mod(G,flag==1?(mod-1)/k:(mod-1)-(mod-1)/k);
        w[0]=1;
        for(int i=1;i<mk;++i)w[i]=1LL*w[i-1]*wn%mod;
        for(int i=0;i<n;i+=k){
            for(int j=0;j<mk;++j){
                int a=i+j,b=i+j+mk;
                int x=A[a],y=1LL*A[b]*w[j]%mod;
                A[a]=(x+y)%mod;
                A[b]=x-y;if(A[b]<0)A[b]+=mod;
            }
        }
    }
    if(flag==-1){
        int c=pow_mod(n,mod-2);
        for(int i=0;i<n;++i)A[i]=1LL*A[i]*c%mod;
    }
}
void Get_inv(int n){
    if(n==1){
        f[0]=pow_mod(g[0],mod-2);
        return;
    }
    Get_inv(n>>1);
    int now=(n<<1);
    for(len=0;(1<<len)<now;++len);
    for(int i=0;i<now;++i)rev[i]=rev[i>>1]>>1|((i&1)<<(len-1));
    static int tmp[maxn];
    for(int i=0;i<n;++i)tmp[i]=g[i];
    FFT(tmp,now,1);FFT(f,now,1);
    for(int i=0;i<now;++i)f[i]=1LL*f[i]*(2LL-1LL*f[i]*tmp[i]%mod+mod)%mod;
    FFT(f,now,-1);
    memset(f+n,0,sizeof(int)*n);
}
void Get_Dao(){
    g[N]=0;
    for(int i=0;i<N;++i)g[i]=1LL*g[i+1]*(i+1)%mod;
}
void Get_Fen(){
    f[0]=0;
    for(int i=N-1;i>=0;--i)f[i]=1LL*f[i-1]*pow_mod(i,mod-2)%mod;
}
 
int main(){
    scanf("%d",&n);
    for(N=1;N<=n;N<<=1);
    jc[0]=1;
    for(int i=1;i<N;++i)jc[i]=1LL*jc[i-1]*i%mod;
    inv[N-1]=pow_mod(jc[N-1],mod-2);
    for(int i=N-2;i>=0;--i)inv[i]=1LL*inv[i+1]*(i+1)%mod;
    for(int i=0;i<N;++i){
        int tmp=(1LL*i*(i-1)/2)%(mod-1);
        g[i]=1LL*pow_mod(2,tmp)*inv[i]%mod;
    }
    Get_inv(N);Get_Dao();N<<=1;
    for(len=0;(1<<len)<N;++len);
    for(int i=0;i<N;++i)rev[i]=rev[i>>1]>>1|((i&1)<<(len-1));
    FFT(g,N,1);FFT(f,N,1);
    for(int i=0;i<N;++i)f[i]=1LL*f[i]*g[i]%mod;
    FFT(f,N,-1);Get_Fen();
    printf("%d\n",1LL*f[n]*jc[n]%mod);
    return 0;
}

cojs 2358 

首先用多项式求逆可以搞出n个点的DAG的个数，这个不再多说

我们不难发现n个点的DAG是由若干个DAG的弱连通图无序组成的

多项式求ln即可

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cstdlib>
#include<algorithm>
#define G 3
using namespace std;

const int maxn=300010;
const int mod=998244353;
int n,N,len,w[maxn];
int jc[maxn],inv[maxn],rev[maxn];
int f[maxn],Inv[maxn],ln[maxn];

int pow_mod(int v,int p){
	int tmp=1;
	while(p){
		if(p&1)tmp=1LL*tmp*v%mod;
		v=1LL*v*v%mod;p>>=1;
	}return tmp;
}
void FFT(int *A,int n,int flag){
	for(int i=0;i<n;++i)if(i<rev[i])swap(A[i],A[rev[i]]);
	for(int k=2;k<=n;k<<=1){
		int mk=(k>>1);
		int wn=pow_mod(G,flag==1?(mod-1)/k:(mod-1)-(mod-1)/k);
		w[0]=1;
		for(int i=1;i<mk;++i)w[i]=1LL*w[i-1]*wn%mod;
		for(int i=0;i<n;i+=k){
			for(int j=0;j<mk;++j){
				int a=i+j,b=i+j+mk;
				int x=A[a],y=1LL*A[b]*w[j]%mod;
				A[a]=(x+y)%mod;
				A[b]=x-y;if(A[b]<0)A[b]+=mod;
			}
		}
	}
	if(flag==-1){
		int c=pow_mod(n,mod-2);
		for(int i=0;i<n;++i)A[i]=1LL*A[i]*c%mod;
	}return;
}
void Get_dao(int n){
	f[n]=0;
	for(int i=0;i<n;++i)f[i]=1LL*f[i+1]*(i+1)%mod;
}
void Get_Fen(int n){
	for(int i=n-1;i>=1;--i)ln[i]=1LL*ln[i-1]*pow_mod(i,mod-2)%mod;
	ln[0]=0;
}
void Get_inv(int n){
	if(n==1){
		Inv[0]=pow_mod(f[0],mod-2);
		return;
	}
	Get_inv(n>>1);
	int now=(n<<1);
	for(len=0;(1<<len)<now;++len);
	for(int i=0;i<now;++i)rev[i]=rev[i>>1]>>1|((i&1)<<(len-1));
	static int tmp[maxn];
	for(int i=0;i<n;++i)tmp[i]=f[i];
	for(int i=n;i<now;++i)tmp[i]=0;
	FFT(Inv,now,1);FFT(tmp,now,1);
	for(int i=0;i<now;++i)Inv[i]=1LL*Inv[i]*(2LL-1LL*Inv[i]*tmp[i]%mod+mod)%mod;
	FFT(Inv,now,-1);
	memset(Inv+n,0,sizeof(int)*n);
}
void Get_ln(int n){
	int now=(n<<1);
	for(int i=0;i<n;++i)Inv[i]=0;
	Get_inv(n);Get_dao(n);
	for(len=0;(1<<len)<now;++len);
	for(int i=0;i<now;++i)rev[i]=rev[i>>1]>>1|((i&1)<<(len-1));
	FFT(Inv,now,1);FFT(f,now,1);
	for(int i=0;i<now;++i)ln[i]=1LL*Inv[i]*f[i]%mod;
	FFT(ln,now,-1);
	memset(ln+n,0,sizeof(int)*n);
	Get_Fen(n);
}

int main(){
	freopen("dagIV.in","r",stdin);
	freopen("dagIV.out","w",stdout);
	scanf("%d",&n);
	for(N=1;N<=n;N<<=1);
	jc[0]=1;
	for(int i=1;i<N;++i)jc[i]=1LL*jc[i-1]*i%mod;
	inv[N-1]=pow_mod(jc[N-1],mod-2);
	for(int i=N-2;i>=0;--i)inv[i]=1LL*inv[i+1]*(i+1)%mod;
	for(int i=0;i<N;++i){
		int tmp=(1LL*i*(i-1)/2)%(mod-1);
		f[i]=1LL*pow_mod((mod+1)>>1,tmp)*inv[i]%mod;
		if(i&1)f[i]=mod-f[i];
	}Get_inv(N);
	for(int i=0;i<=n;++i){
		int tmp=(1LL*i*(i-1)/2)%(mod-1);
		f[i]=1LL*Inv[i]*pow_mod(2,tmp)%mod*jc[i]%mod;
		f[i]=1LL*f[i]*inv[i]%mod;
	}
	for(int i=n+1;i<N;++i)f[i]=0;
	Get_ln(N);
	printf("%d\n",1LL*ln[n]*jc[n]%mod);
	return 0;
}

多项式求exp

用处其实跟多项式求ln是相反的

对于式子g=e^f，假设我们知道f，求g就需要用到多项式求exp了

多项式求exp的方式是利用牛顿迭代法倍增

设F(x)=e^A(x)

我们知道若G(F0(x))=0(mod x^n)

我们有G(F(x))=0=G(F0(x))/0!+G‘(F0(x))/1!*(F(x)-F0(x))(mod x^2n)

转换一下式子我们可以得到F(x)的表达式

我们又知道G(F(x))=lnF(x)-A(x)

带入上面的式子整理之后我们得到F(x)=F0(x)*(1-lnF0(x)+A(x))

无脑倍增就可以了,每次倍增的时候求ln即可，常数巨大，但还是O(nlogn)的

本来应该出到cojs的题目QAQ但是因为自己的权限掉了QAQ

求前n项的Bell数

在之前的文章里提到了Bell数可以利用CDQ+FFT解决

但是我们也提到了Bell数的生成函数e^(e^x-1)

我们对e^x-1做泰勒展开，之后多项式求exp即可

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cstdlib>
#include<algorithm>
#define G 3
using namespace std;

const int mod=998244353;
const int maxn=300010;
int n,N,len,T;
int jc[maxn],inv[maxn],rev[maxn],w[maxn];
int g[maxn],f[maxn],h[maxn],ni[maxn];
int Exp[maxn],ln[maxn],Inv[maxn];
int C[maxn];
struct ASK{
	int n;
}Q[maxn];

int pow_mod(int v,int p){
	int tmp=1;
	while(p){
		if(p&1)tmp=1LL*tmp*v%mod;
		v=1LL*v*v%mod;p>>=1;
	}return tmp;
}
void FFT(int *A,int n,int flag){
	for(int i=0;i<n;++i)if(i<rev[i])swap(A[i],A[rev[i]]);
	for(int k=2;k<=n;k<<=1){
		int mk=(k>>1);
		int wn=pow_mod(G,flag==1?(mod-1)/k:(mod-1)-(mod-1)/k);
		w[0]=1;
		for(int i=1;i<mk;++i)w[i]=1LL*w[i-1]*wn%mod;
		for(int i=0;i<n;i+=k){
			for(int j=0;j<mk;++j){
				int a=i+j,b=i+j+mk;
				int x=A[a],y=1LL*A[b]*w[j]%mod;
				A[a]=(x+y)%mod;
				A[b]=x-y;if(A[b]<0)A[b]+=mod;
			}
		}
	}
	if(flag==-1){
		int c=pow_mod(n,mod-2);
		for(int i=0;i<n;++i)A[i]=1LL*A[i]*c%mod;
	}return;
}
void Get_dao(int n){
	f[n]=0;
	for(int i=0;i<n;++i)f[i]=1LL*f[i+1]*(i+1)%mod; 
}
void Get_Fen(int n){
	for(int i=n-1;i>=1;--i)ln[i]=1LL*ln[i-1]*ni[i]%mod;
	ln[0]=0;
}
void Get_inv(int n){
	if(n==1){
		Inv[0]=pow_mod(f[0],mod-2);
		return;
	}
	Get_inv(n>>1);
	int now=(n<<1);
	for(len=0;(1<<len)<now;++len);
	for(int i=0;i<now;++i)rev[i]=rev[i>>1]>>1|((i&1)<<(len-1));
	static int tmp[maxn];
	for(int i=0;i<n;++i)tmp[i]=f[i];
	for(int i=n;i<now;++i)tmp[i]=0;
	FFT(tmp,now,1);FFT(Inv,now,1);
	for(int i=0;i<now;++i)Inv[i]=1LL*Inv[i]*(2LL-1LL*tmp[i]*Inv[i]%mod+mod)%mod;
	FFT(Inv,now,-1);
	memset(Inv+n,0,sizeof(int)*n);
} 
void Get_ln(int n){
	int now=(n<<1);
	for(int i=0;i<n;++i)Inv[i]=0;
	Get_inv(n);Get_dao(n);
	for(len=0;(1<<len)<now;++len);
	for(int i=0;i<now;++i)rev[i]=rev[i>>1]>>1|((i&1)<<(len-1));
	FFT(Inv,now,1);FFT(f,now,1);
	for(int i=0;i<now;++i)ln[i]=1LL*f[i]*Inv[i]%mod;
	FFT(ln,now,-1);
	memset(ln+n,0,sizeof(int)*n);
	Get_Fen(n);
}
void Get_exp(int n){
	if(n==1){Exp[0]=1;return;}
	Get_exp(n>>1);
	int now=(n<<1);
	for(int i=0;i<n;++i)f[i]=Exp[i],h[i]=g[i];
	Get_ln(n);
	for(int i=0;i<now;++i)ln[i]=((i==0)-ln[i]+h[i]+mod)%mod;
	for(len=0;(1<<len)<now;++len);
	for(int i=0;i<now;++i)rev[i]=rev[i>>1]>>1|((i&1)<<(len-1));
	FFT(Exp,now,1);FFT(ln,now,1);
	for(int i=0;i<now;++i)Exp[i]=1LL*Exp[i]*ln[i]%mod;
	FFT(Exp,now,-1);
	memset(Exp+n,0,sizeof(int)*n);
}

int main(){
	freopen("QAQ_Bell.in","r",stdin);
	freopen("QAQ_Bell.ans","w",stdout);
	scanf("%d",&T);
	for(int i=1;i<=T;++i){
		scanf("%d",&Q[i].n);
		n=max(n,Q[i].n);
	}
	for(N=1;N<=n;N<<=1);
	jc[0]=1;ni[0]=ni[1]=1;
	for(int i=1;i<N;++i)jc[i]=1LL*jc[i-1]*i%mod;
	inv[N-1]=pow_mod(jc[N-1],mod-2);
	for(int i=N-2;i>=0;--i)inv[i]=1LL*inv[i+1]*(i+1)%mod;
	for(int i=2;i<N;++i)ni[i]=(mod-1LL*(mod/i)*ni[mod%i]%mod);
	for(int i=1;i<N;++i)g[i]=inv[i];
	Get_exp(N);
	for(int i=1;i<=T;++i){
		n=Q[i].n;
		printf("%d\n",1LL*Exp[n]*jc[n]%mod);
	}
	return 0;
}

多项式开根

同样是利用倍增

设G^2(x)=F(x)(mod x^n)

我们得到(G^2(x)-F(x))^2=0(mod x^2n)

又可以得到(G^2(x)+F(x))^2=4*G^2(x)*F(x)(mod x^2n)

即((G^2(x)+F(x))/2*G(x))^2=F(x)(mod x^2n)

显然括号里面的东西是我们要求的东西，求解即可

注意到多项式开根的前提是常数项可开根，但如果在模意义下并不要求常数项是完全平方数

只要常数项存在二次剩余就可以了，至于求解的话可以利用原根的性质求解

BZOJ 3625

我们设答案多项式为F(x)

我们得到生成函数F(x)=C(x)F^2(x)+1

得到C(x)F^2(x)-F(x)+1=0

F(x)=(1 +or- sqrt(1-4C(x))/2C(x))=2/(1 +or- sqrt(1-4C(x)))

由于多项式求逆要求常数项存在逆，所以可以确定是加号

之后做多项式开根之后多项式求逆即可

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstdlib>
#include<cstring>
#define G 3
using namespace std;
 
const int maxn=300010;
const int mod=998244353;
int n,m,x,N,len,C;
int f[maxn],rt[maxn],Inv[maxn],rev[maxn];
int w[maxn];
 
void read(int &num){
    num=0;char ch=getchar();
    while(ch<'!')ch=getchar();
    while(ch>='0'&&ch<='9')num=num*10+ch-'0',ch=getchar();
}
int pow_mod(int v,int p){
    int tmp=1;
    while(p){
        if(p&1)tmp=1LL*tmp*v%mod;
        v=1LL*v*v%mod;p>>=1;
    }return tmp;
}
void FFT(int *A,int n,int flag){
    for(int i=0;i<n;++i)if(i<rev[i])swap(A[i],A[rev[i]]);
    for(int k=2;k<=n;k<<=1){
        int mk=(k>>1);
        int wn=pow_mod(G,flag==1?(mod-1)/k:(mod-1)-(mod-1)/k);
        w[0]=1;
        for(int i=1;i<mk;++i)w[i]=1LL*w[i-1]*wn%mod;
        for(int i=0;i<n;i+=k){
            for(int j=0;j<mk;++j){
                int a=i+j,b=i+j+mk;
                int x=A[a],y=1LL*A[b]*w[j]%mod;
                A[a]=(x+y)%mod;
                A[b]=x-y;if(A[b]<0)A[b]+=mod;
            }
        }
    }
    if(flag==-1){
        int c=pow_mod(n,mod-2);
        for(int i=0;i<n;++i)A[i]=1LL*A[i]*c%mod;
    }return;
}
void Get_inv(int n){
    if(n==1){
        Inv[0]=pow_mod(rt[0],mod-2);
        return;
    }
    Get_inv(n>>1);
    int now=(n<<1);
    for(len=0;(1<<len)<now;++len);
    for(int i=0;i<now;++i)rev[i]=rev[i>>1]>>1|((i&1)<<(len-1));
    static int tmp[maxn];
    for(int i=0;i<n;++i)tmp[i]=rt[i];
    for(int i=n;i<now;++i)tmp[i]=0;
    FFT(tmp,now,1);FFT(Inv,now,1);
    for(int i=0;i<now;++i)Inv[i]=1LL*Inv[i]*(2LL-1LL*tmp[i]*Inv[i]%mod+mod)%mod;
    FFT(Inv,now,-1);
    memset(Inv+n,0,sizeof(int)*n);
}
void Get_root(int n){
    if(n==1){rt[0]=1;return;}
    Get_root(n>>1);
    for(int i=0;i<n;++i)Inv[i]=0;
    int now=(n<<1);Get_inv(n);
    static int tmp[maxn];
    for(int i=0;i<n;++i)tmp[i]=f[i];
    for(len=0;(1<<len)<now;++len);
    for(int i=0;i<now;++i)rev[i]=rev[i>>1]>>1|((i&1)<<(len-1));
    FFT(Inv,now,1);FFT(tmp,now,1);
    for(int i=0;i<now;++i)tmp[i]=1LL*tmp[i]*Inv[i]%mod;
    FFT(tmp,now,-1);
    for(int i=0;i<n;++i)rt[i]=1LL*C*(rt[i]+tmp[i])%mod;
}
 
int main(){
    read(n);read(m);C=((mod+1)>>1);
    for(N=1;N<=m;N<<=1);
    for(int i=1;i<=n;++i){
        read(x);
        if(x<=m)f[x]++;
    }
    for(int i=0;i<N;++i)f[i]=((i==0)-4*f[i]+mod)%mod;
    Get_root(N);rt[0]++;
    for(int i=0;i<(N<<1);++i)Inv[i]=0;
    Get_inv(N);
    for(int i=1;i<=m;++i){
        printf("%d\n",(Inv[i]<<1)%mod);
    }return 0;
}

多项式快速幂

即求F^k(x)

我们可以直接做快速幂每次乘完消掉次数界既可以了

但是这样的做法不够高大上

我们知道F^k(x)=exp(ln(F^k(x))

而ln(F^k(x))=k*ln(F(x))

我们做多项式求ln，然后系数乘以k之后exp还原即可

时间复杂度O(nlogn)

 

多项式的学习估计就到此结束啦，以后可能会学一些用FFT优化DP，用FFT优化字符串匹配的题目

当然也会写总结啦

听说多项式还有什么多点插值，扩展欧几里得，拉格朗日反演之类的鬼东西

先暂时弃掉咯

最后附上多项式终极模板，也是zcg他们出的毒瘤题目

cojs 2189 帕秋莉的超级多项式

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<iostream>
#include<cmath>
#define G 3
using namespace std;

const int maxn=300010;
const int mod=998244353;
int n,k,N,C,len;
int rev[maxn],w[maxn];
int f[maxn],rt[maxn];
int Inv[maxn],ln[maxn],Exp[maxn];

int pow_mod(int v,int p){
	int tmp=1;
	while(p){
		if(p&1)tmp=1LL*tmp*v%mod;
		v=1LL*v*v%mod;p>>=1;
	}return tmp;
}
void FFT(int *A,int n,int flag){
	for(int i=0;i<n;++i)if(i<rev[i])swap(A[i],A[rev[i]]);
	for(int k=2;k<=n;k<<=1){
		int mk=(k>>1);
		int wn=pow_mod(G,flag==1?(mod-1)/k:(mod-1)-(mod-1)/k);
		w[0]=1;
		for(int i=1;i<mk;++i)w[i]=1LL*w[i-1]*wn%mod;
		for(int i=0;i<n;i+=k){
			for(int j=0;j<mk;++j){
				int a=i+j,b=i+j+mk;
				int x=A[a],y=1LL*A[b]*w[j]%mod;
				A[a]=(x+y)%mod;
				A[b]=x-y;if(A[b]<0)A[b]+=mod;
			}
		}
	}
	if(flag==-1){
		int c=pow_mod(n,mod-2);
		for(int i=0;i<n;++i)A[i]=1LL*A[i]*c%mod;
	}return;
}
void Get_Fen(int *f,int n){
	for(int i=n-1;i>=1;--i)f[i]=1LL*f[i-1]*pow_mod(i,mod-2)%mod;
	f[0]=0;
}
void Get_dao(int *f,int n){
	f[n]=0;
	for(int i=0;i<n;++i)f[i]=1LL*f[i+1]*(i+1)%mod;
}
void Get_inv(int n,int *h){
	if(n==1){
		Inv[0]=pow_mod(h[0],mod-2);
		return;
	}
	Get_inv(n>>1,h);
	int now=(n<<1);
	static int tmp[maxn];
	for(int i=0;i<n;++i)tmp[i]=h[i];
	for(int i=n;i<now;++i)tmp[i]=0;
	for(len=0;(1<<len)<now;++len);
	for(int i=0;i<now;++i)rev[i]=rev[i>>1]>>1|((i&1)<<(len-1));
	FFT(tmp,now,1);FFT(Inv,now,1);
	for(int i=0;i<now;++i)Inv[i]=1LL*Inv[i]*(2LL-1LL*tmp[i]*Inv[i]%mod+mod)%mod;
	FFT(Inv,now,-1);
	memset(Inv+n,0,sizeof(int)*n);
}
void Get_root(int n){
	if(n==1){
		rt[0]=(int)(sqrt(f[0]));
		return;
	}
	Get_root(n>>1);
	int now=(n<<1);
	for(int i=0;i<now;++i)Inv[i]=0;
	Get_inv(n,rt);
	static int tmp[maxn];
	for(int i=0;i<n;++i)tmp[i]=f[i];
	for(int i=n;i<now;++i)tmp[i]=0;
	for(len=0;(1<<len)<now;++len);
	for(int i=0;i<now;++i)rev[i]=rev[i>>1]>>1|((i&1)<<(len-1));
	FFT(tmp,now,1);FFT(Inv,now,1);
	for(int i=0;i<now;++i)tmp[i]=1LL*tmp[i]*Inv[i]%mod;
	FFT(tmp,now,-1);
	for(int i=0;i<n;++i)rt[i]=1LL*C*(rt[i]+tmp[i])%mod;
}
void Get_ln(int n,int *h){
	int now=(n<<1);
	for(int i=0;i<now;++i)Inv[i]=0;
	Get_inv(n,h);Get_dao(h,n);
	for(len=0;(1<<len)<now;++len);
	for(int i=0;i<now;++i)rev[i]=rev[i>>1]>>1|((i&1)<<(len-1));
	FFT(h,now,1);FFT(Inv,now,1);
	for(int i=0;i<now;++i)ln[i]=1LL*h[i]*Inv[i]%mod;
	FFT(ln,now,-1);
	memset(ln+n,0,sizeof(int)*n);
	Get_Fen(ln,n);
}
void Get_Exp(int n){
	if(n==1){Exp[0]=1;return;}
	Get_Exp(n>>1);
	int now=(n<<1);
	static int g[maxn];
	for(int i=0;i<n;++i)g[i]=Exp[i];
	for(int i=n;i<now;++i)g[i]=0;
	Get_ln(n,g);
	for(int i=0;i<n;++i)ln[i]=((i==0)-ln[i]+f[i]+mod)%mod;
	for(len=0;(1<<len)<now;++len);
	for(int i=0;i<now;++i)rev[i]=rev[i>>1]>>1|((i&1)<<(len-1));
	FFT(ln,now,1);FFT(Exp,now,1);
	for(int i=0;i<now;++i)Exp[i]=1LL*Exp[i]*ln[i]%mod;
	FFT(Exp,now,-1);
	memset(Exp+n,0,sizeof(int)*n);
}

int main(){
	freopen("polynomial.in","r",stdin);
	freopen("polynomial.out","w",stdout);
	scanf("%d",&n);scanf("%d",&k);C=((mod+1)>>1);
	for(int i=0;i<n;++i)scanf("%d",&f[i]);
	for(N=1;N<=n;N<<=1);
	Get_root(N);
	for(int i=0;i<(N<<1);++i)Inv[i]=0;
	Get_inv(N,rt);
	for(int i=0;i<N;++i)f[i]=Inv[i];
	Get_Fen(f,n);
	for(int i=n;i<N;++i)f[i]=0;
	Get_Exp(N);
	for(int i=0;i<N;++i)f[i]=Exp[i];
	for(int i=0;i<(N<<1);++i)Inv[i]=0;
	Get_inv(N,f);
	for(int i=0;i<N;++i)f[i]=Inv[i];
	f[0]++;
	Get_ln(N,f);
	for(int i=0;i<N;++i)f[i]=ln[i];
	f[0]++;
	for(int i=N+1;i<(N<<1);++i)f[i]=0;
	Get_ln(N,f);
	for(int i=0;i<N;++i)f[i]=1LL*k*ln[i]%mod;
	for(int i=0;i<(N<<1);++i)Exp[i]=0;
	Get_Exp(N);
	for(int i=1;i<n;++i)printf("%d ",1LL*Exp[i]*i%mod);
	printf("0\n");
	return 0;
}

感觉生成函数什么的还是很好玩的，非常想找时间学一学

不过貌似NOI用处不大(雾

转载于:https://www.cnblogs.com/joyouth/p/5604618.html

weixin_30737363
关注
多项式exp-ln学习小结
fengqiyuka的博客
03-29 964
前言 其实多项式exp/ln好像考的频率不算很高 不过总得来说只要多做生成函数计数题的话，多多少少还是会遇到的。 ln(x)ln(x)ln(x)的展开式 通过泰勒展开，可以得到： −ln(x)=x1+x22+x33+... -ln(x)=\frac{x}{1}+\frac{x^2}{2}+\frac{x^3}{3}+... −ln(x)=1x​+2x2​+3x3​+... ln(x)ln(x)ln(x)的导数 ln(x)ln(x)ln(x)有一个特殊的性质，就是如果f(x)=ln(x)f(x
牛顿迭代，多项式求逆，除法，开方，exp，ln，求幂
Cyhlnj
05-31 1915
多项式求逆 给定A(x)A(x)A(x)求满足A(x)∗B(x)=1A(x)∗B(x)=1A(x)*B(x)=1的B(x)B(x)B(x) 写成 A(x)∗B(x)=1(mod&amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;nbsp;xn)A(x)∗B(x)=1(mod&amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;nbsp;xn)A(x)*B(x)=1(mod \ x^n) 我们会求A(x)∗B(x)
多项式求ln
gigo_64 sink in my world
08-06 206
求 两边求导 所以对于g’求导，乘g的逆元，然后得到f'，再积分回去。 关于求导： 关于积分： 关于代码 #include<bits/stdc++.h> using namespace std; #define in read() #define int long long int in{ int cnt=0,f=1;char ch=0; while(!isdigit...
BZOJ3456 城市规划 【多项式求ln】
weixin_30268921的博客
06-15 116
题目链接 BZOJ3456 题解 真是一道经典好题，至此已经写了分治\(NTT\)，多项式求逆，多项式求\(ln\)三种写法 我们发现我们要求的是大小为\(n\)无向联通图的数量 而\(n\)个点的无向图是由若干个无向联通图构成的 那么我们设\(F(x)\)为无向联通图数量的指数型生成函数 设\(G(x)\)为无向图数量的指数型生成函数 \(G(x)\)很好求 而 \[G(x) = \frac{F...
指数型生成函数 及 多项式求ln
weixin_30515513的博客
06-15 288
指数型生成函数 我们知道普通型生成函数解决的是组合问题，而指数型生成函数解决的是排列问题 对于数列\(\{a_n\}\)，我们定义其指数型生成函数为 \[G(x) = a_0 + a_1x + a_2\frac{x^2}{2!} + a_3\frac{x^3}{3!} + a_4\frac{x^4}{4!} + \dots = \sum\limits_{i = 0}^{\infty} a_i\f...
多项式小结 （求逆）（ln）（exp）（除法）（多点求值）（快速插值）（任意模数NTT）
FSYo 的博客
11-07 688
多项式求逆 f(x)∗g(x)=1(mod xn)f(x)*g(x)=1(mod\ x^n)f(x)∗g(x)=1(mod xn) 先求出 f(x)∗g(x)′=1(mod xn2)f(x)*g(x)'=1(mod\ x^{\frac{n}{2}})f(x)∗g(x)′=1(mod x2n​) 然后有 f(x)∗g(x)′=f(x)∗g(x)(mod&nbs...
[学习笔记]多项式的各种操作
C202044zxy的博客
12-20 417
多项式求逆 点此看题 可以暴力递推，时间复杂度O(n2)O(n^2)O(n2)。 倍增可以将之优化到O(nlog⁡n)O(n\log n)O(nlogn)，我们先从111开始倍增，假设现在倍增到了ppp，得到了模xpx^pxp的逆元B′(x)B'(x)B′(x)，现在我们要求B(x)B(x)B(x)，显然这时候求出来的B′(x)B'(x)B′(x)是B(x)B(x)B(x)的一部分，有下列关系式：...
多项式小结
m0_49959202的博客
07-23 628
多项式 文章目录多项式一 简介二 拉格朗日插值2.1 知识点2.2 模板2.1 随机给定n个点插值：O(n2logn)O(n^2log_n)O(n2logn​)2.2 给定连续n个点插值:O(n)O(n)O(n)2.3 例题三 快速傅里叶变换3.1 知识点3.2 模板3.3 例题四 快速数论变换4.1 知识点4.2 模板五 多项式求逆六 多项式开方七 多项式除法|取模八 多项式对数函数|指数函数九 牛顿迭代法十 多项式三角函数十二 多项式反三角函数 一 简介 **多项式的度：**对于一个多项式f(x)f(x
【总结】多项式相关的各种算法
远行客
08-22 4564
多项式相关的各种算法
【学习笔记】超简单的多项式三角函数（含全套证明）
繁凡さん的博客
01-09 1650
整理的算法模板合集： ACM模板 目录P5264 多项式三角函数 点我看多项式全家桶(●^∀^●) P5264 多项式三角函数 如果模的是任意p的话可以枚举i，或者使用Cipolla算法。 //#pragma GCC optimize(2) #include <bits/stdc++.h> using namespace std; typedef long long ll; const int N = 3000007; const int p = 998244353, gg = 3, i
并不对劲的多项式求ln,exp
weixin_30851409的博客
03-13 102
ln 解释 设\(g(x)=ln(f(x))\)，两边同时求导，则有：\(g'(x)=ln'(f(x))*f'(x)=f^{-1}(x)*f'(x)\)(1) 因为\(f(x)\)是个多项式，所以设\(f(x)=\sum_{i=0}^{n}a_i*x^i\)，则有\(f'(x)=\sum_{i=0}^{n-1}a_{i+1}*(i+1)*x^i\) 设\(h(x)=f^{-1}(x)*f'(x)\...
多项式对数函数（多项式求ln）小结
a_forever_dream的博客
05-19 1839
多项式对数函数 顾名思义，给出一个多项式 A(x)A(x)A(x)，求出一个多项式 G(x)G(x)G(x)，满足 G(x)=ln⁡A(x)(modxn)G(x)=\ln A(x)\pmod {x^n}G(x)=lnA(x)(modxn)。 设 F(x)=ln⁡(x)F(x)=\ln (x)F(x)=ln(x)，那么有 G(x)=F(A(x))G(x)=F(A(x))G(x)=F(A(x))，对等式两边同时求导，得到 G′(x)=F′(A(x))A′(x) G'(x)=F'(A(x))A'(x) G′(x)
【Luogu4238】多项式取Ln模板题
დ♂Hany01's Blog Space♂ღ
07-28 261
Description 多项式取lnln\ln。 Solution lnf(x)=∫f′(x)f−1(x)ln⁡f(x)=∫f′(x)f−1(x)\ln f(x)=\int f'(x)f^{-1}(x) Code /************************************************ * Au: Hany01 * Date: Jul 2...
[多项式ln][多项式exp][多项式求幂][生成函数][DP][FNT] BZOJ 3684: 大朋友和多叉树
Vectorxj的博客
03-06 822
SolutionSolutionSolution 把DP写成生成函数的形式。f(x)=x+∑d∈Dfd(x)f(x)=x+∑d∈Dfd(x)f(x)=x+\sum_{d\in D}f^d(x)设g(f(x))=xg(f(x))=xg\big(f(x)\big)=x，有g(f(x))g(x)==f(x)−∑d∈Dfd(x)x−∑d∈Dxdg(f(x))=f(x)−∑d∈Dfd(x)g(x)=x−∑...
多项式ln,exp学习小计
Deep_Thinking的博客
06-10 787
主要作为一个模板和总结 多项式ln 求G(x)=ln(F(x))G(x)=ln(F(x))G(x)=ln(F(x)) 两边求导 G′(x)=F′(x)F(x)G'(x)=\frac{F'(x)}{F(x)}G′(x)=F(x)F′(x)​，注意这是复合函数求导。 多项式求逆，求导，乘法即可求出G′(x)G'(x)G′(x) 再积分回去就是G(x)G(x)G(x) 一般会钦定常数项为1，否则ln的时候先求导再积分（这里是不定积分）你会发现常数项丢失了，得到的结果并不是原来的多项式，而是在每一项的基础上除以了
洛谷 4389 付公主的背包——多项式求ln、exp
weixin_33973600的博客
02-18 128
题目：https://www.luogu.org/problemnew/show/P4389 关于泰勒展开： https://blog.csdn.net/SoHardToNamed/article/details/80550935 https://www.cnblogs.com/guo-xiang/p/6662881.html 大概就是：\( f(x) = \sum\limits_{i=0...
多项式全家桶（三）：多项式的ln，exp
周道-Althen的博客
04-03 1150
      &ThinSpace;\ \ \ \ \ \ \,      对数和指数是很重要的东西了，复杂的多项式都和他们有关系，所以说掌握很重要，这里不建议光背板子，因为这两个板子都有致命的限制，而在实际操作的时候，可以通过一些方法绕过这个限制直接求解，这个就很重要了。...
[LOJ#6268] 分拆数 [多项式求ln][多项式exp]
最后一次修改
04-24 262
Link LOJ - https://loj.ac/problem/6268 http://danihao123.is-programmer.com/posts/213392.html
[多项式ln][多项式exp][背包DP][生成函数] LOJ #556. 咱们去烧菜吧
Vectorxj的博客
03-06 894
SolutionSolutionSolution 就是个背包DP咯。 答案长这个样子ans=[xn]∏i=1m1−xaibi+ai1−xaians=[xn]∏i=1m1−xaibi+ai1−xaians=[x^n]\prod_{i=1}^m{1 - x^{a_ib_i+a_i}\over 1-x^{a_i}}考虑取个lnln\lnlnans=∑i=1m(ln(1−xaibi+ai)−ln(1−...
多项式求根与数值计算算法解析
资源摘要信息:"多项式求根" 在数学和计算机科学领域，多项式求根是一个基础而重要的问题，涉及找出一个多项式方程的所有根。根指的是使得多项式等于零的数值。一个n次多项式最多可以有n个根，这些根可能是实数或者...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值