poj 1185 炮兵阵地 3进制+轮廓线

炮兵阵地

Time Limit: 2000MS Memory Limit: 65536KTotal Submissions: 25634 Accepted: 9885

Description

司令部的将军们打算在N*M的网格地图上部署他们的炮兵部队。一个N*M的地图由N行M列组成，地图的每一格可能是山地（用"H" 表示），也可能是平原（用"P"表示），如下图。在每一格平原地形上最多可以布置一支炮兵部队（山地上不能够部署炮兵部队）；一支炮兵部队在地图上的攻击范围如图中黑色区域所示： 

如果在地图中的灰色所标识的平原上部署一支炮兵部队，则图中的黑色的网格表示它能够攻击到的区域：沿横向左右各两格，沿纵向上下各两格。图上其它白色网格均攻击不到。从图上可见炮兵的攻击范围不受地形的影响。 
现在，将军们规划如何部署炮兵部队，在防止误伤的前提下（保证任何两支炮兵部队之间不能互相攻击，即任何一支炮兵部队都不在其他支炮兵部队的攻击范围内），在整个地图区域内最多能够摆放多少我军的炮兵部队。 

Input

第一行包含两个由空格分割开的正整数，分别表示N和M； 
接下来的N行，每一行含有连续的M个字符('P'或者'H')，中间没有空格。按顺序表示地图中每一行的数据。N <= 100；M <= 10。

Output

仅一行，包含一个整数K，表示最多能摆放的炮兵部队的数量。

Sample Input

5 4PHPPPPHHPPPPPHPPPHHP

Sample Output

6

Source

Noi 01

[Submit]   [Go Back]   [Status]   [Discuss]

用三进制表示状态：0表示该位置没有铺设，1表示该位置的上方一格处有铺设，2表示该位置的上方相距为2的格子有铺设，且上方一格没有铺设。

每个位置只需要检查该位置上方距离2以内是否有铺设 还有 左方距离为2以内是否有铺设。

#include<cstdio>#include<string>#include<cstring>#include<iostream>#include<cmath>#include<algorithm>#include<vector>using namespace std;#define all(x) (x).begin(), (x).end()#define for0(a, n) for (int (a) = 0; (a) < (n); (a)++)#define for1(a, n) for (int (a) = 1; (a) <= (n); (a)++)#define mes(a,x,s)  memset(a,x,(s)*sizeof a[0])#define mem(a,x)  memset(a,x,sizeof a)#define ysk(x)  (1<<(x))#define CLEAR(x) { while( x>=sum[m])  x-=sum[m];   }typedef long long ll;typedef pair<int, int> pii;const int INF =0x3f3f3f3f;const int maxn= 100   ;const int maxm= 10   ;const int maxS=  59049   ;char a[maxn+3][maxm+3];int ed,n,m,now;int sum[maxm+3],dp[2][maxS+3];void pre(){    sum[0]=1;    for1(i,maxm)    {        sum[i]=sum[i-1]*3;    }}void update(int &x,int a){    if(a>x)  x=a;}bool no(int s,int pos){    if(s%3==2)    {        return false;    }    if(pos==1)    {        return true;    }    s/=3;    if(s%3==2)    {        return false;    }    return true;}void work(){    ed=sum[m]-1;    now=0;    mes(dp[now],-1,ed+1);    dp[now][0]=0;    int ans=0;    for0(i,n)    {        for0(j,m)        {             now^=1;             mes(dp[now],-1,ed+1);             for(int s=0;s<=ed;s++)  if(~dp[now^1][s])             {                  //不选                  int p=s/sum[m-1];                  p=p%3;                  int p2=p;                  if(p2)  p2--;                  int s2=s*3+p2;                  CLEAR(s2);                  update(dp[now][s2],dp[1^now][s]);                  update(ans,dp[now][s2]);                 //选                 if(a[i][j]=='P'&&!p&&(  !j|| j==1&&no(s,1)||no(s,2)    ))                 {                     s2=s*3+2;                     CLEAR(s2);                     update(dp[now][s2],dp[1^now][s]+1);                     update(ans,dp[now][s2]);                 }             }        }    }    cout<<ans<<endl;}int main(){   pre();   while(~scanf("%d%d",&n,&m))   {       for0(i,n)       {           for0(j,m)           {               scanf(" %c",&a[i][j]);           }       }       work();   }   return 0;}/*1 1P1 2PP1 3PPP1 4PPPP1 5PPPPP*/