真神题,场上只过了2个,硬想了几个小时还是想出来了。
先考虑如何判断是否有解。
我们先观察一个合法的解具备的性质。我们称
P
P
P中前缀的high位置为关键位。对于
X
X
X和
Y
Y
Y分别的high位置,一定会包含所有关键位,反过来,如果我们钦定两个长度相等且不相交的上升子序列分别作为
X
X
X和
Y
Y
Yhigh位置的集合,并且它们包含所有关键位,那么容易构造出
s
s
s。因此,我们只需要考虑钦定两个上升子序列
U
U
U和
V
V
V即可。显然存在一个方案,使得
U
U
U和
V
V
V中至多有一个包含了非关键位。
如果关键位个数
s
s
s为偶数,直接分成两半即可,否则我们必须选至少一个非关键位。我们考虑令关键位的权值为
2
2
2,非关键位的权值为
1
1
1,我们断言有解当且仅当存在一个包含非关键位的上升子序列权值和
>
s
>s
>s。这个是显然的,因为我们可以不断的删去选择的元素直到和为
s
s
s,那么令
U
U
U是此时的子序列,
V
V
V是没选的关键位集合即可。
现在我们还要求求出字典序最小的
s
s
s,考虑按位贪心并判定。我们可以用类似的思路判断当前枚举的位置,具体的,加上
0
0
0后,令前缀中
X
X
X和
Y
Y
Yhigh位置的差为
d
d
d,剩余后缀中关键位数目为
t
t
t。不妨设
d
≥
0
d\geq 0
d≥0,如果
d
d
d和
t
t
t同奇偶,那么我们需要选出一个权值和
≥
d
+
t
\geq d+t
≥d+t的上升子序列接在当前
Y
Y
Y的最大值后面,如果不同奇偶,我们需要选出一个权值和
≥
d
+
t
\geq d+t
≥d+t且包含非关键位的上升子序列接在当前
Y
Y
Y的最大值后面,或者选出一个权值和
≥
d
−
t
\geq d-t
≥d−t且包含非关键位的上升子序列接在当前
X
X
X的最大值后面。
具体可以对每个位置预处理出以它开头,包含与不包含非关键位的上升子序列的最大权值和,然后用可持久化线段树支持查询后缀中大于某个数开头的最大值即可。
时间复杂度
O
(
n
l
o
g
n
)
\mathcal O(nlogn)
O(nlogn)。
#include <bits/stdc++.h>
#define FR first
#define SE second
#define inf 0x3f3f3f3f
using namespace std;
typedef pair<int,int> pr;
namespace SGT {
const int Maxn=5000000;
int ch[Maxn][2],tot;
int max1[Maxn],max2[Maxn];
inline int cpynode(int o) {
tot++;
ch[tot][0]=ch[o][0];ch[tot][1]=ch[o][1];
max1[tot]=max1[o];max2[tot]=max2[o];
return tot;
}
inline void pushup(int o) {
max1[o]=max(max1[ch[o][0]],max1[ch[o][1]]);
max2[o]=max(max2[ch[o][0]],max2[ch[o][1]]);
}
int update(int l,int r,int o,int p,pr q) {
o=cpynode(o);
if (l==r) {
max1[o]=q.FR;max2[o]=q.SE;
}
else {
int m=((l+r)>>1);
if (m>=p) ch[o][0]=update(l,m,ch[o][0],p,q);
else ch[o][1]=update(m+1,r,ch[o][1],p,q);
pushup(o);
}
return o;
}
int query(int l,int r,int o,int p,bool q) {
if (!o) return (!q)?0:-inf;
if (l==r) return (!q)?max1[o]:max2[o];
else {
int m=((l+r)>>1);
if (m<p) return query(m+1,r,ch[o][1],p,q);
else {
int t=query(l,m,ch[o][0],p,q);
return (!q)?max(t,max1[ch[o][1]]):max(t,max2[ch[o][1]]);
}
}
}
}
int root[200005],num[200005];
bool tag[200005];
int pre(int n) {
int s=0,maxn=0;
for(int i=1;i<=n;i++)
if (num[i]>maxn) {
maxn=num[i];
s++;
tag[i]=1;
}
for(int i=n;i>0;i--)
if (tag[i]) {
int t1=SGT::query(1,n+1,root[i+1],num[i]+1,0),t2=SGT::query(1,n+1,root[i+1],num[i]+1,1);
root[i]=SGT::update(1,n+1,root[i+1],num[i],pr(t1+2,t2+2));
}
else {
int t=SGT::query(1,n+1,root[i+1],num[i]+1,0);
root[i]=SGT::update(1,n+1,root[i+1],num[i],pr(t+1,t+1));
}
return s;
}
int main() {
SGT::max1[0]=0;SGT::max2[0]=-inf;
int n;
scanf("%d",&n);
for(int i=1;i<=n;i++) scanf("%d",&num[i]);
int s=pre(n);
if ((s&1)&&SGT::query(1,n+1,root[1],1,1)<s) {
puts("-1");
return 0;
}
int r1=0,r2=0,d=0;
for(int i=1;i<=n;i++) {
s-=tag[i];
int t=d+(num[i]>r1),c=max(r1,num[i]);
bool ok=1;
if ((t^s)&1) {
if (t>=0) {
if (SGT::query(1,n+1,root[i+1],r2+1,1)<s+t&&(t>s||SGT::query(1,n+1,root[i+1],c+1,1)<s-t)) ok=0;
}
else {
if (SGT::query(1,n+1,root[i+1],c+1,1)<s-t&&(-t>s||SGT::query(1,n+1,root[i+1],r2+1,1)<s+t)) ok=0;
}
}
else {
if (t>=0) {
if (SGT::query(1,n+1,root[i+1],r2+1,0)<s+t) ok=0;
}
else {
if (SGT::query(1,n+1,root[i+1],c+1,0)<s-t) ok=0;
}
}
if (ok) {
putchar('0');
d+=(num[i]>r1);
r1=max(r1,num[i]);
}
else {
putchar('1');
d-=(num[i]>r2);
r2=max(r2,num[i]);
}
}
printf("\n");
return 0;
}