E. Replace Sort
题意:
给你长度为
N
<
=
5
e
5
N<=5e5
N<=5e5的序列
A
A
A,和大小为
5
e
5
5e5
5e5的集合
B
B
B,可以用一次操作 将集合
B
B
B 中的一个元素 覆盖 序列
A
A
A上一个位置的元素,即 让
A
[
i
]
=
B
[
j
]
A[i]=B[j]
A[i]=B[j];
问 最小的操作次数,使得 序列
A
A
A变成严格升序。
思路:
考虑一个很显然的
D
P
DP
DP,
d
p
1
[
i
]
dp1[i]
dp1[i]表示,最终序列第
i
i
i 位是
A
[
i
]
A[i]
A[i] 时的最小操作次数;
d
p
2
[
i
]
[
j
]
dp2[i][j]
dp2[i][j] 表示 最终序列第
i
i
i 为 是
B
[
j
]
B[j]
B[j] 时的最小操作次数。
那么状态转移也是显然的:
i
f
A
[
i
−
1
]
<
A
[
i
]
,
d
p
1
[
i
]
=
d
p
1
[
i
−
1
]
if A[i-1] < A[i],dp1[i]=dp1[i-1]
ifA[i−1]<A[i],dp1[i]=dp1[i−1] ;
对
于
所
有
的
j
满
足
B
[
j
]
<
A
[
i
]
,
d
p
1
[
i
]
=
m
i
n
(
d
p
2
[
i
−
1
]
[
j
]
)
对于所有的 j 满足 B[j]<A[i],dp1[i]=min( dp2[i-1][j] )
对于所有的j满足B[j]<A[i],dp1[i]=min(dp2[i−1][j]);
对
于
所
有
的
k
满
足
B
[
k
]
<
B
[
j
]
,
d
p
2
[
i
]
[
j
]
=
m
i
n
(
d
p
2
[
i
−
1
]
[
k
]
+
1
)
对于所有的k 满足 B[k]<B[j],dp2[i][j]=min( dp2[i-1][k]+1 )
对于所有的k满足B[k]<B[j],dp2[i][j]=min(dp2[i−1][k]+1);
基于贪心原理,转移式可以简化为
d
p
2
[
i
]
[
j
]
=
d
p
2
[
i
−
1
]
[
j
−
1
]
+
1
dp2[i][j]=dp2[i-1][j-1]+1
dp2[i][j]=dp2[i−1][j−1]+1;
i f A [ i − 1 ] < B [ j ] , d p 2 [ i ] [ j ] = m i n ( d p 2 [ i ] [ j ] , d p 1 [ i − 1 ] + 1 ) if A[i-1]<B[j],dp2[i][j]=min( dp2[i][j],dp1[i-1]+1 ) ifA[i−1]<B[j],dp2[i][j]=min(dp2[i][j],dp1[i−1]+1);
但是这个状态是
O
(
n
2
)
O(n^2)
O(n2) 的,考虑优化:
对于
d
p
1
dp1
dp1 的转移 只需要支持
d
p
2
dp2
dp2 状态的区间询问最小值;
对于
d
p
2
dp2
dp2 的转移 可以简化为
如果
A
[
i
−
1
]
<
B
[
j
]
A[i-1]<B[j]
A[i−1]<B[j],
d
p
2
[
i
]
[
j
]
=
m
i
n
(
d
p
2
[
i
−
1
]
[
j
−
1
]
+
1
,
d
p
1
[
i
−
1
]
+
1
)
dp2[i][j]=min( dp2[i-1][j-1]+1,dp1[i-1] +1)
dp2[i][j]=min(dp2[i−1][j−1]+1,dp1[i−1]+1)
否则
d
p
2
[
i
−
1
]
[
j
]
=
d
p
2
[
i
−
1
]
[
j
−
1
]
+
1
dp2[i-1][j]=dp2[i-1][j-1]+1
dp2[i−1][j]=dp2[i−1][j−1]+1;
考虑 从
d
p
1
dp1
dp1转移,需要支持 区间取
m
i
n
min
min的操作。
对于
j
j
j的状态值,从
j
−
1
j-1
j−1赋值再加一,可以通过线段树的经典
t
r
i
c
k
trick
trick,维护一个窗口,每次转移 让窗口左移一位,再区间加一 即可。
综上所述,一颗吉老师 就可以了。
具体转移和技巧 看代码~
代码:
#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int INF=1e9;
inline int qread(){
int s=0,w=1;char ch=getchar();
for(;!isdigit(ch);ch=getchar())if(ch=='-')w=-1;
for (;ch>='0'&&ch<='9';ch=getchar())s=(s<<1)+(s<<3)+(ch^48);
return (w==-1?-s:s);}
int n,m;
int st;
int a[500059],b[500050];
int dp[500050];
struct node{
int mx,smx,mi;
int add1,add2;
}tr[4000050];
#define ls p<<1
#define rs p<<1|1
#define mid (l+r>>1)
void pushup(int p){
tr[p].mi=min(tr[ls].mi,tr[rs].mi);
if(tr[ls].mx==tr[rs].mx){
tr[p].mx=tr[ls].mx;
tr[p].smx=max(tr[ls].smx,tr[rs].smx);
}else if(tr[ls].mx>tr[rs].mx){
tr[p].mx=tr[ls].mx;
tr[p].smx=max(tr[ls].smx,tr[rs].mx);
}else {
tr[p].mx=tr[rs].mx;
tr[p].smx=max(tr[ls].mx,tr[rs].smx);
}
}
void update(int p,int k1,int k2){
if(tr[p].mi==tr[p].mx){
tr[p].mi+=k1;tr[p].mx+=k1;
k2=k1;
}else{
tr[p].mx+=k1;tr[p].mi+=k2;tr[p].smx+=k2;
}
tr[p].add1+=k1;tr[p].add2+=k2;
}
void pushdown(int p){
int mx=max(tr[ls].mx,tr[rs].mx);
if(tr[ls].mx==mx)update(ls,tr[p].add1,tr[p].add2);
else update(ls,tr[p].add2,tr[p].add2);
if(tr[rs].mx==mx)update(rs,tr[p].add1,tr[p].add2);
else update(rs,tr[p].add2,tr[p].add2);
tr[p].add1=tr[p].add2=0;
}
void build(int p,int l,int r){
tr[p].add1=tr[p].add2=0;
if(l==r){
if(l<=n)tr[p].mi=tr[p].mx=INF;
else tr[p].mi=tr[p].mx=1;
tr[p].smx=-INF;
return;
}
build(ls,l,mid);
build(rs,mid+1,r);
pushup(p);
}
void update1(int p,int l,int r,int x,int y,int w){
if(x<=l&&r<=y){
update(p,w,w);
return;
}
pushdown(p);
if(x<=mid)update1(ls,l,mid,x,y,w);
if(mid<y)update1(rs,mid+1,r,x,y,w);
pushup(p);
}
void update2(int p,int l,int r,int x,int y,int w){
if(tr[p].mx<=w)return;
if(x<=l&&r<=y&&tr[p].smx<w){
update(p,w-tr[p].mx,0);
return;
}
pushdown(p);
if(x<=mid)update2(ls,l,mid,x,y,w);
if(mid<y)update2(rs,mid+1,r,x,y,w);
pushup(p);
}
int query(int p,int l,int r,int x,int y){
if(x<=l&&r<=y)return tr[p].mi;
pushdown(p);
int ans=INF;
if(x<=mid)ans=min(ans,query(ls,l,mid,x,y));
if(mid<y)ans=min(ans,query(rs,mid+1,r,x,y));
return ans;
}
#undef ls
#undef rs
#undef mid
int main(){
n=qread(),m=qread();
st=n;
for(int i=1;i<=n;i++)a[i]=qread();
for(int i=1;i<=m;i++)b[i]=qread();
sort(b+1,b+1+m);
build(1,1,n+m);
for(int i=2;i<=n;i++){
if(a[i-1]<a[i])dp[i]=dp[i-1];
else dp[i]=INF;
int pos=lower_bound(b+1,b+1+m,a[i])-b;
int res=query(1,1,n+m,st+1,st+pos-1);
pos=lower_bound(b+1,b+1+m,a[i-1])-b;
if(dp[i]>res)dp[i]=res;
st--;
update1(1,1,n+m,st+1,st+m,1);
update2(1,1,n+m,st+pos,st+m,dp[i-1]+1);
}
int ans=dp[n];
int res=query(1,1,n+m,st+1,st+m);
ans=min(ans,res);
if(ans>n)printf("-1\n");
else printf("%d\n",ans);
return 0;
}