问题描述
输入一个正整数n,输出n!的值。
其中n!=123*…*n。
算法描述
n!可能很大,而计算机能表示的整数范围有限,需要使用高精度计算的方法。使用一个数组A来表示一个大整数a,A[0]表示a的个位,A[1]表示a的十位,依次类推。
将a乘以一个整数k变为将数组A的每一个元素都乘以k,请注意处理相应的进位。
首先将a设为1,然后乘2,乘3,当乘到n时,即得到了n!的值。
输入格式
输入包含一个正整数n,n<=1000。
输出格式
输出n!的准确值。
样例输入
10
样例输出
3628800
思路:模拟平常列竖式的算法,将每个数存进数组,按照错位相乘的原则进行计算。
#include <math.h>
#include <stdio.h>
#include <algorithm>
#include <string>
#include <vector>
#include <iostream>
using namespace std;
int a[1001][3001];//前面表示阶数 后面表示位数
int n,len=1;//len为最高位数
void per(int n1)//高精度乘法
{
int b[10],t=0,m=n1,i,j;
while(m) //将数m放进数组里 并求它的数位数t
{ b[t++]=m%10;
m=m/10;
}
for(i=0;i<t;i++) //第二个(即被乘的数)的数位长度
for(j=0;j<len;j++)//第一个数的数位长度
{ a[n1][i+j]+=a[n1-1][j]*b[i];//错位相乘
}
for(i=0;i<len+t;i++)//两个数的乘积后位数一定不超过两个数的位数和
{
if(a[n1][i]>=10) //进位处理
{ a[n1][i+1]+=a[n1][i]/10;
a[n1][i]=a[n1][i]%10;
}
}
int ln1=len+t;
while(!a[n1][ln1])//求最高位数 将后面的多余0删掉
ln1--;
len=ln1+1; //赋给len,因为从0开始,判断的时候以i<len判断,所以+1
}
int main()
{
int i;
cin>>n;
a[1][0]=1;
for(i=2;i<=n;i++)
{
per(i);
}
for(i=len-1;i>=0;i--)
cout<<a[n][i];
return 0;
}