解题思路:
由于n太大需要考虑降幂,利用扩展欧拉定理可以得出a^(p%phi)%mod = a^p%mod,因此我们只需要计算p%phi即可。
讲明两个坑点吧:
1、做了很多矩阵快速幂,做吐了之后忽视了矩阵的转移方式,应该是1,0作为矩阵的初始状态而不是0,1.写出来的话应该不是有这种错误的。
2、再一个是费马小定理推导扩展欧拉的时候是有一个取模phi,这个phi是mod的欧拉函数,因为是素数可以使用费马推到。得到了扩展欧拉的其中一条。如上所示。我们这题的主要思路是由费马得知是指数循环节。因此取模上不是%mod而是%phi。因此利用矩阵快速幂求解指数的值的时候,我们不应该%mod。我一直没转过来弯。
3、扩展欧拉定理有适用范围,尽量不要越界使用,,本题就不可以。
AC代码:
#include <iostream>
#include <vector>
#include <cstring>
#include <string>
#include <cmath>
#include <set>
using namespace std;
typedef long long ll;
typedef vector<ll> vec;
typedef vector<vec> mat;
const int mod = 1000000007;
const int modd = 1000000006;
ll n, a, b;
ll cnt, cns;
ll ZS;
mat mul(mat &A, mat &B) {
mat C(A.size(), vec(B[0].size()));
for (int i = 0; i < A.size(); i++) {
for (int k = 0; k < B.size(); k++) {
for (int j = 0; j < B[0].size(); j++) {
C[i][j] = (C[i][j] + (A[i][k]*B[k][j])%modd+modd)%modd;
}
}
}
return C;
}
mat mat_quick_mod(mat A, ll n) {
mat B(A.size(), vec(A.size()));
for (int i = 0; i < A.size(); i++) {
B[i][i] = 1;
}
while(n) {
if(n&1) {
B = mul(B, A);
}
A = mul(A, A);
n >>= 1;
}
return B;
}
ll quick_mod(ll a, ll n){
ll ans = 1;
while(n) {
if(n&1) ans = ans*a%mod;
a = a*a %mod;
n >>= 1;
}
return ans;
}
ll euler(ll x) {
ll ans = x;
for (int i = 2; i < sqrt(x) + 1; i++) {
if(!(x%i)) {
ans -= ans/i;
while(!(x%i)) x/=i;
}
}
return x>1 ? ans - ans/x : ans;
}
void init() {
ll phi = euler(mod);
n = n%phi;
}
ll solve(ll x, ll y) {
mat A(2, vec(2));
A[0][0] = 1;
A[0][1] = 1;
A[1][0] = 1;
A[1][1] = 0;
A = mat_quick_mod(A, n-2);
mat C(2, vec(1));
C[0][0] = x;
C[1][0] = y;
C = mul(A, C);
return C[0][0];
}
int main() {
while(cin >> a >> b >> n) {
init();
if(n == 0) {
printf("%lld\n",a%mod);
}
else if(n == 1){
printf("%lld\n",b%mod);
}
else {
cnt = solve(1, 0);
cns = solve(1, 1);
printf("%lld\n",(quick_mod(a, cnt)*quick_mod(b, cns))%mod);
}
}
}