8.28 考试第二题(DP + 矩阵快速幂)

原创 2016年08月28日 20:47:05
2.方块游戏(block.c/cpp/pas)(标准题目见段大神的---> http://blog.csdn.net/dy0607/article/details/52347717)
【问题描述】
到了游乐场后,ZQ 为了感谢你,承诺免费请你打俄罗斯方块游戏,在经过你细致的观
察后,发现这个游戏只有如下两种方块:
当然,这两种形状的方块上的每一小块是一样大的。
然而,好动脑筋的 ZQ 想知道在宽为 n,高为 2 的界面内,有多少种消除方法。
例如,当 n=2 时,只有这么 5 种消除方法:
于是作为 ZQ 的伙伴,请你计算出消除方块的方案数 s。
【输入格式】
输入文件 block.in 第一行一个正整数 n。
【输出格式】
输出文件 block.out 一行一个整数 s,所求的方案数对 1000000007 取模后的结果。
【样例输入】
2
【样例输出】
5
【数据规模】
对于 30%的数据,保证 0≤n≤10000。
对于 50%的数据,保证 0≤n≤100000000。

显然本题很容易看出递推式,f[x] =  f[x-1] + 4*f[x-2] + 2*f[x-3]


但是呢,即便用滚动数组,我们也无法完成10^18数据规模的计算。。。

显然。。。我们用矩阵快速幂进行优化

f[x-1] f[x-2] f[x-3]   ---->    f[x]   f[x-1]   f[x-2]

  0         0         0                   0         0        0

  0         0         0                   0         0        0

发现:要用f[x-3], f[x-2], f[x-1]递推f[x], 我们把幂矩阵设成这个样子:

1 1 0

4 0 1

2 0 0

第1列的三个数即为要乘的数字,用来递推出f[x]的,另外的两个1,是用来平移f[x-1]与f[x-2]到前面的,可以自己乘一下试一试。。。


矩阵快速幂我还要多加练习,今天知道是用它却不知道怎么写,

代码上!


#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>

using namespace std;

#define MOD 1000000007
#define LL unsigned long long

struct Matrix{
	LL m[3][3];
	
	Matrix(int d){
		memset(m, 0, sizeof(m));
		if(d == 1){
			m[0][0] = m[0][1] = m[1][2] = 1;
			m[1][0] = 4; m[2][0] = 2;
		}
		if(d == 2){
			m[0][0] = m[1][1] = m[2][2] = 1;
		}
	}
	
	Matrix operator *(const Matrix &x)const{
		Matrix ret(3);
		for(int i = 0; i < 3; i++)
		  for(int j = 0; j < 3; j++)
		    for(int k = 0; k < 3; k++) ret.m[i][j] = (ret.m[i][j] + (m[i][k] * x.m[k][j]) % MOD) % MOD;
		
		return ret;
	}
};

Matrix power(Matrix a, long long num){
	Matrix ret(2);
	if(num == 1) return a;
	
	while(num){
		if(num & 1) ret = ret * a;
		a = a*a;
		num >>= 1;
	}
	
	return ret;
}

int main(){
	freopen("block.in", "r", stdin);
	freopen("block.out", "w", stdout);
	
	long long n;
	scanf("%lld", &n);
	
	if(n == 0 || n == 1){
		printf("1\n");
		return 0;
	}
	if(n == 2){
		printf("5\n");
		return 0;
	}
	
	Matrix g = 1;
	Matrix ans = power(g, n-2);
	
	LL P = ans.m[0][0]*5 + ans.m[1][0] + ans.m[2][0];
	printf("%llu\n", P % MOD);
	
	return 0;
}

更多解释见段大神和松AK的矩阵快速幂加递推。。。

(松AK: blog.csdn.net/frods)

(段大神: blog.csdn.net/dy0607)



I believe I can fly!!!


版权声明:支持的都转走!!!

相关文章推荐

[bzoj 1009] [HNOI2008]GT考试:DP,单串AC自动机,矩阵快速幂

题意:给一个长为M(M<=20)的十进制数字串,求有多少个长为N(N<=10^9)的十进制数字串不包含它,结果对K(K<=10^3)取模。 在黄学长博客的KMP分类中找到这题。并没产生什么好思路……准...

bzoj1009[GT考试] dp+矩阵快速幂

http://www.lydsy.com/JudgeOnline/problem.php?id=1009 题意:构造一个字符串(数字0-9)不含一个子串的个数 我这么蒟蒻肿么可能想到是dp嘛~...

bzoj 1009: [HNOI2008]GT考试(dp+kmp+矩阵快速幂)

1009: [HNOI2008]GT考试 Time Limit: 1 Sec  Memory Limit: 162 MB Submit: 3932  Solved: 2398 [Submit][...

Codeforces 821E Okabe and El Psy Kongroo【Dp+矩阵快速幂】套路题

E. Okabe and El Psy Kongroo time limit per test 2 seconds memory limit per test 256 megabyte...

【POJ 3420】Quad Tiling(dp|递推 +矩阵快速幂)

秋意渐浓,阳光西斜,遍地金黄

poj 2778 AC自动机+DP+矩阵快速幂

#include #include #include #include using namespace std; typedef __int64 type; const int kind=4; //...

矩阵快速幂DP Darth Vader and Tree : CodeForces - 514E

题目:CodeForces - 514Etime limit per test2 seconds memory limit per test256 megabytes When Darth Vad...

HDU 5318 (dp+矩阵快速幂优化)

The Goddess Of The Moon Time Limit: 6000/3000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java...

100道动态规划——28 POJ 2778 DNA Sequence AC自动机+DP+矩阵快速幂 很不错嘛

一眼看出是AC自动机上的DP,因为数据量的问题显然不能记忆化搜索或者是递推,然后就不知道怎么做了         天知道还要联系到矩阵快速幂,图论的东西好久都没写过了=_=         定义一...

POJ 3420 Quad Tiling(状压DP 用矩阵快速幂优化)

如果数据量不那么大那么这题就是一个裸的状压DP。
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)