P4198 楼房重建 | 动态最长上升序列 线段树 区间合并

描述

n个点，初始为0，m次操作，每次修改第$X_i$个数的权值为$X_i/Y_i$，并要求输出此时最长上升序列长度。

输入

1≤$X_i$≤n，1≤$Y_i$≤1e9，1≤n,m≤1e5。

最后结果是整个区间内最长的上升序列长度，维护区间的最大值以及最长上升序列的长度，合并时，左区间的序列长度时一定要统计的，然后用左区间的最大值，在右区间中递归寻找，大于左区间最大值的最长上升序列。

代码

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<set>
#include<map>
#include<list>
#include<stack>
#include<queue>
#include<vector>
#define cl (k<<1)
#define cr (k<<1|1)
#define Mid ((a[k].l+a[k].r)>>1)
#define mem(a,x) memset(a,x,sizeof(a))
using namespace std;
#define bug printf("???\n")

typedef long long LL;

const int Inf=0x3f3f3f3f;
const double eps=1e-7;
const int maxn=1e5+50;

struct N{
    int l,r;
    int len;    //最长上升序列长度
    double h;   //最大值
}a[maxn*4];

void build(int k,int l,int r){
    a[k].l=l, a[k].r=r;
    a[k].h = 0;
    a[k].len = 0;
    if(l==r){
        return;
    }
    build(cl, l, Mid);
    build(cr, Mid+1, r);
}
int query(int k,int l,int r,double lim){
    if(a[k].l==a[k].r){
        return lim < a[k].h;
    }
    if(lim < a[cl].h) return query(cl, l, Mid, lim)+a[k].len-a[cl].len; //若没有这个会超时
    else              return query(cr, Mid+1, r, lim);
}
void update(int k,int i,double x){
    if(a[k].l==a[k].r){
        a[k].h = x;
        a[k].len = (x!=0);
        return;
    }
    if(i<=Mid) update(cl, i, x);
    else       update(cr, i, x);
    a[k].h = max(a[cl].h, a[cr].h);
    a[k].len = a[cl].len + query(cr, a[cr].l, a[cr].r, a[cl].h);
}

int main()
{
    int n,m;
    cin>>n>>m;
    build(1, 1, n);
    for(int i=1; i<=m; i++){
        int a,b;
        scanf("%d%d",&a,&b);
        update(1, a, 1.0*b/a);
        printf("%d\n",query(1, 1, n, 0));
    }
}