这类型的题目其实没什么意思..知道怎么做后,就有固定套路了..而且感觉这东西要出的很难的话,有这种方法解常数会比较大吧..所以一般最多套一些比较简单的直接可以暴力求循环节的题目了..

/** @Date    : 2017-09-26 16:37:05
* @FileName: HDU 3977 斐波那契循环节.cpp
* @Platform: Windows
* @Author : Lweleth (SoungEarlf@gmail.com)
* @Link : https://github.com/
* @Version : $Id$
*/
#include <bits/stdc++.h>
#define LL long long
#define PII pair<int ,int>
#define MP(x, y) make_pair((x),(y))
#define fi first
#define se second
#define PB(x) push_back((x))
#define MMG(x) memset((x), -1,sizeof(x))
#define MMF(x) memset((x),0,sizeof(x))
#define MMI(x) memset((x), INF, sizeof(x))
using namespace std; const int INF = 0x3f3f3f3f;
const double eps = 1e-8; /////////
LL mul(LL x, LL y, LL mod)
{
return (x * y - (LL)(x / (long double)mod * y + 1e-3) * mod + mod) % mod;
} struct Matrix
{
LL m[2][2];
}; Matrix A;
Matrix I = {1, 0, 0, 1}; Matrix multi(Matrix a, Matrix b, LL MOD)
{
Matrix c;
for(int i = 0; i < 2; i++)
{
for(int j = 0; j < 2; j++)
{
c.m[i][j] = 0;
for(int k = 0; k < 2; k++)
c.m[i][j] = (c.m[i][j] % MOD + (a.m[i][k] % MOD) * (b.m[k][j] % MOD) % MOD) % MOD;
c.m[i][j] %= MOD;
}
}
return c;
} Matrix power(Matrix a, LL k, LL MOD)
{
Matrix ans = I, p = a;
while(k)
{
if(k & 1)
{
ans = multi(ans, p, MOD);
k--;
}
k >>= 1;
p = multi(p, p, MOD);
}
return ans;
} LL gcd(LL a, LL b)
{
return b ? gcd(b, a % b) : a;
} const int N = 400005;
const int NN = 5005; LL num[NN], pri[NN];
LL fac[NN];
int cnt, c; bool prime[N];
int p[N];
int k; void isprime()
{
k = 0;
memset(prime, true, sizeof(prime));
for(int i = 2; i < N; i++)
{
if(prime[i])
{
p[k++] = i;
for(int j = i + i; j < N; j += i)
prime[j] = false;
}
}
} LL fpow(LL a, LL b, LL m)
{
LL ans = 1;
a %= m;
while(b)
{
if(b & 1)
{
ans = mul(ans, a , m);
b--;
}
b >>= 1;
a = mul(a , a , m);
}
return ans;
} LL legendre(LL a, LL p)
{
if(fpow(a, (p - 1) >> 1, p) == 1)
return 1;
else return -1;
} void Solve(LL n, LL pri[], LL num[])
{
cnt = 0;
LL t = (LL)sqrt(1.0 * n);
for(int i = 0; p[i] <= t; i++)
{
if(n % p[i] == 0)
{
int a = 0;
pri[cnt] = p[i];
while(n % p[i] == 0)
{
a++;
n /= p[i];
}
num[cnt] = a;
cnt++;
}
}
if(n > 1)
{
pri[cnt] = n;
num[cnt] = 1;
cnt++;
}
} void Work(LL n)
{
c = 0;
LL t = (LL)sqrt(1.0 * n);
for(int i = 1; i <= t; i++)
{
if(n % i == 0)
{
if(i * i == n) fac[c++] = i;
else
{
fac[c++] = i;
fac[c++] = n / i;
}
}
}
} LL get_loop(LL n)
{
Solve(n, pri, num);
LL ans = 1;
for(int i = 0; i < cnt; i++)
{
LL record = 1;
if(pri[i] == 2)
record = 3;
else if(pri[i] == 3)
record = 8;
else if(pri[i] == 5)
record = 20;
else
{
if(legendre(5, pri[i]) == 1)
Work(pri[i] - 1);
else
Work(2 * (pri[i] + 1));
sort(fac, fac + c);
for(int k = 0; k < c; k++)
{
Matrix a = power(A, fac[k] - 1, pri[i]);
LL x = (a.m[0][0] % pri[i] + a.m[0][1] % pri[i]) % pri[i];
LL y = (a.m[1][0] % pri[i] + a.m[1][1] % pri[i]) % pri[i];
if(x == 1 && y == 0)
{
record = fac[k];
break;
}
}
}
for(int k = 1; k < num[i]; k++)
record *= pri[i];
ans = ans / gcd(ans, record) * record;
}
return ans;
} LL fib[5005]; void Init()
{
A.m[0][0] = 1;
A.m[0][1] = 1;
A.m[1][0] = 1;
A.m[1][1] = 0;
fib[0] = 0;
fib[1] = 1;
for(int i = 2; i < 5005; i++)
fib[i] = fib[i - 1] + fib[i - 2];
}
//////////
int main()
{
Init();
isprime();
int T;
cin >> T;
int icas = 0;
while(T--)
{
LL n;
cin >> n;
LL ans = get_loop(n);
printf("Case #%d: %lld\n", ++icas, ans);
}
return 0;
}

