关闭

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

312人阅读 评论(0) 收藏 举报
分类:
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!!!


2
0
查看评论
发表评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场

poj 3070 Fibonacci 【矩阵快速幂 求第N个斐波那契数%1000】

Fibonacci Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 11123   Acc...
  • chenzhenyu123456
  • chenzhenyu123456
  • 2015-08-30 20:01
  • 2193

矩阵快速幂求斐波那契数列(初学整理)

参考文章:                  http://blog.csdn.net/u013795055/article/details/38599321       ...
  • NYIST_TC_LYQ
  • NYIST_TC_LYQ
  • 2016-10-31 13:27
  • 6697

矩阵快速幂 ——(递推表达式)

矩阵快速幂    首先知道矩阵       矩阵(Matrix)是一个按照长方阵列排列的复数或实数集合;      矩阵乘法: 定义:设A为    的矩阵,B为   ...
  • WR_technology
  • WR_technology
  • 2016-05-06 19:44
  • 1360

20161031的考试】模拟,暴力取模+矩阵快速幂,搜索脸的dp

首先……我要大喊三声flaze煞笔 T2……当时刚写完不知哪儿抽了手一抖一个shift+del删源码,解锁成就【mdzz】,然后重写的时候,果断脑残取模优化GG 嗯……T3莫名其妙WA...
  • Flaze_
  • Flaze_
  • 2016-11-01 16:35
  • 277

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

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

BZOJ 1009: [HNOI2008]GT考试(DP+KMP失配+矩阵快速幂)

1009: [HNOI2008]GT考试 Time Limit: 1 Sec Memory Limit: 162 MB Submit: 4181 Solved: 2552 [Submit][St...
  • Ever_glow
  • Ever_glow
  • 2017-10-20 16:52
  • 361

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

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

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

1009: [HNOI2008]GT考试 Time Limit: 1 Sec  Memory Limit: 162 MB Submit: 3932  Solved: 2398 [Submit][...
  • Jaihk662
  • Jaihk662
  • 2017-07-19 02:35
  • 279

BZOJ 1009 GT考试 KMP+DP+矩阵快速幂

KMP+DP+矩阵快速幂(详解KMP算法)
  • neighthorn
  • neighthorn
  • 2016-07-17 20:59
  • 264

[NOIP模拟题][LIS][同余最短路][DP][矩阵快速幂][容斥原理]

小L的二叉树 【题目描述】 勤奋又善于思考的小L接触了信息学竞赛,开始的学习十分顺利。但是,小L对数据结构的掌握实在十分渣渣。 所以,小L当时卡在了二叉树。 在计算机科学中,二叉树是...
  • yznoi2017
  • yznoi2017
  • 2016-11-15 19:58
  • 152
    个人资料
    • 访问:13886次
    • 积分:931
    • 等级:
    • 排名:千里之外
    • 原创:74篇
    • 转载:0篇
    • 译文:0篇
    • 评论:19条
    另一个我
    Free组合~~友链O(∩_∩)O
    文章分类