HDU 5950——Recursive sequence
- Recursive sequence
题意:给定起始的两个数a,b,求第n个数%mod
思路:
第一眼看出来是矩阵快速幂,不过当时没推出来(其实推了一半把自己给否定了)
正解是,根据二项式定理,对于f(n)=f(n-1)+2f(n-2)+n^4可以转换为
f(n)=f(n-1)+2f(n-2)+c(4,0)(n-1)^4+c(4,1)(n-1)^3+c(4,2)(n-1)^2+c(4,3)(n-1)+c(4,4) (#1)
那么,对于下一项:
f(n+1)=f(n)+2f(n-1)+c(4,0)(n)^4+c(4,1)(n)^3+c(4,2)(n)^2+c(4,3)(n)+c(4,4) (#2)
可以发现#2的f(n),f(n-1),n^4,n^3,n^2,n,1都可以通过#1构造出来,那么根据倍数关系构造相关的矩阵即可。
code:
#include <bits/stdc++.h>
using namespace std;typedef long long ll;
const ll mod = 2147493647ll;struct Matirx{int r, c;ll mat[10][10];
};Matirx Mul(Matirx a, Matirx b){Matirx ret;ret.r = a.r;ret.c = b.c;for(int i = 0; i < ret.r; i++){for(int j = 0; j < ret.c; j++){ret.mat[i][j] = 0;for(int k = 0; k < a.c; k++){ret.mat[i][j] += a.mat[i][k]*b.mat[k][j];ret.mat[i][j] %= mod;}}}return ret;
}Matirx Quick_pow(Matirx a, ll n){Matirx ret;ret.r = ret.c = 7;memset(ret.mat, 0, sizeof(ret.mat));for(int i = 0; i < 7; i++) ret.mat[i][i] = 1;while(n){if(n&1) ret = Mul(ret, a);a = Mul(a, a);n >>= 1;}return ret;
}void build(Matirx& t){t.r = t.c = 7;t.mat[0][0] = 1, t.mat[0][1] = 4, t.mat[0][2] = 6, t.mat[0][3] = 4, t.mat[0][4] = 1, t.mat[0][5] = 0, t.mat[0][6] = 0;//1t.mat[1][0] = 0, t.mat[1][1] = 1, t.mat[1][2] = 3, t.mat[1][3] = 3, t.mat[1][4] = 1, t.mat[1][5] = 0, t.mat[1][6] = 0;//2t.mat[2][0] = 0, t.mat[2][1] = 0, t.mat[2][2] = 1, t.mat[2][3] = 2, t.mat[2][4] = 1, t.mat[2][5] = 0, t.mat[2][6] = 0;//2t.mat[3][0] = 0, t.mat[3][1] = 0, t.mat[3][2] = 0, t.mat[3][3] = 1, t.mat[3][4] = 1, t.mat[3][5] = 0, t.mat[3][6] = 0;//4t.mat[4][0] = 0, t.mat[4][1] = 0, t.mat[4][2] = 0, t.mat[4][3] = 0, t.mat[4][4] = 1, t.mat[4][5] = 0, t.mat[4][6] = 0;//5t.mat[5][0] = 0, t.mat[5][1] = 0, t.mat[5][2] = 0, t.mat[5][3] = 0, t.mat[5][4] = 0, t.mat[5][5] = 0, t.mat[5][6] = 1;//6t.mat[6][0] = 1, t.mat[6][1] = 4, t.mat[6][2] = 6, t.mat[6][3] = 4, t.mat[6][4] = 1, t.mat[6][5] = 2, t.mat[6][6] = 1;//7
}int main(){int T;scanf("%d", &T);ll N, a, b;while(T--){scanf("%I64d%I64d%I64d", &N, &a, &b);if(N == 1) printf("%I64d\n", a%mod);else if(N == 2) printf("%I64d\n", b%mod);else{Matirx temp;build(temp);Matirx ret = Quick_pow(temp, N-2);temp.r = 7;temp.c = 1;temp.mat[0][0] = 16,temp.mat[1][0] = 8,temp.mat[2][0] = 4,temp.mat[3][0] = 2,temp.mat[4][0] = 1,temp.mat[5][0] = a%mod, temp.mat[6][0] = b%mod;ret = Mul(ret, temp);printf("%I64d\n", ret.mat[6][0]);}}return 0;
}
HDU 5950——Recursive sequence相关推荐
- HDU 5950 Recursive sequence(矩阵快速幂)
题目链接:Recursive sequence 题意:给出前两项和递推式,求第n项的值. 题解:递推式为:$F[i]=F[i-1]+2*f[i-2]+i^4$ 主要问题是$i^4$处理,容易想到用矩阵 ...
- Recursive sequence HDU - 5950
Recursive sequence HDU - 5950 题意: 给你一个式子:f[n]=2f[n-2]+f[n-1]+n4 给你f[1]和f[2],给你一个n,求f[n] f[1],f[2],n& ...
- HDU 1005 Number Sequence
[题目] Number Sequence Time Limit: 2000/1000 MS (Jav ...
- Codeforces 1106F Lunar New Year and a Recursive Sequence (数学、线性代数、线性递推、数论、BSGS、扩展欧几里得算法)...
Codeforces 1106F Lunar New Year and a Recursive Sequence (数学.线性代数.线性递推.数论.BSGS.扩展欧几里得算法) 哎呀大水题..我写了一 ...
- 数论练习二之BSGS算法——随机数生成器,Matrix,Lunar New Year and a Recursive Sequence,Fermat‘s Last Theorem
[SDOI2013] 随机数生成器 description solution 肯定是非常想找一个通项公式来表示第nnn个数的 依据形式,考虑化成等比数列 xi+1+k=a(xi+k)=a⋅xi+b+t ...
- [矩阵乘法/快速幂专题]Arc of Dream,Recursive sequence,233 Matrix,Training little cats
矩阵快速幂习题 复习矩阵乘法及快速幂模板 乘法模板 快速幂模板 T1:Arc of Dream 题目 题解 code T2:Recursive sequence 题目 题解 code T3:233 M ...
- HDU 5400 Arithmetic Sequence
HDU 5400 Arithmetic Sequence /** HDU 5400 Arithmetic Sequence 直接预处理求解就好了 预处理找出以a[i]结尾最长的subArr长度(满足条 ...
- HDU.1005 Number Sequence
原题 HDU.1005 Number Sequence 分类 杂题 题意 给定一个数列{an}\left\{ a_n \right\}{an}的前两项a1a_1a1.a2a_2a2,以及其递推公 ...
- HDU 1560 DNA sequence(DNA序列)
HDU 1560 DNA sequence(DNA序列) Time Limit: 15000/5000 MS (Java/Others) Memory Limit: 32768/32768 K ...
最新文章
- 送你16个matplotlib绘图实用小技巧(附代码)
- 实现集合类的元素删除和修改的一点实践。。。
- SDN控制器部署相关工作—Vecloud
- python自动化弹框_Python+webdriver自动化脚本弹出框定位
- Python 网络爬虫笔记10 -- Scrapy 使用入门
- php +号在传输参数的过程中被变为空格了_编程难学?web相关知识,跟着淼哥学php全栈之路6...
- Web 2.0 编程思想:16条法则
- mysql修改存储引擎报错_修改mysql存储引擎备份数据库报错及解决方案
- 11g cursor_sharing 参数说明
- Go设计模式之Factory
- 读书笔记 - 《疯狂的站长》
- 医学统计学傻瓜教程(2019版)*
- python下雨动画特效_pygame实现烟雨蒙蒙下彩虹雨
- Axure8元件库.rplib(Iaxure)
- 配置 JAVA_HOME
- Weir Flow Control售予First Reserve事宜完成之后更名为Trillium Flow Technologies
- 用计算机高新证成都入户,关于2021年成都高新区小升初“报名电脑随机录取”的公告发布!...
- Windows 注册表操作 reg 命令详解
- Qt-qss之QSlider滑动条美化
- 搜索门户 Web Server分析