【模板】多项式全家桶

最新推荐文章于 2023-01-18 16:33:47 发布
最新推荐文章于 2023-01-18 16:33:47 发布
阅读量141 收藏
点赞数
原文链接：http://www.cnblogs.com/BAJimH/p/10575044.html
版权

FFT

struct Z
{
    double x,y;
    Z(double _x=0,double _y=0):x(_x),y(_y){};
};
Z operator +(Z x,Z y) {return Z(x.x+y.x,x.y+y.y);}
Z operator -(Z x,Z y) {return Z(x.x-y.x,x.y-y.y);}
Z operator *(Z x,Z y) {return Z(x.x*y.x-x.y*y.y,x.x*y.y+x.y*y.x);}
int bit[M],n,m;
Z wi[M+1],a[M],b[M];
void prp()
{
    fo(i,0,M-1) bit[i]=(bit[i>>1]>>1)|((i&1)<<(L-1));
    fo(i,0,M) wi[i]=Z(cos(2*pi*i/M),sin(2*pi*i/M));
}
void DFT(Z *a,int pd)
{
    fo(i,0,M-1) if(i<bit[i]) swap(a[i],a[bit[i]]);
    int half=1,lim=(M>>1);
    Z v;
    for(int m=2;m<=M;half=m,m<<=1,lim>>=1)
    {
        fo(i,0,half-1)
        {
            Z w=(pd==1)?wi[i*lim]:wi[M-i*lim];
            for(int j=i;j<M;j+=m)
            {
                v=w*a[j+half];
                a[j+half]=a[j]-v;
                a[j]=a[j]+v;
            }
        }
    }
    if(pd<0) fo(i,0,M-1) a[i].x/=M;
}
int main()
{
    cin>>n>>m;
    fo(i,0,n) scanf("%lf",&a[i].x);
    fo(i,0,m) scanf("%lf",&b[i].x);
    prp();
    DFT(a,1),DFT(b,1);
    fo(i,0,M-1) a[i]=a[i]*b[i];
    DFT(a,-1);
    fo(i,0,n+m) 
    {
        int v=(int)(a[i].x+0.5);
        printf("%d ",v);
    }
}

NTT

inline LL md(LL x)
{
    return(x<0)?(x+mo):((x>=mo)?x-mo:x);
}
void NTT(LL *a,int pd,int num)
{
    prp(num);
    fo(i,0,num-1) if(i<bit[i]) swap(a[i],a[bit[i]]);
    int lim=num>>1,half=1;
    LL v;
    for(int m=2;m<=num;half=m,m<<=1,lim>>=1)
    {
        fo(i,0,half-1)
        {
            LL w=(pd==1)?wi[i*lim]:wi[num-i*lim];
            for(int j=i;j<num;j+=m)
            {
                v=(w*a[j+half])%mo;
                a[j+half]=md(a[j]-v);
                a[j]=md(a[j]+v);
            }
        }
    }
    if(pd<0) fo(i,0,num-1) a[i]=a[i]*ny[num]%mo;
}

多项式求逆

