斐波那契数列满足$f(n-1)f(n+1)-f(n)^2=(-1)^n$。
枚举$-1$的符号,根据二次剩余即可求出最多$4$个可能的$f(n+1)$的值。
那么根据$f(n)$和$f(n+1)$,对矩阵做BSGS求出最小的$n$即可。
时间复杂度$O(sqrt{P}log P)$。
#include<cstdio> #include<algorithm> #include<map> using namespace std; typedef long long ll; typedef pair<ll,ll>PI; const int P=1000000009,K=32000; ll X,ans=-1,g[K+5][2][2],w[K+5][2];map<PI,int>T; inline ll po(ll a,ll b,ll P){ll t=1;for(;b;b>>=1,a=a*a%P)if(b&1)t=t*a%P;return t;} ll ToneLLi_Shanks(ll n,ll p){ if(!n)return 0; if(p==2)return n&1?1:-1; if(po(n,p>>1,p)!=1)return -1; if(p&2)return po(n,p+1>>2,p); int s=__builtin_ctzll(p^1); ll q=p>>s,z=2; for(;po(z,p>>1,p)==1;z++); ll c=po(z,q,p); ll r=po(n,q+1>>1,p); ll t=po(n,q,p),tmp; for(int m=s,i;t!=1;){ for(i=0,tmp=t;tmp!=1;i++)tmp=tmp*tmp%p; for(;i<--m;)c=c*c%p; r=r*c%p; c=c*c%p; t=t*c%p; } return r; } inline void mul(ll a[][2],ll b[][2],ll f[][2]){ static ll c[2][2]; c[0][0]=(a[0][0]*b[0][0]+a[0][1]*b[1][0])%P; c[0][1]=(a[0][0]*b[0][1]+a[0][1]*b[1][1])%P; c[1][0]=(a[1][0]*b[0][0]+a[1][1]*b[1][0])%P; c[1][1]=(a[1][0]*b[0][1]+a[1][1]*b[1][1])%P; f[0][0]=c[0][0]; f[0][1]=c[0][1]; f[1][0]=c[1][0]; f[1][1]=c[1][1]; } void work(ll X,ll Y){ T.clear(); for(int i=0;i<K;i++){ ll a=(g[i][0][0]*X+g[i][0][1]*Y)%P; ll b=(g[i][1][0]*X+g[i][1][1]*Y)%P; T[PI(a,b)]=i; } for(int i=1;i<=K;i++){ map<PI,int>::iterator it=T.find(PI(w[i][0],w[i][1])); if(it!=T.end()){ ll now=i*K-it->second; if(ans==-1||ans>now)ans=now; return; } } } void solve(ll A,ll B,ll C){ ll d=((B*B-4*A*C)%P+P)%P; d=ToneLLi_Shanks(d,P); if(d<0)return; d=(d%P+P)%P; ll t=po(A*2%P,P-2,P); work(X,((-B-d)*t%P+P)%P); work(X,((-B+d)*t%P+P)%P); } int main(){ scanf("%lld",&X); if(X>=P)return puts("-1"),0; ll A=P-1,B=X,C=X*X%P; g[0][0][0]=1; g[0][1][1]=1; g[1][0][1]=1; g[1][1][0]=1; g[1][1][1]=1; for(int i=2;i<=K;i++)mul(g[i-1],g[1],g[i]); w[0][1]=1; for(int i=1;i<=K;i++){ w[i][0]=(g[K][0][0]*w[i-1][0]+g[K][0][1]*w[i-1][1])%P; w[i][1]=(g[K][1][0]*w[i-1][0]+g[K][1][1]*w[i-1][1])%P; } solve(A,B,(C-1+P)%P); solve(A,B,(C+1)%P); return printf("%lld",ans),0; }