SSL_1384&&P2704&&ybtoj【动态规划】5章4题【炮兵阵地】

炮兵阵地

题目

P2704

---

解析

本来出题人的本意是写滚动数组，然而合法方案最多只有60个，于是我就没有写
调了一个晚上，问题居然是iostream::sync_with_stdio(0)和cin.tie(0)把输入字符搞错了
考虑像SSL_1383一样找出每行的合法方案

vector<int>s[110],w[110];
inline void dfs(int wcr,int aki,int oiev,int ery,int day)
{
	if(wcr==m){if((aki&ch[day])==aki)s[day].push_back(aki),w[day].push_back(oiev);return;}
	dfs(wcr+1,aki<<1,oiev,ery-1,day);
	if(ery<=0)dfs(wcr+1,aki<<1|1,oiev+1,2,day);
}

这里ch数组代表每行的地形，亦即

for(rr i=1;i<=n;++i)for(rr j=1;j<=m;++j)
 cin>>wcrakioi,ch[i]=ch[i]|((wcrakioi=='P')?(1<<(j-1)):0);

设dp[i][j][k]表示处理到第i行，当前行方案编号为j，上一行方案为k，则显然易见有状态转移方程

for(rr i=0;i<s[1].size();++i)dp[1][i][0]=w[1][i];//第一行
for(rr i=0;i<s[2].size();++i)
 for(rr j=0;j<s[1].size();++j)
  if((!(s[2][i]&s[1][j])))
   dp[2][i][j]=max(dp[2][i][j],dp[1][j][0]+w[2][i]);//第二行
for(rr i=3;i<=n;++i)//枚举行的编号
 for(rr j=0;j<s[i].size();++j)//枚举这一行
  for(rr k=0;k<s[i-1].size();++k)//枚举上一行
   if(!(s[i][j]&s[i-1][k]))//这两行不冲突
    for(rr l=0;l<s[i-2].size();++l)//枚举上两行
     if((!(s[i][j]&s[i-2][l]))&&(!(s[i-1][k]&s[i-2][l])))//上两行与上一行和这一行都不冲突
      dp[i][j][k]=max(dp[i][j][k],dp[i-1][k][l]+w[i][j]);//取max值

易得ans来源：

if(n!=1)for(rr i=0;i<s[n].size();++i)for(rr j=0;j<s[n-1].size();++j)ans=max(ans,dp[n][i][j]);
else for(rr i=0;i<s[1].size();++i)ans=max(ans,dp[1][i][0]);

于是我们就愉快的A掉了这题
代码压行稍显毒瘤，请谨慎食用

code:

#include<cstdio>
#include<vector>
#include<iostream>
#define rr register int
using namespace std;
int n,m,ch[110],dp[110][70][70],ans;
vector<int>s[110],w[110];
inline void dfs(int wcr,int aki,int oiev,int ery,int day)
{
	if(wcr==m){if((aki&ch[day])==aki)s[day].push_back(aki),w[day].push_back(oiev);return;}
	dfs(wcr+1,aki<<1,oiev,ery-1,day);
	if(ery<=0)dfs(wcr+1,aki<<1|1,oiev+1,2,day);
}
char wcrakioi;
int main()
{
	scanf("%d%d",&n,&m);
	for(rr i=1;i<=n;++i)for(rr j=1;j<=m;++j)cin>>wcrakioi,ch[i]=ch[i]|((wcrakioi=='P')?(1<<(j-1)):0);
	for(rr i=1;i<=n;++i)dfs(0,0,0,0,i);
	for(rr i=0;i<s[1].size();++i)dp[1][i][0]=w[1][i];
	for(rr i=0;i<s[2].size();++i)for(rr j=0;j<s[1].size();++j)if((!(s[2][i]&s[1][j])))dp[2][i][j]=max(dp[2][i][j],dp[1][j][0]+w[2][i]);
	for(rr i=3;i<=n;++i)for(rr j=0;j<s[i].size();++j)for(rr k=0;k<s[i-1].size();++k)if(!(s[i][j]&s[i-1][k]))for(rr l=0;l<s[i-2].size();++l)if((!(s[i][j]&s[i-2][l]))&&(!(s[i-1][k]&s[i-2][l])))dp[i][j][k]=max(dp[i][j][k],dp[i-1][k][l]+w[i][j]);
	if(n!=1)for(rr i=0;i<s[n].size();++i)for(rr j=0;j<s[n-1].size();++j)ans=max(ans,dp[n][i][j]);
	else for(rr i=0;i<s[1].size();++i)ans=max(ans,dp[1][i][0]);
	printf("%d",ans);
	return 0;
}