快速幂算法与例题

算法思想：
首先把指数部分化为二进制，比如2^63若是正常计算需要计算63次，那么如果化为 2 ^111111,这样计算的话就是
2^32 * 2^16 * 2^8 *2^4 * 2^2 * 2^1 只需要计算6次，明显比普通的算法快很多。

取模运算的公式：
(a*b)%mod=(a%mod * b%mod )%mod
(a+b)%mod=(a%mod + b%mod) %mod

代码实现：

long long quickpow(long long a, long long b ,long long mod)//a代表底数，b代表指数
{
	long long res = 1;//代表输出结果
	while (b)//当b为0时不进入计算
	{
		if (b & 1) // 位运算符，来判断二进制的b的第一位是否为1，为1则进入计算，为0则不进入计算
			res = res * a % mod;
		a = a * a % mod;
		b >>= 1;//b左移一位
	}
	return res;
}

例题：
已知一个等比数列的某三项分别是a,b,c 且已知第一项是a，求等比数列的第N项最大是多少。

输入：
输入共一行，输出4个整数a,b,c,N

输出：
输出第N项，结果对10^9取模

分析：
结果要对10^9取模，因此需要用快速幂算法。
首先找最简单的情况：a=b或者b=c的时候，说明该数列为常数列。
b/a==c/b的时候，那么公比就为b/a
接下来如果取b/a和c/b中大的为max，小的为min，如果(max/min)%min!=0的话，那么需要对min进行开方，继续判断。

代码实现：

#include<iostream>
using namespace std;
long long quickpow(long long, long long, long long);
int main()
{
	long long  a, b, c, N;
	long long min;
	long long mod = 1e9;//10^9；
	cin >> a >> b >> c >> N;
	if (a == b || b == c || (b / a == c / b))
	{
		cout << a*quickpow(b/a,N-1,mod)%mod;
	}
	else
	{ //先找出两者最小的一个
		if (b / a < c / b)
			min = b / a;
		else
			min = c / b;
		while (((c / b) / min) % min != 0)//如果不为0，开方后继续判断
		{
			min = sqrt(min);
		}
		cout << a * quickpow(min, N - 1,mod) % mod;
	}
	return 0;
}
long long quickpow(long long a, long long b ,long long mod)//快速幂算法
{
	long long res = 1;
	while (b)
	{
		if (b & 1)
			res = res * a % mod;
		a = a * a % mod;
		b >>= 1;
	}
	return res;
}