[BZOJ3594]二维树状数组优化DP

说在前面

刚刚看到了@wzq_qwq 传送门 在消失两年后的一篇博客。还有岛姐的博客传送门，感慨满满啊。
NOIP已经不足一个月了，能不能往后走也只有这一次机会啊…

---

题目

BZOJ3594传送门
题面的话就自己去看吧qwq
已经是十点十分啦，晚自习放课要回寝室了呀…

---

解法&&大常数代码

就是一个dp，然后二维线段树维护
这还是me写的第一个二维线段树代码，没想到居然如此简单…

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std ;

int N , K , a[10005] , Hmax;
class bi_BIT{
    public :
        void Modify( int x , int y , int delta ){
            for( ; x <= K + 1 ; x += x&-x )
                for( int t = y ; t <= Hmax + 1 ; t += t&-t )
                    b[x][t] = max( b[x][t] , delta ) ;
        }
        int Query( int x , int y ){
            int rt = 0 ;
            for( ; x ; x -= x&-x )
                for( int t = y ; t ; t -= t&-t )
                    rt = max( b[x][t] , rt ) ;
            return rt ;
        }
    private:
        int b[505][5505] ;
}B ;

void solve(){
    for( int i = 1 ; i <= N ; i ++ )
        for( int j = K + 1 ; j >= 1 ; j -- )
            B.Modify( j , a[i] + j , B.Query( j , a[i] + j ) + 1 ) ;
    printf( "%d" , B.Query( K + 1 , Hmax + 1 ) ) ;
}

int main(){
    scanf( "%d%d" , &N , &K ) ;
    for( int i = 1 ; i <= N ; i ++ ){
        scanf( "%d" , &a[i] ) ;
        Hmax = max( Hmax , a[i] ) ;
    }
    Hmax += K ;
    solve() ;
}