[已解决]如图问题

奋斗中的鱼

最佳答案

月排行榜 / 总排行榜

liuhongrun2022

2023-9-24 12:24:31

高精度乘法模板加循环也行（和第一个比简单很多）

1. #include<iostream>  
   
2. #include<cstring>  
   
3. #include<algorithm>  
   
4. using namespace std;  
   
5. 
   
6. const int L=110;
   
7. 
   
8. string mul(string a, string b)  
   
9. {  
   
10.     string s;  
    
11.     int na[L], nb[L], nc[L], La = a.size(), Lb = b.size();
    
12.     fill(na,na+L,0); fill(nb,nb+L,0); fill(nc,nc+L,0);
    
13.     for(int i = La - 1; i >= 0; i--) na[La - i] = a[i] - '0';
    
14.     for(int i = Lb - 1; i >= 0; i--) nb[Lb - i] = b[i] - '0';  
    
15.     for(int i = 1; i <= La; i++)  
    
16.         for(int j = 1; j <= Lb; j++)  
    
17.             nc[i + j - 1] += na[i] * nb[j];
    
18.     for(int i = 1; i <= La + Lb; i++)  
    
19.         nc[i + 1] += nc[i] / 10, nc[i] %= 10;
    
20.     if(nc[La + Lb]) s += nc[La + Lb] + '0';
    
21.     for(int i = La + Lb - 1; i >= 1; i--)  
    
22.         s += nc[i] + '0';
    
23.     return s;  
    
24. }  
    
25. 
    
26. int main()  
    
27. {  
    
28.     string a;
    
29.     int exponent;
    
30.     while(cin >> a >> exponent)
    
31.     {
    
32.         string result = "1";
    
33.         for (int i = 0; i < exponent; i++)
    
34.         {
    
35.             result = mul(result, a);
    
36.         }
    
37.         cout << result << endl;  
    
38.     }
    
39.     return 0;  
    
40. }
    

复制代码

跳转到最佳答案楼层

奋斗中的鱼

成功屏蔽GPT

liuhongrun2022

一定要c/c++吗

奋斗中的鱼

liuhongrun2022 发表于 2023-9-24 12:09
一定要c/c++吗

嗯……C++不行吗？

奋斗中的鱼

liuhongrun2022 发表于 2023-9-24 12:09
一定要c/c++吗

想看看能不能实现

liuhongrun2022

一个高精度幂的模板，输入两个数字，然后按回车。

1. #include <iostream>
   
2. #include <cstdio>
   
3. #include <algorithm>
   
4. #include <cstring>
   
5. #include <cmath>
   
6. #include <map>
   
7. #include <queue>
   
8. #include <set>
   
9. #include <vector>
   
10. using namespace std;
    
11. #define L(x) (1 << (x))
    
12. const double PI = acos(-1.0);
    
13. const int Maxn = 133015;
    
14. double ax[Maxn], ay[Maxn], bx[Maxn], by[Maxn];
    
15. char sa[Maxn/2],sb[Maxn/2];
    
16. int sum[Maxn];
    
17. int x1[Maxn],x2[Maxn];
    
18. int revv(int x, int bits)
    
19. {
    
20.     int ret = 0;
    
21.     for (int i = 0; i < bits; i++)
    
22.     {
    
23.         ret <<= 1;
    
24.         ret |= x & 1;
    
25.         x >>= 1;
    
26.     }
    
27.     return ret;
    
28. }
    
29. void fft(double * a, double * b, int n, bool rev)
    
30. {
    
31.     int bits = 0;
    
32.     while (1 << bits < n) ++bits;
    
33.     for (int i = 0; i < n; i++)
    
34.     {
    
35.         int j = revv(i, bits);
    
36.         if (i < j)
    
37.             swap(a[i], a[j]), swap(b[i], b[j]);
    
38.     }
    
39.     for (int len = 2; len <= n; len <<= 1)
    
40.     {
    
41.         int half = len >> 1;
    
42.         double wmx = cos(2 * PI / len), wmy = sin(2 * PI / len);
    
43.         if (rev) wmy = -wmy;
    
44.         for (int i = 0; i < n; i += len)
    
45.         {
    
46.             double wx = 1, wy = 0;
    
47.             for (int j = 0; j < half; j++)
    
48.             {
    
49.                 double cx = a[i + j], cy = b[i + j];
    
50.                 double dx = a[i + j + half], dy = b[i + j + half];
    
51.                 double ex = dx * wx - dy * wy, ey = dx * wy + dy * wx;
    
52.                 a[i + j] = cx + ex, b[i + j] = cy + ey;
    
53.                 a[i + j + half] = cx - ex, b[i + j + half] = cy - ey;
    
54.                 double wnx = wx * wmx - wy * wmy, wny = wx * wmy + wy * wmx;
    
55.                 wx = wnx, wy = wny;
    
56.             }
    
57.         }
    
58.     }
    
59.     if (rev)
    
60.     {
    
61.         for (int i = 0; i < n; i++)
    
62.             a[i] /= n, b[i] /= n;
    
63.     }
    
64. }
    
65. int solve(int a[],int na,int b[],int nb,int ans[])
    
