JZOJ Day 3-B组-T1——学习神技

题目大意：

已知：一个等比数列共有n项，第一项为$A_1$，公比为q
求：$$S_n=\frac{A_1\times (1-q^n)}{1-q}(q̸=1)or{A}_1\times n(q=1)$$

解题思路：

等比数列求和公式$q̸=1$的情况证明：
证：
设$$S=A_1+A_1\times q+A_1\times q^2+...+A_1\times q^n$$ ①
两边同乘$q$，得
$$S\times q=A_1\times q+A_1\times q^2+A_1\times q^3...+A_1\times q^{n+1}$$ ②
用②$-$①得：
$$S\times q-S=A_1\times q^{n+1}-A_1$$
因式分解得：
$$S\times (q-1)=A_1\times (q^{n+1}-1)$$
移项得：
$$S_n=\frac{A_1\times (q^n-1)}{q-1}$$
∵$q̸=1$且$q>1$
所以我们就有：$$\frac{A_1\times (q^n-1)}{q-1}=\frac{A_1\times (1-q^n)}{1-q}$$
证毕。

那啥，上面的证明应该是①$-$②但是我懒得改了
所以在后面我就加上了①$-$②$=$②$-$①的证明
然后正题来了：
快速幂求乘法逆元随便搞

$Acceptedcode$：

#include<cstdio>
#include<iostream>
#define int long long

using namespace std;

const int Mod = 1e9+7;

int A, Q, N, mod;

inline int read() {
	int f = 0; char c = getchar();
	while (!isdigit(c)) c = getchar();
	while (isdigit(c)) f = (f<<1) + (f<<3) + c - 48, c = getchar();
	return f;
}

int ksm(int a, int M) {
	int ans = 1; a %= Mod;
	while (M) {
		if (M & 1) ans = (ans * a) % Mod;
		a = (a * a) % Mod;
		M >>= 1;
	}
	return ans;
}

void write(int x) {
	if (x > 9) write(x/10); putchar(x%10+48); return;
}

void writeln(int x) {
	if (x < 0) x = -x, putchar('-'); write(x); putchar(10); return;
}

signed main() {
	int T = read();
	for (;T--;) {
		A = read(), Q = read(), N = read();
		if (Q == 1) {
			writeln((A * (N % Mod)) % Mod);
		} else {
			mod = ksm(1 - Q, Mod-2);
			writeln(A * (1 - ksm(Q, N)) % Mod * mod % Mod);
		}
	}
	return 0;
}