4358: permu|K-Dtree

最新推荐文章于 2019-07-25 21:08:00 发布

ws_yzy

最新推荐文章于 2019-07-25 21:08:00 发布

阅读量725

点赞数

分类专栏： K-D tree 文章标签： K-DTree

本文链接：https://blog.csdn.net/ws_yzy/article/details/50883566

版权

K-D tree 专栏收录该内容

5 篇文章 0 订阅

订阅专栏

几天之前写了个莫队+线段树T的飞起。。（据说卡常数可以卡过去。。
然后可以KD树做，标记比较难处理。
ymClairs的KD树神之标记http://www.cnblogs.com/clrs97/p/5049090.html

#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<vector>
#include<cstdio>
#include<queue>
#include<cmath>
#include<set>
#include<map>
#define L t[x].l
#define R t[x].r
#define ll long long
#define N 100005
using namespace std;
int sc()
{
    int i=0; char c=getchar();
    while(c>'9'||c<'0')c=getchar();
    while(c>='0'&&c<='9')i=i*10+c-'0',c=getchar();
    return i;
} 
struct W{
    int D[2],mx[2],mn[2],p,l,r;
    int m,d,e,hm,hd,he;
}t[N];
int n,m,D,root,a[N],ans[N];
bool cmp(W a,W b){return a.D[D]<b.D[D];}
void push_up(int x)
{
    t[x].he=t[x].e=-1;
    for(int i=0;i<=1;i++)
    {
        t[x].mn[i]=min(t[x].D[i],min(t[L].mn[i],t[R].mn[i]));
        t[x].mx[i]=max(t[x].D[i],max(t[L].mx[i],t[R].mx[i]));
    }
}
void build(int &x,int l,int r,int dir)
{
    x=l+r>>1;D=dir;
    nth_element(t+l,t+x,t+r+1,cmp);
    if(l<x)build(L,l,x-1,dir^1);
    if(r>x)build(R,x+1,r,dir^1);
    push_up(x);
}
void hdoa(int x,int v)
{
    t[x].hm=max(t[x].hm,t[x].m+v);
    if(t[x].e>-1)
         t[x].he=max(t[x].he,t[x].e+v);
    else t[x].hd=max(t[x].hd,t[x].d+v);
}
void hdoc(int x,int v)
{
    t[x].hm=max(t[x].hm,v);
    t[x].he=max(t[x].he,v);
}
void doc(int x,int v)
{
    t[x].hm=max(t[x].hm,t[x].m=v);
    t[x].he=max(t[x].he,t[x].e=v);
    t[x].d=0;
}
void doa(int x,int v)
{
    t[x].hm=max(t[x].hm,t[x].m+=v);
    if(t[x].e>-1)
         t[x].he=max(t[x].he,t[x].e+=v);
    else t[x].hd=max(t[x].hd,t[x].d+=v);
}
void push_down(int x)
{
    if(t[x].hd)
    {
        if(L)hdoa(L,t[x].hd);
        if(R)hdoa(R,t[x].hd);
        t[x].hd=0;
    }
    if(t[x].he>-1)
    {
        if(L)hdoc(L,t[x].he);
        if(R)hdoc(R,t[x].he);
        t[x].he=-1;
    }
    if(t[x].d)
    {
        if(L)doa(L,t[x].d);
        if(R)doa(R,t[x].d);
        t[x].d=0;
    }else if(t[x].e>-1)
    {
        if(L)doc(L,t[x].e);
        if(R)doc(R,t[x].e);
        t[x].e=-1;
    }
}
void change(int x,int p)
{
    if(t[x].mn[0]>p||t[x].mx[1]<p){doc(x,0);return ;}
    if(t[x].mn[1]>=p&&t[x].mx[0]<=p){doa(x,1);return ;}
    push_down(x);
    if(t[x].D[0]<=p&&t[x].D[1]>=p)
        t[x].hm=max(t[x].hm,++t[x].m);
    else t[x].m=0;
    if(L)change(L,p);
    if(R)change(R,p);
}
void dfs(int x)
{
    if(!x)return ;
    push_down(x);
    ans[t[x].p]=t[x].hm;
    dfs(L);dfs(R);
}
int main()
{
    t[0].mx[0]=t[0].mx[1]=-1e9;
    t[0].mn[0]=t[0].mn[1]=1e9;
    n=sc();m=sc();
    for(int i=1;i<=n;i++)a[sc()]=i;
    for(int i=1;i<=m;i++)
        t[i].D[0]=sc(),t[i].D[1]=sc(),t[i].p=i;
    build(root,1,m,0);
    for(int i=1;i<=n;i++)
        change(root,a[i]);
    dfs(root);
    for(int i=1;i<=m;i++)printf("%d\n",ans[i]);
    return 0;
}

附上莫队+线段树代码（可以尝试一下卡常数。。说不定真能直接卡过去）

#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<vector>
#include<cstdio>
#include<queue>
#include<cmath>
#include<set>
#include<map>
#define N 200105
#define L x<<1
#define R x<<1|1
using namespace std;
int sc()
{
    int i=0,f=1;char c=getchar();
    while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9')i=i*10+c-'0',c=getchar();
    return i*f;
}
struct T{int lmx,rmx,mx,v;}t[N*4];
struct W{int l,r,pos;}a[N];
int bl[N],ans[N],v[N];
int n,m,block,tot;
inline bool cmp1(W a,W b)
{
    return bl[a.l]<bl[b.l]||(bl[a.l]==bl[b.l]&&a.r<b.r);
}
inline void change(int x,int l,int r,int v,int f)
{
    if(l==r)
    {
        t[x].v+=f;
        t[x].mx=t[x].lmx=t[x].rmx=(bool)t[x].v;
        return;
    }
    int mid=l+r>>1;
    if(v<=mid)change(L,l,mid,v,f);else change(R,mid+1,r,v,f);
    t[x].lmx=t[L].lmx; t[x].rmx=t[R].rmx;
    if(t[L].lmx==mid-l+1)t[x].lmx+=t[R].lmx;
    if(t[R].rmx==r-mid)t[x].rmx+=t[L].rmx;
    t[x].mx=max(t[L].rmx+t[R].lmx,max(t[L].mx,t[R].mx));
}
int main()
{
    n=sc(),m=sc();block=sqrt(n);
    for(int i=1;i<=n;i++)
    {
        bl[i]=(i-1)/block+1;
        v[i]=sc();
    }
    for(int i=1;i<=m;i++)
        a[i].l=sc(),a[i].r=sc(),a[i].pos=i;
    sort(a+1,a+m+1,cmp1);
    int l=a[1].l,r=a[1].l-1;
    for(int i=1;i<=m;i++)
    {
        while(l<a[i].l)change(1,1,n,v[l++],-1);
        while(l>a[i].l)change(1,1,n,v[--l],1);
        while(r>a[i].r)change(1,1,n,v[r--],-1);
        while(r<a[i].r)change(1,1,n,v[++r],1);
        ans[a[i].pos]=t[1].mx;
    }
    for(int i=1;i<=m;i++)
        printf("%d\n",ans[i]);
    return 0;
}