这道题的暴力分还是很良心嘛~~~~~

直接刚的话我发现本蒟蒻只会暴力,矩乘根本写不出来,然后让我们找一下规律,我们发现如果我们把这个序列在mod k的意义下摆出,并且在此过程中把值为1的的数减一,我们发现他可以成为一段一段的被0(我们在此只关注减1变为0的点)区间,我们继续分析我们分析出来了这样的性质:如果存在这样的点,那么他右边的点一定是两个重复的数,而且往后是fibonacci数列(重头开始)乘第一个数,那么他之后再出现这样的0,的充要条件是其后存在一个fibonacci数是这段数第一个数的逆元。

我们先介绍一个结论(本蒟蒻并不会证):斐波那契数列模k后一定是0,1,1开头的纯循环,而且这个循环节的长度≤6k。我们开一个数组vis[i]表示第一个在mod k意义下值为i的fibonacci数的位置,这样我们求出某个数的逆元就知道以这个数为首项的数段的长度了(这样我们就可以跳啦)。

现在我说一下到了现在我们做了什么,我们发现在mod k 的意义下,这个数列会是一个由0(再次强调我们在此只关注减1变为0的点)隔开散区间(一整个也算)开头,之后要么成为循环,要么不再有0。先解释成为循环:因为我们在mod k的意义下因此最多不过k个数段就可以找到循环。我们再解释一下不在有0,这个就是当这个数段的第一项在mod k的意义下没有逆元,或者有逆元但是找不到vis[]。

现在我们用k再乘一个不大的数的时间复杂度找出来我们要处理的假fibonacci在哪些位置减一,以及他到底是存在循环节还是到后来没有了0,现在我们就可以想到一些矩阵,进行操作,但是这些矩阵一定要满足可乘。对于循环节我们暴力处理两边和一个循环节,对于最后没有0,我们就在后面直接fibonacci。

坑:I.via[1]!=1!!!!我们要找的是他后面第一个1,而我们前两个数是受法律保护的。

  II.对于处理一段一段的,最后一段可能顶到n,不满

  |||.一定要注意矩阵乘没有交换律

