【输入】

一行两个整数 n P

【输出】

从小到大输出可能的 k,若不存在,输出 None

【样例输入 1】

5 5

【样例输出】

2

【样例解释】

f[0] = 2

f[1] = 2

f[2] = 4

f[3] = 6 mod 5 = 1

f[4] = 5 mod 5 = 0

f[5] = 1

30%的数据保证 n, P ≤ 1000

100%的数据保证 n, P ≤ 10^9


一道算是比较综合的数论题吧,感觉不是很难。

先用矩阵快速幂求出k=1时f[n]的值。

然后解一个k*f[n]+x*p=1的方程。(考到一半惊奇的发现不能快速幂求逆元233)

代码:

#include<bits/stdc++.h>
#define ll long long
#define a1 a[0][0]
#define a2 a[0][1]
#define a3 a[1][0]
#define a4 a[1][1]
using namespace std;
inline ll read(){
    ll ans=0;
    char ch=getchar();
    while(!isdigit(ch))ch=getchar();
    while(isdigit(ch))ans=(ans<<3)+(ans<<1)+(ch^48),ch=getchar();
    return ans;
}
ll n,mod;
struct Node{ll a[2][2];Node(){a1=0,a2=a3=a4=1;}};
inline Node operator*(Node a,Node b){
    Node ret;
    ret.a1=(a.a1*b.a1+a.a2*b.a3)%mod;
    ret.a2=(a.a1*b.a2+a.a2*b.a4)%mod;
    ret.a3=(a.a3*b.a1+a.a4*b.a3)%mod;
    ret.a4=(a.a3*b.a2+a.a4*b.a4)%mod;
    return ret;
}
ll ans=0,a,b;
ll Exgcd(ll a,ll b,ll &x,ll &y){
    ll ret;
    if(!b){
        x=1;
        y=0;
        return a;
    }
    else{
    ret=Exgcd(b,a%b,x,y);
    ll t=x;x=y;y=t-(a/b)*y;
    return ret;
    }
}
ll equa(ll a,ll b,ll c,ll & x,ll y){
    ll d=Exgcd(a,b,x,y);
    if(c%d)return 0;
    ll k=c/d;
    x*=k;y*=k;
    ans=(x%(b/d)+b/d)%(b/d);
    return 1;
}
inline Node ksm(ll p){
    Node x,ret;
    while(p){if(p&1)ret=ret*x;x=x*x,p>>=1;}
    return ret;
}
inline ll gcd(ll a,ll b){while(b){ll t=a;a=b,b=t%a;}return a;}
int main(){
    freopen("kfib.in","r",stdin);
    freopen("kfib.out","w",stdout);
    n=read(),mod=read();
    Node ttmp=ksm(n-2);
    ll tmp=(ttmp.a2+ttmp.a4)%mod;
    if(!tmp){puts("None");return 0;}
    if(gcd(tmp,mod)!=1){puts("None");return 0;}
    ll x,y;
    equa(tmp,mod,1,x,y);
    cout<<ans;
    return 0;
}

