upc 卡德加的兔子线段树 + 矩阵快速幂

卡德加的兔子
时间限制: 5 Sec 内存限制: 128 MB

题目描述
卡德加喜欢养兔⼦。他在达拉然的下⽔道⾥放了N个兔笼（编号1到N），⾥⾯养着他从德拉诺带来的兔⼦。它们的繁殖遵循斐波那契数列的规律：刚开始时，笼⼦⾥有⼀对刚出⽣的兔⼦。每对兔⼦在出⽣第⼆个⽉后，每个⽉都⽣⼀对兔⼦。（第⼀个⽉结束后有1对兔⼦。第⼆个⽉结束后有2对。）
卡德加从苏拉玛的⼤魔导⼠艾利桑德那边学习了先进的扭曲时空法术。
有时候，他会对⼀排连续的兔笼（从第L号到第R号）释放时光流逝法术，让这些兔笼⾥的时间前进K个⽉。另外⼀些时候，他想喂⼀下兔⼦，所以他想知道第L号到第R号兔笼⾥有多少只兔⼦。
（假设这些操作都是在⼀个⽉以内完成的，不需要考虑⾃然时间对兔⼦的影响。）
输入
第⼀⾏两个整数N;M，表⽰兔笼的数量和操作的数量。
接下来M⾏，每⾏包含三个数L;R;K。如果K>0，说明卡德加在使⽤时光流逝，编号L到R的兔笼时间前进K个⽉。如果K=0，说明他只是想喂兔⼦了，输出这些兔笼⾥有多少兔⼦。
输出
对每个喂兔⼦的操作，输出兔⼦的数量。答案模10007。
样例输入 Copy
100 100
26 85 0
2 81 1
67 69 0
22 69 2
27 80 0
79 87 2
57 63 2
76 80 3
95 99 4
45 94 2
27 75 1
22 75 0
12 51 2
25 66 2
11 61 1
83 88 0
27 83 1
14 97 1
71 90 3
10 44 4
73 93 4
1 98 4
23 93 3
56 76 3
13 25 2
57 68 1
18 28 3
19 28 0
21 78 4
10 95 3
93 98 4
73 81 0
28 44 0
20 53 0
59 75 4
53 69 2
9 54 1
55 95 4
14 46 4
22 50 2
35 98 3
34 93 4
91 92 4
43 53 2
45 64 2
58 87 4
64 74 3
36 62 4
16 98 0
71 76 3
39 99 1
6 16 3
12 73 2
37 82 4
56 92 1
81 99 0
9 51 3
27 80 4
13 60 3
24 39 0
13 17 1
46 54 0
81 84 4
45 61 1
21 87 0
12 61 1
52 98 0
25 64 2
3 68 0
61 64 0
24 45 4
16 53 2
23 59 3
68 86 4
55 74 1
68 74 4
22 87 1
5 21 1
3 89 1
5 84 2
18 58 1
47 81 3
74 80 1
54 71 0
9 16 2
29 55 4
12 85 1
49 84 0
30 71 3
50 51 0
51 51 0
2 38 0
14 92 3
5 31 1
5 53 1
34 83 1
29 69 0
34 82 2
50 93 1
83 96 2
样例输出 Copy
60
3
140
607
27
6408
1161
5444
7985
9779
9274
5261
9934
2118
9996
4133
2687
6932
7982
3829
2931
3078
8342
提示

首先要做这个题，需要用矩阵快速幂来求第 n 项的斐波那契数列，因为矩阵乘法是有分配律的，所以当前进 k 天，也就是 { 1 , 1 }{ 1 , 0 } 的 k 次方乘上原来线段树维护的矩阵值，以及懒标记的值，那么这个题就转换成了区间乘法了，唯一需要的一点就是掏出来个好的矩阵板子，以前没有存矩阵的板子，光这个板子就卡了我半天。
注意一下懒标的初始值设为单位矩阵，区间和设位 fib2矩阵({ 1 , 0 }{ 0 , 0 }) ，mp[ 0 ] [ 0 ]这个位置就是斐波那契的第1项，因为题意可知初始的时候都应该是第一项。
pushdown的时候注意哪个矩阵乘上哪个矩阵，矩阵乘法是没有交换律的，乘反了答案肯定是不对的。
套上板子之后，心情就变得愉悦起来了 ~
一开始还想用斐波那契数列模数有循环节打个表来做，结果这样的方法只能支持单点修改，修改一个区间是不现实的。。
改下代码，尽量变得可读性高一点吧，写的太丑了。。

