http://acm.hdu.edu.cn/showproblem.php?pid=1402

题意:

  求$a*b$ 但是$a$和$b$的范围可以达到 $1e50000$

题解:

  显然...用字符串模拟的大数或者压位的大数是无法胜任这种计算的....

  然后,2个大整数相乘,可以理解为卷积,所以就用快速傅里叶变换(FFT)来加速他

  模板题

  简单总结一下对FFT的认知:

FFT用于算卷积,卷积可以理解为两个多项式相乘
显然复杂度是$O(n^2)$的 但是fft可以优化为$O(nlogn)$
如何优化,
考虑2个点确定直线,3个点确定抛物线,$n+1$个点确定了$n$次多项式
我们用$n+1$个点$(x_i,y_i)$表达了原来的式子
如何用点值表达卷积(多项式相乘)?
式子展开可以知道,$x_i$确定的时候,对于点值表达式就对应了点$(x_i,f(x_i)*g(x_i))$
也就是可以$O(1)$得到一个点答案,$n$个点就是$O(n)$
那问题又来了,$n+1$个点的$(x_i,y_i)$,每次求解一次多项式的值,需要$O(n)$
也就是$O(n^2)$复杂度,依然超时
所以,现在考虑如何快速的算多项式的值
因为$x_i$我们可以自己定义,因此利用虚根的特殊性质,我们就能分治处理了
就可以复杂度$O(nlogn)$的得到点值表达式
然后$O(n)$的算点值表达式的值
然后逆运算回到答案

#include <bits/stdc++.h>
#define ll long long
#define usd unsigned
#define IO ios::sync_with_stdio(false);cin.tie(0);cout.tie(0)
#define rep(ii,a,b) for(int ii=a;ii<=b;++i)
#define per(ii,a,b) for(int ii=b;ii>=a;--i)
#define show(x) cout<<#x<<"="<<x<<endl
#define show2(x,y) cout<<#x<<"="<<x<<" "<<#y<<"="<<y<<endl
#define show3(x,y,z) cout<<#x<<"="<<x<<" "<<#y<<"="<<y<<" "<<#z<<"="<<z<<endl
#define show4(w,x,y,z) cout<<#w<<"="<<w<<" "<<#x<<"="<<x<<" "<<#y<<"="<<y<<" "<<#z<<"="<<z<<endl
#define showa(a,b) cout<<#a<<'['<<b<<"]="<<a[b]<<endl
#define pii pair<int,int>
#define cp complex<double>//<complex>头文件里的东西
using namespace std;
const int maxn=2e5+10;
const int maxm=2e6+10;
const int INF=0x3f3f3f3f;
const int mod=1e9+7;
int casn,n,m,k;
const double pi=acos(-1.0);
int lena,lenb,len,pw,ans[maxn],rr[maxn];
cp a[maxn] ,b[maxn];
char sa[maxn],sb[maxn];
void fft(cp *a,int f){//递归分治太慢,用迭代实现
rep(i,0,n-1) if(i<rr[i]) swap(a[i],a[rr[i]]);//分治原来的数组
for(int i=1;i<n;i<<=1){
cp wn(cos(pi/i),sin(f*pi/i)),x,y;//计算虚根
for(int j=0;j<n;j+=(i<<1)){
cp w(1,0);
for(int k=0;k<i;++k,w*=wn){
x=a[j+k],y=w*a[i+j+k];
a[j+k]=x+y;
a[i+j+k]=x-y;
}
}
}
}
int main(){
while(~scanf("%s%s",sa,sb)){
n=1;
len=pw=0;
memset(ans,0,sizeof ans);
lena=strlen(sa);
lenb=strlen(sb);
while(n<=lena+lenb) n<<=1,++pw;
rep(i,0,lena-1) a[i]=sa[lena-i-1]-'0';
rep(i,lena,n) a[i]=0;
rep(i,0,lenb-1) b[i]=sb[lenb-i-1]-'0';
rep(i,lenb,n) b[i]=0;
rep(i,0,n) rr[i]=(rr[i>>1]>>1)|((i&1)<<(pw-1));//分治中的二进制规律,因此预处理参数
fft(a,1);//DFT a
fft(b,1);//DFT b
rep(i,0,n) a[i]*=b[i];
fft(a,-1);//IDFT a*b
rep(i,0,lena+lenb) {
ans[i]+=int(a[i].real()/n+0.5);
ans[i+1]+=ans[i]/10;
ans[i]%=10;
if(ans[i]>0) len=i;
}
per(i,0,len){
putchar (ans[i]+'0');
}
puts("");
}
return 0;
}

  

