Luogu1919 【模板】A*B Problem升级版（FFT快速傅里叶）

原题链接：https://www.luogu.org/problemnew/show/P1919

【模板】A*B Problem升级版（FFT快速傅里叶）

题目描述

给出两个n位10进制整数x和y，你需要计算x*y。

输入输出格式

输入格式：

第一行一个正整数n。第二行描述一个位数为n的正整数x。第三行描述一个位数为n的正整数y。

输出格式：

输出一行，即x*y的结果。（注意判断前导0）

输入输出样例

输入样例#1：

1
3
4

输出样例#1：

说明

数据范围：

n<=60000

来源：bzoj2179

本题数据为洛谷自造数据，使用CYaRon耗时5分钟完成数据制作。

题解

我们把每一位看做多项式的系数，用FFTFFT\mathcal{FFT}做一波多项式乘法即可。

代码

#include<bits/stdc++.h>
#define db double
using namespace std;
const int M=4e6+5;
const db pi=acos(-1.0);
struct cpx{db x,y;}f[M],g[M];
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};}
int n,mx,len,rev[M],ans[M],flag;
char ch[M];
void in()
{scanf("%d",&n);--n;scanf("%s",ch);for(int i=0;i<=n;++i)f[i].x=ch[i]-'0';scanf("%s",ch);for(int i=0;i<=n;++i)g[i].x=ch[i]-'0';
}
void fft(cpx *f,int typ)
{cpx wn,w,x,y;int i,mid,j,k;for(i=0;i<mx;++i)if(i<rev[i])swap(f[i],f[rev[i]]);for(mid=1;mid<mx;mid<<=1)for(j=0,wn=(cpx){cos(pi/mid),typ*sin(pi/mid)};j<mx;j+=mid<<1)for(k=0,w=(cpx){1,0};k<mid;++k,w=w*wn)x=f[j+k],y=w*f[j+mid+k],f[j+k]=x+y,f[j+mid+k]=x-y;
}
void ac()
{for(mx=1;mx<=n+n;mx<<=1,++len);for(int i=0;i<mx;++i)rev[i]=(rev[i>>1]>>1|((i&1)<<(len-1)));fft(f,1);fft(g,1);for(int i=0;i<=mx;++i)f[i]=f[i]*g[i];fft(f,-1);for(int i=0;i<=n+n;++i)ans[i]=int(f[i].x/mx+0.5);for(int i=n+n;i>=1;--i)ans[i-1]+=ans[i]/10,ans[i]%=10;for(int i=0;i<=n+n;++i){if(ans[i])flag=1;if(flag)printf("%d",ans[i]);}
}
int main(){in();ac();}

Luogu1919 【模板】A*B Problem升级版（FFT快速傅里叶）相关推荐

洛谷 P1919 【模板】A*B Problem升级版（FFT快速傅里叶）
题目链接:P1919 [模板]A*B Problem升级版(FFT快速傅里叶) 题意给定两个 \(n\) 位 \(10\) 进制整数 \(x\) 和 \(y\),求 \(x*y\). 思路 FFT ...
[P1919 【模板】A*B Problem升级版（FFT快速傅里叶）（高精乘板子，FFT板子）
P1919 [模板]A*B Problem升级版(FFT快速傅里叶) 分析: 为什么高精乘可以用 FFT 优化??? 众所周知 FFT 是解决多项式之间的乘法运算,将乘数的每一位看成多项式的系 ...
【模板】A*B Problem升级版（FFT快速傅里叶）
总算是弄懂了一点点FFTFFTFFT 但是那个蝴蝶操作还只是背的板子呜呜先放个优化高精度乘法的板子在这里把 P3803 [模板]多项式乘法(FFT) #include <cmath> # ...
洛谷 P1919 【模板】A*B Problem升级版【快速傅里叶变换 FFT】
[洛谷 P1919] [模板]A*B Problem升级版题目大意给你两个超大整数 a,ba,ba,b,问 a∗ba*ba∗b 其中 a,b≤101000000a,b\leq10^{1000000 ...
洛谷 - P3803 【模板】多项式乘法（FFT/NTT）
题目链接:点击查看题目大意:给出两个多项式 F( x ) 和 G( x ) 的系数,求其卷积后的系数题目分析:存一个FFT的模板,原理学不明白,数论和dp都扔给队友了,当个快乐的fw 代码: // ...
FFT(快速傅里叶变化)学习
快速傅里叶变化忽然发现自己FFT的博客都已经鸽了一年了,恰好又讲FFT就补了一篇. 关于多项式多项式是什么垃圾玩意就不用我多说了吧在平常解决各类数学问题中,我们都很容易发现多项式乘法的影子. 我 ...
P3803 【模板】多项式乘法（FFT）
P3803 [模板]多项式乘法(FFT) 题目描述给定一个 n 次多项式 F(x),和一个 m 次多项式 G(x). 请求出 F(x)和 G(x)的卷积. 从低到高输出F(x)*G(x)的系数另一 ...
洛谷 P1919 模板】A*B Problem升级版（FFT快速傅里叶）
https://www.luogu.com.cn/problem/P1919 题目背景本题数据已加强,请使用 FFT/NTT,不要再交 Python 代码浪费评测资源. 题目描述给你两个正整数 a ...
luogu P1919 【模板】A*B Problem升级版（FFT快速傅里叶）
题链:https://www.luogu.com.cn/problem/P1919 思路:将数看成即可. #include <bits/stdc++.h> #define ll long ...
洛谷P1919 【模板】A*B Problem升级版（FFT快速傅里叶）
题目描述给出两个n位10进制整数x和y,你需要计算x*y. 输入输出格式输入格式: 第一行一个正整数n. 第二行描述一个位数为n的正整数x. 第三行描述一个位数为n的正整数y. 输出格式: 输出一 ...