【poj3420】 Quad Tiling
生活随笔
收集整理的這篇文章主要介紹了
【poj3420】 Quad Tiling
小編覺(jué)得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.
http://poj.org/problem?id=3420?(題目鏈接)
題意
給出$n*m$的網(wǎng)格,用$1*2$的方塊覆蓋有多少種方案。
Solution
數(shù)據(jù)很大,不能直接搞了,我們矩乘一下。0表示已放置,1表示未放置。dfs跑出一個(gè)$16*16$的轉(zhuǎn)移矩陣,然后矩乘,最后輸出$ans[0][0]$就可以了。
代碼
// poj3420
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<cmath>
#define LL long long
#define HAS 4001
#define inf 2147483640
#define Pi acos(-1.0)
#define free(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
using namespace std;int n,m;
LL a[16][16],tmp[16][16],t[16][16],ans[16][16];void dfs(int x,int now,int pre) {if (x==4) {++a[pre][now];return;}dfs(x+1,now<<1|1,pre<<1);dfs(x+1,now<<1,pre<<1|1);if (x<3) dfs(x+2,now<<2,pre<<2);
}
void power() {for (int i=0;i<16;i++) {for (int j=0;j<16;j++) t[i][j]=a[i][j],ans[i][j]=0;ans[i][i]=1;}while (n) {if (n&1) {for (int i=0;i<16;i++)for (int j=0;j<16;j++) {tmp[i][j]=0;for (int k=0;k<16;k++) (tmp[i][j]+=ans[i][k]*t[k][j]%m)%=m;}for (int i=0;i<16;i++)for (int j=0;j<16;j++) ans[i][j]=tmp[i][j];}n>>=1;for (int i=0;i<16;i++)for (int j=0;j<16;j++) {tmp[i][j]=0;for (int k=0;k<16;k++) (tmp[i][j]+=t[i][k]*t[k][j]%m)%=m;}for (int i=0;i<16;i++)for (int j=0;j<16;j++) t[i][j]=tmp[i][j];}
}
int main() {dfs(0,0,0);while (scanf("%d%d",&n,&m)!=EOF && n && m) {power();printf("%lld\n",ans[0][0]);}return 0;
}
?
轉(zhuǎn)載于:https://www.cnblogs.com/MashiroSky/p/6417325.html
總結(jié)
以上是生活随笔為你收集整理的【poj3420】 Quad Tiling的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇: JS实战篇
- 下一篇: linux下ndk编译命令行程序及配置