费马小定理·欧拉降幂·矩阵快速幂·M斐波那契数列 HDU - 4549

解题思路：

由于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);
        }
    }
}