2016湖大校赛 L题 The Sequence likes Ladder

#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
const int maxn=10010;
int euler[maxn];
ll a,k,m,n;

int phi()
{
    for(int i=0;i<maxn;i++) euler[i]=i;
    for(int i=2;i<maxn;i++)
    {
        if(euler[i]==i)
        {
            for(int j=i;j<maxn;j+=i)
            {
                euler[j]=euler[j]/i*(i-1);
            }
        }
    }
}

struct Matrix
{
    ll a[2][2];
    Matrix(){memset(a,0,sizeof(a));}
    Matrix operator* (const Matrix &p)
    {
        Matrix res;
        for(int i=0;i<2;i++)
        {
            for(int j=0;j<2;j++)
            {
                for(int k=0;k<2;k++)
                {
                    res.a[i][j]+=(a[i][k]*p.a[k][j]%euler[m]);
                }
                res.a[i][j]%=euler[m];
            }
        }
        return res;
    }
}ans,base;

Matrix quick_pow(Matrix base,ll n)
{
    Matrix res;
    for(int i=0;i<2;i++)
    {
        res.a[i][i]=1;
    }
    while(n)
    {
        if(n&1) res=res*base;
        base=base*base;
        n>>=1;
    }
    return res;
}

ll pow(ll a,ll n)
{
    ll ans=1;
    while(n)
    {
        if(n&1) ans=ans*a%m;
        a=a*a%m;
        n>>=1;
    }
    return ans;
}

void Matrix_init()
{
    ans.a[0][0]=1;
    ans.a[0][1]=1;
    ans.a[1][0]=0;
    ans.a[1][1]=0;
    base.a[0][0]=k;
    base.a[0][1]=0;
    base.a[1][0]=1;
    base.a[1][1]=1;
}

int main()
{
    ll x;
    phi();
    while(~scanf("%lld%lld%lld%lld",&a,&k,&m,&n))
    {
        Matrix_init();
        ans=ans*quick_pow(base,n-1);
        x=ans.a[0][0];
        x=pow(a%m,x)%m;
        printf("%lld\n",x);
    }
    return 0;
}