noi-1733-求10000以内n的阶乘

14:求10000以内n的阶乘

• 
  

总时间限制: 

5000ms 

内存限制: 

655360kB

描述

求10000以内n的阶乘。

输入

只有一行输入，整数n（0<=n<=10000）。

输出

一行，即n!的值。

样例输入

100

样例输出

93326215443944152681699238856266700490715968264381621468592963895217599993229915608941463976156518286253697920827223758251185210916864000000000000000000000000

来源

JP06

//noi_1733
#include
#include
#include
#include
#include

using namespace std;

int c[100005]={0};

void factoriala(int a[],int b){
	for(int i=1;i<=a[0];i++){
		a[i]=a[i]*b;
	}
	for(int i=1;i<=a[0];i++){
		a[i+1]+=a[i]/10;
		a[i]=a[i]%10;
	}
	while(a[a[0]+1]>0){
		a[0]++;
		a[a[0]+1]=a[a[0]]/10;
		a[a[0]]=a[a[0]]%10;
	}
}

void printa(int a[]){
	for(int i=1;i<=a[0];i++){
		cout<>n;
	if(n==0) ;
	else for(int i=1;i<=n;i++){
		factoriala(c,i);
	}
	printa(c);
	return 0;
}

效率较低，还要优化。

万进制优化

//noi_1733_beta1 
#include
#include
#include
#include
#include

using namespace std;

int c[10005]={0};

void factoriala(int a[],int b){
	for(int i=1;i<=a[0];i++){
		a[i]=a[i]*b;
	}
	for(int i=1;i<=a[0];i++){
		a[i+1]=a[i+1]+a[i]/10000;
		a[i]=a[i]%10000;
	}
	while(a[a[0]+1]>0){
		a[0]++;
		a[a[0]+1]=a[a[0]]/10000;
		a[a[0]]=a[a[0]]%10000;
	}
}

void printa(int a[]){
	for(int i=1;i<=a[0];i++){
		if(i==1) cout<>n;
	c[0]=1;
	c[1]=1;
	if(n==0) ;
	else for(int i=1;i<=n;i++){
		factoriala(c,i);
	}
	printa(c);
	return 0;
}