poj 2411 zoj1100

参考周伟的论文：http://download.csdn.net/detail/hearthougan/7045453

第i行的铺设方法，例如(i, j)位置的铺设受到(i-1,j), (i, j-1)影响，因此(i, j)的铺设可能影响到(i-1, j+1), (i, j+1)位置的铺设.首先交代几个变量的含义：

j, curstatus, prestatus,分别表示当前为第j列，当前第i行的状态为curstatus，前一行i-1的状态为prestatus;

每个位置有三种铺法

(1)竖直放置

(2)横向放置

(3)不放

则有：

（1）：j = j+1, curstatus<<1|1, prestatus<<1

（2）：j = j+2, curstatus<<2|3,prestatus<<2|3

（3）：j = j+1,curstatus<<1, prestatus<<1|1

 

curstatus<<1|1即 为把curstatus的二进制表示 后面加上一个1，对于本题来说就是(j+1)列 上放置，
prestatus<<1即为把prestatus的二进制表示后面加上一个0，对于本题来说就是(j+1)列上不放置。
至于curstatus、prestatus的变化，我是纠结了好久啊，不明其因！后来看了很多别人的博客，认为比较合理的应该是：curstatus对应于本行的状态，prestatus对应于上一行的状态，能竖直放置意味着上一行的(j+1)列是空着的，因此此时上一行的状态为prestatus<<1，同时竖置放置了之后，则本行(j+1)列放置了东西，状态于是变为curstatus<<1|1;

注意：

因为算法 复杂度是  n * (m^4) 所以当 m > n的时候，我们交换他们这样适当降低复杂度。
   另外数据到后面比较大，所以使用 long long(不要用__int64!)

下面是不用滚动数组的代码：

 

#include <iostream>
#include <cstring>
#include <cstdio>

using namespace std;

typedef long long LL;
LL dp[15][3000];
int i, n, m;

void DFS(int j, int curstatus, int prestatus)
{
	if(j >= m)
	{
		dp[i][curstatus] += dp[i-1][prestatus];
		return ;
	}

	if(j+1 <= m)
	{
		DFS(j+1, curstatus<<1|1, prestatus<<1);
		DFS(j+1, curstatus<<1, prestatus<<1|1);
	}
	if(j+2 <= m)
		DFS(j+2, curstatus<<2|3, prestatus<<2|3);
}

int main()
{
	while (~scanf("%d %d", &n, &m) && (m+n != 0))
	{
		if(m*n % 2)
		{
			printf("0\n");
			continue ;
		}
		memset(dp, 0, sizeof(dp));

		dp[0][(1<<m) - 1] = 1;
		for (i = 1; i <= n; ++i)
		{
			DFS(0, 0, 0);
		}
		printf("%lld\n", dp[n][(1<<m)-1]);
	}
	return 0;
}

滚动数组的代码：

 

 

#include <iostream>
#include <cstring>
#include <cstdio>

using namespace std;

typedef long long LL;

LL dp[2][3000];
int n, m;
int k;

void Pre(int j, int curstatus)
{
    if(j >= m)
    {
        dp[0][curstatus]++;
        return ;
    }
    if(j+1 <= m)
        Pre(j+1, curstatus<<1);
    if(j+2 <= m)
        Pre(j+2, curstatus<<2|3);
}
void DFS(int j, int curstatus, int prestatus)
{
	if(j == m)
	{
		dp[k][curstatus] += dp[(k+1)%2][prestatus];
		return ;
	}

	if(j+1 <= m)
	{
		DFS(j+1, curstatus<<1|1, prestatus<<1);
		DFS(j+1, curstatus<<1, prestatus<<1|1);
	}

	if(j+2 <= m)
		DFS(j+2, (curstatus<<2)|3, (prestatus<<2)|3);

}

int main()
{
	while(~scanf("%d %d", &n, &m) && (n + m != 0))
	{
		if(n*m % 2)
		{
			printf("0\n");
			continue ;
		}
		if(n < m)
		{
			n ^= m;
			m ^= n;
			n ^= m;
		}

		memset(dp, 0, sizeof(dp));
		//dp[0][(1 << m) - 1] = 1;
		Pre(0, 0);
		for (int i = 2; i <= n; ++i)
		{
		    k = (i+1)%2;
			DFS(0, 0, 0);
			memset(dp[(k+1)%2], 0, sizeof(dp[0]));
		}
		printf("%lld\n",dp[(n+1)%2][(1 << m) - 1]);
	}
	return 0;
}