gym100543E Can't stop playing CERC2014

https://codeforces.com/gym/100543

虽然知道大概是DP,然而我写的假DP又T又WA又RE

对着官网上下载的标程思考了好久,太神了,楼教主诚不欺我,他博客里看到评价波兰选手有很强的灵气,做题非常艺术,这么神的DP现场能过15个队。。

都不像那种数位DP,wqs二分,决策单调性优化等,有时候是讨论情况麻烦,或者说是套路,然而这题是纯巧妙

很显然我们知道在一维版2048合并的过程中,必须满足他是一个A字型的,数字需要先上升再下降。

那么dp[i][j]表示已经用了1-i 这些数字,递减(不包括最高点)的右边的和为 j  

dp[i][j]=-1表示无法得到这种情况,dp[i][j]>=0表示没放a[i]这个数字的时候,也就是1-(i-1)的情况是dp[i][j],也就是从[i-1][dp[i][j]]这个转移过来的

因为是一段上升,一段下降的序列,而且相邻两位不能一样,所以上升的序列之和的二进制表示,每一位是否为1,就表示有没有这个数位,比如15来表示1 2 4 8,这种映射是唯一的,右边下降的序列之和(不包括顶点),同理

那么我们知道,如果a[i]要放在当前序列的左边,则要满足a[i]<=lowbit[sum-j],sum为1- (i-1)之和,而 j 枚举的不包括顶点的下降序列之和,那么sum-j就表示包括顶点的上升序列之和,lowbit[sum-j]就是最左边的这个数字,那么只有a[i]<=最左边的数字,才能放在左边,注意特判一下第一个数字也是可以放在左边的。

如果a[i]需要放在当前序列的右边,由于lowbit[0]=0,那么如果j=0的时候,也就是右侧没有比最高点小的数字时,可以直接放在右边,第二种就是a[i]<=lowbit[j],也就是最右边也就是最小的数字小时,可以放在最右边。

注意此时如果newj=a[i]+j,如果进位了导致highlowbit[newj]>=highlowbit[sum-j],说明右边的位数>=左边的位数了,那么最高点会变高,其实这种情况实际上就是a[i]==lowbit[j],然后导致最高位都进位了,那么此时右边的所有序列相当于合并到了最高点,由于我们的newj是不包含最高点的下降序列值和,所以此时newj=0,dp[i][newj]=j。

dp转移得过程中可能会存在覆盖的情况,但是没关系,因为最后转移链一定是一条可行解。

#include<bits/stdc++.h>
using namespace std;

const int maxl=1010;
const int all=1<<13;

int n,ans,tot;
int a[maxl],lowbit[all+10],highbit[all+10];
int dp[maxl][(1<<13)+5];
char s[maxl][(1<<13)+5];
char b[maxl];

inline void prework()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
		scanf("%d",&a[i]);
}

inline void mainwork()
{
	dp[0][0]=0;
	int sum=0,newj;
	for(int i=1;i<=n;i++)
	{
		for(int j=0;j<=sum+a[i];j++)
			dp[i][j]=-1;
		for(int j=0;j<=sum;j++)
		{
			if(dp[i-1][j]<0) continue;
			if(i==1 || a[i]<=lowbit[sum-j])
			{
				dp[i][j]=j;
				s[i][j]='l';
			}
			if(j==0 || a[i]<=lowbit[j])
			{
				newj=j+a[i];
				if(highbit[newj]>=highbit[sum-j])
					newj=0;
				dp[i][newj]=j;
				s[i][newj]='r';
			}
		}
		sum+=a[i];
	}
	if(dp[n][0]<0 || sum!=lowbit[sum])
	{
		ans=0;
		return;
	}
	ans=1;
	int nowj=0;
	for(int i=n;i>=1;i--)
	{
		b[i]=s[i][nowj];
		nowj=dp[i][nowj];
	}
}

inline void print()
{
	if(!ans)
		puts("no");
	else
	{
		for(int i=1;i<=n;i++)
			printf("%c",b[i]);
		puts("");
	}
}

int main()
{
	for(int i=1;i<=all;i++)
	{
		lowbit[i]=(i&1)?1:lowbit[i/2]*2;
		highbit[i]=(lowbit[i]==i)?i:highbit[i/2]*2;
	}
	int t;
	scanf("%d",&t);
	for(int i=1;i<=t;i++)
	{
		prework();
		mainwork();
		print();
	}
	return 0;
}

 

 

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值