HDU 3977 斐波那契循环节的更多相关文章

  1. HDu4794 斐波那契循环节

    题意:Arnold变换把矩阵(x,y)变成((x+y)%n,(x+2*y)%n),问最小循环节 题解:仔细算前几项能看出是斐波那契数论modn,然后套个斐波那契循环节板子即可 //#pragma GC ...

  2. HDU 2814 斐波那契循环节 欧拉降幂

    一看就是欧拉降幂,问题是怎么求$fib(a^b)$,C给的那么小显然还是要找循环节.数据范围出的很那啥..unsigned long long注意用防爆的乘法 /** @Date : 2017-09- ...

  3. HDU 2855 斐波那契+矩阵快速幂

    http://acm.hdu.edu.cn/showproblem.php?pid=2855 化简这个公式,多写出几组就会发现规律 d[n]=F[2*n] 后面的任务就是矩阵快速幂拍一个斐波那契模板出 ...

  4. HDU 2516 斐波那契博弈

    点这里去看题 n为斐波那契数时,先手败,推断方法见算法讲堂 #include<bits/stdc++.h> using namespace std; int main() { ],i,n, ...

  5. HDU 1021(斐波那契数与因子3 **)

    题意是说在给定的一种满足每一项等于前两项之和的数列中,判断第 n 项的数字是否为 3 的倍数. 斐波那契数在到第四十多位的时候就会超出 int 存储范围,但是题目问的是是否为 3 的倍数,也就是模 3 ...

  6. hdu 5914(斐波拉契数列)

    Triangle Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)Total Su ...

  7. HDU——4549M斐波那契数列(矩阵快速幂+快速幂+费马小定理)

    M斐波那契数列 Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others) Total Su ...

  8. HDU 5686 斐波那契数列、Java求大数

    原题:http://acm.hdu.edu.cn/showproblem.php?pid=5686 当我们要求f[n]时,可以考虑为前n-1个1的情况有加了一个1. 此时有两种情况:当不适用第n个1进 ...

  9. HDU 1021 斐波那契

    参考自:https://www.cnblogs.com/ECJTUACM-873284962/p/6404504.html Fibonacci Again Time Limit: 2000/1000 ...

随机推荐

  1. Android源码项目目录结构

    src: 存放java代码 gen: 存放自动生成文件的. R.java 存放res文件夹下对应资源的id project.properties: 指定当前工程采用的开发工具包的版本 libs: 当前 ...

  2. bat获取当前日期的前一天

    批处理做这样的事情很麻烦,你可以用cscript来实现,比如把下面的内容保存为a.js文件:var d=new Date();d.setTime(d.getTime()-24*3600*1000);v ...

  3. 1014C程序语法树

    程序:冒泡算法C程序 #include <stdio.h> main() { int i,j,temp; int a[10]; for(i=0;i<10;i++) scanf (&q ...

  4. windows下面安装python3遇到的没有添加到环境变量的问题

    windows下面安装python3出现的问题 在官网上面下载最新版的安装包进行安装,并勾选Add Python 3.5 to PATH 安装的过程中可能会出现没有添加到PATH路径的情况 默认的安装 ...

  5. linux 下svn忽略文件

    假设想忽略文件temp 1. cd到temp所在的目录下: 2. svn propedit svn:ignore . 注意:请别漏掉最后的点(.表示当前目录),如果报错请看下面 3. 打开的文件就是忽 ...

  6. asp.netMVC中权限控制论

    这里设想了一个简单的思路,如果用户登录了,也就是session中有值才可以在控制器中操作,这样先添加一个控制器,如下代码: public class AuthController : Controll ...

  7. 【数据库】百万级数据库SQL优化大总结

    网上关于SQL优化的教程很多,但是比较杂乱.近日有空整理了一下,写出来跟大家分享一下,其中有错误和不足的地方,还请大家纠正补充. 这篇文章我花费了大量的时间查找资料.修改.排版,希望大家阅读之后,感觉 ...

  8. SQL入门之多表查询

    如果查询需要针对两个或者更多个表,则在需要涉及到表的连接操作(join). 0.笛卡儿积 最简单的连接方式是直接在from子句中加入两个表,并且用join操作符隔开.形式为Table1 join Ta ...

  9. 51nod 1624 取余最短路(set)

    题意: 佳佳有一个n*m的带权矩阵,她想从(1,1)出发走到(n,m)且只能往右往下移动,她能得到的娱乐值为所经过的位置的权的总和. 有一天,她被下了恶毒的诅咒,这个诅咒的作用是将她的娱乐值变为对p取 ...

  10. 【BZOJ2878】【NOI2012】迷失游乐园(动态规划)

    [BZOJ2878][NOI2012]迷失游乐园(动态规划) 题面 BZOJ 题解 记得以前考试的时候做过这道题目 这题的暴力还是非常显然的,每次\(dfs\)一下就好了. 时间复杂度\(O(n^2) ...