大数阶乘
如何求大数的阶乘,如1000的阶乘,主要是求原理。。
[解决办法]
虽然效率比较低,但是比较适合新手来学习。。。
http://blog.csdn.net/czh3642210/article/details/7350793
[解决办法]
- C/C++ code
#include <iostream>#include <string>#include <algorithm>using namespace std;int COMPARE(string number1, string number2) { int i,j; int length1 = number1.size(); int length2 = number2.size(); if(number1.size() == 0) number1 = "0"; if(number2.size() == 0) number2 = "0"; j = 0; for(i = 0; i < length1; ++i) { if(number1[i] == '0') ++j; else break; } number1 = number1.substr(j); j = 0; for(i = 0; i < length2; ++i) { if(number2[i] == '0') ++j; else break; } number2 = number2.substr(j); length1 = number1.size(); length2 = number2.size(); if(length1 > length2) { return 1; } else if(length1 == length2) { if(number1.compare(number2) > 0) { return 1; } else if(number1.compare(number2) == 0) { return 0; } else { return -1; } } else { return -1; } return 0;}string PLUS(string number1,string number2) { int i; int length1 = number1.size(); int length2 = number2.size(); string result=""; reverse(number1.begin(), number1.end()); reverse(number2.begin(), number2.end()); for(i = 0; i < length1 && i < length2; i++) { char c = (char)(number1[i] + number2[i] - 48); result = result + c; } while(i < length1) { result = result + number1[i]; ++i; } while(i < length2) { result = result + number2[i]; ++i; } int carry = 0; for(i = 0; i < (int)result.size(); ++i) { int value = result[i] - 48 + carry; result[i] = (char)(value % 10 + 48); carry = value / 10; } if(carry !=0 ) { result = result + (char)(carry + 48); } for(i = result.size() - 1; i >= 0; i--) { if(result[i] != '0') break; } result = result.substr(0, i + 1); reverse(result.begin(), result.end()); if(result.length() == 0) result = "0"; return result;}string MINUS(string number1,string number2) { int i; string result = ""; int length1 = number1.size(); int length2 = number2.size(); if(COMPARE(number2,number1) > 0) { return "-" + MINUS(number2, number1); } reverse(number1.begin(),number1.end()); reverse(number2.begin(),number2.end()); for(i = 0; i < length1 && i < length2; i++) { char c = number1[i] - number2[i] + 48; result = result + c; } if(i < length1) { for(; i < length1; i++) { result = result + number1[i]; } } int carry = 0; for(i = 0; i < (int)result.length(); i++) { int value = result[i] - 48 + carry; if(value < 0) { value = value + 10; carry = -1; } else carry = 0; result[i]=(char)(value + 48); } for(i = result.size() - 1; i >= 0; i--) { if(result[i] != '0')break; } result = result.substr(0, i+1); reverse(result.begin(), result.end()); if(result.length()==0) result = "0"; return result;}string MULTIPLY(string number1, string number2) { int i, j; int *iresult; int length1 = number1.size(); int length2 = number2.size(); string result = ""; reverse(number1.begin(), number1.end()); reverse(number2.begin(), number2.end()); iresult = (int*)malloc(sizeof(int) * (length1 + length2 + 1)); memset(iresult, 0, sizeof(int) * (length1 + length2 + 1)); for(i = 0; i < length1; i++) { for(j = 0; j < length2; j++) { iresult[i+j] += ((number1[i] - 48) * (number2[j] - 48)); } } int carry = 0; for(i = 0; i < length1 + length2; i++) { int value = iresult[i] + carry; iresult[i] = value % 10; carry = value / 10; } for(i = length1 + length2 - 1; i >= 0; i--) { if(iresult[i] != 0)break; } for(; i >= 0; i--) { result = result + (char)(iresult[i]+48); } free(iresult); if(result == "") result = "0"; return result;}string factorial(string n) { string temp = "1"; string i; for(i = "1"; COMPARE(i, n) <= 0; i = PLUS(i, "1")) { temp = MULTIPLY(temp, i); } return temp;}int main(void) { cout << factorial("1000") << endl; return 0;}