我居然没看题解瞎搞出来了?

题解:


不难想到找到每次减1的位置,然后减去它对最终答案的贡献。

假设有一个地方是\(x,1(mod~k)\)

那么减了1后就变成了\(x,0\)。

然后可以推到\(x,0,x,x\)。

可以看做以\(x,x\)为开头,做新的序列。

设y为x在mod k意义下的逆元,那么下次1的地方就是原斐波拉契序列中第一次出现y的位置。

如果没有逆元那就结束了。

原斐波拉契序列的非循环长度是\(O(k)\)级别的,这个可以通过随机序列第一次出现相同元素来理解,为在mo意义下应该是可以看做随机的,值域大小是\(O(k^2)\),那么期望长度就是\(O(\sqrt {k^2})\)。

所以斐波拉契序列做个3k左右然后预处理每一个元素第一次出现的位置。

然后x也是有循环的,因为x<k,所以肯定是\(O(k)\)的,那么就可以一起做了。

最后就是个等比矩阵求和,由于不一定有逆,所以要用分治法去求等比矩阵和。

由于循环开始之前和最后结尾的地方都要判,所以有点复杂。

Code:


#include<bits/stdc++.h>
#define fo(i, x, y) for(int i = x, B = y; i <= B; i ++)
#define ff(i, x, y) for(int i = x, B = y; i < B; i ++)
#define fd(i, x, y) for(int i = x, B = y; i >= B; i --)
#define ll long long
#define pp printf
#define hh pp("\n")
using namespace std; ll n; int k, mo; const int N = 3e6 + 5; int f[N], fi[N];
int p[N], p0, us[N]; int gcd(int x, int y) { return (!y ? x : gcd(y, x % y));}
void exgcd(int a, int b, int &x, int &y) {
if(!b) {x = a, y = 0; return;}
exgcd(b, a % b, y, x); y -= (a / b) * x;
}
int qni(int p, int q) {
int x, y;
exgcd(p, q, x, y);
x = (x % q + q) % q;
return x;
} int l = -1, r; void work() {
int x = 1;
while(1) {
if(gcd(x, k) > 1) break;
int y = qni(x, k);
if(fi[y]) {
p0 ++;
p[p0] = p[p0 - 1] + fi[y];
if(us[y]) {
l = us[y]; r = p0;
break;
}
us[y] = p0;
x = (ll) f[fi[y] - 1] * x % k;
} else break;
}
} struct jz {
ll a[2][2];
jz() {
a[0][0] = a[1][1] = 1;
a[0][1] = a[1][0] = 0;
}
}; jz operator * (jz a, jz b) {
jz c;
fo(i, 0, 1) fo(j, 0, 1)
c.a[i][j] = (a.a[i][0] * b.a[0][j] + a.a[i][1] * b.a[1][j]) % mo;
return c;
} jz operator + (jz a, jz b) {
fo(i, 0, 1) fo(j, 0, 1) a.a[i][j] = (a.a[i][j] + b.a[i][j]) % mo;
return a;
} jz operator - (jz a, jz b) {
fo(i, 0, 1) fo(j, 0, 1) a.a[i][j] = (a.a[i][j] - b.a[i][j] + mo) % mo;
return a;
} jz ksm(jz x, ll y) {
jz s;
for(; y; y /= 2, x = x * x)
if(y & 1) s = s * x;
return s;
} jz z; jz ksb(jz x, ll y) {
if(y == 1) return x;
if(y & 1) return ksb(x, y - 1) * x + x;
jz a = ksb(x, y / 2);
return a + a * ksm(x, y / 2);
}
jz calc(jz x, ll y) {
jz s = jz();
if(y > 0) s = s + ksb(x, y);
return s;
} int main() {
scanf("%lld %d %d", &n, &k, &mo);
f[1] = f[2] = 1;
fo(i, 3, 3000000) {
f[i] = (f[i - 2] + f[i - 1]) % k;
if(!fi[f[i]]) fi[f[i]] = i;
}
work();
z.a[0][1] = z.a[1][0] = z.a[1][1] = 1; z.a[0][0] = 0;
ll ans = ksm(z, n).a[1][0];
if(l == -1) {
fo(i, 1, p0) if(p[i] <= n)
ans = (ans - ksm(z, n - p[i]).a[1][1] + mo) % mo;
} else {
fo(i, 1, l) if(p[i] <= n)
ans = (ans - ksm(z, n - p[i]).a[1][1] + mo) % mo;
if(p[l] <= n) {
ll v = (n - p[l]) / (p[r] - p[l]);
fo(i, l + 1, r) {
ll n2 = (ll) p[i] + v * (p[r] - p[l]);
if(n2 <= n) ans = (ans - ksm(z, n - n2).a[1][1] + mo) % mo;
}
if(v > 0) {
ll st = v * (p[r] - p[l]) + p[l];
jz sa; memset(sa.a, 0, sizeof sa.a);
fo(i, l + 1, r) sa = sa + ksm(z, p[r] - p[i]);
jz c = ksm(z, n - st) * sa * calc(ksm(z, p[r] - p[l]), v - 1);
ans = (ans - c.a[1][1] + mo) % mo;
}
}
}
pp("%lld\n", ans);
}

