HDU4565

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4565

题目大意：

　　给出a，b，n，m，求出下式中的 Sn.

                                                    

解题思路：

　　

AC代码：

#include <cstdio>
#include <cstring>

using namespace std;
typedef long long ll;
struct Matrix{
    ll mat[2][2];
};
Matrix Multiply(Matrix x,Matrix y,ll mod){
    Matrix temp;
    memset(temp.mat,0,sizeof(temp.mat));
    for(int i=0;i<2;i++){
        for(int j=0;j<2;j++){
            for(int k=0;k<2;k++){
                temp.mat[i][j]+=(x.mat[i][k]*y.mat[k][j]);
                temp.mat[i][j]%=mod;
            }
        }
    }
    return temp;
}
Matrix Fast_Power(Matrix a,ll m,ll mod){
    Matrix res;
    memset(res.mat,0,sizeof(res.mat));
    for(int i=0;i<2;i++)    res.mat[i][i]=1;
    while(m){
        if(m&1) res=Multiply(res,a,mod);
        m>>=1;
        a=Multiply(a,a,mod);
    }
    return res;
}
int main()
{
    ll a,b,n,m;
    while(scanf("%lld%lld%lld%lld",&a, &b, &n, &m)==4){
        Matrix l={  a,1,
                    0,0 };
        Matrix r={  a,1,
                    b,a };
        r=Fast_Power(r,n-1,m);
        l=Multiply(l,r,m);
        printf("%lld
",(l.mat[0][0]*2)%m);
    }
    return 0;
}