#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define mid (tr[u].l+tr[u].r>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
using namespace std;typedef long long LL;
typedef pair<int,int> PII;const int N=100010,M=2,mod=10007,INF=0x3f3f3f3f;
const double eps=1e-6;struct Mat
{int r,c,mp[5][5];Mat(){ memset(mp,0,sizeof mp); }//构建0矩阵 void build_one()//构建单位矩阵 {memset(mp,0,sizeof mp);mp[0][0]=mp[1][1]=1;}void build_fib1()//构建{1,1}{1,0} {mp[1][0]=mp[0][1]=mp[0][0]=1;mp[1][1]=0;}void build_fib2()//构建{1,0}{0,0} {mp[0][0]=1;mp[1][0]=mp[1][1]=mp[0][1]=0;}Mat operator * (Mat m) const//矩阵乘法 {Mat ans;memset(ans.mp,0,sizeof ans.mp);for(int i=0;i<2;i++)for(int j=0;j<2;j++)for(int k=0;k<2;k++){ans.mp[i][j]+=(mp[i][k]*m.mp[k][j])%mod;ans.mp[i][j]%=mod;}return ans;}Mat operator + (Mat m) const//矩阵加法 {Mat ans;for(int i=0;i<2;i++)for(int j=0;j<2;j++)ans.mp[i][j]=(mp[i][j]+m.mp[i][j])%mod;return ans;}friend Mat operator ^ (Mat m,int b)//矩阵快速幂 {Mat ans;ans.build_one();while(b){if(b&1) ans=ans*m;m=m*m;b>>=1;}return ans;}
};struct Node
{int l,r,flag;//flag表示b是否需要往下传 Mat a,b;//a 是和  b 是懒标
}tr[N<<2];int n,m;void pushup(int u)
{tr[u].a=tr[L].a+tr[R].a;
}void pushdown(int u)
{if(!tr[u].flag) return;tr[L].b=tr[u].b*tr[L].b,tr[R].b=tr[u].b*tr[R].b;//更新儿子的懒标 tr[L].a=tr[u].b*tr[L].a,tr[R].a=tr[u].b*tr[R].a;//更新儿子的和 tr[u].flag=0,tr[L].flag=tr[R].flag=1;//更新标记 tr[u].b.build_one();//将懒标设位单位矩阵
}void build(int u,int l,int r)
{Mat x,y;y.build_one();tr[u]={l,r,0};tr[u].b.build_one();if(l==r){tr[u].a.build_fib2();return;}build(L,l,mid),build(R,mid+1,r);pushup(u);
}void modify(int u,int l,int r,Mat c)
{if(tr[u].l>=l&&tr[u].r<=r){tr[u].flag=1;tr[u].a=c*tr[u].a;tr[u].b=c*tr[u].b;return;}pushdown(u);if(l<=mid) modify(L,l,r,c);if(r>mid) modify(R,l,r,c);pushup(u);
}int query(int u,int l,int r)
{if(tr[u].l>=l&&tr[u].r<=r)return tr[u].a.mp[0][0]%mod;pushdown(u);int x=0;if(l<=mid) x=(x+query(L,l,r))%mod;if(r>mid) x=(x+query(R,l,r))%mod;return x%mod;
}int main()
{//  ios::sync_with_stdio(false);
//  cin.tie(0);Mat temp;temp.build_fib1();cin>>n>>m;build(1,1,n);while(m--){int l,r,k;scanf("%d%d%d",&l,&r,&k);if(k) modify(1,l,r,temp^k);else printf("%d\n",query(1,l,r)%mod);}return 0;
}
/**/

