BNUOJ 52322 Fibonacci of Fibonacci 矩阵乘法+高精度+预处理

题意

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;
}