BZOJ5104 Fib数列 二次剩余、BSGS
发现只有通项公式可以解决考虑通项公式
\(F_n = \frac{1}{\sqrt{5}}((\frac{1+\sqrt{5}}{2})^n - (\frac{1-\sqrt{5}}{2})^n) = a\)
注意到根据二次互反律,在\(\mod 10^9+9\)意义下\(5\)存在二次剩余,所以先把\(\sqrt{5}\)对应的值算出来(实际上是\(383001016\))。
那么原式变为了\((\frac{1+\sqrt{5}}{2})^n - (\frac{1-\sqrt{5}}{2})^n = \sqrt{5}a\),\(x^n-y^n=a\)的方程好像还是不会解TAT
观察上面式子可以发现:\(\frac{1-\sqrt{5}}{2} = -\frac{1}{\frac{1 + \sqrt{5}}{2}}\),所以设\(\frac{1 + \sqrt{5}}{2} = x\),那么原式变为\(x - \frac{(-1)^n}{x} = \sqrt{5}a\)。
按照\(n\)的奇偶性讨论,可以得到一个一元二次方程,使用求根公式求出方程的解。注意到在求根公式中还有一个求\(\sqrt{\Delta}\)的操作,这里还需要计算一次二次剩余,如果这里二次剩余无解,那么方程就是无解的,否则就存在两个解。
然后我们最后的问题就是解出\((\frac{1 + \sqrt{5}}{2})^n = t\)的最小的\(n\),同时需要满足\(n\)是奇数或者是偶数。这里我们同样使用BSGS,只需要在分块的时候把块大小分成偶数,然后对于奇数和偶数分开考虑就可以了。
#include<bits/stdc++.h>
#include<tr1/unordered_map>
using namespace std;
using namespace std::tr1;
const int MOD = 1e9 + 9;
int val , g , sqrt5 , inv2;
int poww(long long a , int b){
int times = 1;
while(b){
if(b & 1) times = times * a % MOD;
a = a * a % MOD; b >>= 1;
}
return times;
}
struct BSGS{
unordered_map < int , int > odd , even , all;
int T , powbs;
void init(int base){
T = sqrt(MOD) + 1; if(T & 1) ++T;
powbs = poww(base , T);
int tms = 1;
for(int i = 0 ; i < T ; ++i , tms = 1ll * tms * base % MOD){
unordered_map < int , int > &now = i & 1 ? odd : even;
now[tms] = i; all[tms] = i;
}
}
int calc(int ans , int val){
unordered_map < int , int > &now = val == 0 ? all : (val == 1 ? odd : even);
int tms = 1ll * powbs * poww(ans , MOD - 2) % MOD;
for(int i = 1 ; i <= T ; ++i , tms = 1ll * tms * powbs % MOD)
if(now.find(tms) != now.end())
return i * T - now[tms];
return 2e9 + 1;
}
}A , B;
int findrt(){
int t = MOD - 1; vector < int > zys;
for(int i = 2 ; i * i <= t ; ++i)
if(t % i == 0){
zys.push_back(MOD / i);
while(t % i == 0) t /= i;
}
if(t - 1) zys.push_back(MOD / t);
for(int i = 2 ; ; ++i){
bool flg = 1;
for(int j = 0 ; j < zys.size() && flg ; ++j)
if(poww(i , zys[j]) == 1) flg = 0;
if(flg) return i;
}
}
int main(){
cin >> val; g = findrt(); A.init(g); sqrt5 = poww(g , A.calc(5 , 0) / 2);
val = 1ll * val * sqrt5 % MOD; inv2 = poww(2 , MOD - 2);
B.init((1ll + sqrt5) * inv2 % MOD);
int t = A.calc((1ll * val * val + MOD - 4) % MOD , 0) , ans = 2e9 + 1;
if(t % 2 == 0){
int sqt = poww(g , t / 2) , ans1 = 1ll * (val + sqt) * inv2 % MOD , ans2 = 1ll * (val - sqt + MOD) * inv2 % MOD;
ans = min(B.calc(ans1 , 1) , B.calc(ans2 , 1));
}
t = A.calc((1ll * val * val + 4) % MOD , 0);
if(t % 2 == 0){
int sqt = poww(g , t / 2) , ans1 = 1ll * (val + sqt) * inv2 % MOD , ans2 = 1ll * (val - sqt + MOD) * inv2 % MOD;
ans = min(ans , min(B.calc(ans1 , 2) , B.calc(ans2 , 2)));
}
if(ans == 2e9 + 1) puts("-1");
else cout << ans;
return 0;
}
BZOJ5104 Fib数列 二次剩余、BSGS的更多相关文章
- bzoj5104 Fib数列(BSGS+二次剩余)
快AFO了才第一次写二次剩余的题…… 显然应该将Fn写成通项公式(具体是什么写起来不方便而且大家也都知道),设t=((1+√5)/2)n,T=√5N,然后可以得到t-(-1)t/t=√5N,两边同时乘 ...
- 【BZOJ5104】Fib数列(BSGS,二次剩余)
[BZOJ5104]Fib数列(BSGS,二次剩余) 题面 BZOJ 题解 首先求出斐波那契数列的通项: 令\(A=\frac{1+\sqrt 5}{2},B=\frac{1-\sqrt 5}{2}\ ...
- bzoj5104: Fib数列
Description Fib数列为1,1,2,3,5,8... 求在Mod10^9+9的意义下,数字N在Fib数列中出现在哪个位置 无解输出-1 Input 一行,一个数字N,N < = 10 ...
- BZOJ5104 Fib数列(二次剩余+BSGS)
5在1e9+9下有二次剩余,那么fib的通项公式就有用了. 已知Fn,求n.注意到[(1+√5)/2]·[(1-√5)/2]=-1,于是换元,设t=[(1+√5)/2]n,原式变为√5·Fn=t-(- ...
- @bzoj - 5104@ Fib数列
目录 @description@ @solution@ @accepted code@ @details@ @description@ Fib数列为1,1,2,3,5,8... 求在Mod10^9+9 ...
- FIB数列
斐波那契级数除以N会出现循环,此周期称为皮萨诺周期. 下面给出证明 必然会出现循环 这是基于下面事实: 1. R(n+2)=F(n+2) mod P=(F(n+1)+F(n)) mod P=(F(n+ ...
- 动态规划之Fib数列类问题应用
一,问题描述 有个小孩上楼梯,共有N阶楼梯,小孩一次可以上1阶,2阶或者3阶.走到N阶楼梯,一共有多少种走法? 二,问题分析 DP之自顶向下分析方式: 爬到第N阶楼梯,一共只有三种情况(全划分,加法原 ...
- UVaLive 3357 Pinary (Fib数列+递归)
题意:求第 k 个不含前导 0 和连续 1 的二进制串. 析:1,10,100,101,1000,...很容易发现长度为 i 的二进制串的个数正好就是Fib数列的第 i 个数,因为第 i 个也有子问题 ...
- 【bzoj5118】Fib数列2 费马小定理+矩阵乘法
题目描述 Fib定义为Fib(0)=0,Fib(1)=1,对于n≥2,Fib(n)=Fib(n-1)+Fib(n-2) 现给出N,求Fib(2^n). 输入 本题有多组数据.第一行一个整数T,表示数据 ...
随机推荐
- 如何在Unity中开发Leap Motion桌面版(Non-VR)APP
最近因需要,翻出几年前的Leapmotion感测器,准备用Unity3D做个互动APP,于是连上官网下载SDK.等下载下来一安装调试,瞬间傻眼,居然要求VR设备.我们Lab倒是不缺VR,有几套VIVE ...
- 如何防范web前端安全攻击
一.对于XSS防御: 1.不要信任任何外部传入的数据,针对用户输入作相关的格式检查.过滤等操作,以及转义字符处理.最普遍的做法就是转义输入输出的内容,对于括号,尖括号,斜杠进行转义 function ...
- C博客作业01——分支,顺序结构
C博客作业01--分支,顺序结构 0.展示PTA总分 1本章学习内容 1.1学习内容总结 1)格式化输出函数printf(),scanf(). 它是什么? 对于初学者而言,一开始了解接触它们,只是被硬 ...
- fluent将出口温度赋值给入口
Fluent版本:Fluent18.2 首先我们启动Fluent 然后按照正常的流程进行模型缩放,材料的设置,边界条件的设置等等,然后初始化. 在完成了算例的初始化以后 (define (OutToI ...
- 软件工程 “校园汇” 个人IDEA竞选分析与总结
IDEA竞选 19/10/8软件工程课上举行了一次IDEA竞选: 我的竞选IDEA是"校友汇",大学生的在线活动中心. 投票结果: 可以看到,校友会(汇)IDEA竞选结果十分惨淡, ...
- SSH框架整合2
===========================================web.xml================================================== ...
- arts lettcode 题目
编写一个函数来查找字符串数组中的最长公共前缀. 如果不存在公共前缀,返回空字符串 "". 示例 1: 输入: ["flower","flow" ...
- 利用Shell命令与HDFS进行交互
以”./bin/dfs dfs”开头的Shell命令方式 1.目录操作 在HDFS中为hadoop用户创建一个用户目录(hadoop用户) 在用户目录下创建一个input目录, HDFS的根目录下创建 ...
- MYSQL 什么时候用单列索引?什么使用用联合索引?
我一个表 students 表,有3个字段 ,id,name,age 我要查询 通过 name 和age,在这两个字段 是创建 联合索引?还是分别在name和age上创建 单列索引呢? 多个字段查询什 ...
- NoSql数据库Redis系列(6)——Redis数据过期策略详解
本文对Redis的过期机制简单的讲解一下 讲解之前我们先抛出一个问题,我们知道很多时候服务器经常会用到redis作为缓存,有很多数据都是临时缓存一下,可能用过之后很久都不会再用到了(比如暂存sessi ...