66. {
    
67.     int len = max(na, nb), ln;
    
68.     for(ln=0; L(ln)<len; ++ln);
    
69.     len=L(++ln);
    
70.     for (int i = 0; i < len ; ++i)
    
71.     {
    
72.         if (i >= na) ax[i] = 0, ay[i] =0;
    
73.         else ax[i] = a[i], ay[i] = 0;
    
74.     }
    
75.     fft(ax, ay, len, 0);
    
76.     for (int i = 0; i < len; ++i)
    
77.     {
    
78.         if (i >= nb) bx[i] = 0, by[i] = 0;
    
79.         else bx[i] = b[i], by[i] = 0;
    
80.     }
    
81.     fft(bx, by, len, 0);
    
82.     for (int i = 0; i < len; ++i)
    
83.     {
    
84.         double cx = ax[i] * bx[i] - ay[i] * by[i];
    
85.         double cy = ax[i] * by[i] + ay[i] * bx[i];
    
86.         ax[i] = cx, ay[i] = cy;
    
87.     }
    
88.     fft(ax, ay, len, 1);
    
89.     for (int i = 0; i < len; ++i)
    
90.         ans[i] = (int)(ax[i] + 0.5);
    
91.     return len;
    
92. }
    
93. string mul(string sa,string sb)
    
94. {
    
95.     int l1,l2,l;
    
96.     int i;
    
97.     string ans;
    
98.     memset(sum, 0, sizeof(sum));
    
99.     l1 = sa.size();
    
100.     l2 = sb.size();
     
101.     for(i = 0; i < l1; i++)
     
102.         x1[i] = sa[l1 - i - 1]-'0';
     
103.     for(i = 0; i < l2; i++)
     
104.         x2[i] = sb[l2-i-1]-'0';
     
105.     l = solve(x1, l1, x2, l2, sum);
     
106.     for(i = 0; i<l || sum[i] >= 10; i++)
     
107.     {
     
108.         sum[i + 1] += sum[i] / 10;
     
109.         sum[i] %= 10;
     
110.     }
     
111.     l = i;
     
112.     while(sum[l] <= 0 && l>0)
     
113.     for(i = l; i >= 0; i--)    ans+=sum[i] + '0';
     
114.     return ans;
     
115. }
     
116. string Pow(string a,int n)
     
117. {
     
118.     if(n==1) return a;
     
119.     if(n&1) return mul(Pow(a,n-1),a);
     
120.     string ans=Pow(a,n/2);
     
121.     return mul(ans,ans);
     
122. }
     
123. int main()
     
124. {
     
125.     cin.sync_with_stdio(false);
     
126.     string a;
     
127.     int b;
     
128.     while(cin>>a>>b) cout<<Pow(a,b)<<endl;
     
129.     return 0;
     
130. }
     

复制代码

liuhongrun2022

高精度乘法模板加循环也行（和第一个比简单很多）

1. #include<iostream>  
   
2. #include<cstring>  
   
3. #include<algorithm>  
   
4. using namespace std;  
   
5. 
   
6. const int L=110;
   
7. 
   
8. string mul(string a, string b)  
   
9. {  
   
10.     string s;  
    
11.     int na[L], nb[L], nc[L], La = a.size(), Lb = b.size();
    
12.     fill(na,na+L,0); fill(nb,nb+L,0); fill(nc,nc+L,0);
    
13.     for(int i = La - 1; i >= 0; i--) na[La - i] = a[i] - '0';
    
14.     for(int i = Lb - 1; i >= 0; i--) nb[Lb - i] = b[i] - '0';  
    
15.     for(int i = 1; i <= La; i++)  
    
16.         for(int j = 1; j <= Lb; j++)  
    
17.             nc[i + j - 1] += na[i] * nb[j];
    
18.     for(int i = 1; i <= La + Lb; i++)  
    
19.         nc[i + 1] += nc[i] / 10, nc[i] %= 10;
    
20.     if(nc[La + Lb]) s += nc[La + Lb] + '0';
    
21.     for(int i = La + Lb - 1; i >= 1; i--)  
    
22.         s += nc[i] + '0';
    
23.     return s;  
    
24. }  
    
25. 
    
26. int main()  
    
27. {  
    
28.     string a;
    
29.     int exponent;
    
30.     while(cin >> a >> exponent)
    
31.     {
    
32.         string result = "1";
    
33.         for (int i = 0; i < exponent; i++)
    
34.         {
    
35.             result = mul(result, a);
    
36.         }
    
37.         cout << result << endl;  
    
38.     }
    
39.     return 0;  
    
40. }
    

复制代码

liuhongrun2022

另外，这两个例子都是C++的

liuhongrun2022

我要最佳答案！

人造人

1. sh-5.1$ cat main.cpp
   
2. #include <iostream>
   
3. #include <boost/multiprecision/cpp_int.hpp>
   
4. 
   
5. using std::cin, std::cout, std::endl;
   
6. using boost::multiprecision::cpp_int;
   
7. using boost::multiprecision::pow;
   
8. 
   
9. int main() {
   
10.     cout << pow(cpp_int(2), 100) << endl;
    
11.     return 0;
    
12. }
    
13. sh-5.1$ ./main
    
14. 1267650600228229401496703205376
    
15. sh-5.1$

复制代码

歌者文明清理员

isdkz用的机器人好像可以看图片