hdu 5451 Best Solver (特征方程求通项+广义Fibonacci数列找循环节)

题目大意： y=(5+26√)1+2x 。给出x和M，求这个式子的整数部分对M取模的值

恶补了一下线性代数。。。

记λ1=5+2√6，λ2=5-2√6。则λ1λ2=1,λ1+λ2=10

则λ1、λ2可以看作特征方程x^2-10x+1=0的两根。

该特征方程对应的数列递推公式为a[n]=10*a[n-1]-a[n-2]。

由定理3，得到通项a[n]=(5+2√6)^n+(5-2√6)^n

观察通项，发现(5+2√6)^n>1，(5-2√6)^n<1，因此，(5+2√6)^n的整数部分的值即为a[n] - 1。

到这里，可以利用递推公式a[n]=10*a[n-1]-a[n-2]，构造矩阵来找循环节，相关方法：点击打开链接

对于该题，发现n=2^x指数非常大，而模数M<=46337，打表发现对于每一个模数，数列a[n]都具有循环节，并且都比较小。

因此，也可以直接暴力找一下循环节。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<map>
#include<set>
#include<queue>
#include<stack>
using namespace std;
typedef double db;
typedef __int64 LL;
#define mod 1000000007
#define maxn 50005
LL n;
int r[maxn],a[maxn],M;
int POW(int a,LL b,int p)
{
    int ans=1;
    while(b)
    {
        if(b&1) ans=ans*a%p;
        a=a*a%p;
        b>>=1;
    }
    return ans;
}
int main()
{
    LL T;
    int ca=1,i;
    cin>>T;
    while(T--)
    {
        scanf("%I64d%d",&n,&M);
        printf("Case #%d: ",ca++);
        a[0]=2%M,a[1]=10%M;
        if(!r[M]){
            for(i=2;;++i)
            {
                a[i]=(10*a[i-1]%M-a[i-2]+M)%M;
                if(a[i-1]==a[0]&&a[i]==a[1]) {r[M]=i-1;break;}
            }
        }
        int p=(POW(2,n,r[M])+1)%r[M];
        for(i=2;i<=p;++i) a[i]=(10*a[i-1]%M-a[i-2]+M)%M;
        printf("%d\n",(a[p]-1+M)%M);
    }
    return 0;
}