bzoj4553【TJOI2016&HEOI2016】序列

标签：

Description

Input

Output

Sample Input

Sample Output

动态规划+CDQ分治

每种变化只有一个只发生变化，所以每个位置i只有最大值mx[i]和最小值mn[i]有用。

DP方程：f[i]=max{f[j]}+1，其中j<i，a[j]≤mn[i]且mx[j]≤a[i]。

用CDQ分治优化DP。最初的时候所有编号是递增的，然后每次计算左边对右边影响的时候，让左面的a值有序，右面的mn值有序，这样就满足了前两个条件，第三个条件mx[j]≤a[i]用一个树状数组维护。每次计算完之后按照a值归并排序。

注意，递归到每一层的顺序是：处理左边→计算左边对右边的影响→处理右边。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<algorithm>
#define F(i,j,n) for(int i=j;i<=n;i++)
#define D(i,j,n) for(int i=j;i>=n;i--)
#define ll long long
#define N 100005
using namespace std;
int n,m,ans,f[N],s[N];
struct data{int v,mn,mx,id;}a[N],t[N];
inline int read()
{
	int x=0,f=1;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;
}
inline bool cmpmn(data a,data b){return a.mn<b.mn;}
inline bool cmpid(data a,data b){return a.id<b.id;}
inline int query(int x)
{
	int ret=0;
	for(;x;x-=(x&(-x))) ret=max(ret,s[x]);
	return ret;
}
inline void change(int x,int y){for(;x<=100000;x+=(x&(-x))) s[x]=max(s[x],y);}
inline void clear(int x){for(;x<=100000;x+=(x&(-x))) s[x]=0;}
void solve(int l,int r)
{
	if (l>=r) return;
	int mid=(l+r)>>1;
	solve(l,mid);
	sort(a+mid+1,a+r+1,cmpmn);
	int tmp=l;
	F(i,mid+1,r)
	{
		while (tmp<=mid&&a[tmp].v<=a[i].mn) change(a[tmp].mx,f[a[tmp].id]),tmp++;
		f[a[i].id]=max(f[a[i].id],query(a[i].v)+1);
	}
	F(i,l,tmp-1) clear(a[i].mx);
	sort(a+mid+1,a+r+1,cmpid);
	solve(mid+1,r);
	int l1=l,l2=mid+1;tmp=l;
	while (l1<=mid&&l2<=r) t[tmp++]=a[l1].v<a[l2].v?a[l1++]:a[l2++];
	while (l1<=mid) t[tmp++]=a[l1++];
	while (l2<=r) t[tmp++]=a[l2++];
	F(i,l,r) a[i]=t[i];
}
int main()
{
	n=read();m=read();
	F(i,1,n) a[i].v=a[i].mn=a[i].mx=read(),a[i].id=i,f[i]=1;
	F(i,1,m){int x=read(),y=read();a[x].mn=min(a[x].mn,y);a[x].mx=max(a[x].mx,y);}
	solve(1,n);
	F(i,1,n) ans=max(ans,f[i]);
	printf("%d\n",ans);
	return 0;
}

bzoj4553【TJOI2016&HEOI2016】序列

标签：

原文地址：http://blog.csdn.net/aarongzk/article/details/51676056