POJ 3420 Quad Tiling 状压DP+矩阵快速幂

26 篇文章 1 订阅
7 篇文章 0 订阅

链接:http://poj.org/problem?id=3420

题意:给一个4*N(1 ≤ N ≤ 1e9)的矩形空间,并且给不限块数的1*2的多米诺骨牌,问是由多少种方式能把这个矩形空间填满。

思路:看到这种问题果断想到状压,虽然是在看矩阵的时候看到的这道题。dp[i][j]表示在第i行状态为j的情况下的填满方式数,j的二进制表示中0表示对应位置上一行的骨牌是竖放,或者对应位置的骨牌是横放,1则表示该行该位置的骨牌是竖放。由于N最大1e9所以O(n)的DP绝对超时,用矩阵快速幂来加速DP递推。

代码:

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <map>
#include <cstdlib>
#include <queue>
#include <stack>
#include <vector>
#include <ctype.h>
#include <algorithm>
#include <string>
#include <set>
#define PI acos(-1.0)
#define INF 0x7fffffff
#define eps 1e-8
typedef long long LL;
typedef unsigned long long ULL;
using namespace std;
#define maxn 20
#define maxm 20
using namespace std;
int a,MOD;
struct Matrix
{
    int n,m;
    LL a[maxn][maxm];
    void clear()
    {
        n=m=0;
        memset(a,0,sizeof(a));
    }
    Matrix operator +(const Matrix &b) const
    {
        Matrix tmp;
        tmp.n=n;
        tmp.m=m;
        for(int i=0; i<n; i++)
            for(int j=0; j<m; j++)
                tmp.a[i][j]=a[i][j]+b.a[i][j];
        return tmp;
    }
    Matrix operator -(const Matrix &b) const
    {
        Matrix tmp;
        tmp.n=n;
        tmp.m=m;
        for(int i=0; i<n; i++)
            for(int j=0; j<m; j++)
                tmp.a[i][j]=a[i][j]-b.a[i][j];
        return tmp;
    }
    Matrix operator *(const Matrix &b) const
    {
        Matrix tmp;
        tmp.clear();
        tmp.n=n;
        tmp.m=b.m;
        for(int i=0; i<n; i++)
            for(int j=0; j<b.m; j++)
                for(int k=0; k<m; k++)
                    tmp.a[i][j]=(tmp.a[i][j]+(a[i][k]*b.a[k][j])%MOD)%MOD;
        return tmp;
    }
};
Matrix M_quick_pow(Matrix &m,int k)
{
    Matrix tmp;
    tmp.n=m.n;
    tmp.m=m.m;
    for(int i=0; i<tmp.n; i++)
        for(int j=0; j<tmp.n; j++)
            if(i==j)
                tmp.a[i][j]=1;
            else tmp.a[i][j]=0;
    while(k)
    {
        if(k&1)
            tmp=tmp*m;
        k>>=1;
        m=m*m;
    }
    return tmp;
}
int flag (int i,int j)
{
    if(((i%3==0)&&(j%3==0)&&(i!=6)&&(j!=6))||((i==6)&&(j==9))||((i==9)&&(j==6)))
    {
        if((i&j)==0)
        return 1;
    }
    return 0;
}
int main()
{
    Matrix M,N,ans;
    while(scanf("%d%d",&a,&MOD)&&a&&MOD)
    {
        memset(N.a,0,sizeof(N.a));
        memset(M.a,0,sizeof(M.a));
        for(int i=0; i<16; i++)
        {
            for(int j=0; j<16; j++)
            {
                if(flag(i,j))
                    M.a[i][j]=1;
            }
        }
        M.m=M.n=16;
        ans=M_quick_pow(M,a);
        N.m=1;
        N.n=16;
        N.a[0][0]=1;
        ans=ans*N;
        printf("%lld\n",ans.a[0][0]%MOD);
    }
    return 0;
}

  • 2
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值