1. 程式人生 > 其它 >luogu P6564 [POI2007] 堆積木KLO

luogu P6564 [POI2007] 堆積木KLO

題面傳送門
這個東西顯然dp可以搞。
\(f_i\)\(i\)處符合的最大符合數。
考慮一個\(j\)能轉移到\(i\)的條件。
\(j<i\)\(a_i<a_j\),\(i-j\geq a_i-a_j\)
把第三個變形一下就是\(i-a_i\geq j-a_j\)
好的我們可以cdq三維偏序了
注意觀察一下,後面兩個式子加起來就是第一個式子。
所以我們只要搞二維偏序即可。
樹狀陣列隨便維護。時間複雜度\(O(nlogn)\)
code:

#include <vector>
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<cmath>
#include<algorithm>
#include<bitset>
#include<set>
#include<map>
#define I inline
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define abs(x) ((x)>0?(x):-(x))
#define l(x) x<<1
#define r(x) x<<1|1
#define re register
#define ll long long
#define db double
#define N 100000
#define eps (1e-5)
#define mod (1<<31)
#define U unsigned int
using namespace std;
int n,x,Maxn,G[N+5<<5],now,ans;vector<int> F[N+5<<4];
I void get(int x,int y){while(x<=Maxn) G[x]=max(G[x],y),x+=x&-x;}
I int find(int x){int ans=0;while(x) ans=max(ans,G[x]),x-=x&-x;return ans;}
int main(){
	freopen("1.in","r",stdin);
	re int i,j;scanf("%d",&n);for(i=1;i<=n;i++) scanf("%d",&x),i>=x&&(F[i-x].push_back(x),Maxn=max(Maxn,x));
	for(i=0;i<=n;i++){
		for(j=0;j<F[i].size();j++)now=F[i][j],get(now,find(now-1)+1);
	}printf("%d\n",find(Maxn));
}