题意
f(n)表示斐波那契数列的第n项,f(0)=0,f(1)=1,现在给出一个n,求f(f(n)) mod 20160519
n<=10^100(原题是10^9,但我的程序可以跑过10^100)
分析
一开始想的是费马小定理,后来发现是不行的,于是就去问了一发,发现了一个很脑洞的想法:找寻环节。也就是找到在mod MOD的情况下最小的k使得f(k) mod MOD=0,f(k+1) mod MOD=1且k!=0,那么我们要求的就是f(f(n) mod k),同理只要我们再找一波在mod k的情况下的循环节就好了。
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define MOD1 20160519
#define MOD2 26880696
#define MOD3 746688
#define ll long long
using namespace std;
struct arr
{
int a[3][3];
}w,c;
char str[205];
int divi(int p)
{
ll s=0;int len=strlen(str);
for (int i=0;i<len;i++)
{
s=s*10+str[i]-'0';
s=s%p;
}
return s;
}
void clear()
{
w.a[1][1]=0;w.a[1][2]=w.a[2][1]=w.a[2][2]=1;
c.a[1][1]=0;c.a[1][2]=c.a[2][1]=c.a[2][2]=1;
}
void mul(arr x,arr y,int p)
{
memset(c.a,0,sizeof(c.a));
for (int i=1;i<=2;i++)
for (int j=1;j<=2;j++)
for (int k=1;k<=2;k++)
c.a[i][j]=(c.a[i][j]+(ll)x.a[i][k]*y.a[k][j]%p)%p;
}
void ksm(int x,int y)
{
if (x==1) return;
ksm(x/2,y);
mul(c,c,y);
if (x%2==1) mul(c,w,y);
}
int main()
{
//freopen("na.in","r",stdin);
//freopen("na.out","w",stdout);
int t;
scanf("%d",&t);
while (t--)
{
clear();
scanf("%s",str);
int x=divi(MOD3);
if (!x)
{
printf("0\n");
continue;
}
if (x==1)
{
printf("1\n");
continue;
}
ksm(x,MOD2);
x=c.a[2][1];
if (!x)
{
printf("0\n");
continue;
}
if (x==1)
{
printf("1\n");
continue;
}
clear();
ksm(x,MOD1);
printf("%d\n",c.a[2][1]);
}
return 0;
}