题目大意:$FFT$,给你两个多项式,请输出乘起来后的多项式。

题解:$FFT$,由于给的$n$不是很大,也可以用$NTT$做

卡点:

C++ Code:

 FFT:

#include <cstdio>
#include <cmath>
using namespace std;
const double Pi = acos(-1);
int n, m;
struct complex {
double r, i;
complex (double a = 0, double b = 0) {r = a, i = b;}
complex operator + (complex a) {return (complex) {r + a.r, i + a.i};}
complex operator - (complex a) {return (complex) {r - a.r, i - a.i};}
complex operator /= (int a) {r /= a, i /= a;}
complex operator * (complex a) {return (complex) {r * a.r - i * a.i, r * a.i + i * a.r};}
} a[500000], b[500000];
int rev[500000], dig, l;
void swap(complex &a, complex &b) {complex t = a; a = b; b = t;}
void FFT(complex *a, int op) {
for (int i = 0; i < l; i++) if (i < rev[i]) swap(a[i], a[rev[i]]);
for (int mid = 1; mid < l; mid <<= 1 ) {
complex Wn(cos(Pi / mid), op * sin(Pi / mid));
for (int i = 0; i < l; i += (mid << 1)) {
complex W(1, 0);
for (int j = 0; j < mid; j++, W = W * Wn) {
complex X = a[i + j], Y = W * a[i + j + mid];
a[i + j] = X + Y;
a[i + j + mid] = X - Y;
}
}
}
if (op == -1) for (int i = 0; i <= l; i++) a[i] /= l;
}
int main() {
scanf("%d%d", &n, &m);
for (int i = 0; i <= n; i++) scanf("%lf", &a[i].r);
for (int i = 0; i <= m; i++) scanf("%lf", &b[i].r);
l = 1; while (l <= (n + m)) l <<= 1, dig++;
for (int i = 0; i < l; i++) rev[i] = (rev[i >> 1] >> 1) | ((i & 1) << (dig - 1));
FFT(a, 1), FFT(b, 1);
for (int i = 0; i < l; i++) a[i] = a[i] * b[i];
FFT(a, -1);
for (int i = 0; i <= n + m; i++) printf("%d ", int(a[i].r + 0.5));
return 0;
}

NTT:

#include <cstdio>
#define int long long
using namespace std;
const int maxn = 2100010;
const int mod = 998244353;
const int P = 3, invP = 332748118;
int n, m;
int a[maxn], b[maxn], rev[maxn], l, dig;
int Inv[2040826], invl;
inline void swap(int &a, int &b) {a ^= b ^= a ^= b;}
int inv(int i) {
if (i < 2040826) {
if (Inv[i]) return Inv[i];
return (Inv[i] = inv(mod % i) * (mod - mod / i) % mod);
}else return inv(mod % i) * (mod - mod / i) % mod;
}
inline int pw(int base, int p) {
int ans = 1;
for (p <<= 1; p >>= 1; (base *= base) %= mod) if (p & 1) (ans *= base) %= mod;
return ans;
}
void NTT(int *a, int op) {
int Yx;
if (op == 1) Yx = P; else Yx = invP;
for (int i = 0; i < l; i++) if (i < rev[i]) swap(a[i], a[rev[i]]);
for (int mid = 1; mid < l; mid <<= 1) {
int Wn = pw(Yx, (mod - 1) / (mid << 1));
for (int i = 0; i < l; i += (mid << 1)) {
int W = 1;
for (int j = 0; j < mid; j++, W = W * Wn % mod) {
int X = a[i + j], Y = W * a[i + j + mid] % mod;
a[i + j] = (X + Y) % mod;
a[i + j + mid] = (X - Y + mod) % mod;
}
}
}
if (op == -1) for (int i = 0; i < l; i++) a[i] = (a[i] * invl) % mod;
}
signed main() {
Inv[0] = Inv[1] = 1;
scanf("%lld%lld", &n, &m);
for (int i = 0; i <= n; i++) scanf("%lld", &a[i]);
for (int i = 0; i <= m; i++) scanf("%lld", &b[i]);
l = 1; while (l <= n + m) l <<= 1, dig++; invl = inv(l);
for (int i = 1; i < l; i++) rev[i] = (rev[i >> 1] >> 1) | ((i & 1) << (dig - 1));
NTT(a, 1), NTT(b, 1);
for (int i = 0; i < l; i++) (a[i] *= b[i]) %= mod;
NTT(a, -1);
for (int i = 0; i <= n + m; i++) printf("%lld ", a[i]);
return 0;
}

  

  

  

