P5175 数列（矩阵快速幂）

P5175 数列

anb=(x×an−1+y×an−2)2x2×an−12+y2×an−22+2×x×y×an−1an−2x2×an−12+y2×an−22+2×x×y×an−2(x×an−2+y×an−3)x2×an−12+y2×an−22+2×x×y×(x×an−22+y×an−2×an−3)a_n ^ b = \left(x \times a_{n - 1} + y \times a_{n - 2}\right) ^ 2\\ x ^ 2 \times a_{n - 1} ^ 2 + y ^ 2 \times a_{n - 2} ^ 2 + 2 \times x \times y \times a_{n - 1} a_{n - 2}\\ x ^ 2 \times a_{n - 1} ^ 2 + y ^ 2 \times a_{n - 2} ^ 2 + 2 \times x \times y \times a_{n - 2} \left(x \times a_{n - 2} + y \times a_{n - 3} \right)\\ x ^ 2 \times a_{n - 1} ^ 2 + y ^ 2 \times a_{n - 2} ^ 2 + 2 \times x \times y \times \left(x \times a_{n - 2} ^ 2 + y \times a_{n - 2} \times a_{n - 3} \right)\\ anb=(x×an−1+y×an−2)2x2×an−12+y2×an−22+2×x×y×an−1an−2x2×an−12+y2×an−22+2×x×y×an−2(x×an−2+y×an−3)x2×an−12+y2×an−22+2×x×y×(x×an−22+y×an−2×an−3)
Sn+1=Sn+an+12S_{n + 1} = S_{n} + a_{n + 1} ^ 2Sn+1=Sn+an+12,

an+12=x2×an+y2×an−1+2×x×y×an×an−1a_{n + 1} ^ 2 = x ^ 2 \times a_n + y ^ 2 \times a_{n - 1} + 2 \times x \times y \times a_n \times a_{n - 1}an+12=x2×an+y2×an−1+2×x×y×an×an−1,

an+1an=x×an2+y×an×an−1a_{n + 1} a_n = x \times a_n ^ 2 + y \times a_{n} \times a_{n - 1}an+1an=x×an2+y×an×an−1,

[Snan2an−12anan−1]\begin{bmatrix}S_n & a_n ^ 2 & a_{n - 1} ^ 2 & a_n a_{n - 1}\end{bmatrix}[Snan2an−12anan−1] ×\times× [1000x2x21xy2y2002×x×y2×x×y0y]\begin{bmatrix} 1 & 0 & 0 & 0\\ x ^ 2 & x ^ 2 & 1 & x\\ y ^ 2 & y ^ 2 & 0 & 0\\ 2 \times x \times y & 2 \times x \times y & 0 & y\end{bmatrix}⎣⎢⎢⎡1x2y22×x×y0x2y22×x×y01000x0y⎦⎥⎥⎤

#include <bits/stdc++.h>using namespace std;const int mod = 1e9 + 7;inline int add(int x, int y) {return x + y < mod ? x + y : x + y - mod;
}struct Matrix {int a[4][4];
};Matrix operator * (Matrix a, Matrix b) {Matrix ans;for (int i = 0; i < 4; i++) {for (int j = 0; j < 4; j++) {ans.a[i][j] = 0;for (int k = 0; k < 4; k++) {ans.a[i][j] = add(ans.a[i][j], 1ll * a.a[i][k] * b.a[k][j] % mod);}}}return ans;
}int main() {//  freopen("in.txt", "r", stdin);// freopen("out.txt", "w", stdout);int T;scanf("%d", &T);while (T--) {long long n;int a1, a2, x, y;scanf("%lld %d %d %d %d", &n, &a1, &a2, &x, &y);if (n == 1) {printf("%d\n", 1ll * a1 * a1 % mod);}else if (n == 2) {printf("%d\n", add(1ll * a1 * a1 % mod, 1ll * a2 * a2 % mod));}else {Matrix ans = {add(1ll * a1 * a1 % mod, 1ll * a2 * a2 % mod), int(1ll * a2 * a2 % mod), int(1ll * a1 * a1 % mod), int(1ll * a1 * a2 % mod),0, 0, 0, 0,0, 0, 0, 0,0, 0, 0, 0,};Matrix a = {1, 0, 0, 0,int(1ll * x * x % mod), int(1ll * x * x % mod), 1, x,int(1ll * y * y % mod), int(1ll * y * y % mod), 0, 0,int(2ll * x * y % mod), int(2ll * x * y % mod), 0, y};n -= 2;while (n) {if (n & 1) {ans = ans * a;}a = a * a;n >>= 1;}printf("%d\n", ans.a[0][0]);}}return 0;
}

P5175 数列（矩阵快速幂）相关推荐

HDU4549 M斐波那契数列(矩阵快速幂+费马小定理)
Problem Description M斐波那契数列F[n]是一种整数数列,它的定义如下: F[0] = a F[1] = b F[n] = F[n-1] * F[n-2] ( n > 1 ) ...
UVA10689 Yet another Number Sequence【数列+矩阵快速幂】
Let's define another number sequence, given by the following function: f(0) = a f(1) = b f(n) = f(n ...
HDU 4549 M斐波那契数列(矩阵快速幂费马小定理)
ps:今天和战友聊到矩阵快速幂,想到前几天学长推荐去刷矩阵专题,挑了其中唯一一道中文题,没想到越过山却被河挡住去路... 题目链接:[kuangbin带你飞]专题十九矩阵 R - M斐波那契数列 T ...
斐波那契数列矩阵快速幂
题意用矩阵乘法求fibonacci数列的第n项. Solution 矩乘入门题啊,题目把题解已经说的很清楚里= =. 矩乘其实很简单,通过自己YY或者是搜索对于一个递推公式求出它所对应的矩阵,然后套 ...
Fibonacci数列矩阵快速幂
点击打开链接问题 A: 一道简单的Fibonacci 时间限制: 1 Sec 内存限制: 32 MB 题目描述有如下数列:F(0) = 7, F(1) = 11, F(n) = F(n - 1) ...
HDU4549 M斐波那契数列 —— 斐波那契、费马小定理、矩阵快速幂
题目链接:https://vjudge.net/problem/HDU-4549 M斐波那契数列 Time Limit: 3000/1000 MS (Java/Others) Memory Li ...
hdu 4549 M斐波那契数列（费马小定理 + 二分快速幂 + 矩阵快速幂）
M斐波那契数列 Time Limit: 3000/1 ...
矩阵快速幂及斐波那契数列模板
本篇博客先给出矩阵快速幂以及利用矩阵快速幂求斐波那契数列的模板,讲解待更新-- const int N=10; int tmp[N][N]; void multi(int a[][N],int b[] ...
牛客18987 粉嘤花之恋（矩阵快速幂、斐波那契数列）
链接:https://ac.nowcoder.com/acm/problem/18987 来源:牛客网时间限制:C/C++ 1秒,其他语言2秒空间限制:C/C++ 32768K,其他语言65536 ...
矩阵快速幂求斐波那契数列（初学整理）
参考文章: http://blog.csdn.net/u013795055/article/details/38599321 http://blog.csdn.net/g_congratulation ...

P5175 数列（矩阵快速幂）

P5175 数列

P5175 数列（矩阵快速幂）相关推荐

最新文章

热门文章