참고 :
제가 직접 짠 큰 수 팩토리얼을 구하는 소스코드입니다. 최적화나 어떤 알고리즘을 쓴것은 아니고 막코딩입니다.
10000팩토리얼의 자리수가 3만자리정도 되기때문에 결과값이 저장될 벡터의 크기는 40000으로 잡았습니다.
더 큰수의 팩토리얼을 구하려면 자리수를 늘려주면 됩니다.
하지만 초보자분들이나 안짜보신분들은 꼭 먼저 짜보시고, 성공하신 다음 소스코드를 참조하시기 바랍니다.
그것이 실력향상의 지름길이니까요...^^
#include <iostream>
#include <vector>
using namespace std;
int main()
{
vector<int>* vecData = new vector<int>;
int length = 40000; // 결과값의 자리수
vecData->reserve(length);
vecData->assign(length, -1);
int n = 1000;
int carry = 0;
(*vecData)[0] = 1;
int offset = 1;
for(int i=2;i <= n;i++){
for(int j=0;j < offset;j++){
(*vecData)[j] = (*vecData)[j] * i + carry;
carry = 0;
if((*vecData)[j] > 9){
if((*vecData)[j+1] == -1){
offset++;
(*vecData)[j+1] = 0;
}
carry = (*vecData)[j] / 10;
(*vecData)[j] = (*vecData)[j] % 10;
}
}
}
for(int i=0;i < offset;i++){
cout << (*vecData)[offset-i-1];
}
delete vecData;
return 0;
}
#include <vector>
using namespace std;
int main()
{
vector<int>* vecData = new vector<int>;
int length = 40000; // 결과값의 자리수
vecData->reserve(length);
vecData->assign(length, -1);
int n = 1000;
int carry = 0;
(*vecData)[0] = 1;
int offset = 1;
for(int i=2;i <= n;i++){
for(int j=0;j < offset;j++){
(*vecData)[j] = (*vecData)[j] * i + carry;
carry = 0;
if((*vecData)[j] > 9){
if((*vecData)[j+1] == -1){
offset++;
(*vecData)[j+1] = 0;
}
carry = (*vecData)[j] / 10;
(*vecData)[j] = (*vecData)[j] % 10;
}
}
}
for(int i=0;i < offset;i++){
cout << (*vecData)[offset-i-1];
}
delete vecData;
return 0;
}
Google의 C++ 라이브러리
Thread Programming
