标签:
Time Limit: 2000MS | Memory Limit: 65536K | |
Total Submissions: 20940 | Accepted: 8104 |
Description
Input
Output
Sample Input
5 4 PHPP PPHH PPPP PHPP PHHP
Sample Output
6
思路:状压dp。按行进行递推,每一行的状态只与上一行和上上一行有关,而每一行的结果只是由上一行推出,满足无后效性,可dp,将每一行的状态用二进制整数存,筛掉明显不符合条件的,进行dp。
dp(i,j,k)表示第i行的决策做出后,第i行的状态为j,第i-1行的状态为k的炮兵数。
dp(i,j,k)=max(dp(i,j,k),dp(i-1,k,t)+sum[j]) (i,j,k状态存在且不互相冲突)
注意初始化筛掉不存在或不符合条件的状态
#include<iostream> #include<cstdio> #include<cstring> #include<cstdlib> #include<algorithm> using namespace std; const int maxn=2200; const int INF=(1<<29); int N,M; char s[maxn]; int Map[maxn]; int sum[maxn*300]; int dp[2][maxn][maxn]; int state[maxn*300],cnt; void InitSum() { memset(sum,0,sizeof(sum)); for(int i=0;i<cnt;i++){ for(int j=0;j<M;j++){ if(state[i]&(1<<j)) sum[i]++; } } } bool jud(int n) { if((n&(n<<1))||(n&(n<<2))) return 0; return 1; } void Init() { cnt=0; for(int i=0;i<(1<<M);i++){ if(jud(i)) state[cnt++]=i; } } int main() { cin>>N>>M; memset(Map,0,sizeof(Map)); for(int i=1;i<=N;i++){ scanf("%s",s); for(int j=0;j<strlen(s);j++){ if(s[j]==‘H‘) Map[i]|=(1<<j); } } Init(); InitSum(); memset(dp,0,sizeof(dp)); for(int i=1;i<=N;i++){ for(int j=0;j<cnt;j++){ if(state[j]&Map[i]) continue; for(int k=0;k<cnt;k++){ if(i==1) dp[i%2][j][k]=sum[j]; else if((state[k]&Map[i-1])||(state[j]&state[k])) continue; else if(i==2) dp[i%2][j][k]=max(dp[i%2][j][k],dp[(i+1)%2][k][0]+sum[j]); else{ for(int t=0;t<cnt;t++){ if(!(state[t]&Map[i-2])&&!(state[j]&state[t])&&!(state[k]&state[t])) dp[i%2][j][k]=max(dp[i%2][j][k],dp[(i+1)%2][k][t]+sum[j]); } } } } } int ans=0; for(int j=0;j<cnt;j++){ for(int k=0;k<cnt;k++){ if(dp[N%2][j][k]>ans) ans=dp[N%2][j][k]; } } cout<<ans<<endl; return 0; }
标签:
原文地址:http://www.cnblogs.com/--560/p/4504897.html