题35.C++练习-3-3快速幂 等比数列求和

---

题35.C++练习-3-3快速幂 等比数列求和

---

一、题目

二、题解

本题若使用等比数列求和公式+快速幂当然会因为幂结果数据过大导致求余gg然后拉闸了。。。所以只能在用原公式计算过程中进行快速幂+求余操作并必须要让这个过程尽可能少。那就想想能不能来波分治吧，对原公式进行改造，分成左半段和右半段，但这样由于power奇偶会导致分的段有点不同，所以我们再分一波奇偶，处理后这些只是一堆递归而已，而且每次都将过程减少了一半，效率大大提高了。

#include <bits/stdc++.h>
#define Mod 1000000007
using namespace std;

typedef long long ll;
ll fastPower(ll base, ll power, ll mod) {
	ll result = 1;
	while (power > 0) {
		if (power & 1) result = result * base % mod;
		power >>= 1;
		base = base * base % mod;
	}
	return result;
}

/*
当power为偶数时，项数为奇，
sum
=1+base ^ 1+base ^ 2+...+base ^ power
=1+base ^ 2 +…+base ^ (power/2-1)+base ^ (power/2)+...+base ^ (power-1)+base ^ (power)
= (1+base ^ (power/2)) * (1+base ^ 1 +…+base ^ (power/2-1))+base ^ power。
当power为奇数时，项数为偶，同理，
sum
=(1+base ^ ((power+1)/2)) * (1+base ^ 1+base ^ 2+...+base ^ ((power-1)/2))。
*/
ll sum(ll base, ll power, ll mod) {
	if (power == 0)
		return 1;
	else if (power & 1) {
		return ((1 + fastPower(base, (power + 1) / 2, mod)) * sum(base, (power - 1) / 2, mod)) % mod;
	}
	else
		return ((1 + fastPower(base, power / 2, mod)) * sum(base, power / 2 - 1, mod) + fastPower(base, power, mod)) % mod;
}

int main()
{
    ll a,n;
    cin>>a>>n;
    cout<<sum(a,n,Mod);
}

取模”运算的运算法则：
1.(a + b) % p = (a % p + b % p) % p
2.(a - b) % p = (a % p - b % p ) % p
3.(a * b) % p = (a % p * b % p) % p
4.(a * b * c) % p = (a % p * b % p * c % p) % p …
（利用3，4…只需要在循环乘积的每一步都提前进行“取模”运算可得到整个乘积结果的模，注意%和*，/是同一优先级）

快速幂：
每一步都把指数分成两半，而相应的底数做平方运算（而碰到奇数底的时候本质其实是抽出了一个1，然后剩下的指数自然是偶数然后才/2），到指数为0为止。（点我，康康牛人的快速幂文章）

//快速幂
long long fastPower(long long base, long long power) {
    long long result = 1;
    while (power > 0) {
        if (power & 1)//if(power%2==1) 
        {
            result = result * base;
        }
        power >>= 1;//power=power/2
        base = (base * base);
    }
    return result;
}

---