//
// main.cpp
// 高精度阶乘
//
// Created by 张嘉韬 on 16/1/28.
// Copyright © 2016年 张嘉韬. All rights reserved.
//
#include <iostream>
#include <cstring>
using namespace std;
const int maxn=3000;
int main(int argc, const char * argv[]) {
int n,f[3010],c;
cin>>n;
memset(f,0,sizeof(f));
f[1]=1,c=0;
for(int i=2;i<=n;i++)
{
for(int j=1;j<=maxn;j++)
{
int temp=f[j]*i+c;
f[j]=temp%10;
c=temp/10;
}
}
int lenc=maxn;
while(f[lenc]==0) lenc--;
for(int i=lenc;i>=1;i--) cout<<f[i];
cout<<endl;
return 0;
}
总结
1.储存
进行高进度运算,首先我们要分析出我们需要多大的数组来储存我们的数据,我们的高精度是数组中的一个元素来储存大数的一位,所以假如我们要储存的数据是1*10^45次方那么我们只需用45位大小的数组就可以储存的下。
2.乘法运算
从这道题目中我们可以对乘法运算有一个不同于以往直接的高精度乘法的运算方式因为我们的高进度乘法往往是两个大大数相乘,所以只能将将两个数都拆开来运算,但是阶乘运算一般都是一个大数乘以一个小数,所以我们只需要将大数进行拆分和小数直接进行相乘即可,,如547*25 我们可以让5 4 7 分别乘以25,也就是理解成有500个25和40个25再加上7个25,有了这样的方法后我们发现,我们的运算从对需要相乘的两个数都需要每一位进行遍历来运算变为了把一个数做为整体遍历另一个数,这样一来运算就简单了很多。
547
* 25
-------
3.进位处理
乘法运算中每次乘积产生的进位是一个很麻烦的问题,这里提供了一个很好的解决方案,来看这段代码
for(int j=1;j<=maxn;j++)
{
int temp=f[j]*i+c;
f[j]=temp%10;
c=temp/10;
}
我们首先将该位运算得到的乘积和之前的c加起来的到一个和,再对这个和取模后的到了该位的正确结果,然而我们并不着急把所有的进位直接安排到位,而是将temp/10后的到的结果进到下一位的运算中来逐步的将进位分配到位。