NYOJ 88 汉诺塔（一）（快速幂）

NYOJ 88 汉诺塔（一）

解题思路：根据汉诺塔的函数表达式，得出放不同的金片数，所需对应的移动步数。函数表达式，我上网搜了下，函数推到如下：

计算移动次数：
如果要计算一共移动了多少次，找出规律即可。
假设移动n个盘子需要移动f(n)次，所以把n-1个盘子移动到b柱子上，需要移动f(n-1)次，然后把第n个盘子移动到c柱子上，需要移动1次，最后把n-1个盘子移动到c柱子上，需要移动f(n-1)次，综上所述，一共移动了
                       f(n) = 2 f(n-1) + 1
  而：
  f(1) = 1;
  f(2) = 2*1+ 1;
  f(3) = 2(2*1+ 1)+ 1;
  .......

  f(n) = 2^n -1.

但有了公式，不是就万事大吉了。当你愉快的写完，编译连接运行，输入测试数据，全对；于是快速提交，然而显示你超时了。。。为神马呢？？？

原因是你用乘法遍历的时候，后台的测试数据n最大的有10亿，如果算下来，是非常慢的，所以也就很容易超时。具体我们来看个例子，你就懂了。

设A为矩阵，求A的N次方，N很大，1000000左右吧。。。
先看小一点的，A的9次方
A^9
= A*A*A*A*A*A*A*A*A  【一个一个乘，要乘9次】
= A*(A*A)*(A*A)*(A*A)*(A*A)【保持格式的上下统一，所以加上这句】
 = A*(A^2)^4 【A平方后，再四次方，还要乘上剩下的一个A，要乘6次】
= A*((A^2)^2)^2【A平方后，再平方，再平方，还要乘上剩下的一个A，要乘4次】

也算是一种二分思想的应用吧，1000000次幂，暴力要乘1000000次，快速幂就只要（log2底1000000的对数） 次，大约20次。

因此，我们需要对乘法进行优化，于是我们就用到了快速幂。这里给大家一个百度百科关于快速幂的链接：快速幂算法。

于是，我们对代码进行修改，此时，又发现当输入111时，尽然输出-1，输入1111更是如此。把中间值base输出来，看看，发现到6万后，全部都输出0了。因为中间变量base是进行乘法运算的，很容易超出整型，所以要用双长整型才行（长整型最大为10亿，但中间变量最大可取100万，两个一百万相乘，结果肯定超过10亿）。此处改完后，就顺利被AC了。

#include<stdio.h>
#define MOD 1000000
int main(){
	long long int N,m;
	scanf("%d",&N);
	while(N--){
		scanf("%lld",&m);
		long long int sum = 1,base = 2;
		while(m != 0){
			if(m%2)  sum = sum*base%MOD;
			base *= base;
			base %= MOD;
//			printf("%lld.....\n",base);
			m /= 2;
		}
		printf("%lld\n",sum-1);
	}
	return 0;
}

或者把快速幂写成函数形式：

#include<stdio.h>
#define MOD 1000000
int fastpow(long int a,long int b){
	long long int result = 1,base = a;
	while(b != 0){
		if(b%2) result = result*base%MOD;
		base *= base;
		base %= MOD;
		b /= 2;
	}
	return result;
}

int main(){
	long int N,m;
	scanf("%ld",&N);
	while(N--){
		scanf("%ld",&m);
		printf("%lld\n",fastpow(2,m)-1);
	}
	return 0;
}