大数阶乘 c++

就是求n！

（一）

//为了明白演示，将max改为4，并以n=4为样例 //由最后一句可以看出b是用来存储结果的，由高位到低位；b[max-1]为个位，b[max-2]为十位，依此类推 //a和b存储的方式一样，但是用来存储中间结果的 //作者没有使用乘法，显然是要以加法来替代的。 #include "stdafx.h" #include <iostream> #define max 1000//实际应用时，这个要改成比较大的数，如1000；那就表示能表达多达1000位的数 using namespace std; void Print(int array[])//打印 { for(int i=0;i<max;i++) cout<<array[i]<<" "; cout<<endl; } long NFai(long n)//求n!，作对比用 { int res=1; for(long i=2;i<=n;i++) { res=res*i; } return res; } void main() { int a[max],b[max]; int i,j,k,l=0,n; cout <<"input a number:" <<endl; n = 13; //以输入4为样例。 //初始化a，b for(i=0;i <max-1;i++) { a[i]=0;b[i]=0; } a[max-1]=1;b[max-1]=1; for(i=1;i <=n;i++) //对于1到n每个数 { for(j=max-1;j>=0;j--) //由低位到高位 { for(k=1;k <i;k++) //假如i=3，那此时a就是存储以前的结果3=1*2 { int f=0;//表示进几次位 while(1) { if((b[j]+a[j]-f*10) <10) //如果不够进位 { b[j]=b[j]+a[j]-f*10;//求本位的值 break; } else //如果够进位 { b[j-1]=b[j-1]+1; //高位加1 f++; } } } //cout<<"a="; //Print(a); //cout<<"b="; //Print(b); a[j]=b[j]; //将本位的值保存到a } cout<<"数"<<i<<"处理完毕"<<endl; } while(b[l]==0)l++; //越过高位为0的部分 cout<<"结果："; while(l <max) //这看出b是用来存储结果的，由高位到低位；b[max-1]为个位，b[max-2]为十位，依此类推 {cout<<b[l];l++;} cout<<endl; cout<<NFai(n)<<endl; }

对于
for(k=1;k <i;k++) 这句，
假如此时i=3，j=max-1，
那新个位 就是 3个以前结果的个位相加并进位后剩下的个位。（循环开始时a[max-1],b[max-1]都是以前结果的个位）

（二）

转自http://www.programfan.com/blog/article.asp?id=568

#include<iostream> #include<stack> #include<fstream> using namespace std; struct Node { int i; Node *next; }; Node *createNode() { Node *node=new Node; node->i=0; node->next=NULL; return node; } Node* MulLink(Node *head,int i) { Node *p=head; head->i=1; p=head; int temp=0; for(int j=1;j<=i;j++) { while(p) { p->i*=j; p->i+=temp; if(p->i>9) { temp=p->i/10; if(!p->next) p->next=createNode(); p->i%=10; p=p->next; } else { temp=0; p=p->next; } } p=head; } return head; } void displayLink(Node *head) { fstream outfile; stack<int> display; outfile.open("example.txt",ios::out|ios::ate); while(head) { display.push(head->i); head=head->next; } while(!display.empty()) { cout<<display.top(); outfile<<display.top(); display.pop(); } cout<<endl; } void freeLink(Node *head) { Node *p=head; while(p) { head=p->next; delete p; p=head; } } int main(int argc,char *argv[]) { int n; cout<<"Please input the number"<<endl; cin>>n; Node *head=createNode(); MulLink(head,n); displayLink(head); freeLink(head); system("PAUSE"); return 0; }

（三）

转自http://download.csdn.net/source/629944

#include <stdlib.h> #include <time.h> #include <stdio.h> #include <string.h> #define MAXSIZE 100000 #define RADIX 10000 #define LOG10_RADIX 4 void Change(long a[] , int &m) ; int Factorial(long, long a[]); int main() { long a[MAXSIZE] = {0} ; int n; printf("Input a number: "); scanf("%d", &n); long t; t = clock(); int m = Factorial(n, a); t = clock()-t; printf("%d", a[m]); for(int l=m-1 ; l>=0 ; l--) //输出结果 { printf("%04d", a[l]) ; if((l*LOG10_RADIX) % 10000 == 0) system("pause"); } printf("/n"); printf("Time cost: %dms/n", t); system("pause") ; return 0; } void Change(long a[] , int &m) { int b = 0 ; //b表示低位向高位的进位 for(int i=0 ; i<m ; i++) { a[i] += b ; if(a[i] < RADIX) //当a[i]<9时，说明这里不做进位处理 // <= { b = 0 ; } else if(a[i] >= RADIX) //当a[i]不是最高位时，满足条件只进位 { b = a[i] / RADIX ; a[i] -= b* RADIX ; } } a[i] += b; while(a[i] >= RADIX) { b = a[i] / RADIX ; a[i] -= b* RADIX ; i++ ; a[i] = b ; ++m; // 增加一个结点 } } // 求n的阶乘，结果存放到a数组中 int Factorial(long n, long a[]) { memset(a, 0, MAXSIZE*sizeof(long)); a[0] = 1 ; int m =0 , b=0 ; for(int i=2 ; i<=n ; i++) { for(int j=0 ; j <= m ; j++) //把每一个元素分别乘以i，并保存在原来的位置 { a[j] = a[j] * i ; } Change(a, m) ; //对上面的每个元素分别做调整 ,每乘一次做一次调整 } return m; }

（四）

转自http://hi.baidu.com/peng3409/blog/item/87917cd997b747ec38012f6d.html

HDU 1042—题目要求计算10000！大约有35000~36000位，用上篇介绍的方法是满足不了需求的。分析，我们可以知道计算时所需要的时间耗费在进位处理上，由于我们使用10进制方式，结果导致计算长度len很大。为了减少计算时间，做一点改动，将进位改为10^10，这样进位长度大大缩减，由于10^10超过int的范围，我们使用__int64（注意，两个下划线）。提交到杭州电子科技大学的ACM判题系统显示耗时1720MS，如果想进一步减少时间可以通过扩大基数来实现，不过原理就类似了。

#include <iostream> #include <cmath> using namespace std; const int max_len=36000;// const __int64 base5=100000; const __int64 base=base5*base5; __int64 x[max_len]; int len; void fact(int n) { __int64 carry; int i,j; for(i=0;i<max_len;++i) x[i]=0; x[0]=1; len=1; for(i=1;i<=n;++i) { carry=0; for(j=0;j<len;++j) { x[j]=carry+x[j]*i; carry=x[j]/base; x[j]=x[j]%base; } while(carry) { x[j++]=carry%base; carry=carry/base; len++; } } } int main() { int n; while(scanf("%d",&n)!=EOF&&n>=0) { fact(n); int i; printf("%I64d",x[len-1]); for(i=len-2;i>=0;i--) {//Èôx[i]²»×ã9λÊý£¬Ç°Ãæ²¹Áã¡£ if(x[i]<10) printf("000000000"); else if(x[i]<100)printf("00000000"); else if(x[i]<1000)printf("0000000"); else if(x[i]<10000)printf("000000"); else if(x[i]<100000)printf("00000"); else if(x[i]<1000000)printf("0000"); else if(x[i]<10000000)printf("000"); else if(x[i]<100000000)printf("00"); else if(x[i]<1000000000)printf("0"); printf("%I64d",x[i]); } printf("/n"); } return 0; }