洛谷 P4512 【模板】多项式除法 ntt

本文介绍了一种使用快速傅里叶变换(FFT)进行多项式除法的方法,通过模998244353下的运算实现了多项式的高效相除,并给出了完整的C++代码实现。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题目描述

给定一个 nn 次多项式 F(x) 和一个 mm 次多项式 G(x) ,请求出多项式 Q(x)Q(x)R(x)R(x) ,满足以下条件:

Q(x)Q(x)次数为 nmn−mR(x)R(x)次数小于mm
F(x)=Q(x)G(x)+R(x)
所有的运算在模 998244353998244353 意义下进行。

输入输出格式

输入格式:
第一行两个整数 nnm,意义如上。
第二行 n+1n+1 个整数,从低到高表示 F(x)F(x) 的各个系数。
第三行 m+1m+1 个整数,从低到高表示 G(x)G(x) 的各个系数。

输出格式:
第一行 nm+1n−m+1 个整数,从低到高表示 Q(x)Q(x) 的各个系数。
第二行 mm 个整数,从低到高表示 R(x) 的各个系数。
如果 R(x)R(x) 不足 m1m−1 次,多余的项系数补 00

输入输出样例

输入样例#1:
5 1
1 9 2 6 0 8
1 7
输出样例#1:
237340659 335104102 649004347 448191342 855638018
760903695
说明
对于所有数据,1m<n105,给出的系数均属于[0,998244353)Z[0,998244353)∩Z

分析:
我们设degA=ndegA=ndegB=mdegB=mdegC=nm+1degC=n−m+1degD<mdegD<m
设一个nn次多项式A的变换为xnA(1x)xnA(1x),其实就是把AA翻转。
因为A(x)=B(x)C(x)+D(x),有

A(1x)=B(1x)C(1x)+D(1x)A(1x)=B(1x)∗C(1x)+D(1x)

两边同时乘xnxn
xnA(1x)=xmB(1x)xnmC(1x)+xnm+1xm1D(1x)xnA(1x)=xmB(1x)∗xn−mC(1x)+xn−m+1∗xm−1D(1x)

那么如果我们在mod xnm+1mod xn−m+1意义下跑,那么DD就没掉了,即
xnA(1x)xmB(1x)xnmC(1x) (mod xnm+1)

xnmC(1x)A(x)B(x) (mod xnm+1)xn−mC(1x)≡A′(x)B′(x) (mod xn−m+1)

一个多项式求逆后就可以得到B1(x)B′−1(x),卷上A(x)A′(x)就可以得到C(x)C′(x)
D(x)=A(x)B(x)C(x)D(x)=A(x)−B(x)∗C(x)。就完成了。

代码:

// luogu-judger-enable-o2
#include <iostream>
#include <cmath>
#include <cstdio>
#include <algorithm>
#define LL long long

const LL G=3;
const int maxn=3e5+7;
const LL mod=998244353;

using namespace std;

int n,m,len;
LL f[maxn],g[maxn],inv[maxn],c[maxn],r[maxn],x[maxn],y[maxn],w[maxn];

LL power(LL x,LL y)
{
    if (y==1) return x;
    LL c=power(x,y/2);
    c=(c*c)%mod;
    if (y%2) c=(c*x)%mod;
    return c;
}

void ntt(LL *a,LL f)
{
    for (LL i=0;i<len;i++)
    {
        if (i<r[i]) swap(a[i],a[r[i]]);
    }   
    w[0]=1;
    for (LL i=2;i<=len;i*=2)
    {
        LL wn;
        if (f==1) wn=power(G,(mod-1)/i);
             else wn=power(G,(mod-1)-(mod-1)/i);
        for (LL j=i/2;j>=0;j-=2) w[j]=w[j/2];
        for (LL j=1;j<i/2;j+=2) w[j]=(w[j-1]*wn)%mod;
        for (LL j=0;j<len;j+=i)
        {
            for (LL k=0;k<i/2;k++)
            {
                LL u=a[j+k],v=(a[j+k+i/2]*w[k])%mod;
                a[j+k]=(u+v)%mod;
                a[j+k+i/2]=(u-v+mod)%mod;
            }
        }
    }
    if (f==-1)
    {
        LL inv=power(len,mod-2);
        for (LL i=0;i<len;i++) a[i]=(a[i]*inv)%mod;
    }
}

void init(LL len)
{
    LL k=trunc(log(len+0.5)/log(2));
    for (LL i=0;i<len;i++)
    {
        r[i]=(r[i>>1]>>1)|((i&1)<<(k-1));
    }
}


void NTT(LL *a,LL *b,LL *c,LL n,LL m)
{   
    len=1;
    while (len<=(n+m)) len*=2;
    init(len);
    for (LL i=0;i<len;i++)
    {
        if (i<n) x[i]=a[i]; else x[i]=0;
        if (i<m) y[i]=b[i]; else y[i]=0;
    }   
    ntt(x,1); ntt(y,1);
    for (LL i=0;i<len;i++) c[i]=x[i]*y[i]%mod;
    ntt(c,-1);
}

void getinv(LL *a,LL *b,int deg)
{
    if (deg==1)
    {       
        b[0]=power(a[0],mod-2);
        return;     
    }
    LL d=(deg+1)/2;
    getinv(a,b,d);   
    NTT(a,b,c,m+1,d);
    c[0]=(2+mod-c[0])%mod;
    for (int i=1;i<=m+d;i++) c[i]=(mod-c[i])%mod;
    NTT(c,b,b,m+d+1,d);
    for (int i=deg;i<len;i++) b[i]=0;
}


int main()
{
    scanf("%d%d",&n,&m);
    for (int i=0;i<=n;i++) scanf("%lld",&f[i]);
    for (int i=0;i<=m;i++) scanf("%lld",&g[i]);
    reverse(f,f+n+1);
    reverse(g,g+m+1);   
    getinv(g,inv,n-m+1);
    NTT(f,inv,c,n+1,n-m+1);
    reverse(c,c+n-m+1);
    for (int i=0;i<=n-m;i++) printf("%lld ",c[i]);
    printf("\n");
    reverse(f,f+n+1);
    reverse(g,g+m+1);
    NTT(g,c,c,m+1,n-m+1);
    for (int i=0;i<m;i++) printf("%lld ",(f[i]-c[i]+mod)%mod);
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值