题解 P5110 【块速递推】

一道不错的题，教会我好多新套路。

首先得到递推式的特征方程

$x^2-233x-666=0$

易得$x=\frac{233\pm 13\sqrt{337}}{2}$

设通项式$A_n=a{x}_1^n+b{x}_2^n$

因为$A_0=0,A_1=1$

所以

$$\{\begin{array}{l}\begin{array}{c}a+b=0\\ a{x}_1+b{x}_2=1\end{array}\end{array}$$

解得

$$\{\begin{array}{l}\begin{array}{c}a=\frac{1}{13\sqrt{337}}\\ b=-\frac{1}{13\sqrt{337}}\end{array}\end{array}$$

于是得通项式$A_n=\frac{1}{13\sqrt{337}}((\frac{233+13\sqrt{337}}{2}{)}^n-(\frac{233-13\sqrt{337}}{2}{)}^n)$

用$BSGS$求二次剩余，不难解得在$\mathrm{m}\mathrm{o}\mathrm{d}1000000007$ 时，通项式为

$A_n\equiv 233230706(94153035^n-905847205^n)(\mathrm{m}\mathrm{o}\mathrm{d}1000000007)$

考虑到快速幂带$log$，无法通过本题。。。

然而因为模数为奇质数

于是根据费马小定理，显然任意数的幂都有长度为$1000000006$的循环节，因此指数被缩小到$1000000006$的范围

又因为底数与模数确定

所以我们可以选取一个与$\sqrt{1000000006}$同阶的数$m$，预处理出$x^0,x^1,...,x^{m-1}$和$x^0,x^m,...,x^{m^2}$

然后对于每一个询问$x^k$，$O(1)$求$x^{\lfloor k/m\rfloor }\times x^{k\mathrm{m}\mathrm{o}\mathrm{d}m}$即可

于是此题得解。。。

#include<cstdio>
const int P=1000000007;
const int a=233230706,b1=94153035,b2=905847205;
const int BLOCK_SIZE=32768;

int hbp1[BLOCK_SIZE],lbp1[BLOCK_SIZE];
int hbp2[BLOCK_SIZE],lbp2[BLOCK_SIZE];

void init()
{
    lbp1[0]=1;
    for(int i=1;i<BLOCK_SIZE;++i)lbp1[i]=(long long)lbp1[i-1]*b1%P;
    hbp1[0]=1;
    hbp1[1]=(long long)lbp1[BLOCK_SIZE-1]*b1%P;
    for(int i=2;i<BLOCK_SIZE;++i)hbp1[i]=(long long)hbp1[i-1]*hbp1[1]%P;
    lbp2[0]=1;
    for(int i=1;i<BLOCK_SIZE;++i)lbp2[i]=(long long)lbp2[i-1]*b2%P;
    hbp2[0]=1;
    hbp2[1]=(long long)lbp2[BLOCK_SIZE-1]*b2%P;
    for(int i=2;i<BLOCK_SIZE;++i)hbp2[i]=(long long)hbp2[i-1]*hbp2[1]%P;
}
int power1(int k)
{
    return (long long)hbp1[k/BLOCK_SIZE]*lbp1[k%BLOCK_SIZE]%P;
}
int power2(int k)
{
    return (long long)hbp2[k/BLOCK_SIZE]*lbp2[k%BLOCK_SIZE]%P;
}

namespace Mker
{
    unsigned long long SA,SB,SC;
    void init(){scanf("%llu%llu%llu",&SA,&SB,&SC);}
    unsigned long long rand()
    {
        SA^=SA<<32,SA^=SA>>13,SA^=SA<<1;
        unsigned long long t=SA;
        SA=SB,SB=SC,SC^=t^SA;return SC;
    }
}

int main()
{
    int T;
    scanf("%d",&T);
    init();
    Mker::init();
    int res=0;
    while(T--)
    {
        int n=Mker::rand()%1000000006;
        res^=(long long)a*(power1(n)-power2(n)+P)%P;
    }
    printf("%d",res);
}