「Codechef FN」Fibonacci Number

传送门

---

problem

设 $F_i$ 表示斐波那契数列：

$$\{\begin{array}{ll}{F}_n=n,& n\le 1\\ F_n=F_{n-1}+F_{n-2},& \mathrm{o}\mathrm{t}\mathrm{h}\mathrm{e}\mathrm{r}\mathrm{w}\mathrm{i}\mathrm{s}\mathrm{e}\end{array}$$

给定 $c$ 和 $p$，求满足 $F_n\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}p)$ 的最小的 $n$。

保证 $11\le p\le 2\times 10^9$，$p$ 是素数，且 $p\mathrm{m}\mathrm{o}\mathrm{d}10$ 是一个完全平方数。

---

solution

首先，斐波那契数列是有通项公式的：

$$F_n=\frac{1}{\sqrt{5}}\left((\frac{1+\sqrt{5}}{2}{)}^n-(\frac{1-\sqrt{5}}{2}{)}^n\right)$$

那么设 $\phi =\frac{1+\sqrt{5}}{2}$，那么 $-\frac{1}{\phi }=\frac{1-\sqrt{5}}{2}$，那么相当于是解：

$$\begin{array}{rl}\frac{1}{\sqrt{5}}({\phi }^n-(-\frac{1}{\phi }{)}^n)& \equiv c(\mathrm{m}\mathrm{o}\mathrm{d}p)\\ {\phi }^n-\frac{(-1{)}^n}{{\phi }^n}& \equiv \sqrt{5}c(\mathrm{m}\mathrm{o}\mathrm{d}p)\end{array}$$

又令 $\varphi ={\phi }^n$，继续化简，有：

$$\varphi -\frac{(-1{)}^n}{\varphi }\equiv \sqrt{5}c(\mathrm{m}\mathrm{o}\mathrm{d}p)$$

然后用个求根公式就可以解出 $\varphi$：

$$\varphi \equiv \frac{\sqrt{5}c\pm \sqrt{5c^2+4(-1{)}^n}}{2}(\mathrm{m}\mathrm{o}\mathrm{d}p)$$

那么用二次剩余就可以解 $\varphi$ 了，又由于 $\varphi ={\phi }^n$，再套个 BSGS 解 $n$ 即可。注意解得时候要特判 $n$ 的奇偶性来解。

然后 $5$ 是模 $p$ 意义下的二次剩余，直接解就可以了。

由于 $p$ 是素数，且 $p\mathrm{m}\mathrm{o}\mathrm{d}10$ 是完全平方数，所以 $p=\pm 1(\mathrm{m}\mathrm{o}\mathrm{d}5)$，那么 $p^{\frac{5-1}{2}}=1$，即 $p$ 是模 $5$ 意义下的二次剩余。
根据二次互反律，设 $p,q$ 是不同的奇素数，有：

$$\left(\frac{p}{q}\right)\left(\frac{q}{p}\right)=(-1{)}^{\frac{p-1}{2}\frac{q-1}{2}}$$

其中 $\left(\frac{p}{q}\right)$ 是勒让德符号。
那么把 $5$ 和 $p$ 代进上式，就不难发现 $5$ 是模 $p$ 意义下的二次剩余。

---

code

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
int T,C,P,s5,U,V;
int add(int x,int y)  {return (ll)x+y>=P?(ll)x+y-P:x+y;}
int dec(int x,int y)  {return (ll)x-y< 0?(ll)x-y+P:x-y;}
int mul(int x,int y)  {return (ll)x*y>=P?(ll)x*y%P:x*y;}
int power(int a,int b){
	int ans=1;
	for(;b;b>>=1,a=mul(a,a))  if(b&1)  ans=mul(ans,a);
	return ans;
}
int Inv(int x)  {return power(x,P-2);}
void Max(int &x,int y)  {if(x<y)x=y;}
void Min(int &x,int y)  {if(x>y)x=y;}
namespace Cipolla{
	int I,val;
	struct num{
		int x,y;
		num(int x=0,int y=0):x(x),y(y){}
		friend num operator*(const num &a,const num &b){
			return num(add(mul(a.x,b.x),mul(val,mul(a.y,b.y))),add(mul(a.y,b.x),mul(a.x,b.y)));
		}
		friend num operator^(num a,int b){
			num ans(1,0);
			for(;b;b>>=1,a=a*a)  if(b&1)  ans=ans*a;
			return ans;
		}
	};
	int Sqrt(int n){
		if(!n)  return 0;
		if(power(n,(P-1)>>1)==P-1)  return -1;
		while(1){
			I=rand(),val=dec(mul(I,I),n);
			if(power(val,(P-1)>>1)==P-1)  break;
		}
		return (num(I,1)^((P+1)>>1)).x;
	}
}
using Cipolla::Sqrt;
unordered_map<int,int>Hash;
int BSGS(int a,int b,int type){
	Hash.clear();
	int now=b,t=ceil(sqrt(P));
	for(int i=0;i<t;++i)  Hash[now]=i,now=mul(now,a);
	now=1,a=power(a,t);
	if(!a)  return b?P+1:1;
	for(int i=0,pos;i<=t;++i){
		pos=(Hash.find(now)==Hash.end())?-1:Hash[now],now=mul(now,a);
		if(pos>=0&&i*t-pos>=0&&((i*t-pos)&1)==type)  return i*t-pos;
	}
	return P+1;
}
int main(){
	scanf("%d",&T);
	while(T--){
		scanf("%d%d",&C,&P);
		s5=Sqrt(5),U=mul(add(1,s5),Inv(2)),V=mul(C,s5);
		int d1=add(mul(V,V),4),d2=dec(mul(V,V),4),ans=P+1;
		if(!d1||power(d1,(P-1)>>1)==1){
			int tmp=Sqrt(d1),val;
			val=mul(add(V,tmp),Inv(2)),Min(ans,BSGS(U,val,0));
			val=mul(dec(V,tmp),Inv(2)),Min(ans,BSGS(U,val,0));
		}
		if(!d2||power(d2,(P-1)>>1)==1){
			int tmp=Sqrt(d2),val;
			val=mul(add(V,tmp),Inv(2)),Min(ans,BSGS(U,val,1));
			val=mul(dec(V,tmp),Inv(2)),Min(ans,BSGS(U,val,1));
		}
		printf("%d\n",ans==P+1?-1:ans);
	}
	return 0;
}