P4093 [TJOI2016]序列 [CDQ+DP]

传送门

我们发现

f[i] = max(f[j])+1(j<=i,Max[j]<=a[i],a[j]<=Min[i])

三维偏序问题用CDQ分治实现

对于[L-R]区间的DP的值, 我们先递归处理[L-Mid]

然后用[L-Mid]的值更新 [Mid+1-R], 最后递归处理[Mid+1 - R]

中途左边按a[i]排序, 右边按Min[i]排序, 然后两个指针扫就可以了

#include<bits/stdc++.h>
#define N 100050
using namespace std;
int n,m,c[N];
struct Node{
	int val, Max, Min, id, f;
}q[N];
bool cmpMin(Node a, Node b){ return a.Min < b.Min;}
bool cmpval(Node a, Node b){ return a.val < b.val;}
bool cmpid(Node a, Node b){ return a.id < b.id;}
int read(){
	int cnt = 0; char ch = 0;
	while(!isdigit(ch)) ch = getchar();
	while(isdigit(ch)) cnt = cnt * 10 + (ch-'0'), ch = getchar();
	return cnt;
}
void add(int x,int val){ for(;x<=N-50;x+=x&-x) c[x] = max(c[x],val); }
int quary(int x){int ans=0; for(;x;x-=x&-x) ans = max(ans, c[x]); return ans;}
void clear(int x){ for(;x<=N-50;x+=x&-x) c[x] = 0;}
void CDQ(int l,int r){
	if(l==r) return;
	int mid = (l+r) >> 1;
	CDQ(l,mid);
	sort(q+l, q+mid+1, cmpval);
	sort(q+mid+1, q+r+1, cmpMin);
	
	int i = l, j = mid+1;
	for(;j<=r;j++){
		while(i<=mid && q[j].Min >= q[i].val) add(q[i].Max, q[i].f), i++;
		q[j].f = max(q[j].f, quary(q[j].val) + 1);
	}
	for(i=l;i<=mid;i++) clear(q[i].Max);
	sort(q+l,q+r+1,cmpid);
	CDQ(mid+1,r);
}
int main(){
	n = read(), m = read();
	for(int i=1;i<=n;i++){
		int x = read();
		q[i] = (Node){x, x, x, i, 1};
	}
	for(int i=1;i<=m;i++){
		int x = read(), val = read();
		q[x].Min = min(q[x].Min, val);
		q[x].Max = max(q[x].Max, val);
	}
	CDQ(1,n); int ans = 0;
	for(int i=1;i<=n;i++) ans = max(ans, q[i].f);
	printf("%d",ans); return 0;
}

 

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

FSYo

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值