[bzoj 2957] 楼房重建

传送门

Description

动态维护\(a_i=\frac{w_i}{i}\)严格上升序列的长度

Solution 

用线段树维护每个区间的答案，以及每个区间的最大值

合并的时候，显然最多只会往一边传下去，所以复杂度为\(O(\log len)\)

总复杂度是\(O(m\log^2 n)\)

Code 

#include<bits/stdc++.h>
#define db double
#define ll long long
#define mid (T[x].l+T[x].r>>1)
using namespace std;
inline int read()
{
    int x=0;char ch=getchar();
    while(ch<'0'||ch>'9')ch=getchar();
    while(ch>='0'&&ch<='9')x=(x<<1)+(x<<3)+ch-'0',ch=getchar();
    return x;
}

const int MN=1e5+5;

struct{int l,r,ans;db mx;}T[MN<<2];

int G(int x,db v)
{
    if(T[x].l==T[x].r) return T[x].mx>v;
    if(T[x<<1].mx<=v) return G(x<<1|1,v);
	return G(x<<1,v)+T[x].ans-T[x<<1].ans;
}

void up(int x)
{
    int l=x<<1,r=x<<1|1;
    T[x].ans=T[l].ans;T[x].mx=max(T[l].mx,T[r].mx);
    if(T[r].mx>T[l].mx)T[x].ans+=G(r,T[l].mx);
}

void Build(int x,int l,int r)
{
    T[x].mx=0;T[x].ans=1;
    if((T[x].l=l)==(T[x].r=r))return;
	Build(x<<1,l,mid),Build(x<<1|1,mid+1,r);
}

void Modify(int x,int a,db v)
{
    if(T[x].l==T[x].r)return(void)(T[x].mx=v);
    Modify(x<<1|(a>mid),a,v);up(x);
}

int main()
{
    int n=read(),m=read(),x;Build(1,0,n);
    while(m--) x=read(),Modify(1,x,(db)read()/x),printf("%d\n",T[1].ans-1);
    return 0;
}

---

Blog来自PaperCloud，未经允许，请勿转载，TKS！