Project Euler16.求2^1000的各位数之和

我们知道 2^10=1024,它的各位数之和为1+0+2+4=7。
对于2^1000这样一个庞大的数据,恐怕是我们普通计算机所不能承受的,因此我们不能直接球出2^1000的结果。显然现在我们只有找其他的途径来解决。

之前我们见过这样一个题目:给你100个50位数,求它们和的前十位 。直接对一个50位的数进行加法运算显然也是行不通的,于是我想到是将这100个50位的数存储在一个

a[100][50]的二维数组中,然后在对每一位进行相应的处理。如今我们又碰到了同样的问题,是否能用相同的方法来解决这样的问题呢?答案是肯定的.

定义一个二维数组a[n][m],前面的n代表阶数,后面的m表示的是2^n的位数。m的取值是有一定技巧的,我们知道2^10=1024,约为1000,如是1000^100的位数不超过400。

对于本题而言,我们只须定义一个a[1000][400]的二维数组。这样表述可能还有点抽象,我们选取一个例子:a[1000][400]中a[100][0]表示的具体含义是2^100的个位。

算法实现:

  • 初始化条件:a[n][m] = {2};
  • 递推关系:a[i][j] = 2*a[i-1][j-1]/10 + 2*a[i-1][j]%10,其中1 <= i < n,0 <= j < m。

然而我们现在面临这着一个很严重的问题,a[n][m]中你n,m的取值范围,本台计算m能取到最大值为258(一个32位的机器上,定义的局部变量数组分配内存空间一般不超过1M,作为全局变量而言可分配的内存空间可达到2G,也有说法是3G)。所以定一个a[1000][400]的二维数组是有误的,并且存在大量内存的浪费,如是我们需要另辟蹊径。有人马上会想到STL中容器的概念,当然这种方法行得通,但我们现在使用另外一种方法:定义两个一维数组a[n],b[n]。这样定义的好处是:

  • 解决了内存浪费的现象;
  • 解决了上述下标取值不可达的现象。

下面我们将算法跟新如下:

  1. 初始化条件:a[n]={2};
  2. 递推关系:设0 <= i < 400,0 <= j < 1000;
    • 外部循环条件:0< j < 1000
      • 内部循环条件1:0 <= i < 400
        • b[i]=a[i];
      • 内部循环条件2:0 <= i < 400
        • 若i==0时,a[i]=2 * b[i]  % 10;否则,a[i] = 2 * b[i]  % 10 +  2 * b[i-1]  / 10
    • 程序终止

序代码如下:

#include <IOSTREAM>
#include <TIME.H>
using namespace std;
#define X 1000
#define Y 400
int Fun()
{	
	int a[Y]={2},b[Y];
	int sum=0;
	for (int i=1;i<X;++i)
	{
		for (int j=0;j<Y;++j)
			b[j]=a[j];
		for (int n=0;n<Y;++n)
			if (n==0)
				a[n]=2*b[n]%10;
			else
			    a[n]=2*b[n-1]/10+2*b[n]%10;
	}
	for(int m=0;m<Y;m++)
			sum+=a[m];	
	return sum;
}
int main()
{

	time_t start_time,end_time;
    start_time=clock();
	cout<<Fun()<<endl;
	end_time=clock();
	cout<<"Run Time:"<<(end_time-start_time)<<" ms"<<endl;
	return 0;
}


运行结果:1366 ,Run Time: 5ms




评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值