【DP+树状数组】波动序列

题目

思路

用树状数组优化即可

代码

#include<bits/stdc++.h>
using namespace std;
const int N=1e5+77;
int n,m,aii,a[4][N],f[4][N],yjy[4][N<<2],b[N<<2];
int main()
{
	freopen("sequence.in","r",stdin); freopen("sequence.out","w",stdout);
	scanf("%d",&n);
	for(int j=0; j<3; j++) for(int i=1; i<=n; i++)
	{
		scanf("%d",&a[j][i]); b[j*n+i]=a[j][i];
	}
	sort(b+1,b+3*n+1); m=unique(b+1,b+3*n+1)-b;
	for(int i=1; i<=n; i++)
	{
		a[0][i]=lower_bound(b+1,b+m,a[0][i])-b,
		a[1][i]=m-(lower_bound(b+1,b+m,a[1][i])-b),
		a[2][i]=lower_bound(b+1,b+m,a[2][i])-b,
		a[3][i]=m-a[2][i];
	}
	for(int i=1; i<=n; i++)
	{
		for(int j=0; j<4; j++) 
		{
			for(int k=a[j][i]; k; k-=k&-k) f[j][i]=max(f[j][i],yjy[j][k]); 
			f[j][i]++;
		}
		for(int k=a[0][i]; k<m; k+=k&-k) yjy[0][k]=max(yjy[0][k],f[0][i]);
		for(int k=a[2][i]; k<m; k+=k&-k) yjy[0][k]=max(yjy[0][k],f[2][i]);
		for(int k=m-a[1][i]; k<m; k+=k&-k) yjy[0][k]=max(yjy[0][k],f[1][i]);
		for(int k=m-a[3][i]; k<m; k+=k&-k) yjy[0][k]=max(yjy[0][k],f[3][i]);
		
		for(int k=m-a[0][i]; k<m; k+=k&-k) yjy[1][k]=max(yjy[1][k],f[0][i]);
		for(int k=m-a[2][i]; k<m; k+=k&-k) yjy[1][k]=max(yjy[1][k],f[2][i]);
		for(int k=a[1][i]; k<m; k+=k&-k) yjy[1][k]=max(yjy[1][k],f[1][i]);
		for(int k=a[3][i]; k<m; k+=k&-k) yjy[1][k]=max(yjy[1][k],f[3][i]);
		
		for(int k=a[0][i]; k<m; k+=k&-k) yjy[2][k]=max(yjy[2][k],f[0][i]);
		for(int k=a[2][i]; k<m; k+=k&-k) yjy[2][k]=max(yjy[2][k],f[2][i]);
		for(int k=m-a[1][i]; k<m; k+=k&-k) yjy[2][k]=max(yjy[2][k],f[1][i]);
		
		for(int k=m-a[0][i]; k<m; k+=k&-k) yjy[3][k]=max(yjy[3][k],f[0][i]);
		for(int k=a[1][i]; k<m; k+=k&-k) yjy[3][k]=max(yjy[3][k],f[1][i]);
		for(int k=a[3][i]; k<m; k+=k&-k) yjy[3][k]=max(yjy[3][k],f[3][i]);
		
		aii=max(aii,max(max(f[0][i],f[1][i]),max(f[2][i],f[3][i])));
	}
	printf("%d",aii); return 0;

}