#include <cstdio>
#include <cmath>
#include <cstring>
#include <iostream>
using namespace std;
typedef long long LL;
const LL N=;
LL n,k,p,vis[N],f[N*],Ola,Stop,L,R,pos[N];
LL a[][],b[][],temp_ab[][],F[],temp_F[],S[][];
bool huzhi[N];
LL GCD(LL x,LL y){
return x==?y:GCD(y%x,x);
}
inline LL Min(LL x,LL y){
return x<y?x:y;
}
inline void get_Ola(){
LL lim=(LL)sqrt(k+0.5);
LL x=k;Ola=k;
for(LL i=;i<=lim;i++)
if(x%i==){
Ola=Ola/i*(i-);
while(x%i==)x/=i;
}
if(x!=){
Ola=Ola/x*(x-);
}
}
//*******************GCD&&Ola*********************//
inline void Multi_One(){
memset(temp_F,,sizeof(temp_F));
for(int i=;i<=;i++)
for(int j=;j<=;j++)
temp_F[i]=(temp_F[i]+F[j]*a[j][i]%p+p)%p;
memcpy(F,temp_F,sizeof(F));
}
inline void Multi_Two(){
memset(temp_ab,,sizeof(temp_ab));
for(int i=;i<=;i++)
for(int j=;j<=;j++)
for(int l=;l<=;l++)
temp_ab[i][j]=(temp_ab[i][j]+a[i][l]*a[l][j]%p+p)%p;
memcpy(a,temp_ab,sizeof(a));
}
inline void Multi_Three(){
memset(temp_F,,sizeof(temp_F));
for(int i=;i<=;i++)
for(int j=;j<=;j++)
temp_F[i]=((temp_F[i]+F[j]*b[j][i]%p)%p+p)%p;
memcpy(F,temp_F,sizeof(F));
}
inline void Multi_Four(){
memset(temp_ab,,sizeof(temp_ab));
for(int i=;i<=;i++)
for(int j=;j<=;j++)
for(int l=;l<=;l++)
temp_ab[i][j]=(temp_ab[i][j]+S[i][l]*a[l][j]%p+p)%p;
memcpy(S,temp_ab,sizeof(S));
}
inline void Multi_Five(){
memset(temp_ab,,sizeof(temp_ab));
for(int i=;i<=;i++)
for(int j=;j<=;j++)
for(int l=;l<=;l++)
temp_ab[i][j]=(temp_ab[i][j]+a[i][l]*a[l][j]%p)%p;
memcpy(a,temp_ab,sizeof(a));
}
inline void Multi_Six(){
memset(temp_ab,,sizeof(temp_ab));
for(int i=;i<=;i++)
for(int j=;j<=;j++)
for(int l=;l<=;l++)
temp_ab[i][j]=(temp_ab[i][j]+S[i][l]*b[l][j]%p+p)%p;
memcpy(S,temp_ab,sizeof(S));
}
inline void Multi_Seven(){
memset(temp_F,,sizeof(temp_F));
for(int i=;i<=;i++)
for(int j=;j<=;j++)
temp_F[i]=(temp_F[i]+F[j]*S[j][i]%p+p)%p;
memcpy(F,temp_F,sizeof(F));
}
inline void Multi_E(){
memset(temp_ab,,sizeof(temp_ab));
for(int i=;i<=;i++)
for(int j=;j<=;j++)
for(int l=;l<=;l++)
temp_ab[i][j]=(temp_ab[i][j]+S[i][l]*S[l][j]%p+p)%p;
memcpy(S,temp_ab,sizeof(S));
}
//*********************Multi**********************//
inline LL Pow(LL x,LL y,LL P){
LL ans=;
while(y){
if(y&)ans=ans*x%P;
y>>=,x=x*x%P;
}
return ans;
}
inline void POW(int step){
while(step){
if(step&)Multi_Seven();
step>>=,Multi_E();
}
}
//**********************POW***********************//
inline void Stop_Forever(){
LL now=,step=;
while(){
if(step>=n){
Stop=n;
break;
}
if(pos[now]){
L=pos[now],R=step-;
break;
}
pos[now]=step;
if(huzhi[now]==){
Stop=step-;break;
}
LL Now=Pow(now,Ola-,k);
if(vis[Now]==){
Stop=step-;break;
}
step+=vis[Now];
now=f[vis[Now]-]*now%k;
}
}
//*********************Judge*********************//
inline void F_H(){
f[]=,f[]=;
for(int i=;i<=*k;i++)
f[i]=(f[i-]+f[i-])%k,vis[f[i]]=vis[f[i]]==?i:vis[f[i]];
for(int i=;i<k;i++)
if(GCD(i,k)==)huzhi[i]=;
}
inline void Pre(){
b[][]=;
b[][]=;
b[][]=-,b[][]=;
F[]=,F[]=F[]=;
}
inline void Init(){
scanf("%lld%lld%lld",&n,&k,&p);
F_H();
get_Ola();
Stop_Forever();
Pre();
}
//***********************Pre**********************//
inline void GO(LL step){
memset(a,,sizeof(a));
a[][]=;
a[][]=,a[][]=;
a[][]=;
while(step){
if(step&)Multi_One();
step>>=,Multi_Two();
}
}
inline void GO_ON(int step){
memset(a,,sizeof(a));
a[][]=;
a[][]=,a[][]=;
a[][]=;
while(step){
if(step&)Multi_Four();
step>>=,Multi_Five();
}
}
inline void Get_Boss(LL &now,LL &step){
S[][]=;
S[][]=;
S[][]=;
while(step<R){
LL Now=Pow(now,Ola-,k);
LL y=vis[Now];
if(!step)y--;
GO_ON(y);
if(!step)y++;
step+=y;
Multi_Six();
now=f[vis[Now]-]*now%k;
}
}
//****************Carry On******************//
inline void Work_Stop(){
LL now=,step=;
while(step<Stop){
LL Now=Pow(now,Ola-,k);
LL y=Min(n-step,vis[Now]);
if(!step)y--;
GO(y);
if(!step)y++;
step+=y;
if(y==vis[Now])Multi_Three();
now=f[vis[Now]-]*now%k;
}
GO(n-step);
printf("%lld",F[]);
}
inline void Work_Forever(){
LL now=,step=;
while(step<L-){
LL Now=Pow(now,Ola-,k);
LL y=vis[Now];
if(!step)y--;
GO(y);
if(!step)y++;
step+=y;
Multi_Three();
now=f[vis[Now]-]*now%k;
}
Get_Boss(now,step);
POW((n-L+)/(R-L+));
step=(n-L+)/(R-L+)*(R-L+)+L-;
while(step<n){
LL Now=Pow(now,Ola-,k);
LL y=Min(n-step,vis[Now]);
GO(y);
step+=y;
if(y==vis[Now])Multi_Three();
now=f[vis[Now]-]*now%k;
}
printf("%lld",F[]);
}
//*********************War*****************//
int main(){
Init();
if(Stop)Work_Stop();
else Work_Forever();
return ;
}