#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <cmath>
#include <algorithm>
#include <cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fod(i,a,b) for(int i=a;i>=b;i--)
#define M 262144
#define mo 998244353
#define LL long long
#define N 100005
using namespace std;
LL a[M+1],b[M+1],c[M+1],wi[M+1],wg[M+1],ny2,d[M+1];
int bit[M+1],n;
LL ksm(LL k,LL n)
{
    LL s=1;
    for(;n;n>>=1,k=k*k%mo) if(n&1) s=s*k%mo;
    return s;
}
void prp(int num,int L)
{
    ny2=ksm(num,mo-2);
    fo(i,0,num) wi[i]=wg[i*(M/num)];
    fo(i,0,num-1) bit[i]=(bit[i>>1]>>1)|((i&1)<<(L-1));
}
void NTT(LL *a,bool pd,int num)
{
    fo(i,0,num-1) if(bit[i]<i) swap(a[bit[i]],a[i]);
    int half=1,lim=num>>1;
    LL v;
    for(int m=2;m<=num;half=m,m<<=1,lim>>=1)
    {
        for(register int i=0;i<half;i++)
        {
            LL w=(!pd)?wi[i*lim]:wi[num-i*lim];
            for(register int j=i;j<num;j+=m)
            {
                v=a[j+half]*w%mo;
                a[j+half]=(a[j]-v+mo)%mo;
                a[j]=(a[j]+v)%mo;
            }
        }
    }   
    if(pd) fo(i,0,num-1) a[i]=a[i]*ny2%mo;
}
void make(int l,LL *a,LL *b)
{
    b[0]=ksm(a[0],mo-2);
    for(int m=1,t=2,num=4,cnt=2;m<l;m=t,t=num,num<<=1,cnt++)
    {
        prp(num,cnt);
        fo(i,0,m-1) c[i]=a[i],d[i]=b[i];
        fo(i,m,t-1) c[i]=a[i];
        fo(i,t,num-1) c[i]=0;
        NTT(c,0,num),NTT(b,0,num);
        fo(i,0,num-1) b[i]=b[i]*b[i]%mo*c[i]%mo;
        NTT(b,1,num);
        fo(i,0,t-1) b[i]=((LL)2*d[i]-b[i]+mo)%mo;
        fo(i,t,num-1) b[i]=0;
    }       
}
int main()
{
    cin>>n;
    fo(i,0,n-1) scanf("%lld",&a[i]);
    wg[0]=1;
    LL c=ksm(3,(mo-1)/M);
    fo(i,1,M) wg[i]=wg[i-1]*c%mo;
    make(n,a,b);
    fo(i,0,n-1) printf("%lld ",b[i]);
}

多项式除法（多项式取模）

#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <cmath>
#include <algorithm>
#include <cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fod(i,a,b) for(int i=a;i>=b;i--)
#define M 524288
#define mo 998244353
#define LL long long
#define N 100005
using namespace std;
LL a[M+1],b[M+1],c[M+1],wi[M+1],wg[M+1],ny2,d[M+1],r[M+1],q[M+1],l2[M+1],cf[20],f[M+1];
int bit[M+1],n,m;
LL ksm(LL k,LL n)
{
    LL s=1;
    for(;n;n>>=1,k=k*k%mo) if(n&1) s=s*k%mo;
    return s;
}
void prp(int num)
{
    ny2=ksm(num,mo-2);
    fo(i,0,num) wi[i]=wg[i*(M/num)];
    fo(i,0,num-1) bit[i]=(bit[i>>1]>>1)|((i&1)<<(l2[num]-1));
}
void NTT(LL *a,bool pd,int num)
{
    fo(i,0,num-1) if(bit[i]<i) swap(a[bit[i]],a[i]);
    int half=1,lim=num>>1;
    LL v;
    for(int m=2;m<=num;half=m,m<<=1,lim>>=1)
    {
        for(register int i=0;i<half;i++)
        {
            LL w=(!pd)?wi[i*lim]:wi[num-i*lim];
            for(register int j=i;j<num;j+=m)
            {
                v=a[j+half]*w%mo;
                a[j+half]=(a[j]-v+mo)%mo;
                a[j]=(a[j]+v)%mo;
            }
        }
    }   
    if(pd) fo(i,0,num-1) a[i]=a[i]*ny2%mo;
}
void make(int l,LL *a,LL *b)
{
    b[0]=ksm(a[0],mo-2);
    for(int m=1,t=2,num=4;m<l;m=t,t=num,num<<=1)
    {
        prp(num);
        fo(i,0,m-1) c[i]=a[i],d[i]=b[i];
        fo(i,m,t-1) c[i]=a[i];
        fo(i,t,num-1) c[i]=0;
        NTT(c,0,num),NTT(b,0,num);
        fo(i,0,num-1) b[i]=b[i]*b[i]%mo*c[i]%mo;
        NTT(b,1,num);
        fo(i,0,t-1) b[i]=((LL)2*d[i]-b[i]+mo)%mo;
        fo(i,t,num-1) b[i]=0;
    }       
}
void rev(int num,LL *a,LL *b)
{
    fo(i,0,num-1) b[i]=a[num-i-1];
}
void div(LL *a,LL *b,LL *d,LL *r)
{
    rev(m+1,b,r);
    make(n-m+1,r,d);
    int num=cf[l2[n+m+1]];
    prp(num);
    fo(i,0,n) f[i]=a[n-i];
    fo(i,n-m+1,num-1) f[i]=0,d[i]=0;
    NTT(f,0,num),NTT(d,0,num);
    fo(i,0,num-1) d[i]=d[i]*f[i]%mo;
    NTT(d,1,num);
    fo(i,n-m+1,num-1) d[i]=0;
    fo(i,0,(n-m)>>1) swap(d[i],d[n-m-i]);
    fo(i,0,num-1) r[i]=d[i],f[i]=b[i];
    num=cf[l2[n+1]];
    prp(num);
    NTT(f,0,num),NTT(r,0,num);
    fo(i,0,num-1) r[i]=r[i]*f[i]%mo;
    NTT(r,1,num);
    fo(i,0,num-1) r[i]=(a[i]-r[i]+mo)%mo;
}
int main()
{
    cf[0]=1;
    fo(i,1,19) cf[i]=cf[i-1]*2,l2[cf[i]]=i;
    cin>>n>>m;
    fo(i,0,n) scanf("%lld",&a[i]);
    fo(i,0,m) scanf("%lld",&b[i]);
    wg[0]=1;
    LL c=ksm(3,(mo-1)/M);
    fod(i,M-1,1) if(!l2[i]) l2[i]=l2[i+1];
    fo(i,1,M) wg[i]=wg[i-1]*c%mo;
    div(a,b,q,r);
    fo(i,0,n-m) printf("%lld ",q[i]);
    printf("\n");
    fo(i,0,m-1) printf("%lld ",r[i]);
}

