long long pow(long long a,int k,long long mod)//快速幂取模
{
    long long b=1;
    while(k)
    {
        if(k&1) b=b*a%mod;
        a=a*a%mod;
        k>>=1;
    }
    return b;
}

struct Mat
{
    long long m[3][3];
    Mat(){memset(m,0,sizeof(m));}
    void init(int c,int d,int x)
    {
        memset(m,0,sizeof(m));
        m[0][0]=d;m[0][1]=c;m[0][2]=x;m[1][0]=m[2][2]=1;
    }
};
Mat multi(const Mat &a,const Mat &b)
{
    Mat c;
    for(int i=0;i<3;i++)
    for(int j=0;j<3;j++)
    for(int k=0;k<3;k++)
    c.m[i][j]=(c.m[i][j]+a.m[i][k]*b.m[k][j]%m)%m;
    return c;
}
Mat pow(Mat &a,int k)
{
    Mat b;
    for(int i=0;i<3;i++) b.m[i][i]=1;
    while(k)
    {
        if(k&1) b=multi(b,a);
        a=multi(a,a);
        k>>=1;
    }
    return b;
}