[数学-构造矩阵]NEFU 1113
然而以下居然不能继续推到sn的公式!!!
!
这道题考察的就是求随意数列的前n项和,在sn的递推公式不太明显的时候。用矩阵解决。
设矩阵A=。矩阵F0=
" />
那么设矩阵S=(A+A2+A3…. + An)*F0
终于答案就是矩阵S内两个元素之和。
那么怎么求A+A2+A3…. + An ?
能够继续构造例如以下的分块矩阵,当中 I 是单位矩阵
设R=。则有: R2=,R3=
能够发现右上角即为 I + A + A^2 + ... + A^n,多一个 I 后面给减掉就能够了
能够用高速幂求出R^n;
然而上面的方法对于此题仍然tle,看了标码发现。能通过推导进一步缩小矩阵的阶数,我这里的R是四阶。而标码里的运算仅仅有三阶。
继续思考:
看看能不能直接推导得到S的通项公式。看解说:
T[i] = dp[i][0]+dp[i][1]
=6*dp[i-1][1]+5*dp[i-1][0]
=6*T[i-1]-dp[i-1][0]
=6*T[i-1]-T[i-2]
依据S[i]=S[i-1]+T[i]能够计算出:
S[i]=S[i-1]+ 6*T[i-1]-T[i-2]
则有公式:
设R= ,搞定!
总结:矩阵的应用,细致学习上面的构造矩阵和推导过程。第一种构造分块矩阵的方法非常实用,它对sn公式不好直接构造矩阵的时候适用。
但假设像上面S[i]=S[i-1]+ 6*T[i-1]-T[i-2]这种公式能够推导出sn的递推矩阵。能够减少复杂度。
#include<stdio.h>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
using namespace std;
const long long mod = 1000000007;
struct Ma
{
long long m[3][3];
};
Ma operator * (Ma a,Ma b)
{
Ma c;
for(int i=0; i<3; i++)
for(int j=0; j<3; j++)
{
c.m[i][j] = 0;
for (int k = 0; k < 3; k++)
{
c.m[i][j]+=(a.m[i][k]*b.m[k][j]);
if(c.m[i][j]>=mod||c.m[i][j]<=-mod) c.m[i][j]%=mod;
//c.m[i][j]%=mod;
}
}
return c;
}
Ma mm[65];
long long cal(long long n)
{
int cur=0;
Ma ans= {1,6,-1,
0,6,-1,
0,1,0
};
while(n)
{
if(n&1)
{
ans = ans*mm[cur];
}
cur++;
n>>=1;
}
long long tmp=41*ans.m[0][0]%mod+35*ans.m[0][1]%mod+6*ans.m[0][2]%mod+mod;
tmp%=mod;
while(tmp<0) tmp+=mod;
printf("%lld\n",tmp);
return tmp;
}
//long long ans[10000005];
int main()
{
Ma tmp= {1,6,-1,
0,6,-1,
0,1,0
};
mm[0] = tmp;
for(int i=1; i<64; i++) mm[i]=mm[i-1]*mm[i-1];
long long n;
while(scanf("%lld",&n)!=EOF)
{
if(n==1) puts("6");
else if(n==2) puts("41");
else cal(n-3);
}
return 0;
}
[数学-构造矩阵]NEFU 1113的更多相关文章
- BZOJ-2326 数学作业 矩阵乘法快速幂+快速乘
2326: [HNOI2011]数学作业 Time Limit: 10 Sec Memory Limit: 128 MB Submit: 1564 Solved: 910 [Submit][Statu ...
- POJ 3233 Matrix Power Series(构造矩阵求等比)
Description Given a n × n matrix A and a positive integer k, find the sum S = A + A2 + A3 + … + Ak. ...
- Number Sequence(HDU 1005 构造矩阵 )
Number Sequence Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)T ...
- hdu 5015 233 Matrix(构造矩阵)
http://acm.hdu.edu.cn/showproblem.php?pid=5015 由于是个二维的递推式,当时没有想到能够这样构造矩阵.从列上看,当前这一列都是由前一列递推得到.依据这一点来 ...
- poj 3735 Training little cats(构造矩阵)
http://poj.org/problem?id=3735 大致题意: 有n仅仅猫,開始时每仅仅猫有花生0颗,现有一组操作,由以下三个中的k个操作组成: 1. g i 给i仅仅猫一颗花生米 2. e ...
- HDU 3306 Another kind of Fibonacci ---构造矩阵***
Another kind of Fibonacci Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Jav ...
- 【CodeForces】708 B. Recover the String 数学构造
[题目]B. Recover the String [题意]找到一个串s,满足其中子序列{0,0}{0,1}{1,0}{1,1}的数量分别满足给定的数a1~a4,或判断不存在.数字<=10^9, ...
- 构造矩阵解决这个问题 【nyoj299 Matrix Power Series】
矩阵的又一个新使用方法,构造矩阵进行高速幂. 比方拿 nyoj299 Matrix Power Series 来说 给出这样一个递推式: S = A + A2 + A3 + - + Ak. 让你求s. ...
- UVa 11149 Power of Matrix (矩阵快速幂,倍增法或构造矩阵)
题意:求A + A^2 + A^3 + ... + A^m. 析:主要是两种方式,第一种是倍增法,把A + A^2 + A^3 + ... + A^m,拆成两部分,一部分是(E + A^(m/2))( ...
随机推荐
- css ::selection 的妙用
1.选中页面文字和元素时的背景颜色 ::selection { background: #25b864; color: #fff; } 2.不能选择页面内容(但可以拖拽内容进行复制.挺好玩的) ::s ...
- 转 RabbitMQ 基础概念及 Spring 的配置和使用 推荐好文 举例讲解
从不知道到了解—RabbitMQ 基础概念及 Spring 的配置和使用 原理同上 请求地址:http://localhost:8080/home?type=3&routing_key=myO ...
- go包管理工具glide
一.命令 安装 $ go get github.com/Masterminds/glide $ go install github.com/Masterminds/glide 初始化项目并在当前项目下 ...
- Flowable 的event介绍
1 事件分为两种,一种是抛出:当流程执行到这时,抛出一个事件:另一种是捕获:当流程执行到这时,他就会等待一个事件的发生. 一个事件需要有事件定义,否则不会做任何“特殊”的事.对于一个流程实例,一个开始 ...
- Unity创建一个简易的弹簧(弹动)效果
参考文章:http://www.cnblogs.com/SkyD/archive/2008/09/05/1284778.html 主要依据胡克公式F=-k·x.这里k是倔度系数,可以理解为k值越大弹性 ...
- Atitti html5 h5 新特性attilax总结
Atitti html5 h5 新特性attilax总结 Attilax觉得不错的新特性 3.语义Header和Footer (The Semantic Header and Footer) 8.占位 ...
- 【Unity】第11章 第三人称角色控制器和碰撞体
分类:Unity.C#.VS2015 创建日期:2016-05-02 一.简介 第三人称视角控制器涉及的相关概念有: 1.刚体(Rigidbody). 2.碰撞体(Collider).包括球体碰撞体( ...
- java 字符串格式化
转:http://blog.sina.com.cn/s/blog_af26e3330101988v.html 1.对整数进行格式化:%[index$][标识][最小宽度]转换方式 我们可以看到,格式化 ...
- IE报错:The given path's format is not supported
在使用FileUpload控件进行上传EXCEL文件时,本地调试上传无问题,但是发布之后报地址无效错误 一.出现这个错误的主要原因是,在本地上传图片的时候HttpPostedFileBase对象里面保 ...
- # file Python-3.4.7.tar.xz Python-3.4.7.tar.xz: XZ compressed data
# file Python-3.4.7.tar.xz Python-3.4.7.tar.xz: XZ compressed data # xz -d Python-3.4.7.tar.xz # ls ...