[洛谷P3803] 【模板】多项式乘法(FFT, NTT)的更多相关文章

  1. 洛谷.3803.[模板]多项式乘法(FFT)

    题目链接:洛谷.LOJ. FFT相关:快速傅里叶变换(FFT)详解.FFT总结.从多项式乘法到快速傅里叶变换. 5.4 又看了一遍,这个也不错. 2019.3.7 叕看了一遍,推荐这个. #inclu ...

  2. 洛谷.3803.[模板]多项式乘法(NTT)

    题目链接:洛谷.LOJ. 为什么和那些差那么多啊.. 在这里记一下原根 Definition 阶 若\(a,p\)互质,且\(p>1\),我们称使\(a^n\equiv 1\ (mod\ p)\ ...

  3. 洛谷.4238.[模板]多项式求逆(NTT)

    题目链接 设多项式\(f(x)\)在模\(x^n\)下的逆元为\(g(x)\) \[f(x)g(x)\equiv 1\ (mod\ x^n)\] \[f(x)g(x)-1\equiv 0\ (mod\ ...

  4. P3803 [模板] 多项式乘法 (FFT)

    Rt 注意len要为2的幂 #include <bits/stdc++.h> using namespace std; const double PI = acos(-1.0); inli ...

  5. 洛谷.4512.[模板]多项式除法(NTT)

    题目链接 多项式除法 & 取模 很神奇,记录一下. 只是主要部分,更详细的和其它内容看这吧. 给定一个\(n\)次多项式\(A(x)\)和\(m\)次多项式\(D(x)\),求\(deg(Q) ...

  6. 洛谷 P4512 [模板] 多项式除法

    题目:https://www.luogu.org/problemnew/show/P4512 看博客:https://www.cnblogs.com/owenyu/p/6724611.html htt ...

  7. 洛谷 P4238 [模板] 多项式求逆

    题目:https://www.luogu.org/problemnew/show/P4238 看博客:https://www.cnblogs.com/xiefengze1/p/9107752.html ...

  8. FFT/NTT总结+洛谷P3803 【模板】多项式乘法(FFT)(FFT/NTT)

    前言 众所周知,这两个东西都是用来算多项式乘法的. 对于这种常人思维难以理解的东西,就少些理解,多背板子吧! 因此只总结一下思路和代码,什么概念和推式子就靠巨佬们吧 推荐自为风月马前卒巨佬的概念和定理 ...

  9. 洛谷p3803 FFT入门

    洛谷p3803 FFT入门 ps:花了我一天的时间弄懂fft的原理,感觉fft的折半很神奇! 大致谈一谈FFT的基本原理: 对于两个多项式的卷积,可以O(n^2)求出来(妥妥的暴力) 显然一个多项式可 ...

  10. 【Uoj34】多项式乘法(NTT,FFT)

    [Uoj34]多项式乘法(NTT,FFT) 题面 uoj 题解 首先多项式乘法用\(FFT\)是一个很久很久以前就写过的东西 直接贴一下代码吧.. #include<iostream> # ...

随机推荐

  1. Windows使用Node.js自动生成Vue.js模版环境部署步骤-----记录

    node.js官网下载并安装node 进入node文档目录下,运行cmd 输入 node -v 查看node版本 出现表示安装完成 输入 npm -v 显示npm版本信息 安装cnpm 输入 npm ...

  2. POJ1985 树的直径(BFS

    Cow Marathon   Description After hearing about the epidemic of obesity in the USA, Farmer John wants ...

  3. JavaSE基础复习---2---2018/9/28

    目录: 1.数据类型 2.变量 3.数组 1.数据类型 谈到java的数据类型,必须知道java是强类型语言.首先,每个变量有类型,每个表达式有类型,而且每种类型是严格定义的.其次,所有的数值传递,不 ...

  4. P1049 装箱问题

    装箱问题 题目描述 有一个箱子容量为V(正整数,0<=V<=20000),同时有n个物品(0<n<=30,每个物品有一个体积(正整数). 要求n个物品中,任取若干个装入箱内,使 ...

  5. c++ map的使用方法

    1.头文件:#include<map> 2.定义:map<typename1,typename2> mp 注:字符串数组只能用string而不能使用char[] 3.访问方式: ...

  6. 笔记-flask-原理及请求处理流程

    笔记-flask-原理及请求处理流程 1.      服务器声明及运行 最基本的flask项目代码如下 from flask import Flask app = Flask(__name__) @a ...

  7. JAVA大作业汇总1

    JAVA大作业 代码 ``` package thegreatwork; import javafx.application.; import javafx.scene.control.; impor ...

  8. 基于Ubuntu Server 16.04 LTS版本安装和部署Django之(三):设置上传文件夹权限(这里测试用完全共享)

    基于Ubuntu Server 16.04 LTS版本安装和部署Django之(一):安装Python3-pip和Django 基于Ubuntu Server 16.04 LTS版本安装和部署Djan ...

  9. uber司机已经激活了,就是还没有上传头

    滴快车单单2.5倍,注册地址:http://www.udache.com/ 如何注册Uber司机(全国版最新最详细注册流程)/月入2万/不用抢单:http://www.cnblogs.com/mfry ...

  10. 2648: SJY摆棋子

    2648: SJY摆棋子 https://www.lydsy.com/JudgeOnline/problem.php?id=2648 分析: k-d tree 模板题. 代码: #include< ...