快速幂算法

时间限制：1.0s 内存限制：256.0MB

问题描述：
　　将一个数N分为多个正整数之和，即N=a1+a2+a3+…+ak，定义M=a1a2a3*…*ak为N的潜能。
　　给定N，求它的潜能M。
　　由于M可能过大，只需求M对5218取模的余数。
　　

输入格式：
　　输入共一行，为一个正整数N。

输出格式
　　输出共一行，为N的潜能M对5218取模的余数。

数据规模和约定
　　1<=N<10^18

用这个题目举例，我们来分析一下快速幂算法在题目中的使用。

首先我们来分析一下这个题目，要使M最大，我们就不能把N分成很大的数相乘，4以上就不行了，5=2+3，5<2*3，6=3+3，6<3*3......所以我们要把N分成尽可能多的小于等于4的数字。继续分析，我们发现我们需要分成最多的3和4，2，分成1是没用的，1乘以任何数都等于那个数本身。

按照上述的方法分解N，N可能是全部由3组成，这是最简单的一种情况，那么我们要求的M就是3的n/3次方，再对5218取模。但是我们还可能遇到N%3==1的情况，这个时候我们就可以把N分成若干个3跟一个4，这样分解N的方法会使M最大，这个时候我们需要计算3的n/3-1次方再乘上4，然后对5218取模，输出结果。同理，遇到N%3==2的情况，我们就把N分成n/3个3和一个2，计算3的n/3次方再乘上2，最后对5218取模。这样我们就分情况讨论了N的不同取值的M不同计算方法。

代码如下：

#include<iostream>
#include<cmath>
using namespace std;
long long int n;
int main()
{
	cin >> n;
	long long int t = 1, m = n / 3, q = n % 3;
	if (q == 1) m--, t *= 4;
	else if (q == 2) t * 2;
	t *= pow(3, m);
	cout << t % 5218;
	return 0;
}

但是我们在编写代码的时候发现，3的次方这个指数函数的增长是爆炸型的，我们无法直接计算出如此大的天文数字，就会让程序算出来的结果是错误的，这个时候我们就需要去想办法解决这个问题。怎么解决呢？

面对这个题目我给出两种解决方案：

第一种我们需要使用取模的公式，(a*b)%c==(a%c*b%c)%c

这就意味着我们不需要在最后对结果取模，而是再计算过程中只要数值过大我们就取模，为了使程序更简单，我们干脆就每次计算时都取一次模，这样的话就能保持数值一直在long long int的范围内了。代码如下：

#include<iostream>
using namespace std;
long long int n;
int main()
{
	cin >> n;   //n是输入的N，t是M，m是N中3的个数，q是N分完3后剩下的数
	long long int t = 1, m = n / 3, q = n % 3;
	if (q == 1) m--, t *= 4;        //这个时候就把N分为若干个3和一个4
	else if (q == 2) t * 2;         //这个时候就把N分为若干个3和一个2
	while (m--) {
		t *= 3;
		t = t % 5218;        //每次t乘以3都做一次取模，保证t的数值一直小于5218
	}
	cout << t;
    return 0;
}

但是我们分析这种算法的时间复杂度，发现是O(n)，那我们能不能降低这种算法的时间复杂度呢？

这个时候我们想到计算幂较大的数时可以使用快速幂，快速乘算法来简化计算，我们来看看怎么在这个题目中学习并运用这种算法。

快速幂顾名思义就是快速计算出幂值较大时的值，我们知道3^4=81，9^2=81，81^1=81，同样的结果，用循环来写，第一种需要循环4次，第二种需要循环2次，第三种只需要循环1次，每次让底数平方，让指数除以2，这样减少了循环的次数，是不是就可以提升我们算法的性能了呢？对于上面的while循环我们只需要logn次就可以算出结果，那么我们的时间复杂度是不是就降低了呢？

代码如下：

long long int a = 3;
while (m != 0) {
	 if (m % 2) t = (t*a)%5218;  //如果m是奇数，我们将t乘以一遍a；
	 a = (a*a)% 5218;            //每次都让平方一遍（为了防止a的值溢出，我们每次都取一次模）
	 m /= 2;
	}

当然我们在这里也可以用位运算判断m的大小和除以2的操作

代码如下：

long long int a = 3;
while (m != 0) {
	 if (m&1) t = (t*a)%5218;  //m&1可以判断m的奇偶性
	 a = (a*a)% 5218;
	 m>>=1;    //每次二进制数向右移动一位就等同于除以2，余数会被省略
	}

这样的话我们就解决了数值过大溢出的问题，也学习到了快速幂的用法。