【NOI2011】兔农(循环节)的更多相关文章

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

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

  2. BZOJ2432 [Noi2011]兔农

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

  3. 2432: [Noi2011]兔农 - BZOJ

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

  4. 【BZOJ 2432】 [Noi2011]兔农 矩乘+数论

    这道题的暴力分还是很良心嘛~~~~~ 直接刚的话我发现本蒟蒻只会暴力,矩乘根本写不出来,然后让我们找一下规律,我们发现如果我们把这个序列在mod k的意义下摆出,并且在此过程中把值为1的的数减一,我们 ...

  5. NOI2011 兔农

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

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

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

  7. BZOJ 2432 兔农

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

  8. HDU 5895 Mathematician QSC(矩阵乘法+循环节降幂+除法取模小技巧+快速幂)

    传送门:HDU 5895 Mathematician QSC 这是一篇很好的题解,我想讲的他基本都讲了http://blog.csdn.net/queuelovestack/article/detai ...

  9. hdu 2837 Calculation 指数循环节套路题

    Calculation Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total ...

  10. HDU 3746 (KMP求最小循环节) Cyclic Nacklace

    题意: 给出一个字符串,要求在后面添加最少的字符是的新串是循环的,且至少有两个循环节.输出最少需要添加字符的个数. 分析: 假设所给字符串为p[0...l-1],其长度为l 有这样一个结论: 这个串的 ...

随机推荐

  1. nginx支持webSocket ws请求

    服务端webSocket的java配置文件: @Override public void registerStompEndpoints(StompEndpointRegistry registry) ...

  2. Vue学习笔记【32】——Vue路由(watch、computed和methods之间的对比)

    computed属性的结果会被缓存,除非依赖的响应式属性变化才会重新计算.主要当作属性来使用: methods方法表示一个具体的操作,主要书写业务逻辑: watch一个对象,键是需要观察的表达式,值是 ...

  3. 【InnoDB】缓冲池

    索引目录 INNODB的体系结构 缓冲池 缓存中页的定位: checkpoint技术 INNODB的关键特性 插入缓冲 change buffer 两次写 以下的资料总结自:官方文档和<MySQ ...

  4. linux基础知识-常用命令

    ifconfig :查看当前ip hostname:查看主机名 vim /etc/hosts:修改地址映射 service iptables status : 查看防火墙状态 chkconfig ip ...

  5. OC学习篇之---类目的概念和使用

    上一篇文章介绍了OC中的@class关键字的使用http://blog.csdn.net/jiangwei0910410003/article/details/41774747,这一篇我们介绍一下,O ...

  6. SCP-bzoj-1019

    项目编号:bzoj-1019 项目等级:Safe 项目描述: 戳这里 特殊收容措施: 对于一个hanoi,知道了各种移动操作的优先级,也就确定了方案.可以证明对于盘子数为N的hanoi,任意移动方案都 ...

  7. 2019牛客多校第四场D-triples I 贪心

    D-triples 题意 给你一个\(n\),问至少有几个数或运算起来可以等于\(n\),并且输出数量和这个几个数.题目说明给的\(n\)一定符合条件(不会输出\(n= 1\) 之类不存在情况). 思 ...

  8. Service系统服务(一):安装一个KVM服务器、KVM平台构建及简单管理、virsh基本管理操作、xml配置文件的应用、为虚拟机制作快照备份、快建新虚拟机

    一.安装一个KVM服务器 目标: 本例要求准备一台 RHEL7.2 服务器,将其搭建为KVM平台,主要完成下列操作: 1> 关闭本机的SELinux保护.防火墙服务   2> 挂载RHEL ...

  9. linux系统的文件保护

    一些文件在Linux下看上去可能一切正常,但当您尝试删除的时候,居然也会报错,就像下边一样: [root@linux236 root]# ls -l 1.txt-rw-r--r-- 1 root ro ...

  10. express框架总结

    1.express教程及api : http://www.runoob.com/nodejs/nodejs-express-framework.html 2.nodejs的express自动生成项目框 ...