HDU 5434 Peace small elephant 状压dp+矩阵快速幂
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了HDU 5434 Peace small elephant 状压dp+矩阵快速幂相关的知识,希望对你有一定的参考价值。
题目链接:
http://acm.hdu.edu.cn/showproblem.php?pid=5434
Peace small elephant
Accepts: 38
Time Limit: 10000/5000 MS (Java/Others)
Memory Limit: 65536/65536 K (Java/Others)
问题描述
小明很喜欢国际象棋,尤其喜欢国际象棋里面的大象(只要无阻挡能够斜着走任意格),但是他觉得国际象棋里的大象太凶残了,于是他想到了小象,
小象就没有大象那么凶残,它的攻击范围是它当前格子直角所斜对的格子。现在小明要在棋盘上放很多个小象,有趣的是,当两个小象所在格子有公共边时,
它们将合体变成合体象,多个小象满足条件也会合体,合体象的攻击范围也是它所覆盖格子区域直角所斜对的格子,现在要求任何一个象的攻击范围上是空的(即不摆放棋子),
小明的棋盘很特殊,有m*nm∗n个格子,求满足条件的摆放的方案数,由于方案数太大,需要对10000000071000000007取模。
下面给出几种形状下的象的攻击范围图,叉号表示攻击范围。
输入描述
输入有多组数据(最多55组),每组数据有两个整数n,mn,m含义如题目描述。
1 \leq m \leq 7,1 \leq n \leq 10000000001≤m≤7,1≤n≤1000000000
输出描述
每组数据对应输出一行包含一个整数,表示满足条件的摆放的方案数。
输入样例
1 1 2 3
输出样例
2 50
题解:
状压dp+矩阵快速幂。
做了好几天,终于敲出来了!。。
详细题解,明天补充,先睡为敬>
代码:
1 #include<iostream> 2 #include<cstdio> 3 #include<cstring> 4 using namespace std; 5 6 const int maxn = 155; 7 const int mod = 1000000007; 8 9 typedef long long LL; 10 11 struct Matrix { 12 int n, m; 13 int val[maxn][maxn]; 14 Matrix(int n,int m) :n(n),m(m) {} 15 Matrix() {} 16 void init(int n, int m) { this->n = n; this->m = m; } 17 friend Matrix operator * (const Matrix& mat1, const Matrix& mat2) { 18 Matrix ret(mat1.n, mat2.m); 19 for (int i = 0; i < ret.n; i++) { 20 for (int j = 0; j < ret.m; j++) { 21 ret.val[i][j] = 0; 22 for (int k = 0; k < mat1.m; k++) { 23 ret.val[i][j] += (LL)mat1.val[i][k] * mat2.val[k][j]%mod; 24 ret.val[i][j] %= mod; 25 } 26 } 27 } 28 return ret; 29 } 30 }; 31 32 void power(Matrix& mat, int n, Matrix& ans) { 33 while (n > 0) { 34 if (n % 2) ans = mat*ans; 35 mat = mat*mat; 36 n /= 2; 37 } 38 } 39 40 int _n, m; 41 42 bool isOk(int s1, int s2) { 43 for (int i = 0; i<m; i++) { 44 if ((s1&(1 << i)) && !(s2&(1 << i))) { 45 int j; 46 j = i - 1; 47 if (j >= 0) { 48 if ((s2&(1 << j)) && !(s1&(1 << j))) return false; 49 } 50 j = i + 1; 51 if (j<m) { 52 if ((s2&(1 << j)) && !(s1&(1 << j))) return false; 53 } 54 } 55 } 56 return true; 57 } 58 59 Matrix mat, ans; 60 61 void init() { 62 mat.init(1 << m, 1 << m); 63 for (int i = 0; i<mat.n; i++) { 64 for (int j = 0; j<mat.m; j++) { 65 if (isOk(i, j)) mat.val[i][j] = 1; 66 else mat.val[i][j] = 0; 67 } 68 } 69 70 ans.init(1 << m, 1); 71 for (int i = 0; i < ans.n; i++) ans.val[i][0] = 1; 72 } 73 74 int main() { 75 while (scanf("%d%d", &_n, &m) == 2 && _n) { 76 init(); 77 power(mat, _n - 1, ans); 78 int res = 0; 79 for (int i = 0; i < ans.n; i++) { 80 res += ans.val[i][0]; 81 res %= mod; 82 } 83 printf("%d\n", res); 84 } 85 return 0; 86 }
以上是关于HDU 5434 Peace small elephant 状压dp+矩阵快速幂的主要内容,如果未能解决你的问题,请参考以下文章
UVA326 LA5434 POJ1538 Extrapolation Using a Difference Table二项式