Quad Tiling

poj 3420

题目大意

在一个4×n的棋盘上，用1×2的多米诺骨牌把他填满，问有多少种方法

输入样例

输出样例

1
11
95

数据范围

1⩽N⩽1091 \leqslant N \leqslant 10^91⩽N⩽109
0<M⩽1050 < M \leqslant 10^50<M⩽105

解题思路

对于放多米诺骨牌可以用状压DP
但直接状压会TLE
考虑矩阵乘法
可以先把状态之间的关系预处理出来
然后快速幂即可

代码

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define ll long long
#define wyc mod
using namespace std;
ll n, mod;
struct matrix
{ll n, m, a[20][20];matrix operator *(const matrix &b) const{matrix c;c.n = n;c.m = b.m;for (int i = 1; i <= c.n; ++i)for (int j = 1; j <= c.m; ++j)c.a[i][j] = 0;for (int i = 1; i <= c.n; ++i)for (int k = 1; k <= m; ++k)for (int j = 1; j <= c.m; ++j)c.a[i][j] = (c.a[i][j] + a[i][k] * b.a[k][j] % wyc) % wyc;return c;}
}A, B;
ll ggg(ll x, ll y, ll z, ll g)
{return x + (y<<1) + (z<<2) + (g<<3) + 1;
}
void pp(ll x, ll y, ll z, ll g)
{ll xx = x^1, yy = y^1, zz = z^1, gg = g^1;B.a[ggg(xx, yy, zz, gg)][ggg(x, y, z, g)] = 1;if (!xx && !yy){B.a[ggg(1, 1, zz, gg)][ggg(x, y, z, g)] = 1;if (!zz && !gg)B.a[ggg(1, 1, 1, 1)][ggg(x, y, z, g)] = 1;}if (!yy && !zz)B.a[ggg(xx, 1, 1, gg)][ggg(x, y, z, g)] = 1;if (!zz && !gg)B.a[ggg(xx, yy, 1, 1)][ggg(x, y, z, g)] = 1;return;
}
void Counting(ll x)
{while(x){if (x&1) A = A * B;B = B * B;x>>=1; }return;
}
int main()
{while(1){scanf("%lld%lld", &n, &mod);if (!n) break;for (int i = 1; i <= 16; ++i){A.a[1][i] = 0;for (int j = 1; j <= 16; ++j)B.a[i][j] = 0;}A.n = 1;A.m = B.n = B.m = 16;A.a[1][1] = 1;//初始状态A.a[1][4] = 1;A.a[1][7] = 1;A.a[1][13] = 1;A.a[1][16] = 1;for (int i = 0; i <= 1; ++i)for (int j = 0; j <= 1; ++j)for (int k = 0; k <= 1; ++k)for (int c = 0; c <= 1; ++c)pp(i, j, k, c);//预处理矩阵Counting(n - 1);//快速幂printf("%lld\n", A.a[1][16] % wyc);}return 0;
}

【矩阵乘法】Quad Tiling（poj 3420）相关推荐

线性代数三之状压DP的矩阵加速——Quad Tiling，Bus公交线路
状压与矩阵加速的藕断丝连 Quad Tiling description solution code [Hnoi2010]Bus 公交线路 description solution code Quad ...
POJ3420 Quad Tiling【矩阵快速幂】
Quad Tiling Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 5008 Accepted: 2269 Descripti ...
【矩阵乘法】Matrix Power Series（poj 3233）
Matrix Power Series poj 3233 题目大意给你一个矩阵A,让你求S=A+A2+A3+-+AkS = A + A^2 + A^3 + - + A^kS=A+A2+A3+-+Ak ...
POJ3420 Quad Tiling（模板+矩阵快速幂）
Quad Tiling Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 4107 Accepted: 1878 Descripti ...
poj 3233 矩阵乘法（分块矩阵）
POJ 3233 题解:Sn为所求矩阵, 则这样, 此题就变成了求矩阵幂和矩阵乘法, 分块矩阵乘法和普通矩阵一样的. code: /* adrui's submission Language : C ...
如何在CPU上优化GEMM矩阵乘法
如何在CPU上优化GEMM矩阵乘法 How to optimize GEMM on CPU (TL;DR) TVM 提供抽象接口,允许用户分别描述算法和算法的实现组织(所谓的调度).通常,在高性能调度 ...
HDU 2865 Birthday Toy [Polya 矩阵乘法]
传送门题意: 相邻珠子不能相同,旋转等价.$n$个珠子$k$中颜色,求方案数首先中间珠子$k$种选择,$k--$ 如果没有相邻不同的限制,就和$POJ\ 2154$一样了 $|C(f)|=k^{\ ...
2.6 矩阵乘法重要特例
重要特例 AmnA_{mn}Amn , BnpB_{np}Bnp 两个矩阵相乘,当 m,n,pm,n,pm,n,p 有的为 111 时,是极其重要的特例,具有重要意义. 定义列向量列数为 11 ...
信息学奥赛一本通提高篇第5章矩阵乘法
例1 矩阵AXB 信息学奥赛一本通(C++版)在线评测系统 [矩阵乘法]矩阵A×B_Uletay-CSDN博客矩阵乘法--矩阵A×B_vina的博客-CSDN博客一本通1641[例 1]矩阵 A× ...

【矩阵乘法】Quad Tiling（poj 3420）

Quad Tiling

poj 3420

题目大意

输入样例

输出样例

数据范围

解题思路

代码

【矩阵乘法】Quad Tiling（poj 3420）相关推荐

最新文章

热门文章