【CodeVS 3123】 高精度练习之超大整数乘法
来源:互联网 发布:补偿网络 编辑:程序博客网 时间:2024/05/19 20:44
RE了2发 init竟然会开小.......
1 #include <cstdio>
2 #include <algorithm>
3 #include <complex>
4 #include <iostream>
5 #include <cmath>
6 using namespace std;
7 const double PI= acos(-1);
8 int init[100000*4+10];
9 #define C complex<double>
10 C a[100000*4+1],b[100000*4+1];
11 int n,m,rev[100000*4];
12 void FFT(C a[],int type)
13 {
14 for(int i=0;i<n;i++) if(rev[i]<i) swap(a[rev[i]],a[i]);
15 for(int i=2;i<=n;i<<=1)
16 {
17 C wn(cos(2*PI/i),type*sin(2*PI/i));
18 for(int j=0;j<n;j+=i)
19 {
20 C w(1,0),x,y;
21 for(int k=0;k<i/2;k++)
22 {
23 x=a[j+k];
24 y=a[j+k+i/2]*w;
25 a[j+k]=x+y;
26 a[j+k+i/2]=x-y;
27 w=w*wn;
28 }
29 }
30 }
31 }
32 int main()
33 {
34 //freop0en("a.in","r",stdin);
35 char c=getchar();
36 while(c!=' ') init[++n]=c-'0',c=getchar();
37 for(int i=n;i>=1;i--) a[n-i]=init[i];
38 c=getchar();
39 while(c!='\n') init[++m]=c-'0',c=getchar();
40 for(int i=m;i>=1;i--) b[m-i]=init[i];
41
42 n--,m--;
43 // for(int i=0;i<=n;i++) cout<<a[i]<<' ';cout<<endl;
44 // for(int i=0;i<=m;i++) cout<<b[i]<<' ';cout<<endl;
45 m=n+m;int L=0;
46 for(n=1;n<=m;n<<=1)L++; L--;
47 for(int i=0;i<n;i++) rev[i]=((rev[i>>1]>>1)|((i&1)<<L));
48
49 FFT(a,1);FFT(b,1);
50 for(int i=0;i<n;i++) a[i]*=b[i];
51 FFT(a,-1);
52 for(int i=0;i<=m;i++) init[i]=(int)(a[i].real()/n+0.5);//,cout<<init[i]<<' ';cout<<endl;
53 init[m+1]=0;
54 for(int i=0;i<=m;i++) init[i+1]+=init[i]/10,init[i]%=10;//,cout<<init[i+1]<<' ';cout<<endl;
55 if(init[m+1]!=0) m++;
56 for(int i=m;i>=0;i--) printf("%d",init[i]);
57 return 0;
58 }
2 #include <algorithm>
3 #include <complex>
4 #include <iostream>
5 #include <cmath>
6 using namespace std;
7 const double PI= acos(-1);
8 int init[100000*4+10];
9 #define C complex<double>
10 C a[100000*4+1],b[100000*4+1];
11 int n,m,rev[100000*4];
12 void FFT(C a[],int type)
13 {
14 for(int i=0;i<n;i++) if(rev[i]<i) swap(a[rev[i]],a[i]);
15 for(int i=2;i<=n;i<<=1)
16 {
17 C wn(cos(2*PI/i),type*sin(2*PI/i));
18 for(int j=0;j<n;j+=i)
19 {
20 C w(1,0),x,y;
21 for(int k=0;k<i/2;k++)
22 {
23 x=a[j+k];
24 y=a[j+k+i/2]*w;
25 a[j+k]=x+y;
26 a[j+k+i/2]=x-y;
27 w=w*wn;
28 }
29 }
30 }
31 }
32 int main()
33 {
34 //freop0en("a.in","r",stdin);
35 char c=getchar();
36 while(c!=' ') init[++n]=c-'0',c=getchar();
37 for(int i=n;i>=1;i--) a[n-i]=init[i];
38 c=getchar();
39 while(c!='\n') init[++m]=c-'0',c=getchar();
40 for(int i=m;i>=1;i--) b[m-i]=init[i];
41
42 n--,m--;
43 // for(int i=0;i<=n;i++) cout<<a[i]<<' ';cout<<endl;
44 // for(int i=0;i<=m;i++) cout<<b[i]<<' ';cout<<endl;
45 m=n+m;int L=0;
46 for(n=1;n<=m;n<<=1)L++; L--;
47 for(int i=0;i<n;i++) rev[i]=((rev[i>>1]>>1)|((i&1)<<L));
48
49 FFT(a,1);FFT(b,1);
50 for(int i=0;i<n;i++) a[i]*=b[i];
51 FFT(a,-1);
52 for(int i=0;i<=m;i++) init[i]=(int)(a[i].real()/n+0.5);//,cout<<init[i]<<' ';cout<<endl;
53 init[m+1]=0;
54 for(int i=0;i<=m;i++) init[i+1]+=init[i]/10,init[i]%=10;//,cout<<init[i+1]<<' ';cout<<endl;
55 if(init[m+1]!=0) m++;
56 for(int i=m;i>=0;i--) printf("%d",init[i]);
57 return 0;
58 }
0 0
- 【CodeVS 3123】 高精度练习之超大整数乘法
- CODEVS 3123 高精度练习之超大整数乘法
- wikioi3123 高精度练习之超大整数乘法
- cv3123 高精度练习之超大整数乘法(FFT)
- 高精度-codevs-3117高精度练习之乘法
- 高精度练习之超大整数开根
- 【快速傅里叶变换】【FFT】【WikiOI】【P3132】【高精度练习之超大整数乘法】
- Codevs 3117 高精度练习之乘法
- codevs 3117 高精度练习之乘法
- 【codevs 3115 3116 3117】高精度练习之加/减/乘法
- codevs 3119 高精度练习之大整数开根
- 超大整数乘法模板(高精度乘以低精度)
- 高精度-codevs-3115高精度练习之减法
- 高精度-codevs-3116高精度练习之加法
- 3117 高精度练习之乘法
- 3117 高精度练习之乘法
- p3117 高精度练习之乘法
- 3117 高精度练习之乘法
- 【BZOJ 4016】[FJOI2014]最短路径树问题
- 【BZOJ 1095】 [ZJOI2007]Hide 捉迷藏
- 【BZOJ 3144】 [Hnoi2013]切糕
- 【UOJ 34】多项式乘法
- java数据结构 深搜+栈
- 【CodeVS 3123】 高精度练习之超大整数乘法
- 【BZOJ 1415】 [Noi2005]聪聪和可可
- 【BZOJ 2301】 单选错位
- Java String的构造和getByte
- 【BZOJ 1877】 [SDOI2009]晨跑
- KMP模板
- 【BZOJ 1355】 [Baltic2009]Radio Transmission
- 博客开更!!
- 1036: [ZJOI2008]树的统计Count