多点求值（洛谷 P5050）

#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <cmath>
#include <cstring>
#define fo(i,a,b) for(int i=a;i<=b;++i)
#define fod(i,a,b) for(int i=a;i>=b;--i)
#define M 262144
#define mo 998244353
#define N 100005
#define LL long long
using namespace std;
LL a[M+1],b[M+1],c[M+1],u1[M+1],wg[M+1],u2[M+1],wi[M+1],r[M+1],q[M+1],cf[19],l2[M+1],ny[M+1],a1[20][M+1],ans[N],s1[20][M+1];
int bit[M+1],n,m,fi[20][N+1],le[20][N+1],lf;
LL ksm(LL k,LL n)
{
    LL s=1;
    for(;n;n>>=1,k=k*k%mo) if(n&1) s=s*k%mo;
    return s;
}
void prp(int num)
{
    fo(i,0,num) wi[i]=wg[i*(M/num)];
    fo(i,0,num-1) bit[i]=(bit[i>>1]>>1)|((i&1)<<(l2[num]-1));
}
void NTT(LL *a,bool pd,int num)
{
    fo(i,0,num-1) if(i<bit[i]) swap(a[i],a[bit[i]]);
    LL w,v;
    for(int lim=num>>1,m=2,half=1;m<=num;half=m,lim>>=1,m<<=1)
    {
        fo(i,0,half-1)
        {
            w=(!pd)?wi[i*lim]:wi[num-i*lim];
            for(int j=i;j<num;j+=m)
            {
                v=w*a[j+half]%mo;
                a[j+half]=(a[j]-v+mo)%mo;
                a[j]=(a[j]+v)%mo;
            }
        }
    }
    if(pd) fo(i,0,num-1) a[i]=a[i]*ny[num]%mo;
}
void make(int l,LL *a,LL *b)
{
    b[0]=ksm(a[0],mo-2);
    for(int m=1,t=2,num=4;m<=l;m=t,t=num,num<<=1)
    {
        prp(num);
        fo(i,0,m-1) c[i]=a[i],u2[i]=b[i];
        fo(i,m,t-1) c[i]=a[i],u2[i]=b[i]=0;
        fo(i,t,num-1) c[i]=u2[i]=b[i]=0;
        NTT(c,0,num),NTT(b,0,num);
        fo(i,0,num-1) b[i]=b[i]*b[i]%mo*c[i]%mo;
        NTT(b,1,num);
        fo(i,0,t-1) b[i]=((LL)2*u2[i]-b[i]+mo)%mo;
        fo(i,t,num-1) b[i]=0;
    }
}
void rev(int num,LL *a,LL *b)
{
    fo(i,0,num-1) b[i]=a[num-i-1];
}
void div(int n,int m,LL *a,LL *b,LL *d,LL *r)
{
    rev(m,b,r);
    fo(i,m,n-m+1) r[i]=0;
    make(n-m+1,r,d);
    int num=cf[l2[n]]*2;
    prp(num),rev(n,a,u2);
    fo(i,n-m+1,num-1) d[i]=u2[i]=0;
    
    NTT(d,0,num),NTT(u2,0,num);
    fo(i,0,num-1) d[i]=d[i]*u2[i]%mo;
    NTT(d,1,num);
    fo(i,n-m+1,num-1) d[i]=0;
    fo(i,0,(n-m)>>1) swap(d[i],d[n-m-i]);

    num=cf[l2[n]],prp(num);
    fo(i,0,m-1) r[i]=d[i],u2[i]=b[i];
    fo(i,m,num-1) r[i]=d[i],u2[i]=0;
    NTT(r,0,num),NTT(u2,0,num);
    fo(i,0,num-1) r[i]=r[i]*u2[i]%mo;
    NTT(r,1,num);
    fo(i,0,m-1) r[i]=(a[i]-r[i]+mo)%mo;
    fo(i,m,num-1) r[i]=0;
}
void prd(int t,int l,int r)
{
    if(l==r)
    {
        fi[t][l]=lf+1;
        le[t][l]=2;
        a1[t][fi[t][l]]=-b[l],a1[t][fi[t][l]+1]=1;
        lf+=2;
        return;
    }
    int mid=(l+r)>>1;
    prd(t+1,l,mid),prd(t+1,mid+1,r);
}
void doit(int t,int l,int r)
{
    if(l==r) return;
    int mid=(l+r)>>1;
    doit(t+1,l,mid),doit(t+1,mid+1,r);
    fi[t][l]=fi[t+1][l];
    int num=cf[l2[le[t+1][l]+le[t+1][mid+1]-1]];
    prp(num);
    fo(i,0,num-1)
    {
        u1[i]=(i<le[t+1][l])?a1[t+1][fi[t+1][l]+i]:0;
        u2[i]=(i<le[t+1][mid+1])?a1[t+1][fi[t+1][mid+1]+i]:0;
    }
    NTT(u1,0,num),NTT(u2,0,num);
    fo(i,0,num-1) u1[i]=u1[i]*u2[i]%mo;
    NTT(u1,1,num);
    le[t][l]=le[t+1][l]+le[t+1][mid+1]-1;
    fo(i,0,le[t][l]) a1[t][fi[t][l]+i]=u1[i];
}
void query(int t,int l,int r,int n)
{
    if(l==r) ans[l]=s1[t][0];
    else
    {
        int mid=(l+r)>>1;
        div(n,le[t+1][l],s1[t],a1[t+1]+fi[t+1][l],q,s1[t+1]);
        query(t+1,l,mid,le[t+1][l]-1);
        div(n,le[t+1][mid+1],s1[t],a1[t+1]+fi[t+1][mid+1],q,s1[t+1]);
        query(t+1,mid+1,r,le[t+1][mid+1]-1);
    }
}
int main()
{
    cf[0]=1;
    fo(i,1,18) cf[i]=cf[i-1]<<1,l2[cf[i]]=i;
    fod(i,M-1,1) if(!l2[i]) l2[i]=l2[i+1];
    ny[1]=1;
    fo(i,2,M) ny[i]=(-ny[mo%i]*(LL)(mo/i)%mo+mo)%mo;
    wg[0]=1;
    LL c=ksm(3,(mo-1)/M);
    fo(i,1,M) wg[i]=wg[i-1]*c%mo;
    cin>>n>>m;
    fo(i,0,n) scanf("%lld",&a[i]);
    fo(i,0,m-1) scanf("%lld",&b[i]);
    prd(0,0,m-1);
    doit(0,0,m-1);
    fo(i,0,n) s1[0][i]=a[i];
    query(0,0,m-1,n+1);
    fo(i,0,m-1) printf("%lld\n",ans[i]);
}

