任意模数多项式乘法MTT(可拆系数FFT)详解

更好的阅读体验

任意模数多项式乘法

前言:
在教练讲的时候脑子并不清醒,所以没听懂。后来自己看博客学会了,但目前只学了一种方法:可拆系数FFT。为了方便日后复习,决定先写下这个的笔记,关于三模数NTT下次再补。

建议:准备好演算纸和笔,本篇含有大量推算部分。

为什么不直接使用NTT/FFT

此处的模数是任意的,所以我们使用NTT时,有局限性。只有当模数满足下列情况时才可使用NTT。

模数为 P P P P = r × 2 k + 1 P=r\times 2^k + 1 P=r×2k+1,其中 k k k足够大,满足 2 k ≥ N 2^k \geq N 2kN,其中 N N N为多项式扩充后的项数,多项式乘法都需要将项数扩充到 2 2 2的幂。

如果不满足上述条件,就不能直接使用NTT。

那为什么不使用FFT带模运算?
首先不考虑模数的情况下,多项式结果的系数不能超过 d o u b l e double double能表示的精度(一般在 1 0 16 10^{16} 1016)。超过后, d o u b l e double double所表示的结果将不再精确。

那怎么办?目前可以给出两种解决方法

  1. 可拆系数FFT
  2. 三模数NTT (可是我还没学懂QAQ)

所以,向可拆系数FFT进发!

可拆系数FFT

怎么用可拆系数FFT

首先我们还是先假设一个情景:

求: c 1 c_1 c1 c 2 c_2 c2的卷积

我们可以把一个数拆成 a × 2 15 + b a\times 2^{15}+b a×215+b的形式(不一定是 2 15 2^{15} 215,大概在 值域 \sqrt{值域} 值域 内)。

c 1 = a 1 × 2 15 + a 2 c 2 = b 1 × 2 15 + b 2 c_1=a_1\times 2^{15}+a_2 \\ c_2=b_1\times 2^{15}+b_2 c1=a1×215+a2c2=b1×215+b2

那这俩的积为

c 1 × c 2 = ( a 1 × 2 15 + a 2 ) × ( b 1 × 2 15 + b 2 ) = a 1 b 1 × 2 30 + ( a 1 b 2 + a 2 b 1 ) × 2 15 + a 2 b 2 \begin{aligned} c_1\times c_2&=(a_1\times 2^{15}+a_2)\times(b_1\times 2^{15}+b_2)\\&=a_1b_1\times 2^{30}+(a_1b_2+a_2b_1)\times 2^{15}+a_2b_2 \end{aligned} c1

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值