2018.08.30 NOIP模拟 kfib(矩阵快速幂+exgcd)的更多相关文章

  1. 2018.08.30 NOIP模拟 wall(模拟)

    [问题描述] 万里长城是中国强大的标志,长城在古代的用途主要用于快速传递军事消息和抵御 外敌,在长城上的烽火台即可以作为藏兵的堡垒有可以来点燃狼烟传递消息. 现在有一段 万里长城,一共有 N 个烽火台 ...

  2. 2018.08.30 NOIP模拟 graph(dfs序/树剖+线段树)

    [描述] 给你一个图,一共有 N 个点,2*N-2 条有向边. 边目录按两部分给出 1. 开始的 n-1 条边描述了一颗以 1 号点为根的生成树,即每个点都可以由 1 号点 到达. 2. 接下来的 N ...

  3. 2018.09.25 poj3070 Fibonacci(矩阵快速幂)

    传送门 矩阵快速幂板题,写一道来练练手. 这一次在poj做题总算没忘了改万能库. 代码: #include<iostream> #include<cstdio> #define ...

  4. CodeChef February Challenge 2018 Broken Clock (三角函数推导 + 矩阵快速幂)

    题目链接  Broken Clock   中文题面链接 令$cos(xα) = f(x)$ 根据三角函数变换公式有 $f(x) = \frac{2d}{l} f(x-1) - f(x-2)$ 我们现在 ...

  5. 2018.10.19 NOIP训练 桌子(快速幂优化dp)

    传送门 勉强算一道dp好题. 显然第kkk列和第k+nk+nk+n列放的棋子数是相同的. 因此只需要统计出前nnn列的选法数. 对于前mmm%nnn列,一共有(m−1)/n+1(m-1)/n+1(m− ...

  6. 2018.08.29 NOIP模拟 pmatrix(线性筛)

    [问题描述] 根据哥德巴赫猜想(每个不小于 6 的偶数都可以表示为两个奇素数之和),定义 哥德巴赫矩阵 A 如下:对于正整数对(i,j),若 i+j 为偶数且 i,j 均为奇素数,则 Ai,j = 1 ...

  7. 2018.08.22 NOIP模拟 string(模拟)

    string [描述] 给定两个字符串 s,t,其中 s 只包含小写字母以及*,t 只包含小写字母. 你可以进行任意多次操作,每次选择 s 中的一个*,将它修改为任意多个(可以是 0 个)它的前一个字 ...

  8. 2018.08.29 NOIP模拟 movie(状压dp/随机化贪心)

    [描述] 小石头喜欢看电影,选择有 N 部电影可供选择,每一部电影会在一天的不同时段播 放.他希望连续看 L 分钟的电影.因为电影院是他家开的,所以他可以在一部电影播放过程中任何时间进入或退出,当然他 ...

  9. 2018.08.22 NOIP模拟 or(线段树)

    or [描述] 构造一个长度为 n 的非负整数序列 x,满足 m 个条件,第 i 个条件为x[li] | x[li+1] | - | x[ri]=pi. [输入] 第一行两个整数 n,m.接下来 m ...

随机推荐

  1. xe7 控件升级

    rm.ehlib.synedit OK SynSQLSyn1->TableNames 为NULL,导致添加数据失败,XE6正常 放在按钮里也不正常,就不说初始化次序引起的.

  2. MySQL优化order by导致的 using filesort

    using filesort 一般出现在 使用了 order by 语句当中. using filesort不一定引起mysql的性能问题.但是如果查询次数非常多,那么每次在mysql中进行排序,还是 ...

  3. linux c++下载http文件并显示进度<转>

    #include <stdio.h> #include <sys/types.h> #include <sys/socket.h> #include <net ...

  4. ABAP-年月期间搜索帮助

    selection-screen begin of block block1 with frame title text-. parameters:p_mon1 like s031-spmon, p_ ...

  5. Apache配置本地域名

    打开Apache的安装目录,找到httpd.conf文件,分别去掉下面两行文字前面的#号. LoadModule vhost_alias_module modules/mod_vhost_alias. ...

  6. ccttp图片处理过程

    1.python输出404的图片 #!/usr/bin/python # coding: utf-8 import psycopg2 import sys from datetime import * ...

  7. Object-c 构造、析构函数

    一.构造函数 在OC中凡是已init开头的函数我们都称之为构造函数,在声明构造函数的时候,不带参数的一般直接声明为“-(id)init”,带参数的一般声明为“-(id)initWith...”. @i ...

  8. ant使用备忘

    ant是一个脚本构建工具,可能就是持续集成里面所需要的构建工具. 如果使用eclipse,里面会自带有ant工具,不需要再安装了,创建一个build.xml(或者其他的名字都可以),使用ant来运行就 ...

  9. spring 中c3p0的优化配置

    jdbc.properties driverClass=com.mysql.jdbc.Driver jdbcUrl=jdbc:mysql://localhost:3306/testdb user=ro ...

  10. linux下面得小数计算

    可以通过命令行向awk中传递参数. 这样子传递进去的参数,在awk命令中可以访问.每一项都必须作为单一的一个参数来解释.所以,等号之间不能有空格. 比如说我们传递进去一个学生名字,想查这个学生的得分. ...