upc 卡德加的兔子线段树 + 矩阵快速幂相关推荐

2018.10.09 ZYH的斐波那契数列（线段树+矩阵快速幂）
描述 ZYH最近研究数列研究得入迷啦! 现在有一个斐波拉契数列(f[1]=f[2]=1,对于n>2有f[n]=f[n-1]+f[n-2]), 但是斐波拉契数列太简单啦,于是ZYH把它改成了斐波拉 ...
NOIP模拟题 [线段树][矩阵快速幂]
有一定难度,要深入挖掘问题特性. T1: 题意: 给定一个序列,每次操作把操作位置及其后面比它小的数按顺序排列(整体上仍在原来的位置),求每次操作后的逆序对数. 分析: 每个数对逆序对数都有一个贡献( ...
AtCoder abc256全题解（区间合并模板、矩阵快速幂优化dp、线段树……）
文章目录 A B C-枚举 D-区间合并模板 E-图论建模,函数图的性质题意思路代码 F-树状数组题意思路代码 G-矩阵快速幂优化dp H-线段树思路实现传送门本文CSDN 本文j ...
ZYH的斐波那契数列【线段树动态开点+矩阵快速幂求斐波那契】
描述 ZYH最近研究数列研究得入迷啦! 现在有一个斐波拉契数列(f[1]=f[2]=1,对于n>2有f[n]=f[n-1]+f[n-2]), 但是斐波拉契数列太简单啦,于是ZYH把它改成了斐波拉 ...
牛客 - 求函数(线段树+区间合并/线段树+矩阵维护)
题目链接:点击查看题目大意:现在有 n 个函数,每个函数都是诸如 f( x ) = k * x + b 的形式,只是每个函数的 k 和 b 都是相互独立的,现在给出两个操作: 1 pos k b:将 ...
[选拔赛1]花园（矩阵快速幂），JM的月亮神树（最短路），保护出题人（斜率优化）
多年不考试,一夜回到解放前 T1:花园 title solution code T2:月亮神树 title solution code T3:保护出题人 title solution code T1: ...
浅谈矩阵矩阵快速幂动态dp 矩阵求逆
文章目录 1 矩阵 1.1 矩阵乘法 1.2 单位矩阵 1.3 矩阵乘法的结合律 1.4 矩阵快速幂 1.5 矩阵快速幂的应用 1.5.1 1.5.2 1.5.3 1.5.4 1.5.5 1.5.6 ...
POJ 2778 DNA Sequence [AC自动机 + 矩阵快速幂]
http://poj.org/problem?id=2778 题意:给一些只由ACGT组成的模式串,问有多少种长度为n且不含有给出的模式串的DNA序列. 自动机的状态转换可以看成一个有向图(有重边的) ...
【BZOJ】2553: [BeiJing2011]禁忌 AC自动机+期望+矩阵快速幂
[题意]给定n个禁忌字符串和字符集大小alphabet,保证所有字符在集合内.一个字符串的禁忌伤害定义为分割能匹配到最多的禁忌字符串数量(一个可以匹配多次),求由字符集构成的长度为Len的字符串的期望 ...
[转]快速矩阵快速幂
大学生程序代写矩阵快速幂出处:http://www.cnblogs.com/yan-boy/archive/2012/11/29/2795294.html 矩阵的快速幂是用来高效地计算矩阵的高次 ...

upc 卡德加的兔子线段树 + 矩阵快速幂

upc 卡德加的兔子线段树 + 矩阵快速幂相关推荐

最新文章

热门文章

upc 卡德加的兔子 线段树 + 矩阵快速幂

upc 卡德加的兔子 线段树 + 矩阵快速幂相关推荐

最新文章

热门文章

upc 卡德加的兔子线段树 + 矩阵快速幂

upc 卡德加的兔子线段树 + 矩阵快速幂相关推荐