二元数列 - 矩阵快速幂

题目地址：http://www.51cpc.com/web/problem.php?id=4267

 

Knowledge Point：https://www.cnblogs.com/liubilan/p/9450568.html

 

Summarize：

　　1. 矩阵乘法 a*b 与 b*a 不同，注意顺序；

　　2. 初始其中一个是乘数应该为对角矩阵，且对角线上元素全部为1，作用类似于整数1；

　　3. 关于函数返回值，由于试了很久我也没能直接返回二维数组，所以在这里用结构体封装二维数组，函数返回结构体；

 

附代码：

#include<iostream>
#include<cstring>
using namespace std;

#define LL long long
#define MOD 20181225
LL x0, y0, a, b, c, d, n;
struct Matrix {
    LL w[2][2];
}matrix;

Matrix calculate(Matrix x, Matrix y)
{    
    Matrix tmp;
    memset(tmp.w, 0, sizeof(tmp.w));
    for(int i=0; i<2; i++)
        for(int j=0; j<2; j++) {
            for(int k=0; k<2; k++)
                tmp.w[i][j] += x.w[i][k]*y.w[k][j];
            tmp.w[i][j]%=MOD;
        }
    return tmp;
}

void Quick_Pow()
{
    Matrix tmp;
    tmp.w[0][0]=tmp.w[1][1]=1;
    tmp.w[0][1]=tmp.w[1][0]=0;
    
    while(n) {
        if(n&1) {
            tmp = calculate(matrix, tmp);
        }
        matrix = calculate(matrix, matrix);
        n>>=1;
    }
    
    LL xn = (tmp.w[0][0]*x0%MOD + tmp.w[0][1]*y0%MOD)%MOD;
    LL yn = (tmp.w[1][0]*x0%MOD + tmp.w[1][1]*y0%MOD)%MOD;
    cout<<xn<<' '<<yn<<endl;
}

int main()
{
    ios::sync_with_stdio(false);
    while(cin>>x0>>y0>>a>>b>>c>>d>>n) {
        matrix.w[0][0]=a, matrix.w[0][1]=b;
        matrix.w[1][0]=c, matrix.w[1][1]=d;
        Quick_Pow();
    }
    return 0;
}