P3803 【模板】多项式乘法（FFT）

题目描述

给定一个 n 次多项式 F(x)，和一个 m 次多项式 G(x)。
请求出 F(x)和 G(x)的卷积。
从低到高输出F(x)*G(x)的系数
另一种问法：
如果有两个无限序列a和b，那么它们卷积的结果是：
求出yn值

题解：

模板题
建议背过模板
讲的非常不错的博客

代码：

#include <cmath>
#include <cctype>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;const int maxn=2e6+5;
const double pi=3.1415926535898;
int t, n, m, len=1, l, r[maxn*2];struct Cpx{  //复数double x, y;Cpx (double t1=0, double t2=0){ x=t1, y=t2; }
}A[maxn*2], B[maxn*2], C[maxn*2];
Cpx operator +(Cpx a, Cpx b){ return Cpx(a.x+b.x, a.y+b.y); }
Cpx operator -(Cpx a, Cpx b){ return Cpx(a.x-b.x, a.y-b.y); }
Cpx operator *(Cpx a, Cpx b){ return Cpx(a.x*b.x-a.y*b.y, a.x*b.y+a.y*b.x); }void fdft(Cpx *a, int n, int flag){  //快速将当前多项式从系数表达转换为点值表达for (int i=0; i<n; ++i) if (i<r[i]) swap(a[i], a[r[i]]);for (int mid=1; mid<n; mid<<=1){  //当前区间长度的一半Cpx w1(cos(pi/mid), flag*sin(pi/mid)), x, y;for (int j=0; j<n; j+=(mid<<1)){  //j:区间起始点Cpx w(1, 0);for (int k=0; k<mid; ++k, w=w*w1){  //系数转点值x=a[j+k], y=w*a[j+mid+k];a[j+k]=x+y; a[j+mid+k]=x-y;}}}
}inline int getint(int &x){char c; int flag=0;for (c=getchar(); !isdigit(c); c=getchar())if (c=='-') flag=1;for (x=c-48; c=getchar(), isdigit(c);)x=(x<<3)+(x<<1)+c-48;return flag?x:-x;
}int main(){getint(n); getint(m); int x;for (int i=0; i<=n; ++i) getint(x), A[i].x=x;for (int i=0; i<=m; ++i) getint(x), B[i].x=x;while (len<=n+m) len<<=1, ++l;  //idft需要至少l1+l2个点值for (int i=0; i<len; ++i)  //编号的字节长度为lr[i]=(r[i>>1]>>1)|((i&1)<<(l-1));fdft(A, len, 1); fdft(B, len, 1);for (int i=0; i<len; ++i) C[i]=A[i]*B[i];fdft(C, len, -1);  //idftfor (int i=0; i<=n+m; ++i) printf("%d ", int(C[i].x/len+0.5));return 0;
}

P3803 【模板】多项式乘法（FFT）相关推荐

多项式乘法 FFT模板
题目传送门在下只是来存个板子,,(板子还是洛谷找的2333) 证明的话,太(wo)难(bu)写(hui),就先留个坑吧,,, NOIP后,如果没退役,我会回来填坑的,, #include<bi ...
解题报告（二）多项式问题（多项式乘法及其各种运算）（ACM/ OI）超高质量题解
整理的算法模板合集: ACM模板点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划繁凡出品的全新系列:解题报告系列 -- 超高质量算法题单,配套我写的超高质量的题解和代码,题目难度不一 ...
洛谷 - P3803 【模板】多项式乘法（FFT/NTT）
题目链接:点击查看题目大意:给出两个多项式 F( x ) 和 G( x ) 的系数,求其卷积后的系数题目分析:存一个FFT的模板,原理学不明白,数论和dp都扔给队友了,当个快乐的fw 代码: // ...
【luogu P3803】【模板】多项式乘法（NTT）
[模板]多项式乘法(NTT) 题目链接:luogu P3803 题目大意给你两个多项式,要你求它们的卷积. 思路这次我们写 NTT 的做法. 它的优点就是它可以取模,而且不会有精度问题,而且会比 ...
luogu P4725 多项式对数函数 (模板题、FFT、多项式求逆、求导和积分)
luogu P4725 多项式对数函数 (模板题.FFT.多项式求逆.求导和积分) 手动博客搬家: 本文发表于20181125 13:25:03, 原地址https://blog.csdn.net/s ...
P4245 【模板】任意模数多项式乘法
P4245 [模板]任意模数多项式乘法 https://www.luogu.com.cn/blog/AzusaCat/solution-p4245 首先这类问题指的是对于一个非NTT模数,我们如何计算 ...
P4245 【模板】任意模数多项式乘法(NTT)
题意: P4245 [模板]任意模数多项式乘法题解: NTT模板,记录一下代码: #include <bits/stdc++.h>using namespace std;#define ...
luogu P4512 多项式除法 (模板题、FFT、多项式求逆)
luogu P4512 多项式除法 (模板题.FFT.多项式求逆) 手动博客搬家: 本文发表于20181206 14:42:53, 原地址https://blog.csdn.net/suncongbo ...
luogu P4238 多项式求逆 (模板题、FFT)
luogu P4238 多项式求逆 (模板题.FFT) 手动博客搬家: 本文发表于20181125 13:21:46, 原地址https://blog.csdn.net/suncongbo/artic ...

P3803 【模板】多项式乘法（FFT）

题目描述

题解：

代码：

P3803 【模板】多项式乘法（FFT）相关推荐

最新文章

热门文章