2019牛客多校第五场 B generator 1

传送门：https://ac.nowcoder.com/acm/contest/885/B

据说这题byf在他自己讲课的时候讲过这个方法，然后其他队都做出来了，byf在巨大压力下脑补出了这个10次幂的矩阵快速幂，最后两分钟过了。然后一堆水题没看

时限卡得比较死，我最后一起写，每次自己乘10次，超时了，还是得用2次幂的快速幂算这10次才能过

#include
using namespace std;
typedef long long LL;
const int size=1e6+5;
char n[size];
int x0,x1,a,b;
typedef long long LL;
LL mod;
struct mat
{
    int m[3][3];
    mat(){for(int i=0;i<3;i++)for(int j=0;j<3;j++) m[i][j]=0;}
    friend mat operator*(mat x,mat y)
    {
        mat ans;
        for(int i=1;i<=2;i++)
        {
            for(int j=1;j<=2;j++)
            {
                for(int k=1;k<=2;k++)
                {
                    ans.m[i][k] = (ans.m[i][k]+1LL*x.m[i][j]*y.m[j][k])%mod;
                }
            }
        }
        return ans;
    }
};
inline mat quick_pow(mat a,int b)
{
    mat ans;
    ans.m[1][1]=1;ans.m[2][2]=1;
    while(b)
    {
        if(b&1) ans=ans*a;
        a=a*a;
        b>>=1;
    }
    return ans;
}
int main()
{
    scanf("%d%d%d%d",&x0,&x1,&a,&b);
    scanf("%s%lld",n,&mod);
    mat ori;
    int len =strlen(n);
    ori.m[1][2]=1,ori.m[2][1]=b;ori.m[2][2]=a;
    mat ans;
    ans.m[1][1]=1;ans.m[2][2]=1;
    for(int i=len-1;i>=0;i--)
    {
        ans=ans*quick_pow(ori,n[i]-'0');
        ori=quick_pow(ori,10);
    }
    printf("%lld\n",(1LL*x0*ans.m[1][1]+1LL*x1*ans.m[1][2])%mod);
}