[AGC028-E][树状数组]High Elements

本文介绍了一种解决字典序最小问题的方法,通过按位确定思想和分析序列的性质,利用前缀最大值的概念,确保序列拆分后的合法性。文章详细解释了如何通过枚举和数学方程确定序列的拆分方式,并使用了BIT树进行最大值查询,以实现高效处理。

翻译

在这里插入图片描述

题解

菜逼选手又来报到啦!
对于字典序最小的问题,我们显然是用按位确定的思想
定义几个变量方便使用
cnt0cnt0cnt0表示AAA序列当前有多少个前缀最大值,cnt1cnt1cnt1表示BBB序列当前有多少个前缀最大值
mx0mx0mx0表示AAA序列当前的最大值,mx1mx1mx1表示BBB序列当前的最大值
对于第iii位的确定工作,先分析能对序列大小做出贡献的序列的性质
假设我们后面能分成两个前缀最大值序列CCCDDD,满足siz(A)+siz(C)=siz(B)+siz(D)siz(A)+siz(C)=siz(B)+siz(D)siz(A)+siz(C)=siz(B)+siz(D),且满足ACACAC合并与BDBDBD合并后前缀最大值仍是前缀最大值,显然剩下的数总有方法填入的,故合法
分析序列CCCDDD的性质
先找出原序列的所有前缀最大值,不妨设这个序列为EEE
给出结论就是C,DC,DC,D之中至少有一个序列满足是EEE的完全子集,证明考虑反证
C,DC,DC,D中均存在至少一个数非原序列前缀最大值,则交换两数后C,DC,DC,D的前缀最大值个数均减少111,因为能制裁他们的都在对面的序列。连续做如上过程,我们总能使得其中至少一个序列成为EEE的完全子集
不妨枚举哪个序列成为了完全子集,以CCC为例
如果DDD序列中存在EEExxx个元素,存在非EEE序列中的mmm个元素,且第iii位之后存在TTTEEE中的元素,我们可以列出方程
cnt0+T−x=cnt1+x+mcnt0+T-x=cnt1+x+mcnt0+Tx=cnt1+x+m
移项可知是
cnt0−cnt1+T=2∗x+mcnt0-cnt1+T=2*x+mcnt0cnt1+T=2x+m
左边是常数,我们仅需要知道右边能构成这样的序列即可
将是EEE序列中的元素的权设为222,其他的元素权设为111,则我们需要知道一个上升序列且其权能够组成2∗x+m2*x+m2x+m
分奇偶讨论,我们发现如果一个数uuu能够组成,则u−2u-2u2也能够组成
维护两个bitbitbit,支持最大值查询即可

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<queue>
#include<vector>
#include<ctime>
#include<map>
#include<bitset>
#include<set>
#define LL long long
#define mp(x,y) make_pair(x,y)
#define pll pair<long long,long long>
#define pii pair<int,int>
using namespace std;
inline int read()
{
	int f=1,x=0;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
int stack[20];
inline void write(int x)
{
	if(x<0){putchar('-');x=-x;}
    if(!x){putchar('0');return;}
    int top=0;
    while(x)stack[++top]=x%10,x/=10;
    while(top)putchar(stack[top--]+'0');
}
inline void pr1(int x){write(x);putchar(' ');}
inline void pr2(int x){write(x);putchar('\n');}
const int MAXN=200005;
int n;
int lowbit(int x){return x&-x;}
struct bittree
{
	int s[MAXN];
	void modify(int x,int c){for(;x<=n;x+=lowbit(x))s[x]=max(s[x],c);}
	int qry(int x){int ret=0;for(;x>=1;x-=lowbit(x))ret=max(ret,s[x]);return ret;}
}bi[2];
int L[MAXN],R[MAXN],a[MAXN];
int is[MAXN],sum;
int re[2][MAXN][20];
void rebuild(int i)
{
	int u=n-a[i]+1;
	for(int x=u,cnt=1;x<=n;x+=lowbit(x),cnt++)bi[0].s[x]=re[0][i][cnt];
	for(int x=u,cnt=1;x<=n;x+=lowbit(x),cnt++)bi[1].s[x]=re[1][i][cnt];
}
void getit(int i)
{
	int u=n-a[i]+1;
	for(int x=u,cnt=1;x<=n;x+=lowbit(x),cnt++)re[0][i][cnt]=bi[0].s[x];
	for(int x=u,cnt=1;x<=n;x+=lowbit(x),cnt++)re[1][i][cnt]=bi[1].s[x];
}
void init()
{
	int mx=0;
	for(int i=1;i<=n;i++)if(a[i]>mx)is[i]=1,mx=a[i],sum++;
	for(int i=n;i>=1;i--)
	{
		int u=n-a[i]+1;
		int u1=bi[0].qry(u-1),u2=bi[1].qry(u-1);
		if(is[i])u1+=2,u2+=2;
		else u1+=1,u2+=1;
		getit(i);
		bi[u1&1].modify(u,u1);bi[u2&1].modify(u,u2);
	}
}
int ans[MAXN];
int mx[2],cnt[2];
bool ok(int i,int o)
{
	cnt[o]+=(a[i]>mx[o]);
	int u=n-mx[o^1]+1,c1=bi[0].qry(u-1),c2=bi[1].qry(u-1);
	int cal=cnt[o]-cnt[o^1]+sum;
	cnt[o]-=(a[i]>mx[o]);
	if(cal>=0)
	{
		if((cal&1)&&cal<=c2)return true;
		if((!(cal&1))&&cal<=c1)return true;
	}
	int ls=mx[o];
	cnt[o]+=(a[i]>mx[o]);mx[o]=max(mx[o],a[i]);
	u=n-mx[o]+1,c1=bi[0].qry(u-1),c2=bi[1].qry(u-1);
	cal=cnt[o^1]-cnt[o]+sum;
	mx[o]=ls;cnt[o]-=(a[i]>mx[o]);
	if(cal>=0)
	{
		if((cal&1)&&cal<=c2)return true;
		if((!(cal&1))&&cal<=c1)return true;
	}
	return false;
}
int main()
{
	n=read();
	for(int i=1;i<=n;i++)a[i]=read();
	init();
	for(int i=1;i<=n;i++)
	{
		rebuild(i);sum-=is[i];
		if(ok(i,0))ans[i]=0,cnt[0]+=(a[i]>mx[0]),mx[0]=max(mx[0],a[i]);
		else if(ok(i,1))ans[i]=1,cnt[1]+=(a[i]>mx[1]),mx[1]=max(mx[1],a[i]);
		else return puts("-1"),0;
	}
	for(int i=1;i<=n;i++)putchar(ans[i]+'0');
	puts("");
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值