POJ 3254，POJ 1185（状态压缩DP）

题目链接： POJ-3254 Corn Fields

思路：

我们可以用二进制枚举第一行状态，用1代表在这种植，0代表不种植，第一行确定了，第二行就可以根据第一行不断累加。dp[row][state]代表第row行用state状态时的方法数。如果这个state合法，则dp[row][state]=dp[row-1][state1]+dp[row-1][state2]+…

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
const int mod=100000000;
int tot,state[1001],n,m,dp[13][1001],row[13];
void unit()
{
    int i;
    tot=0;
    for(i=0;i<(1<<m);i++)
    {
        if(i&(i<<1))
            continue;
        state[tot++]=i;
    }
}
int main()
{
    int i,j,k,a;
    cin>>n>>m;
    memset(row,0,sizeof(row));
    memset(dp,0,sizeof(dp));
    unit();
    for(i=1;i<=n;i++)
        for(j=1;j<=m;j++)
    {
        cin>>a;
        if(a==1)
            row[i]+=1<<(m-j);
    }
    dp[0][0]=1;
    for(i=1;i<=n;i++)
    {
        for(j=0;j<tot;j++)          //列举本行状态。
        {
            if((row[i]&state[j])!=state[j])
                continue;
            for(k=0;k<tot;k++)              //列举上一行状态
            {
                if(state[k]&state[j])continue;
                dp[i][j]=(dp[i][j]+dp[i-1][k])%mod;
            }
        }
    }
    int ans=0;
    for(i=0;i<tot;i++)
    {
        ans+=dp[n][i];
        ans=ans%mod;
    }
    cout<<ans<<endl;

}

题目链接： POJ-1185 炮兵阵地

思路：

上面的题是求方法总数，这个题是求最大炮兵数目，也就是求某一种方法中的最大炮兵放置数目，因此我们还需要一个函数用来求炮兵放置的数目，也就是getx()函数。用row[i]存储第i行原始状态，P平原代表1,然后用二进制不断枚举第一行各种状态，并存储起来，因为这个题是牵涉到影响上下两行，所以需要三维dp,dp[row][i][j]=max(dp[row-1][j][k1],dp[row-1][j][k2]…)+number[i].
dp[row][i][j]代表第row行的状态为i且 第row-1行状态为j时的最大炮兵放置数目。

#include <bits/stdc++.h>
using namespace std;
int tot,n,m,state[1001],number[1001],dp[101][101][101],row[101];
int getx(int x)
{
    int num=0,i;
    for(i=0;i<15;i++)
    {
        if(x & (1<<i))
            num++;
    }
    return num;
}
void unit()
{
    int i;
    tot=0;
    for(i=0;i<(1<<m);i++)
    {
        if( !(i&(i<<1)) && !(i&(i<<2)))
        {
        state[tot]=i;                  //state数组用来存储各个状态
        number[tot++]=getx(i) ;          //number数组用来存储该状态下有几座炮兵部队
        }
    }
}
int main()
{
    int t,i,j,k;
    char c;
    memset(row,0,sizeof(row));
    memset(dp,0,sizeof(dp));
    cin>>n>>m;
    for(i=1;i<=n;i++)
    for(j=1;j<=m;j++)
      {
        cin>>c;
        if(c=='P')
        row[i]+=1<<(m-j);             // row[]用来存储每一行的原始状态。
      }
    unit();
    for(i=0;i<tot;i++)                     //dp初始化处理
        {
        if((state[i]&row[1])!=state[i])continue;
        dp[1][i][0]=number[i];
        }
    for(i=2;i<=n;i++)
        {
        for(j=0;j<tot;j++)               //枚举本行状态。
           {
        if((row[i]&state[j])!=state[j])continue;
             for(k=0;k<tot;k++)               //枚举上一行状态
             {
               if((state[k]&row[i-1])!=state[k])continue;
               if(state[k]&state[j])continue;
               for(t=0;t<tot;t++)             //枚举上上行状态
               {
                  if(state[t]&state[j])continue;
                  if(state[t]&state[k])continue;
                  dp[i][j][k]=max(dp[i][j][k],dp[i-1][k][t]+number[j]);
               }
             }
           }
        }
    int sum=0;
    for(i=0;i<tot;i++)
        for(j=0;j<tot;j++)
        {
        sum=max(sum,dp[n][i][j]);
        }
    cout<<sum<<endl;
}