求1000的阶乘

引用 50 楼 runffer_yang 的回复:

引用 46 楼 groovy2007 的回复:
不是有BigInteger吗，只要内存够用，再大的数也能算


如果你去微软或者谷歌面试,那么很遗憾,面试官不会让你用BigInteger或者BitDecimal,
事实上这个题目不是考用什么类型,而是考高精度计算的算法+数据结构.

------
由于 int 类型数据长度的限制，采用每位数字与乘数相乘的做法我认为并不好，
也就是说使用该方法做出的阶乘计算受到 int 类型最大数字的限制，并不能称
为大数计算。

所谓的大数计算不管是 200 亿的阶乘还是 2 的阶乘在计算上除了内存空间的限
制之外，不应受到编程语言数值类型最大数的限制。
------
牛淫
------
#include<stdio.h> 
#define N 1000 //要算的阶乘数 N!
long s[N]={1,1},n=N,t=2,a=1,b=0; 
int main()
{ 
  for(; a <= *s || (++t <= n ? (b = 0, a = 1) : 0); (*s == a++ && b) ? (*s)++ : 0) 
  s[a]=(b+=s[a]*t)%10000,b/=10000; 
  for(printf("%d",s[*s]);--*s>0;)
printf("%04d",s[*s]); 
printf("\n");
  return 0; 
} 


这是一个北大大二的女生写的，我看不懂。可以很快求出N的阶乘！
------
可以运用指针或文本来输出，自己去学吧，其实不难！
------
给你一个关于求PI的任意高阶精度值得算法，你看看里面的方法！
#include <iostream>
#include <fstream>
#include <ctime>
#include "windows.h"
using namespace std;

void arctg(int * arctgx, int x, int len)
{
int n = 1, s = 1,c = 4,d = 0;
int i, *t, *a, eps = 0;
t = new int[len+1];
a = new int[len+1];
for(i=0;i<len;i++)
{
t[i] = c / x;
d = c % x;
c = d*10;
}
while(eps<len-1)
{
d = 0;
for(i=0;i<len;i++)
{
c = t[i]+d*10;
a[i] = c / n;
d = c % n;
}
for(i=1;i<len;i++)
{
if (a[i])
{
eps = i;
break;
}
else eps++;
}
if(s==1)
{
for(i=len-1;i>0;i--)
{
c = arctgx[i] + a[i];
arctgx[i] = c%10;
arctgx[i-1]+= c/10;
}
}
else
{
for(i=len-1;i>0;i--)
{

arctgx[i]-=a[i];
if (arctgx[i]<0)
{
arctgx[i] += 10;
arctgx[i-1]-=1;
}
}
}
n = n+2;
s = -s;
d = 0;
for(i=0;i<len;i++)
{
c = t[i]+d*10;
t[i] = c / x;
d = c % x;
}
d = 0;
for(i=0;i<len;i++)
{
c = t[i]+d*10;
t[i] = c / x;
d = c % x;
}
}
delete []t;
delete []a;
}


int main()
{
DWORD dwStart,dwStop; //调用系统功能来计时
dwStart=GetTickCount(); //开始计时
int c, d, len;
int decimal, *arctg5, *arctg239, *pi;
cout << "请输入所要计算的π的精度值:";
cin >> decimal;
len = decimal+10;
arctg5 = new int[len+1], arctg239 = new int[len+1], pi = new int[len+1];
cout << "计算进行中，请稍等……" << endl;
for(int i=0;i<len;i++)
{
arctg5[i] = 0;
arctg239[i] = 0;
}
arctg(arctg5, 5, len);
arctg(arctg239, 239, len);
d = 0;
for(i=len-1;i>=0;i--)
{
c = arctg5[i]*4 + d;
pi[i] = c % 10;
d = c / 10;
}
delete []arctg5;
for(i=len-1;i>0;i--)
{
pi[i]-=arctg239[i];
if (pi[i]<0)
{
pi[i] += 10;
pi[i-1]-=1;
}
}
delete []arctg239;
cout<< "π= " << pi[0] << '.';
  for(i=1;i<=decimal;i++)
{
cout<<pi[i];
}
delete []pi;
dwStop=GetTickCount(); //计时结束
cout<<endl;
cout<<"计算用时:"<<(double)(dwStop-dwStart)/1000<<"秒"<<endl;
return 0;
}

------
for循环啊！
不过要注意数据的类型！
double int 可能可以
------