poj1185 状态压缩的dp

NOI炮兵阵地，经典题！

【题意】略

【解答】

我一开始的思路是基于三进制的，0表示这个方格不受控制，1表示这个方格收到距离为1的方格的控制，2表示这个方格收到距离为2的方格的控制。后来发现有点烦！所以换了一个。令f[i][j][k]表示第i行，i-1行状态为k，i-2行的状态为j。j，k均为01串，1表示有炮。

则状态转移就很方便了，这里留作思考。

【优化】

上述算法的时间复杂度还是很高的，可以发现，满足条件的01串远远小于1024种（只有60个），所以可以预处理一下。这样就快多了。

还可以使用滚动数组优化空间。

但时间还是有点紧，建议程序写的细致些，控制常数。

#include <iostream>
using namespace std;

int f[105][62][62],b[1050];
bool a[105][12];
int i,j,k,t,l,m,n,tot=0,ans;

int isok(int row,int x,int y,int z)
{
    int i,ans=0;
    x=b[x];y=b[y];z=b[z];
    for (i=0;i<m;i++)
        if (z>>i&1)
        {
           if (!a[row][i] || x>>i&1 || y>>i&1) 
              return -1;
           ans++;
        }
    return ans;
}
    
int main()
{
    freopen("pin.txt","r",stdin);
    freopen("pou.txt","w",stdout);
    cin >> n >> m;
    getchar();
    for (i=0;i<n;i++)
    {
        for (j=0;j<m;j++)
            a[i][j]=(getchar()=='P');
        getchar();
    }
    for (i=0;i<(1<<m);i++)
    {
        int v[12]={0};
        for (j=0;j<m;j++)
        if ((i>>j)&1)
           for ((j-2>=0)?k=j-2:k=0;k<=j+2 && k<m;k++)
               v[k]++;
        bool ok=1;
        for (j=0;j<m;j++)
            if (i>>j&1 && v[j]>1) ok=0;
        if (!ok) continue;
        b[++tot]=i;
    }
    for (i=0;i<=n;i++)
        for (j=1;j<=tot;j++)
            for (k=1;k<=tot;k++)
                f[i][j][k]=-1;
    f[0][1][1]=0;        
    for (i=0;i<n;i++)
        for (j=1;j<=tot;j++)
            for (k=1;k<=tot;k++)
            if (f[i][j][k]!=-1)
                for (l=1;l<=tot;l++)
                {
                    t=isok(i,j,k,l);
                    if (t>=0 && f[i][j][k]+t>f[i+1][k][l])
                       f[i+1][k][l]=t+f[i][j][k];
                }
    ans=0;
    for (i=1;i<=tot;i++)
        for (j=1;j<=tot;j++)
            if (f[n][i][j]>ans)
               ans=f[n][i][j];
    cout << ans;
    return 0;
}