HDU - 1402 A * B Problem Plus FFT裸题的更多相关文章

  1. HDU 1402 A * B Problem Plus (FFT模板题)

    FFT模板题,求A*B. 用次FFT模板需要注意的是,N应为2的幂次,不然二进制平摊反转置换会出现死循环. 取出结果值时注意精度,要加上eps才能A. #include <cstdio> ...

  2. hdu 1402 A * B Problem Plus FFT

    /* hdu 1402 A * B Problem Plus FFT 这是我的第二道FFT的题 第一题是完全照着别人的代码敲出来的,也不明白是什么意思 这个代码是在前一题的基础上改的 做完这个题,我才 ...

  3. HDU 1402 A * B Problem Plus (FFT求高精度乘法)

    A * B Problem Plus Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Other ...

  4. HDU - 1402 A * B Problem Plus (FFT实现高精度乘法)

    题意:计算A*B,A,B均为长度小于50000的整数. 这是FFT在大整数相乘中的一个应用,我本来想用NTT做的,但NTT由于取模很可能取炸,所以base必须设得很小,而且效率也比不上FFT. A和B ...

  5. HDU 5832 A water problem(某水题)

    p.MsoNormal { margin: 0pt; margin-bottom: .0001pt; text-align: justify; font-family: Calibri; font-s ...

  6. HDU 1402 A * B Problem Plus ——(大数乘法,FFT)

    因为刚学fft,想拿这题练练手,结果WA了个爽= =. 总结几点犯的错误: 1.要注意处理前导零的问题. 2.一定要注意数组大小的问题.(前一个fft的题因为没用到b数组,所以b就没管,这里使用了b数 ...

  7. hdu 2993 MAX Average Problem(斜率DP入门题)

    题目链接:hdu 2993 MAX Average Problem 题意: 给一个长度为 n 的序列,找出长度 >= k 的平均值最大的连续子序列. 题解: 这题是论文的原题,请参照2004集训 ...

  8. HDU 1402 A * B Problem Plus 快速傅里叶变换 FFT 多项式

    http://acm.hdu.edu.cn/showproblem.php?pid=1402 快速傅里叶变换优化的高精度乘法. https://blog.csdn.net/ggn_2015/artic ...

  9. hdu 1402 A * B Problem Plus (FFT模板)

    A * B Problem Plus Problem Description Calculate A * B. Input Each line will contain two integers A ...

随机推荐

  1. Linux记录- Linux下限制rm -rf /

    操作说明: 为了防止在linux下执行操作的时候误操作rm -rf /,或者rm -rf 一些比较重要的目录,我们做以下操作来限制rm的删除 1.下载源码安装包 https://raw.githubu ...

  2. CSS3笔记3

    1.CSS的层叠性 <!DOCTYPE html> <html lang="en"> <head> <meta charset=" ...

  3. 【leetcode-74】搜索二维矩阵

    (较简单,但犯错太多) 编写一个高效的算法来判断 m x n 矩阵中,是否存在一个目标值.该矩阵具有如下特性: 每行中的整数从左到右按升序排列. 每行的第一个整数大于前一行的最后一个整数. 示例 1: ...

  4. 微信小程序,错误{"errMsg":"request:fail 小程序要求的 TLS 版本必须大于等于 1.2"}

    解决方法一: 开发环境,项目--->勾选不校验即可 解决办法二: 在 PowerShell中运行以下内容, 然后重启服务器 # Enables TLS R2 and Windows # Thes ...

  5. jsp查询修改

    也是只帖代码....不讲解. 1.search.jsp <%@ page language="java" import="java.util.*" pag ...

  6. write命令帮助文档(ubuntu 18.04)

    WRITE() BSD General Commands Manual WRITE() NAME write — send a message to another user SYNOPSIS wri ...

  7. tensorflow can not find libcusolver.so.8.0

    ImportError: libcusolver.so.8.0: cannot open shared object file: No such file or directory solution: ...

  8. HTML字体自动换行第二行缩进一格

    p {text-indent: -2em; margin-left: 2em;}设置text-indent: -2em;以后p标签中第一行文字向左偏移,这样第二行开始的文字就等于缩进了,但是这样设置会 ...

  9. [C++]3-1 得分(Score ACM-ICPC Seoul 2005,UVa1585)

    Question 习题3-1 得分(Score ACM-ICPC Seoul 2005,UVa1585) 题目:给出一个由O和X组成的串(长度为1~80),统计得分. 每个O的分数为目前连续出现的O的 ...

  10. Coursera, Machine Learning, notes

      Basic theory (i) Supervised learning (parametric/non-parametric algorithms, support vector machine ...