很容易看出来这道题是求模n意义下fib数列的最小循环节

对于fib数列的最小循环节的求法,我们可以这样:

1、令n=p1^m1 * p2^m2 * p3^m3……

2、分别计算fib数列在模p1^m1,p2^m2……意义下的最小循环节

3、模n意义下的最小循环节为2步骤各循环节的LCM

首先步骤三是非常容易证明的,CRT直接可以看出来

我们考虑如何计算p1^m1的循环节

这里有一个定理:G(p1^m1)=G(p1)*p1^(m-1)

我们只需要计算G(p)就可以了

我们知道对于Fib数列,满足f(i)=1.0/sqrt(5) *( (1+sqrt(5))/2)^i - ((1-sqrt(5))/2)^i)

那么在模p意义下存在以下两种情况:

1、5是模p意义下的二次剩余,即满足欧拉判别 5^((p-1)/2)=1

那么原式中的sqrt(5)在模意义下可以用整数x替代

由费马小定理很容易知道F(p-1)=0

又知道

((1+x)/2)^p=(1+x)/2

((1-x)/2)^p=(1-x)/2

两式做差等于x

所以F(p)=1

因此我们知道p-1是Fib数列的一个循环节,那么其最小循环节一定是p-1的约数

枚举+矩阵乘法判断即可

2、5不是模p意义下的二次剩余,则存在5^((p-1)/2)=-1

那么由二项式定理得

(1+sqrt(5))^p=1+(sqrt(5))^p

(1-sqrt(5))^p=1-(sqrt(5))^p

(因为其系数的组合数里含有p这个因子的时候值为0)

两式做差得2*sqrt(5)^p

进而得F(p)=(sqrt(5)/2)^(p-1)

我们知道在模p意义下2^(p-1)的逆元为1

所以F(p)=(sqrt(5))^(p-1)=5^((p-1)/2)=-1

同理,我们可以得到F(p+1)=0

设Fib数列在模意义下第一个F(i)=0且i不等于0的位置为i

设F(i-1)=a 易证(a,p)=1

则一定有i是p+1的约数

且F(p)=a^j*F(p-i*j)

显然F(p-i*j)=a,所以F(p)=a^(j+1)

进而知道F(2*p+1)=a^(2*j+2)=(-1)^2=1

F(2*p+2)=0

所以2*(p+1)是Fib数列的一个循环节,那么其最小循环节一定是2*(p+1)的约数

枚举+矩阵乘法判断即可

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
using namespace std; const int maxn=20010;
typedef long long LL;
LL n,mod,now;
LL a,b,ans;
LL st[maxn],top=0;
struct Matrix{
LL a[2][2];
void clear(){memset(a,0,sizeof(a));}
void print(){
for(int i=0;i<2;++i){
for(int j=0;j<2;++j){
printf("%lld ",a[i][j]);
}printf("\n");
}printf("\n");
}
}A;
LL mul(LL a,LL b){
LL s=0;
while(b){
if(b&1)s=(s+a)%mod;
a=(a<<1)%mod;b>>=1;
}return s;
}
LL GCD(LL a,LL b){return b==0?a:GCD(b,a%b);}
LL LCM(LL a,LL b){return a*b/GCD(a,b);}
LL Get_pow(LL v,LL p){
LL tmp=1;
while(p){
if(p&1)tmp=mul(tmp,v);
v=mul(v,v);p>>=1;
}return tmp;
}
Matrix operator *(const Matrix &A,const Matrix &B){
Matrix C;C.clear();
for(int i=0;i<2;++i){
for(int j=0;j<2;++j){
for(int k=0;k<2;++k){
C.a[i][j]=C.a[i][j]+mul(A.a[i][k],B.a[k][j]);
if(C.a[i][j]>=mod)C.a[i][j]-=mod;
}
}
}return C;
}
Matrix pow_mod(Matrix v,LL p){
Matrix tmp;tmp.clear();
tmp.a[0][0]=tmp.a[1][1]=1;
while(p){
if(p&1)tmp=tmp*v;
v=v*v;p>>=1;
}return tmp;
}
void F(LL p){
A.a[0][0]=1;A.a[0][1]=1;
A.a[1][0]=1;A.a[1][1]=0;
A=pow_mod(A,p-1);
b=(A.a[0][0]+A.a[0][1])%mod;
a=(A.a[1][0]+A.a[1][1])%mod;
}
LL Get_loop(LL k){
if(k==2)return 3;
else if(k==3)return 8;
else if(k==5)return 20;
LL p;mod=k;top=0;
if(Get_pow(5LL,(mod-1)>>1)==1)p=mod-1;
else p=2*(mod+1);
for(LL i=1;i*i<=p;++i){
if(p%i==0){
LL x=i,y=p/i;
F(x);
if(a==0&&b==1)return x;
if(y!=x)st[++top]=y;
}
}
while(top){
F(st[top]);
if(a==0&&b==1)return st[top];
top--;
}return 0;
} int main(){
while(scanf("%lld",&n)==1){
now=n;ans=1;
for(LL i=2;i*i<=now;++i){
if(now%i==0){
LL S=1;
while(now%i==0)now/=i,S*=i;
S/=i;S=S*Get_loop(i);
ans=LCM(ans,S);
}
}
if(now>1){
LL S=Get_loop(now);
ans=LCM(ans,S);
}
if(ans%2==0)ans/=2;
printf("%lld\n",ans);
}return 0;
}

  