多项式开根（洛谷P5205）

#include <bits/stdc++.h>
#define fo(i,a,b) for(int i=a;i<=b;++i)
#define fod(i,a,b) for(int i=a;i>=b;--i)
#define M 524288
#define LL long long
#define mo 998244353
using namespace std;
int l2[M+1],cf[22],n,t,bit[M+1];
LL wi[M+1],wg[M+1],b[M+1],u1[M+1],u2[M+1],u3[M+1],ny[M+1],c[M+1];
void prp(int num)
{
    fo(i,0,num) bit[i]=(bit[i>>1]>>1)|((i&1)<<(l2[num]-1)),wi[i]=wg[i*(M/num)];
}
LL ksm(LL k,LL n)
{
    LL s=1;
    for(;n;n>>=1,k=k*k%mo) if(n&1) s=s*k%mo;
    return s;
}
void NTT(LL *a,bool pd,int num)
{
    fo(i,0,num-1) if(i<bit[i]) swap(a[i],a[bit[i]]);
    LL v;   
    for(register int lim=num>>1,h=1,m=2;m<=num;h=m,lim>>=1,m<<=1)
    {
        for(register int j=0;j<num;j+=m)
        {
            fo(i,0,h-1)
            {
                v=((!pd)?wi[i*lim]:wi[num-i*lim])*a[i+j+h]%mo;  
                a[i+j+h]=(a[i+j]-v+mo)%mo;
                a[i+j]=(a[i+j]+v)%mo;
            }
        }
    }
    if(pd) fo(i,0,num-1) a[i]=a[i]*ny[num]%mo;
}
void inv(int l,LL *a,LL *b)
{
    b[0]=ksm(a[0],mo-2);
    for(int m=1,t=2,num=4;m<l;m=t,t=num,num<<=1)
    {
        prp(num);
        fo(i,0,m-1) u1[i]=a[i],u2[i]=b[i];
        fo(i,m,t-1) u1[i]=a[i],u2[i]=b[i]=0;
        fo(i,t,num-1) u1[i]=0,u2[i]=b[i]=0;
        NTT(u1,0,num),NTT(b,0,num);
        fo(i,0,num-1) b[i]=b[i]*b[i]%mo*u1[i]%mo,u1[i]=0;
        NTT(b,1,num);
        fo(i,0,t-1) b[i]=((LL)2*u2[i]%mo-b[i]+mo)%mo,u2[i]=0;
        fo(i,t,num-1) b[i]=0;
    }       
}
void make(int l,LL *a,LL *b)
{
    b[0]=1;
    for(int m=1,t=2,num=4;m<l;m=t,t=num,num<<=1)
    {
        inv(t,b,u3);
        fo(i,0,t-1) u1[i]=a[i];
        fo(i,t,num-1) u1[i]=u3[i]=0;
        prp(num);
        NTT(u1,0,num),NTT(u3,0,num);
        fo(i,0,num-1) u1[i]=u1[i]*u3[i]%mo,u3[i]=0;
        NTT(u1,1,num);
        fo(i,0,t-1) b[i]=(u1[i]+b[i]+mo)%mo*ny[2]%mo,u1[i]=0;
        fo(i,t,num-1) b[i]=0,u1[i]=0;
    }
}
int main()
{
    cf[0]=1;
    fo(i,1,19) l2[cf[i]=cf[i-1]<<1]=i;
    fod(i,M-1,2) if(!l2[i]) l2[i]=l2[i+1];
    wg[0]=1,wg[1]=ksm(3,(mo-1)/M);
    fo(i,2,M) wg[i]=wg[i-1]*wg[1]%mo;
    ny[1]=1;
    fo(i,2,M) ny[i]=(-ny[mo%i]*(LL)(mo/i)%mo+mo)%mo;
    
    int n;
    cin>>n;
    fo(i,0,n-1) scanf("%lld",&b[i]);
    make(cf[l2[n]],b,c);
    fo(i,0,n-1) printf("%lld ",c[i]);
}

