CodeForces - 385E Bear in the Field(矩阵快速幂)

题目描述：

点击打开链接

题意给你六个数n,sx,sy,dx,dy,t。n表示田地的大小,sx,sy是熊的初始位置，dx,dy为熊每秒的速度增加数,t为熊走的时间。接着还有条件是这片田上长着berry初始数量为x+y为每个格子的两坐标值和，并且没过一秒会长一个berry。熊每到一个格子上会吃掉这个格子上的所有berry并且速度上升这个值。然后要求的是t秒之后熊的位置。这题很明显是个递推的。Xn=Xn-1+dxn-1+Xn-1+Yn-1+Tn-1,Yn=Yn-1+dyn-1+Xn-1+Yn-1+Tn-1,dxn=Xn-1+Yn-1+dxn-1+Tn-1,dyn=Xn-1+Yn-1+dyn-1+Tn-1,Tn=Tn-1+1；这是最基本递推式，但是为了方便处理，我们最好把坐标系从1-n移到0-n-1。因为这道题涉及到对n取模的过程，如果用1-n的坐标系的话，那么一旦取模的话原本处于n的位置可能会被取为0这样答案就出错了。所以移动坐标系，那么初始位置sx-1,sy-1然后在递推式中除时间外各项+2，最后算完答案后都+1就好了。算答案的过程就直接矩阵快速幂就好。

AC代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<string>
#include<stack>
#include<set>
#include<queue>
#include<algorithm>
using namespace std;
struct MAT
{
    long long s[6][6];
    MAT () {
        memset(s,0,sizeof(s));
    }
};

long long n,sx,sy,dx,dy,t;

MAT MAT_MUL(MAT m1,MAT m2)
{
    MAT m3;
    for (int i=0;i<6;i++) {
        for (int j=0;j<6;j++) {
            for (int k=0;k<6;k++)
            m3.s[i][j]=(m3.s[i][j]+m1.s[i][k]*m2.s[k][j]%n+n)%n;
        }
    }
    return m3;
}
int main()
{
    scanf("%lld%lld%lld%lld%lld%lld",&n,&sx,&sy,&dx,&dy,&t);
    MAT m1;
    m1.s[0][0]=2;m1.s[0][1]=1;m1.s[0][2]=1;m1.s[0][3]=0;m1.s[0][4]=1;m1.s[0][5]=2;
    m1.s[1][0]=1;m1.s[1][1]=2;m1.s[1][2]=0;m1.s[1][3]=1;m1.s[1][4]=1;m1.s[1][5]=2;
    m1.s[2][0]=1;m1.s[2][1]=1;m1.s[2][2]=1;m1.s[2][3]=0;m1.s[2][4]=1;m1.s[2][5]=2;
    m1.s[3][0]=1;m1.s[3][1]=1;m1.s[3][2]=0;m1.s[3][3]=1;m1.s[3][4]=1;m1.s[3][5]=2;
    m1.s[4][0]=0;m1.s[4][1]=0;m1.s[4][2]=0;m1.s[4][3]=0;m1.s[4][4]=1;m1.s[4][5]=1;
    m1.s[5][0]=0;m1.s[5][1]=0;m1.s[5][2]=0;m1.s[5][3]=0;m1.s[5][4]=0;m1.s[5][5]=1;
    MAT m2;
    m2.s[0][0]=sx-1;m2.s[1][0]=sy-1;m2.s[2][0]=dx;m2.s[3][0]=dy;m2.s[4][0]=0;m2.s[5][0]=1;

    while(t) {
        if (t&1) m2=MAT_MUL(m1,m2);
        m1=MAT_MUL(m1,m1);
        t=t>>1;
    }

    long long ansx=m2.s[0][0];
    long long ansy=m2.s[1][0];
    printf("%lld %lld\n",ansx+1,ansy+1);
    return 0;
}