免費論壇 繁體 | 簡體
Sclub交友聊天~加入聊天室當版主
分享
新浪微博
QQ空间
Facebook
Google+
Plurk
Twitter
Line
Board logo

標題: c++大數階層計算/large factorial [打印本頁]
作者: Smallcat    時間: 2019-12-5 17:16     標題: c++大數階層計算/large factorial

  1. #include<iostream>
  2. using namespace std;

  4. //目前極限是1142!(2998個位數)

  6. #define MAX 1000

  8. int multiply(int x, int res[], int res2[], int res3[], int res_size);

  10. void factorial(int n)
  11. {
  12.         int res[MAX],res2[MAX],res3[MAX];
  13.         string s;
  14.        
  15.         res[0] = 1;
  16.         int size = 1;

  18.         for (int x=2; x<=n; x++)
  19.                 size = multiply(x, res, res2, res3, size);
  20.         if(size>2000){
  21.                 for (int i=size-1-2000; i>=0; i--)
  22.                 s += res3[i]+'0';
  23.                 for (int i=2000-1-1000; i>=0; i--)
  24.                 s += res2[i]+'0';
  25.                 for (int i=1000-1; i>=0; i--)
  26.                 s += res[i]+'0';
  27.         }
  28.         if(size>1000 && size<=2000){
  29.                 for (int i=size-1-1000; i>=0; i--)
  30.                 cout << res2[i];
  31.                 for (int i=1000-1; i>=0; i--)
  32.                 cout << res[i];
  33.         }
  34.         if(size<=1000){
  35.                 for (int i=size-1; i>=0; i--)
  36.                 s += res[i]+'0';
  37.         }
  38.         cout << "size:" << size << endl;
  39.         cout << s << endl;
  40. }

  42. int multiply(int x, int res[], int res2[], int res3[], int res_size)
  43. {
  44.         int carry = 0;
  45.        
  46.         if(res_size<=1000){
  47.                 for (int i=0; i<res_size; i++)
  48.             {
  49.                 int prod = res[i] * x + carry;

  51.                 res[i] = prod % 10;   

  53.                 carry  = prod/10;     
  54.             }
  55.         }else if(res_size<2000){
  56.                 for (int i=0; i<1000; i++)
  57.             {
  58.                 int prod = res[i] * x + carry;

  60.                 res[i] = prod % 10;   

  62.                 carry  = prod/10;     
  63.                
  64.                    }
  65.                    for (int i=0; i<res_size-1000; i++)
  66.                 {                       
  67.                         int prod = res2[i] * x + carry;

  69.                         res2[i] = prod % 10;
  70.                                
  71.                         carry = prod/10;       
  72.                 }
  73.         }else{
  74.                 for (int i=0; i<1000; i++)
  75.                 {
  76.                         int prod = res[i] * x + carry;

  78.                         res[i] = prod % 10;
  79.        
  80.                         carry = prod/10;         
  81.                 }
  82.                 for (int i=0; i<1000; i++)
  83.                 {
  84.                         int prod = res2[i] * x + carry;

  86.                         res2[i] = prod % 10;
  87.        
  88.                         carry = prod/10;       
  89.                 }
  90.                 for (int i=0; i<res_size-2000; i++)
  91.                 {
  92.                         int prod = res3[i] * x + carry;

  94.                         res3[i] = prod % 10;
  95.        
  96.                         carry = prod/10;       
  97.                 }
  98.         }
  99.                 while (carry)
  100.                 {
  101.                         if(res_size<1000){
  102.                                 res[res_size] = carry%10;
  103.                                 carry = carry/10;
  104.                         }else if(res_size<2000){
  105.                                 res2[res_size-1000] = carry%10;
  106.                                 carry = carry/10;
  107.                         }else{
  108.                                 res3[res_size-2000] = carry%10;
  109.                                 carry = carry/10;
  110.                         }
  111.                         res_size++;
  112.                 }
  113.                 return res_size;
  114.        
  115. }

  117. int main()
  118. {
  119.         int a;
  120.         for(int i=0;i<10;i++){
  121.                 for(int j=0;j<10;j++){
  122.                         cout<<j;
  123.                 }
  124.         }
  125.         cout<<endl;
  126.         while(1){
  127.                 cin>>a;
  128.                 factorial(a);
  129.         }
  130.         return 0;
  131. }
複製代碼




歡迎光臨 小貓貓大聯盟!幻與想的境界\(0w0)/ (http://smallcat.utmall.com/) Powered by Discuz! 7.2