2021 ICPC Southeastern Europe Regional Contest E. Replace Sort

E. Replace Sort

题意：
给你长度为$N<=5e5$的序列$A$，和大小为$5e5$的集合$B$，可以用一次操作 将集合$B$ 中的一个元素 覆盖 序列 $A$上一个位置的元素，即 让 $A[i]=B[j]$；
问 最小的操作次数，使得 序列$A$变成严格升序。

思路：
考虑一个很显然的 $DP$，$dp1[i]$表示，最终序列第 $i$ 位是 $A[i]$ 时的最小操作次数；$dp2[i][j]$ 表示 最终序列第 $i$ 为 是$B[j]$ 时的最小操作次数。

那么状态转移也是显然的：
$ifA[i-1]<A[i]，dp1[i]=dp1[i-1]$ ;
$对于所有的j满足B[j]<A[i]，dp1[i]=min(dp2[i-1][j])$；

$对于所有的k满足B[k]<B[j]，dp2[i][j]=min(dp2[i-1][k]+1)$；
基于贪心原理，转移式可以简化为 $dp2[i][j]=dp2[i-1][j-1]+1$；

$ifA[i-1]<B[j]，dp2[i][j]=min(dp2[i][j]，dp1[i-1]+1)$；

但是这个状态是 $O(n^2)$ 的，考虑优化：
对于 $dp1$ 的转移 只需要支持 $dp2$ 状态的区间询问最小值；
对于 $dp2$ 的转移 可以简化为
如果 $A[i-1]<B[j]$， $dp2[i][j]=min(dp2[i-1][j-1]+1，dp1[i-1]+1)$
否则 $dp2[i-1][j]=dp2[i-1][j-1]+1$；
考虑 从$dp1$转移，需要支持 区间取$min$的操作。
对于 $j$的状态值，从$j-1$赋值再加一，可以通过线段树的经典$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;
}