【HDU5950】Recursive sequence(矩阵快速幂)
BUPT2017 wintertraining(15) #6F
题意
\(f(1)=a,f(2)=b,f(i)=2*(f(i-2)+f(i-1)+i^4)\)
给定n,a,b ,\(N,a,b < 2^{31}\),求f(n)% 2147493647。
题解
i^4=(i-1)^4+4*(i-1)^3+6*(i-1)^2+4*(i-1)+1
\]
我们可以构造出矩阵乘法
\begin{matrix}
f_{i}\\
f_{i-1}\\
i^4\\
i^3\\
i^2\\
i\\
1\\
\end{matrix}
\right]
=
\left[
\begin{matrix}
1&2&1&4&6&4&1\\
1&0&0&0&0&0&0\\
0&0&1&4&6&4&1\\
0&0&0&1&3&3&1\\
0&0&0&0&1&2&1\\
0&0&0&0&0&1&1\\
0&0&0&0&0&0&1\\
\end{matrix}
\right]
*
\left[
\begin{matrix}
f_{i-1}\\
f_{i-2}\\
(i-1)^4\\
(i-1)^3\\
(i-1)^2\\
i-1\\
1\\
\end{matrix}
\right]
\]
B为\([f_2,f_1,2^4,2^3,2^2,2,1]^T\)于是\(f(n)=A^{n-2}*B\)的第一项。
有了递推关系,再用矩阵快速幂解决就好了。
代码
#include <cstdio>
#include <cstring>
#include <algorithm>
#define ll long long
#include <iostream>
using namespace std;
const ll mod=2147493647;
struct Mat{
int r,c;
ll a[10][10];
Mat(int _r,int _c){
r=_r;c=_c;
memset(a,0,sizeof a);
}
Mat operator *(const Mat &b)const{
Mat c(r,b.c);
for(int i=0;i<r;i++)
for(int j=0;j<b.c;j++)
for(int k=0;k<b.r;k++){
c.a[i][j]=(c.a[i][j]+a[i][k]*b.a[k][j]%mod)%mod;
}
return c;
}
}A(7,7),B(7,1);
Mat qpow(Mat a,int b){
Mat c(a.r,a.c);
for(int i=0;i<a.r;i++)c.a[i][i]=1;
while(b){
if(b&1)c=c*a;
b>>=1;
a=a*a;
}
return c;
}
int main() {
int at[10][10]={{1,2,1,4,6,4,1},
{1,0,0,0,0,0,0},
{0,0,1,4,6,4,1},
{0,0,0,1,3,3,1},
{0,0,0,0,1,2,1},
{0,0,0,0,0,1,1},
{0,0,0,0,0,0,1}};
for(int i=0;i<7;i++)for(int j=0;j<7;j++)A.a[i][j]=at[i][j];
int t,n,a,b;
cin>>t;
while(t--){
scanf("%d%d%d",&n,&a,&b);
B.a[0][0]=b;B.a[1][0]=a;
B.a[6][0]=1;
for(int i=5;i>1;i--)B.a[i][0]=B.a[i+1][0]*2;
if(n==1){
printf("%d\n",a);
}else if(n==2){
printf("%d\n",b);
}else{
Mat C=qpow(A,n-2)*B;
printf("%lld\n",C.a[0][0]);
}
}
return 0;
}
【HDU5950】Recursive sequence(矩阵快速幂)的更多相关文章
- HDU5950 Recursive sequence —— 矩阵快速幂
题目链接:https://vjudge.net/problem/HDU-5950 Recursive sequence Time Limit: 2000/1000 MS (Java/Others) ...
- HDU5950 Recursive sequence (矩阵快速幂加速递推) (2016ACM/ICPC亚洲赛区沈阳站 Problem C)
题目链接:传送门 题目: Recursive sequence Time Limit: / MS (Java/Others) Memory Limit: / K (Java/Others) Total ...
- HDU5950 Recursive sequence (矩阵快速幂)
Recursive sequence Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Other ...
- HDU 5950 - Recursive sequence - [矩阵快速幂加速递推][2016ACM/ICPC亚洲区沈阳站 Problem C]
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5950 Farmer John likes to play mathematics games with ...
- 5950 Recursive sequence (矩阵快速幂)
题意:递推公式 Fn = Fn-1 + 2 * Fn-2 + n*n,让求 Fn; 析:很明显的矩阵快速幂,因为这个很像Fibonacci数列,所以我们考虑是矩阵,然后我们进行推公式,因为这样我们是无 ...
- CF1106F Lunar New Year and a Recursive Sequence——矩阵快速幂&&bsgs
题意 设 $$f_i = \left\{\begin{matrix}1 , \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ i < k\\ ...
- hdu 5950 Recursive sequence 矩阵快速幂
Recursive sequence Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Other ...
- hdu-5667 Sequence(矩阵快速幂+费马小定理+快速幂)
题目链接: Sequence Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others) ...
- UVA - 10689 Yet another Number Sequence 矩阵快速幂
Yet another Number Sequence Let’s define another number sequence, given by the foll ...
- Yet Another Number Sequence——[矩阵快速幂]
Description Everyone knows what the Fibonacci sequence is. This sequence can be defined by the recur ...
随机推荐
- matplotlib 入门之Pyplot tutorial
文章目录 pyplot 介绍 修饰你的图案 格式字符串 [color][marker][line] Colors Markers Line Styles 利用关键字作图(大概是数据映射到属性吧) 传入 ...
- 安装pandas时出现环境错误
在安装pandas时出现Could not install packages due to an EnvironmentErrorConsider using the `--user` option ...
- Django的模板语言
Django模板系统 官方文档 常用语法 只需要记两种特殊符号: {{ }}和 {% %} 变量相关的用{{}},逻辑相关的用{%%}. 变量 {{ 变量名 }} 变量名由字母数字和下划线组成. ...
- 多线程系列之四:Guarded Suspension 模式
一,什么是Guarded Suspension模式如果执行现在的处理会造成问题,就让执行处理的线程等待.这种模式通过让线程等待来保证实例的安全性 二,实现一个简单的线程间通信的例子 一个线程(Clie ...
- 阿里云ECS服务器云监控(cloudmonitor)Go语言版本插件安装卸载与维护
云监控Go语言版本插件安装_主机监控_用户指南_云监控-阿里云https://help.aliyun.com/document_detail/97929.html 云监控cloudmonitor 1. ...
- LR 两种html与url录制
一直在使用LR,对于Html_based script和Url-based script 两种录制方式之间,要如何选择,仍是一知半解.最近测试时遇到同样的业务功能,两种录制方式的脚本,单次执行时间差别 ...
- <c:forEach varStatus="status">中 varStatus的作用
varStatus是<c:forEach>jstl循环标签的一个属性,varStatus属性. varStatus=“status”事实上定义了一个status名的对象作为varStatu ...
- Satis搭建composer私有库(自定义下载目录)
在我们的日常php开发中需要使用大量的第三方包和类库, 怎么管理是一个问题, 我们用的Yii2框架, 但是并没有把composer用起来, 由于最近更换为docker部署项目, 于是想起来用compo ...
- java依赖的斗争:依赖倒置、控制反转和依赖注入
控制反转(Inversion Of Controller)的一个著名的同义原则是由Robert C.Martin提出的依赖倒置原则(Dependency Inversion Principle),它的 ...
- [转帖]Centos7 yum安装Chrome浏览器
Centos7 yum安装Chrome浏览器 https://www.cnblogs.com/ianduin/p/8727333.html以及https://blog.csdn.net/libaine ...