鱼头最近查刷题有点紧,事情又多,红了。不得已敲了几发水题,罪恶感爆棚。
不过敲水题过程中发现了这样一个题:
POJ2663:问3*n的方格用1*2的domino填充有多少种方法。
瞎几把YY了之后得出递推公式如下
a[i]=3*a[i-1]+2*b[i-1];
b[i]=a[i-1]+b[i-1];
其中a[i]是n=2*i的答案啦(显然n为奇数的结果为0)
b[i]是这种图形的答案
。。。111
。。。111
。。。011
然后就xjb搞完了,代码:
#include <cstdio>
#include <cstring>
using namespace std;
typedef long long ll;
ll a[30];
ll b[30];
int main()
{
memset(a,0,sizeof(a));
memset(b,0,sizeof(b));
a[1]=3;b[1]=1;
for(int i=2;i<=20;i++)
{
a[i]=3*a[i-1]+2*b[i-1];
b[i]=a[i-1]+b[i-1];
}
int n;
while(scanf("%d",&n),n!=-1)
{
if(!n){printf("1\n");continue;}
else if(n%2){printf("0\n");continue;}
else printf("%lld\n",a[n/2]);
}
return 0;
}
然后考虑m*n的情况,找到POJ2411,记起来宝哥报告讲的状压dp什么的,然而当时只是理论ac了,遂开敲。
题意:就把上面的3换成m,还是问domino填充种数。
做法就是把每行保存成一个2进制数,0表示空,1表示放,注意:竖着放的第一行先空着,视为0。
举个例子:
1–11—
1–11–1
–1—–1
–1——
对应
01100111
11111110
11011111
11111111
状态转移就是dp[i][j]=sum(dp[i-1][j1]),其中j1为与j不冲突的状态。
因为懒,没高兴想通过j构造j1,直接枚举所有j1判断是否和j冲突了,所以1300ms过的(菜如我);
代码:
#include <cstdio>
#include <cstring>
using namespace std;
typedef long long ll;
int m,n;
ll dp[12][4000];
bool judge(int p,int q)
{
int pos=0;
for(int i=0;i<n;)
{
if(p&(1<<i))
{
if(q&(1<<i))
if(i==n-1||!(p&(1<<(i+1)))||!(q&(1<<(i+1))))return false;
else i+=2;
else i++;
}
else
{
if(q&(1<<i))i++;
else return false;
}
}
return true;
}
ll dfs(int f,int s)
{
if(dp[f][s]!=-1)return dp[f][s];
ll ans=0;
for(int i=0;i<(1<<n);i++)
if(judge(s,i))ans+=dfs(f-1,i);
return dp[f][s]=ans;
}
bool jj1(int num)
{
for(int i=0;i<n;)
{
if(num&(1<<i))
if(i==n-1||!(num&(1<<(i+1))))return false;
else i+=2;
else i++;
}
return true;
}
int main()
{
while(~scanf("%d %d",&m,&n))
{
if(!m&&!n)break;
memset(dp,-1,sizeof(dp));
for(int i=0;i<(1<<n);i++)
if(jj1(i))dp[1][i]=1;
else dp[1][i]=0;
ll ans=dfs(m,(1<<n)-1);
printf("%lld\n",ans);
//printf("%lld\n",dp[2][1]);
}
return 0;
}
然后是POJ1185,就是那个炮兵阵地的题,题是中文的,略题意。
与上题不同的是要保存两行的状态,dp[i][now][past]=max(dp[i-1][past][pastspast]),其中pastspast和now和past和地图都不冲突,bingo!
因为每行每个之间至少距离2,所以根本没有1024种,只有六十几吧,先预处理一下就ok了。
还要讨论只有一行的状况。
代码:
#include <cstdio>
#include <cstring>
using namespace std;
typedef long long ll;
int m,n;
ll dp[12][4000];
bool judge(int p,int q)
{
int pos=0;
for(int i=0;i<n;)
{
if(p&(1<<i))
{
if(q&(1<<i))
if(i==n-1||!(p&(1<<(i+1)))||!(q&(1<<(i+1))))return false;
else i+=2;
else i++;
}
else
{
if(q&(1<<i))i++;
else return false;
}
}
return true;
}
ll dfs(int f,int s)
{
if(dp[f][s]!=-1)return dp[f][s];
ll ans=0;
for(int i=0;i<(1<<n);i++)
if(judge(s,i))ans+=dfs(f-1,i);
return dp[f][s]=ans;
}
bool jj1(int num)
{
for(int i=0;i<n;)
{
if(num&(1<<i))
if(i==n-1||!(num&(1<<(i+1))))return false;
else i+=2;
else i++;
}
return true;
}
int main()
{
while(~scanf("%d %d",&m,&n))
{
if(!m&&!n)break;
memset(dp,-1,sizeof(dp));
for(int i=0;i<(1<<n);i++)
if(jj1(i))dp[1][i]=1;
else dp[1][i]=0;
ll ans=dfs(m,(1<<n)-1);
printf("%lld\n",ans);
//printf("%lld\n",dp[2][1]);
}
return 0;
}