题目描述

a[1]=a[2]=a[3]=1

a[x]=a[x-3]+a[x-1] (x>3)

求a数列的第n项对1000000007（10^9+7）取余的值。

输入输出格式

输入格式：

第一行一个整数T，表示询问个数。

以下T行，每行一个正整数n。

输出格式：

每行输出一个非负整数表示答案。

输入输出样例

输入样例#1：

输出样例#1：

4
9
19

说明

对于30%的数据 n<=100；

对于60%的数据 n<=2*10^7；

对于100%的数据 T<=100，n<=2*10^9；

题解：

其实这篇本来想写日记的，但是老师突然来到了我的身边，我就迅速把这道刚做完的题粘到了这里。

有始有终，就写吧。

谁可以告诉我，为什么我的暴力没有分！

众人OS：你错了呗。

#include <cstdio>
#include <cstring>
using namespace std;typedef long long LL;
const int mod=1e9+7;
int T, n;
LL tmp[3][3]={{0,0,1},{1,0,0},{0,1,1}};struct Matrix33{LL mat[3][3];Matrix33 operator *(Matrix33 b){Matrix33 m;for (int i=0; i<3; ++i) for (int j=0; j<3; ++j){m.mat[i][j]=0;for (int k=0; k<3; ++k)m.mat[i][j]=(m.mat[i][j]+(mat[i][k]*b.mat[k][j]%mod))%mod;}return m;}
}beg, unit, plus;Matrix33 get_mat(int n){memcpy(plus.mat, tmp, sizeof(tmp));Matrix33 ans=unit;while (n){if (n&1) ans=ans*plus;plus=plus*plus;n>>=1;}return ans;
}int main(){unit.mat[0][0]=unit.mat[1][1]=unit.mat[2][2]=1;beg.mat[0][0]=beg.mat[0][1]=beg.mat[0][2]=1;scanf("%d", &T);for (int tt=0; tt<T; ++tt){scanf("%d", &n);if (n<4) printf("1\n");else printf("%lld\n", (beg*get_mat(n-3)).mat[0][2]);}return 0;
}

一世安宁

转载于:https://www.cnblogs.com/GTBA/p/9441336.html

洛谷 P1939 【模板】矩阵加速（数列）相关推荐

信息学奥赛一本通 1937：【06NOIP普及组】数列 | 洛谷 P1062 [NOIP2006 普及组] 数列
[题目链接] ybt 1937:[06NOIP普及组]数列洛谷 P1062 [NOIP2006 普及组] 数列 [题目考点] 1. 数制 [解题思路] 如果k为2,那么这个数列第1项为202^02 ...
洛谷P1129 [ZJOI2007] 矩阵游戏题解
洛谷P1129 [ZJOI2007] 矩阵游戏题解题目链接:P1129 [ZJOI2007] 矩阵游戏题意:给定一张有黑白棋子的正方形棋盘,问存不存在解法使得经过若干次交换行或列的操作后,左上角 ...
【bzoj3240 洛谷P1397】矩阵游戏[NOI2013]（矩阵乘法+卡常）
题目传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=3240 这道题其实有普通快速幂+费马小定理的解法--然而我太弱了,一开始只想到了矩阵乘法的 ...
专题·树链剖分【including 洛谷·【模板】树链剖分
初见安~~~终于学会了树剖~~~ [兴奋]当初机房的大佬在学树剖的时候我反复强调过:"学树剖没有前途的!!!" 恩.真香. 一.重链与重儿子所谓树剖--树链剖分,就是赋予一个链的 ...
洛谷·【模板】点分树 | 震波【including 点分树
初见安-这里是传送门:洛谷P6329 [模板]点分树 | 震波一.点分树其实你会点分治的话,点分树就是把点分治时的重心提出来重新连城一棵树. 比如当前点是u,求出子树v的重心root后将root与 ...
洛谷P1939 【模板】矩阵加速（数列）
题目描述 a[1]=a[2]=a[3]=1 a[x]=a[x-3]+a[x-1] (x>3) 求a数列的第n项对1000000007(10^9+7)取余的值. 输入格式第一行一个整数T,表示询 ...
【洛谷P3390】矩阵快速幂（模板）
贴一下矩阵快速幂的模板 #include<iostream> #include<cstdio> #include<cstring> #include<stri ...
【洛谷 P5550】 Chino的数列【矩阵乘法】
解题思路矩阵乘法呀考虑构造一个操作矩阵 CC[1234]\begin{bmatrix} 1&2&3&4\\ \end{bmatrix}[1234] 首先把特殊情况s ...
洛谷 p3372 模板-线段树 1
题目描述如题,已知一个数列,你需要进行下面两种操作: 1.将某区间每一个数加上x 2.求出某区间每一个数的和输入输出格式输入格式: 第一行包含两个整数N.M,分别表示该数列数字的个数和操作的总个 ...

洛谷 P1939 【模板】矩阵加速（数列）