[bzoj 2957] 楼房重建

最新推荐文章于 2019-08-30 20:43:00 发布

weixin_34344403

最新推荐文章于 2019-08-30 20:43:00 发布

阅读量77

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/PaperCloud/p/10919912.html

博客围绕动态维护\(a_i=\frac{w_i}{i}\)严格上升序列的长度展开。采用线段树维护每个区间的答案和最大值，合并时最多往一边传，复杂度为\(O(\log len)\)，总复杂度是\(O(m\log^2 n)\)，还给出代码来源。

传送门

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！

转载于:https://www.cnblogs.com/PaperCloud/p/10919912.html