hdu-5667（矩阵快速幂）

前两天就已经学习了矩阵快速幂，当时的我就学了个矩阵快速幂的模板，而且我还以为我自己学会了，实在是太搞笑了，矩阵快速幂最难的地方就是构造矩阵，而我根本一点都不会构造矩阵，没办法，从新学呗，

hdu-5667

这个题，我个人感觉也挺难的，因为矩阵好构造，这就是这个题的难点吧，

思路：发现a^b，和f[i-1]^c之类的东西，我们很明显吧这个幂变成乘，很自然的想到对数。问题是对什么取对数，最后发现对a取对数是合适的。

loga（fi）=loga（a^b*f[i-1]^c*f[i-2]）=loga(a^b)+loga(f[i-1]^c)+loga(f[i-2])，我们设k[i]=loga(fi)，所以k[i]=b+c*k[i-1]+k[i-2]。我们可以通过矩阵快速幂算出k[n]，然后a^k[n]=f[n]，可以直接用快速幂算出。

这里注意一下，由于k[n]非常大，所以为了使得a^k[n]%p==(a^(k[n]%y))%p，根据费马小定理如何gcd(a,p)=1，那么 a^(p-1)≡1（mod p），所以a^(p-1)%p=a^0%p，所以循环节为p-1，所以a^(k[n]%(p-1))%p。

这里还有注意一下如果a%p==0，a^(k[n]%(p-1))%p当k[n]=p-1的时候，a^(k[n]%(p-1))%p=1，但是实际上a^(k[n]%(p-1))%p=0，会造成错误，所以需要特判a%p==0的情况，不过这个题好像后台题目有点水，不必判断也能过貌似，所以代码里面没有体现

我们可以得出一个结论，那就是，a^(k[n]%(p-1))%p==(a^(k[n]%y))%p 前提p为质数，

废话不多说了下面直接给出代码吧，注意当n=1和n=2的时候要记得分类讨论，

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define met(Q,QQ) memset(Q,QQ,sizeof(Q))
ll n,a,b,c,p,mod;
struct matrix
{
    ll mat[5][5];
};
matrix cheng(matrix A,matrix B)
{
    matrix C;
    met(C.mat,0);
    for(int i=1;i<=3;i++)
        for(int j=1;j<=3;j++)
            for(int k=1;k<=3;k++)
                C.mat[i][j]=(C.mat[i][j]+A.mat[i][k]*B.mat[k][j])%mod;
    return C;

}
matrix P(matrix A,ll n)
{
    matrix C;
    met(C.mat,0);
    for(int i=1;i<=3;i++)
        C.mat[i][i]=1;
    while(n)
    {
        if(n&1) C=cheng(A,C);
        A=cheng(A,A);
        n=n/2;
    }
    return C;
}
ll PP(ll a,ll n)
{
    ll sum=1;
    a=a%p;
    while(n)
    {
        if(n&1) sum=sum*a%p;
        a=a*a%p;
        n=n/2;
    }
    return sum;
}

int main()
{
    ll t;
    scanf("%lld",&t);
    while(t--)
    {
        scanf("%lld %lld %lld %lld %lld",&n,&a,&b,&c,&p);
        mod=p-1;
        matrix A;
        met(A.mat,0);
        A.mat[1][1]=c;A.mat[1][2]=1;
        A.mat[1][3]=b;A.mat[2][1]=1;
        A.mat[3][3]=1;
        if(n==1)
        {
            printf("1\n");
            continue;
        }
        if(n==2)
        {
            ll sum=PP(a,b);
            printf("%lld\n",sum);
            continue ;
        }
        A=P(A,n-2);
        ll sum=A.mat[1][1]*b+A.mat[1][3];
        sum=PP(a,sum);
        printf("%lld\n",sum);


    }
    return 0;
}