hdoj1588(2)

void getSum(){
    n--;
    int cnt=0;
    bool dig[10];
    while(n){
        dig[cnt++]=n&1;
        n>>=1;
    }
    //S=0;C=E
    memset(S,0,sizeof(S));
    memset(C,0,sizeof(C));
    C[1][1]=C[2][2]=1;
    while(cnt--){
        //T=S*(C+E)
        D[1][1]=C[1][1]+1;D[1][2]=C[1][2];//D=C+E
        D[2][1]=C[2][1];D[2][2]=C[2][2]+1;
        memset(T,0,sizeof(T));
        for(int i=1;i<=2;i++)//T=S*D
          for(int j=1;j<=2;j++)
            for(int k=1;k<=2;k++){
                T[i][j]+=S[i][k]*D[k][j];
                T[i][j]%=M;
            }
        //D=C^2
        memset(D,0,sizeof(D));
        for(int i=1;i<=2;i++)
          for(int j=1;j<=2;j++)
            for(int k=1;k<=2;k++){
                D[i][j]+=C[i][k]*C[k][j];
                D[i][j]%=M;
            }
        if(dig[cnt]){
            //C=D*A=C^2*A
            memset(C,0,sizeof(C));
            for(int i=1;i<=2;i++)
              for(int j=1;j<=2;j++)
                for(int k=1;k<=2;k++){
                    C[i][j]+=D[i][k]*K[k][j];
                    C[i][j]%=M;
                }
            //S=T+C
            S[1][1]=(T[1][1]+C[1][1])%M;
            S[1][2]=(T[1][2]+C[1][2])%M;
            S[2][1]=(T[2][1]+C[2][1])%M;
            S[2][2]=(T[2][2]+C[2][2])%M;
        }
        else{
            //C=D;
            C[1][1]=D[1][1];C[1][2]=D[1][2];
            C[2][1]=D[2][1];C[2][2]=D[2][2];
            //S=T;
            S[1][1]=T[1][1];S[1][2]=T[1][2];
            S[2][1]=T[2][1];S[2][2]=T[2][2];
        }
    }
    S[1][1]++;
    S[2][2]++;
}

void getAns(){
    //T=S*B
    memset(T,0,sizeof(T));
    for(int i=1;i<=2;i++)
      for(int j=1;j<=2;j++)
        for(int k=1;k<=2;k++){
            T[i][j]+=S[i][k]*B[k][j];
            T[i][j]%=M;
        }
    //ans
    ans=T[1][2];
}
int main(){
    //freopen("input","r",stdin);
    A[1][1]=0;A[1][2]=A[2][1]=A[2][2]=1;
    while(scanf("%d%d%d%d",&k,&b,&n,&M)!=EOF){
        //求B=A^b
        getB();
        //求K=A^k
        getK();
        //求K^0+K^1+..+K^(n-1)
        getSum();
        //求ans
        getAns();
        printf("%lld\n",ans);
    }
}