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;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值