c語言之大數階乘

大數階乘

時間限制：3000 ms  |  內存限制：65535 KB

難度：3

描述

我們都知道如何計算一個數的階乘，可是，如果這個數很大呢，我們該如何去計算它並輸出它？

輸入

輸入一個整數m(0<m<=5000)

輸出

輸出m的階乘，並在輸出結束之後輸入一個換行符

樣例輸入

50

樣例輸出

30414093201713378043612608166064768844377641568960512000000000000

思路：

第一步：輸入數字

第二步：計算大概存儲多少位

第三步：根據估算的位數創建數組，分配存儲空間

第四步：從2開始到n（n爲輸入的要求的階乘數）進行循環，計算每個數組元素中要存的數 [在此需要一個進位函數]

{

進位函數：

主要是對每一位進行判斷

if(該位上的數字小於等於9)

不進位

else if(該位大於9並且該位不是最高位)

保存該位的進位數

保存該位

else if(該位大於9並且是最高位)

while（只要該位大於9）

{

保存進位；

保存該位；

下一位的賦值；

}

}

第五步輸出：（輸出之前先去確定位數[確定位數的方法，從估算位開始向下尋找，只要滿足該位不是0，即課確定位數]）從高向低一次輸出即可

具體代碼如下：

#include<stdio.h>
#include<malloc.h>
#include<math.h>

void arry(int fact[],int pos)
{
	int i,array = 0;
	for(i=0;i<=pos;i++)
	{
		fact[i] += array;
		if(fact[i] <=9)
			array = 0;
		else if(fact[i] >9 && i<pos)
		{
			array = fact[i] / 10;
			fact[i] = fact[i] % 10;
		}
		else if(fact[i]>9 && i>=pos)
		{
			while(fact[i] > 9)
			{
				array = fact[i] / 10;
				fact[i] = fact[i] % 10;
				i++;
				fact[i] = array;
			}
		}
	}
}

int main()
{
	int num,digit,pos;
	double sum=0;
	int * fact;
	int i,j;

//	printf("請輸入要計算的數字：");
	scanf("%d",&num);

	for(i=1;i<=num;i++)
		sum += log10(i);
	
	digit = (int)sum +1;

//	printf("%d",digit);

	fact = (int *)malloc(sizeof(int) * (digit+1));

	for(i=0;i<=digit;i++)
		fact[i] = 0;
	fact[0] = 1;

	for(i=2;i<=num;i++)
	{
		for(j=digit;j>=0;j--)
			if(fact[j] != 0)
			{
				pos = j;
				break;
			}

		for(j=0;j<=pos;j++)
			fact[j] *= i;
		arry(fact,pos);
	}

	for(i=digit;i>=0;i--)
	{
		if(fact[i] != 0)
		{
			pos = i;
			break;
		}
	}


	for(i=pos;i>=0;i--)
	{
		printf("%d",fact[i]);
	}

	printf("\n");

	return 0;
}