多项式取ln（洛谷P4725）

#include <bits/stdc++.h>
#define fo(i,a,b) for(int i=a;i<=b;++i)
#define fod(i,a,b) for(int i=a;i>=b;--i)
#define mo 998244353
#define LL long long
#define M 262144
#define N 100005
#define L 18
using namespace std;
LL wi[M+1],wg[M+1],a[M+1],b[M+1],ny[M+1];
int bit[M+1],n,cf[L+1],l2[M+1];
LL ksm(LL k,LL n)
{
    LL s=1;
    for(;n;n>>=1,k=k*k%mo) if(n&1) s=s*k%mo;
    return s;
}
void prp(int num)
{
    fo(i,0,num)
    {
        bit[i]=(bit[i>>1]>>1)|((i&1)<<(l2[num]-1));
        wi[i]=wg[i*(M/num)];
    }
    ny[num]=ksm(num,mo-2);
}
void NTT(LL *a,bool pd,int num)
{
    fo(i,0,num-1) if(bit[i]<i) swap(a[i],a[bit[i]]);
    LL v;
    for(int h=1,m=2,lim=num>>1;m<=num;h=m,m<<=1,lim>>=1)
    {
        for(int j=0;j<num;j+=m)
        {
            fo(i,0,h-1)
            {
                v=((!pd)?wi[i*lim]:wi[num-i*lim])*a[i+j+h]%mo;
                a[i+j+h]=(a[i+j]-v+mo)%mo;
                a[i+j]=(a[i+j]+v)%mo;
            }
        }
    }
    if(pd) fo(i,0,num-1) a[i]=a[i]*ny[num]%mo;
}
void inv(int n,LL *a,LL *b)
{
    static LL u1[M+1],u2[M+1];
    fo(i,0,n-1) b[i]=0;
    b[0]=ksm(a[0],mo-2);
    for(int m=1,t=2,num=4;m<n;m=t,t=num,num<<=1)
    {
        prp(num);
        fo(i,0,num-1) u1[i]=u2[i]=0;
        fo(i,0,m-1) u1[i]=b[i];
        fo(i,0,t-1) u2[i]=a[i];
        NTT(u1,0,num),NTT(u2,0,num);
        fo(i,0,num-1) u1[i]=u1[i]*u1[i]%mo*u2[i]%mo;
        NTT(u1,1,num);
        fo(i,0,t-1) b[i]=((LL)2*b[i]-u1[i]+mo)%mo;
    }
}
void getln(int n,LL *a,LL *b)
{
    static LL u3[M+1],u4[M+1];
    int num=cf[l2[2*n+1]];
    inv(n+1,a,u3);
    fo(i,n+1,num) u3[i]=0;
    prp(num);
    fo(i,0,num-1) u4[i]=0;
    fo(i,0,n-1) u4[i]=a[i+1]*(LL)(i+1)%mo;
    NTT(u4,0,num),NTT(u3,0,num);
    fo(i,0,num-1) u4[i]=u3[i]*u4[i]%mo;
    NTT(u4,1,num);
    b[0]=0;
    fo(i,1,n-1) b[i]=u4[i-1]*ny[i]%mo;
}
int main()
{
    cin>>n;
    cf[0]=1;
    fo(i,1,L) l2[cf[i]=cf[i-1]<<1]=i;
    fod(i,M-1,2) if(!l2[i]) l2[i]=l2[i+1];
    ny[1]=1;
    fo(i,2,M) ny[i]=(-ny[mo%i]*(LL)(mo/i)%mo+mo)%mo; 
    wg[0]=1,wg[1]=ksm(3,(mo-1)/M);
    fo(i,2,M) wg[i]=wg[i-1]*wg[1]%mo;
    fo(i,0,n-1) scanf("%lld",&a[i]);
    getln(n,a,b);
    fo(i,0,n-1) printf("%lld ",b[i]);
}