【BZOJ 2432】 [Noi2011]兔农 矩乘+数论的更多相关文章

  1. 2432: [Noi2011]兔农 - BZOJ

    Description 农夫栋栋近年收入不景气,正在他发愁如何能多赚点钱时,他听到隔壁的小朋友在讨论兔子繁殖的问题. 问题是这样的:第一个月初有一对刚出生的小兔子,经过两个月长大后,这对兔子从第三个月 ...

  2. [BZOJ2432][Noi2011]兔农 矩阵乘法+exgcd

    2432: [Noi2011]兔农 Time Limit: 10 Sec  Memory Limit: 256 MB Description 农夫栋栋近年收入不景气,正在他发愁如何能多赚点钱时,他听到 ...

  3. bzoj 2437[Noi2011]兔兔与蛋蛋 黑白染色二分图+博弈+匈牙利新姿势

    noi2011 兔兔与蛋蛋 题目大意 直接看原题吧 就是\(n*m\)的格子上有一些白棋和一些黑棋和唯一一个空格 兔兔先手,蛋蛋后手 兔兔要把与空格相邻的其中一个白棋移到空格里 蛋蛋要把与空格相邻的其 ...

  4. BZOJ2432 [Noi2011]兔农

    本文版权归ljh2000和博客园共有,欢迎转载,但须保留此声明,并给出原文链接,谢谢合作. 本文作者:ljh2000作者博客:http://www.cnblogs.com/ljh2000-jump/转 ...

  5. NOI2011 兔农

    http://www.lydsy.com/JudgeOnline/problem.php?id=2432 感觉是day1中最难的一题,还好出题人很良心,给了75分部分分. 还是跪拜策爷吧~Orz ht ...

  6. bzoj 2437: [Noi2011]兔兔与蛋蛋

    Description Solution 考虑犯错误的条件:之前是处于必胜状态,该操作之后就变成了必败状态. 我们可以把这个过程看成两人对网格图进行黑白染色,变成了一个二分图模型,即当前位置向相邻不同 ...

  7. BZOJ.2437.[NOI2011]兔兔与蛋蛋游戏(二分图博弈 匈牙利)

    题目链接 首先空格的移动等价于棋子在黑白格交替移动(设起点移向白格就是黑色),且不会走到到起点距离为奇数的黑格.到起点距离为偶数的白格(删掉就行了),且不会重复走一个格子. (然后策略就同上题了,只不 ...

  8. 【BZOJ2432】【NOI2011】兔农(数论,矩阵快速幂)

    [BZOJ2432][NOI2011]兔农(数论,矩阵快速幂) 题面 BZOJ 题解 这题\(75\)分就是送的,我什么都不想写. 先手玩一下,发现每次每次出现\(mod\ K=1\)的数之后 把它减 ...

  9. 【BZOJ 2437】 2437: [Noi2011]兔兔与蛋蛋 (博弈+二分图匹配**)

    未经博主同意不得转载 2437: [Noi2011]兔兔与蛋蛋 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 693  Solved: 442 Des ...

随机推荐

  1. Python接受流式输入

    随笔记录——Python接受终端入若干行输入 Python接受终端的若干行输入时,比较常用的input()不再好用. 1. 导入sys模块: import sys 2. for循环接受输入: for ...

  2. CRC校验8

    什么是CRC校验? CRC即循环冗余校验码:是数据通信领域中最常用的一种查错校验码,其特征是信息字段和校验字段的长度可以任意选定.循环冗余检查(CRC)是一种数据传输检错功能,对数据进行多项式计算,并 ...

  3. ubuntu 14.04离线安装docker和docker compose

    准备安装包 1.下载docker docker支持以下版本的ubuntu系统Artful 17.10 (Docker CE 17.11 Edge)Zesty 17.04Xenial 16.04 (LT ...

  4. discuzX3.2 X3.4网站漏洞修复 SQL注入与请求伪造攻击利用与修复

    2018年12月9日,国内某安全组织,对discuz X3.2 X3.4版本的漏洞进行了公开,这次漏洞影响范围较大,具体漏洞是discuz 的用户前段SQL注入与请求伪造漏洞,也俗称SSRF漏洞,漏洞 ...

  5. Leecode刷题之旅-C语言/python-100相同的树

    /* * @lc app=leetcode.cn id=100 lang=c * * [100] 相同的树 * * https://leetcode-cn.com/problems/same-tree ...

  6. C# FTP上传文件时出现"应 PASV 命令的请求,服务器返回了一个与 FTP 连接地址不同的地址。"的错误

    FTP上传文件时出现"应 PASV 命令的请求,服务器返回了一个与 FTP 连接地址不同的地址."的错误 解决方法是在原代码上增加这句话 reqFTP.UsePassive = f ...

  7. 3155: Preprefix sum

    3155: Preprefix sum https://www.lydsy.com/JudgeOnline/problem.php?id=3155 分析: 区间修改,区间查询,线段树就好了. 然后,这 ...

  8. linux redhat 打开防火墙中的某个端口

    服务器成功监听了一个端口(如 5500),但是外面连接不进来,telnet其端口不通,解决办法如下(在root用户下): $ /sbin/iptables -I INPUT -p tcp --dpor ...

  9. 我所认识的XPath

    实例demo 测试demo所需要xml测试数据 <?xml version="1.0" encoding="iso-8859-1"?> <bo ...

  10. 21天学习caffe(一)

    ubuntu环境安装caffe1 安装依赖 apt-get install libatlas-base-dev apt-get install python-dev apt-get install l ...