• [模板][P4238]多项式求逆


    NTT多项式求逆模板,详见代码

    #include <map>
    #include <set>
    #include <stack>
    #include <cmath>
    #include <queue>
    #include <cstdio>
    #include <cstring>
    #include <cstdlib>
    #include <iostream>
    #include <algorithm>
    #define ls p<<1 
    #define rs p<<1|1
    using namespace std;
    typedef long long ll;
    const int mxn=5e5+5,mod=998244353,G=3,Gi=332748118;
    int n,l,lim=1,a[mxn],b[mxn],c[mxn],r[mxn];
    inline int read() {
        char c=getchar(); int x=0,f=1;
        while(c>'9'||c<'0') {if(c=='-') f=-1;c=getchar();}
        while(c<='9'&&c>='0') {x=(x<<3)+(x<<1)+(c&15);c=getchar();}
        return x*f;
    }
    inline int chkmax(int &x,int y) {if(x<y) x=y;}
    inline int chkmin(int &x,int y) {if(x>y) x=y;}
    
    struct ed {
        int to,nxt;
    }t[mxn<<1];
    
    int qpow(int a,int b) 
    {
        int res=1,base=a;
        while(b) {
            if(b&1) res=1ll*res*base%mod;
            base=1ll*base*base%mod;
            b>>=1;
        }
        return res;
    }
    
    void NTT(int *p,int opt)
    {
        for(int i=0;i<lim;++i)
            if(i<r[i]) swap(p[i],p[r[i]]);
        for(int mid=1;mid<lim;mid<<=1) {
            int wn=qpow(opt==1?G:Gi,(mod-1)/(mid<<1));
            for(int len=mid<<1,j=0;j<lim;j+=len) {
                int w=1;
                for(int k=0;k<mid;++k,w=1ll*w*wn%mod) {
                    int x=p[j+k],y=1ll*w*p[j+mid+k]%mod;
                    p[j+k]=(x+y)%mod; p[j+mid+k]=(x-y+mod)%mod;
                }
            }
        }	
    }
    
    void solve(int deg,int *a,int *b)
    {
        if(deg==1) {b[0]=qpow(a[0],mod-2); return ;}
        solve((deg+1)>>1,a,b); lim=1,l=0;
        while(lim<=2*deg) lim<<=1,++l;
        for(int i=0;i<lim;++i) 
            r[i]=(r[i>>1]>>1)|((i&1)<<(l-1));
        for(int i=0;i<deg;++i) c[i]=a[i];	
        for(int i=deg;i<lim;++i) c[i]=0;
        NTT(b,1); NTT(c,1); 
        for(int i=0;i<lim;++i)
            b[i]=1ll*(2ll-1ll*c[i]*b[i]%mod+mod)%mod*b[i]%mod;
        NTT(b,-1); int inv=qpow(lim,mod-2);
        for(int i=0;i<lim;++i) b[i]=1ll*b[i]*inv%mod;
        for(int i=deg;i<lim;++i) b[i]=0;
    }
    
    int main()
    {
        n=read();
        for(int i=0;i<n;++i) a[i]=read();
        solve(n,a,b);
        for(int i=0;i<n;++i) printf("%d ",b[i]);
        return 0;
    }
    
    
  • 相关阅读:
    Oracle数据库编程
    使用JDBC处理MySQL大文本和大数据
    phpcms之带图片的登录信息(带cookie版)(由于cookie和PHPCMS的原因,这个暂时无法使用,看新的)
    phpcms之调用导航栏
    phpcms之修改默认显示文字
    phpcms之创建自己的路径
    dedecms 动态tab写法
    关于dedecms的操作
    12/23
    12/21
  • 原文地址:https://www.cnblogs.com/list1/p/10504738.html
Copyright © 2020-2023  润新知