转载于:https://www.cnblogs.com/BAJimH/p/10575044.html

weixin_33912638
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
多项式全家桶
Recalling_Clouds的博客
01-23 720
下面的所有知识均属多项式范畴。 前言 本文略去多项式的定义，毕竟初中数学都讲过。 本文会出现一些可能没有见过的表达： 对于一个关于 xxx 的多项式 f(x)f(x)f(x)，有时也称其最高次项次数为 f(x)f(x)f(x) 的度数。 多项式同样定义了四则运算，包括常数与多项式的四则运算和多项式之间的四则运算，同时也衍生出了同余符号 ≡\equiv≡、对数函数 ln⁡(f(x))\ln(f(x))ln(f(x))，指数函数 exp⁡(f(x))\exp(f(x))exp(f(x)) 等等。 下文一般采
NTT笔记和多项式全家桶
qq_54026781的博客
11-21 828
NTT和多项式全家桶模板
[学习笔记] 多项式全家桶 学习笔记
STJqwq的博客
08-20 2596
写在前面：建议初一就学懂这个东西，这是多项式题目的基础。 0.复数类的实现 定义一个复数为 a+bia+bia+bi 。 加法： a+bi+c+di=(a+c)+(b+d)ia+bi+c+di=(a+c)+(b+d)ia+bi+c+di=(a+c)+(b+d)i 减法：同上，系数改成 −1-1−1 乘法：(a+bi)(c+di)=ac+adi+bci+bdi2=(ac−bd)+(ad+bc)i(a+bi)(c+di)=ac+adi+bci+bdi^2=(ac-bd)+(ad+bc)i(a+bi)(c+di)
【学习笔记】多项式全家桶
cqbzlydd的博客
01-18 435
终于明白为什么叫多项式全家桶了
多项式全家桶 基础知识+各种操作+例题（细节证明）
DD(XYX)的硫氰化氘博客
02-19 6191
只讲了一些基本操作，过于深奥的多项式三角函数、快速插值可以出门看 OneInDark 的博客
c++模板实现多项式相加
03-03
本话题将详细讲解如何利用C++模板实现链表表示的多项式相加。 首先，我们需要了解多项式的基本概念。一个多项式是由常数、变量及其指数组合而成的数学表达式，例如2x^3 + 5x^2 - 3x + 7。在C++中，我们可以使用链表...
多项式拟合平滑.docx
最新发布
06-27
多项式拟合平滑是一种数据分析技术，用于揭示数据集中的潜在趋势和模式。该方法通过构建一个多项式函数，使这个函数尽可能地贴近原始数据点，从而实现对数据的平滑表示。这种方法在处理含有噪声或者不规则变化的数据...
多项式拟合c++程序
01-18
对数据进行多项式拟合，可得拟合多项式的系数，法方程的系数矩阵，平方误差。
多项式拟合_多项式拟合_VBa_
10-01
在Excel VBA（Visual Basic for Applications）环境中，我们可以利用VBA强大的编程能力来实现复杂的数学运算，例如多项式拟合。多项式拟合是一种通过拟合数据点来创建一个多项式函数的过程，常用于数据分析、预测和...
多项式.zip_多项式
07-15
用链表实现多项式的加减乘除以及多项式的打印显示
进军多项式（一）：多项式全家桶模板
jwg2732的专栏
02-09 4067
多项式乘法(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}
多项式知识小全+全家桶学习小结
a_forever_dream的博客
06-03 938
废话 大概用了两三个星期来搞多项式和生成函数，人都要被虐傻了qwq…… 下面是学到的所有知识的小结，应该是很适合初学者食用的教程，不过一定要按顺序看。 多项式全家桶 快速傅利叶变换FFT 快速数论变换NTT 多项式求逆 学接下来的东西需要的前置知识：微积分的一点初步知识，建议在这里学，看前 444 集大概就够了？泰勒展开与牛顿迭代，以及关于 eee 的一些知识。 多项式ln 多项式exp 多项式开根 多项式除法 多项式多点求值 任意模数NTT 在学的过程中同时也在膜拜师兄的优秀板子，所以每一篇的码风不尽相同
【学习笔记】多项式全家桶（包含全套证明）
热门推荐
繁凡さん的博客
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.多项式反...
C语言课程设计一元多项式.doc
12-03
首先提示用户输入两个多项式A和B，输入的形式是先输入多项式A的项数，然后回车，接着输入多项式A的第一项的系数，每一项按照先输入系数，然后输入空格或回车，再输入多项式A的第一项的指数，然后继续输入多项式A的...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值