hdu 4794 FIb求循环节的更多相关文章

  1. Cyclic Nacklace - HDU 3746(next求循环节)

    题目大意:给你一些串,问如果想让这个串里面的循环节至少循环两次,需要添加几个字符(只能在最前面或者最后面添加).比如ababc 需要添加5个就是添加ababc. 分析:其实字符串的长度len-next ...

  2. hdu 3746 kmp求循环节

    题意就是将所给的字符串变成多个完整的循环(至少两个),然后给出最少需要添加的字符数.

  3. hdu 4291(矩阵+暴力求循环节)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4291 思路:首先保留求出循环节,然后就是矩阵求幂了. #include<iostream> ...

  4. hdu 3374 String Problem (字符串最小最大表示 + KMP求循环节)

    Problem - 3374   KMP求循环节. http://www.cnblogs.com/wuyiqi/archive/2012/01/06/2314078.html   循环节推导的证明相当 ...

  5. hdu 4291 矩阵幂 循环节

    http://acm.hdu.edu.cn/showproblem.php?pid=4291 凡是取模的都有循环节-----常数有,矩阵也有,并且矩阵的更奇妙: g(g(g(n))) mod 109  ...

  6. 区间DP+next求循环节 uva 6876

    // 区间DP+next求循环节 uva 6876 // 题意:化简字符串 并表示出来 // 思路:dp[i][j]表示 i到j的最小长度 // 分成两部分 再求一个循环节 #include < ...

  7. HDU 4259 - Double Dealing(求循环节)

    首先将扑克牌进行一次置换,然后分解出所有的循环节,所有循环节的扑克牌个数的最小公倍数即为答案 #include <stdio.h> #include <string.h> #i ...

  8. 51nod 1126 求递推序列的第N项 思路:递推模拟,求循环节。详细注释

    题目: 看起来比较难,范围10^9 O(n)都过不了,但是仅仅是看起来.(虽然我WA了7次 TLE了3次,被自己蠢哭) 我们观察到 0 <= f[i] <= 6 就简单了,就像小学初中学的 ...

  9. UVA 12012 Detection of Extraterrestrial(KMP求循环节)

    题目描述 E.T. Inc. employs Maryanna as alien signal researcher. To identify possible alien signals and b ...

随机推荐

  1. Objective-C 【动态类型检测&响应方法】

    ------------------------------------------- 动态类型检测 代码: #import <Foundation/Foundation.h> @inte ...

  2. HTML5的Web SQL Database

    本文将介绍 Web SQL Database 规范中定义的三个核心方法: openDatabase:这个方法使用现有数据库或新建数据库来创建数据库对象 transaction:这个方法允许我们根据情况 ...

  3. .net设计模式之观察者模式

    摘要     在今天的设计模式系列中我给大家带来了观察者模式,首先我会以一个生动的故事引入观察者模式的应用的场景,然后描述这个场景中出现的问题,最后我们提出观察者模式的解决方案,并给出C#语言实现的代 ...

  4. 基于MRG_MyISAM引擎的Mysql分表

    正常情况下的分表,都是直接创建多个相同结构的表,比如table_1.table_2...最近碰到一个特殊需求,需要创建一个主表,所有分表的数据增删改查,全部自动实时更新到主表,这个时候可以使用MRG_ ...

  5. 3.IP转发

    1. "vim  /usr/lib/sysctl.d/00-system.conf"在#Disable netfilter on bridges.栏下面添加行:"net. ...

  6. DTCMS展示一级栏目并展示各自栏目下的二级栏目

    c#代码中 <!--C#代码--> <%csharp%> string parent_id=DTRequest.GetQueryString("parent_id&q ...

  7. 在Tomcat中部署Java Web应用程序

    在Tomcat中部署Java Web应用程序有两种方式:静态部署和动态部署.在下文中$CATALINA_HOME指的是Tomcat根目录. 一.静态部署      静态部署指的是我们在服务器启动之前部 ...

  8. [分享]好用的Markdown编辑器

  9. sublime 修改Xdebug插件快捷键

    最近在用Xdebug插件 感觉挺好用 但是快捷键不太舒服,特别是调试下一步的时候,比较麻烦,按键较多: 所以想DIY下 但是preferences->package setting ->X ...

  10. php取整的几种方法

    php取整的几种方式. floor 舍去法取整 语法格式:float floor ( float value )返回不大于value 的下一个整数,将value 的小数部分舍去取整.floor() 返 ...