【BZOJ4358】permu

版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/CreationAugust/article/details/51231980

Description

给出一个长度为n的排列P(P1,P2,…Pn),以及m个询问。每次询问某个区间[l,r]中,最长的值域
连续段长度。
Input

第一行两个整数n,m。
接下来一行n个整数,描述P。
接下来m行,每行两个整数l,r,描述一组询问。
Output

对于每组询问,输出一行一个整数,描述答案。
Sample Input

8 3

3 1 7 2 5 8 6 4

1 4

5 8

1 7
Sample Output

3

3

4
HINT

对于询问[1,4],P2,P4,P1组成最长的值域连续段[1,3];

对于询问[5,8],P8,P5,P7组成最长的值域连续段[4,6];

对于询问[1,7],P5,P7,P3,P6组成最长的值域连续段[5,8]。

1<=n,m<=50000

Source

By sumix173

跟Kangroo一样..改成把给定的n个排列数排序然后插入就行了
现在发现那种[L,R]->(L,R)的方式就是KDtree处理一系列使用树套树的序列问题的通法,甚至可以使用DFS序扩展到树上.

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MAXN 50010
#define Dnum 2
#define GET (ch>='0'&&ch<='9')
using namespace std;
inline void in(int &x)
{
    char ch=getchar();x=0;
    while (!GET)    ch=getchar();
    while (GET) x=x*10+ch-'0',ch=getchar();
}
int n,m,root;
bool cmp_d;
int id[MAXN];
struct KDtree
{
    int ch[2],f,d[Dnum],minn[Dnum],maxn[Dnum],maxv[Dnum],maxt[Dnum],mint;
    inline void init()  {   for (int i=0;i<Dnum;++i)    minn[i]=maxn[i]=d[i],maxv[i]=maxt[i]=0; }
    inline bool operator < (const KDtree& a)const   {   return d[cmp_d]<a.d[cmp_d]; }
}tree[MAXN];
struct Query
{
    int p,id;
    inline bool operator < (const Query& a)const    {   return p<a.p;   }
}q[MAXN];
inline void add(int rt,int x,int y)
{
    if (~x)
    {
        if (~tree[rt].mint)    tree[rt].maxv[0]=max(tree[rt].maxv[0],x-tree[rt].maxt[0]-1);
        else    tree[rt].mint=x;
        tree[rt].maxv[1]=max(tree[rt].maxv[1],x-tree[rt].maxt[1]-1);
        tree[rt].maxt[0]=tree[rt].maxt[1]=y;
    }
}
inline void push_up(int rt)
{
    for (int x=0,i=0;i<2;++i)
        if ((x=tree[rt].ch[i]))
            for (int j=0;j<Dnum;++j)
                tree[rt].minn[j]=min(tree[rt].minn[j],tree[x].minn[j]),
                tree[rt].maxn[j]=max(tree[rt].maxn[j],tree[x].maxn[j]);
}
inline void push_down(int rt)
{
    int x=tree[rt].mint,y=tree[rt].maxt[0];
    for (int i=0,j=0;i<2;++i)   if ((j=tree[rt].ch[i])) add(j,x,y);
    tree[rt].mint=-1;
}
int rebuild(int l=1,int r=m,bool d=0,int f=0)
{
    cmp_d=d;int mid=(l+r)>>1;nth_element(tree+l,tree+mid,tree+r+1);
    id[tree[mid].f]=mid;tree[mid].f=f;tree[mid].init();
    if (l!=mid) tree[mid].ch[0]=rebuild(l,mid-1,d^1,mid);
    if (r!=mid) tree[mid].ch[1]=rebuild(mid+1,r,d^1,mid);
    return push_up(mid),mid;
}
inline int check(int x1,int y1,int x2,int y2,int x,int y)
{
    int ret=0;
    ret+=(x1<=x&&y1>=y);ret+=(x1<=x&&y2>=y);ret+=(x2<=x&&y1>=y);ret+=(x2<=x&&y2>=y);
    return ret?(ret==4?1:2):0;
}
void modify(int rt,int x,int y,int id)
{
    int flag=check(tree[rt].minn[0],tree[rt].minn[1],tree[rt].maxn[0],tree[rt].maxn[1],x,y);
    if (!flag)  {   add(rt,id,id);return;   }
    if (flag==1)    return;
    if (tree[rt].d[0]>x||tree[rt].d[1]<y)   tree[rt].maxv[1]=max(tree[rt].maxv[1],id-tree[rt].maxt[1]-1),tree[rt].maxt[1]=id;
    push_down(rt);
    for (int i=0;i<2;++i)   if (tree[rt].ch[i]) modify(tree[rt].ch[i],x,y,id);
}
void update(int rt,int x)
{
    push_down(rt);
    tree[rt].maxv[1]=max(tree[rt].maxv[1],n-tree[rt].maxt[1]);
    tree[rt].maxv[1]=max(tree[rt].maxv[1],x);
    tree[rt].maxv[0]=max(tree[rt].maxv[0],x);
    tree[rt].maxv[1]=max(tree[rt].maxv[1],tree[rt].maxv[0]);
    for (int i=0;i<2;++i)   if (tree[rt].ch[i]) update(tree[rt].ch[i],tree[rt].maxv[0]);
}
int main()
{
    in(n);in(m);
    for (int i=1;i<=n;++i)  in(q[i].p),q[i].id=i;
    for (int i=1;i<=m;++i)  in(tree[i].d[0]),in(tree[i].d[1]),tree[i].f=i;
    root=rebuild();sort(q+1,q+n+1);
    for (int i=1;i<=n;++i)  modify(root,q[i].id,q[i].id,i);
    update(root,0);
    for (int i=1;i<=m;++i)  printf("%d\n",tree[id[i]].maxv[1]);
}
阅读更多
想对作者说点什么?
相关热词

博主推荐

换一批

没有更多推荐了,返回首页