Luogu P3263 [JLOI2015]有意义的字符串

Link
设(e=frac{b+sqrt d}2,i=frac{b-sqrt d}2)。
显然(f_n=e^n+i^n)是一个整数，且(f_n=(e+i)f_{n-1}+eif_{n-2})。
递推式中的(e+i=b,ei=frac{b^2-d}4)，根据题目条件这两个也是整数。
因此我们可以利用矩阵快速幂求出(f_nmod P)。
显然(iin(-1,0])，因此(i^nin(-1,0])，因此(lfloor e^n floorequiv f_n-lceil |i^n| ceilpmod P)，直接判断即可。

#include<cstdio>
#include<cstring>
using u64=unsigned long long;
const u64 P=7528443412579576937ll;
u64 inc(u64 a,u64 b){return a+=b,a>=P? a-P:a;}
u64 mul(u64 a,u64 b){u64 r=0;for(;b;b>>=1,a=inc(a,a))if(b&1)r=inc(a,r);return r;}
struct matrix{u64 a[3][3];matrix(){memset(a,0,72);}u64*operator[](int x){return a[x];}}E,I;
matrix operator*(matrix a,matrix b)
{
    matrix c;
    for(int i=1;i<=2;++i) for(int j=1;j<=2;++j) for(int k=1;k<=2;++k) c[i][j]=inc(c[i][j],mul(a[i][k],b[k][j]));
    return c;
}
int main()
{
    u64 b,d,n;int f;
    scanf("%llu%llu%llu",&b,&d,&n);
    if(!n) return !printf("1");
    f=--n&1,E[1][2]=(d-b*b)/4,E[2][1]=1,E[2][2]=b,I[1][1]=I[2][2]=1;
    for(;n;n>>=1,E=E*E) if(n&1) I=E*I;
    printf("%llu",inc(mul(2,I[1][2]),mul(b,I[2][2]))-(d*d^b&&f));
}