Fibonacc 数列模n的循环节

最新推荐文章于 2021-04-22 16:56:48 发布

Lifel

最新推荐文章于 2021-04-22 16:56:48 发布

阅读量532

点赞数

本文链接：https://blog.csdn.net/Lifel/article/details/75948618

版权

本文介绍了如何计算斐波那契数列对任意模n的循环节，通过分解n为质因数的幂次，逐个求解每个质因数的循环节，并结合中国剩余定理确定最终的循环节。重点讨论了二次剩余的概念和如何使用矩阵乘法判断循环节长度。

摘要由CSDN通过智能技术生成

只谈解法,不说(hui)证明

步骤:
对于一个斐波拉契数列求模n的循环节:
1:任意一个n都可以分解为 priem[i]^time[i]
2:对于每一个prime[i]求出对它的循环节 g[i],则任意prime[i]^time[i]循环节的大小为:
k[i]= g[i]^(time[i]-1)
3:则对于模n的循环节为lcm(k[i]);

求解g[i]
二次剩余:对于一元二次同余方程有无的情况的讨论
判定对于a%p 来说
if(pow(a,(p-1)/2)mod p == 1 return 1//则a是p的二次剩余 else return -1//不是 —— 称为 legendre符号
ps:二次剩余只对奇素数成立,对2需要特殊判定

而对于prime[i]来说, 如果5是prime[i]的二次剩余则循环节为 priem[i]-1中的因子其中一个
否则为 (prime[i]+1)*2 中因子的一个用矩阵乘法判断即可:当 x==0&&y==1 时为循环节(x,y含义看代码)

Code:

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<vector>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fod(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
typedef long long ll;
const int N = 400005; 
struct Marix 
{
    int val[2][2];
};

Marix mutli(Marix a,Marix b,ll mod) 
{
    Marix c;
    for(int i=0;i<2;i+