汉诺双塔问题python_【汉诺塔进阶】递归与高精度（下）

上一个专栏提到，我们应该用高精度算法来存储大数，那么什么是高精度算法呢？精度算法(High Accuracy Algorithm)就是处理大数字的数学计算方法。在一般的科学计算中，会经常算到小数点后几百位或者更多，当然也可能是几千亿几百亿的大数字。一般这类数字我们统称为高精度数，高精度算法是用计算机对于超大数据的一种模拟加，减，乘，除，乘方，阶乘，开方等运算。具体的实现过程就是定义一个数组，数组中每一个数存储一位，这样你的数组的长度就是你能存储的最大的位数。这样一来就大大提高了数据存储的精度。

回到“汉诺双塔问题”，我们上一节用了递推的算法使得速度大大提升，大家可以回顾我上一篇文章。

虽然我们解决了时间问题，但是数据依然会很快溢出。所以这就是典型的需要运用高精加和高精乘算法来实现的例子。

下面贴出我的代码，具体的解释在注释中。

#include

#include

using namespace std;

int a[201]={0},j,len=1;//定义201位数。

void gjj()//高精加

{

a[1]+=2;

j=1;

while(a[j]>=10&&j

{

a[j]-=10;

a[j+1]++;

j++;

}

if(a[len]>=10)//如果最高位进位长度加1

{

a[len]-=10;

a[++len]++;

}

}

void gjc()//高精乘

{

for(j=1; j<=len; j++)//把每一位都乘2

{

a[j]*=2;

}

for(j=1; j

{

if(a[j]>=10)

{

a[j]-=10;

a[j+1]++;

}

}

if(a[len]>=10)//如果最高位进位长度+1

{

a[len]-=10;

a[++len]++;

}

}

int main()

{

int i,n;

cin>>n;

for(i=1; i<=n; i++)

{

gjc();

gjj();

}

for(j=len;j>=1;j--)

{

cout<

}

return 0;

}

下面是当n=200时的运算结果：运算结果

由此可见我们成功的存储了一位高精度数并输出了出来。

下面总结利用高精度算法的要点：估计最长数据的位数len_MAX，并初始化一个数组a[len_MAX]={0}；，初始化最初的长度len=0；

从最低位开始对每一位分别运算。

从最低位到最高位的前一位，对超过10的位数进行进位。

最后判断最高位是否超过10.若超过则进行进位，并将len++；

以上就是高精度算法的要点，下面练习一道进阶例题。

P1591 阶乘数码

题目描述

求 n! 中某个数码出现的次数。

输入格式

第一行为 t(t≤10)，表示数据组数。接下来 t 行，每行一个正整数 n(n≤1000) 和数码 a。

输出格式

对于每组数据，输出一个整数，表示 n! 中 a 出现的次数。

输入输出样例

输入 #1                              输出 #1

2                                         1

5 2                                       2

7 0

源代码如下：(附加了求阶乘功能)

#include

#include

using namespace std;

int a[11][10001]= {0},len[11]= {0},flag[11]= {0},j;

void gjc(int n,int m)//高精乘

{

for(j=1; j<=len[m]; j++)

{

a[m][j]*=n;

}

for(j=1; j

{

if(a[m][j]>=10)

{

while((a[m][j]>=10))

{

int b=a[m][j]/10;

a[m][j]%=10;

a[m][j+1]+=b;

}

}

}

if(a[m][len[m]]>=10)//如果最高位进位长度+1

{

while(a[m][len[m]]>=10)

{

int b=a[m][len[m]]/10;

a[m][len[m]]%=10;

len[m]++;

a[m][len[m]]+=b;

}

}

}

void fact(int n,int m)//求阶乘

{

a[m][1]=1;

while(n)

{

gjc(n,m);

n--;

}

}

void find_num(int k,int m)//找数

{

for(j=1; j<=len[m]; j++)

{

if(k==a[m][j])

{

flag[m]++;

}

}

}

int main()

{

int n[11],k[11],m,i;

cin>>m;

for(i=1; i<=m; i++)

{

cin>>n[i]>>k[i];

}

for(i=1; i<=m; i++)

{

len[i]=1;

fact(n[i],i);

find_num(k[i],i);

cout<< flag[i]<< endl;

}

for(int l=1; l<=m; l++)

{

cout<

for(i=len[m]; i>0; i--) //附加的验证阶乘的函数

{

cout<