bzoj 4552: [Tjoi2016&Heoi2016]排序

4552: [Tjoi2016&Heoi2016]排序

Description

在2016年,佳媛姐姐喜欢上了数字序列。因而他经常研究关于序列的一些奇奇怪怪的问题,现在他在研究一个难题
,需要你来帮助他。这个难题是这样子的:给出一个1到n的全排列,现在对这个全排列序列进行m次局部排序,排
序分为两种:1:(0,l,r)表示将区间[l,r]的数字升序排序2:(1,l,r)表示将区间[l,r]的数字降序排序最后询问第q
位置上的数字。

Input

输入数据的第一行为两个整数n和m。n表示序列的长度,m表示局部排序的次数。1 <= n, m <= 10^5第二行为n个整
数,表示1到n的一个全排列。接下来输入m行,每一行有三个整数op, l, r, op为0代表升序排序,op为1代表降序
排序, l, r 表示排序的区间。最后输入一个整数q,q表示排序完之后询问的位置, 1 <= q <= n。1 <= n <= 10^5
,1 <= m <= 10^5

Output

 输出数据仅有一行,一个整数,表示按照顺序将全部的部分排序结束后第q位置上的数字。

Sample Input

6 3
1 6 2 5 3 4
0 1 4
1 3 6
0 2 4
3

Sample Output

5
题解:
我们二分这个数mid,判断就是把序列中>=mid的数设为1,<mid的设为
如果是0操作,就把0全部放到区间的最前面,1放到最后面,用线段树实现好了。
1操作的话也差不多。
最后看看第q个是1还是0.
#include<stdio.h>
#include<iostream>
#include<algorithm>
using namespace std;
#define p1 (p<<1)
#define p2 (p<<1|1)
const int N=100005;
struct node
{
    int a,b,c;
}p[N];
int n,m,i,ans,k,a[N],b[N],t[N<<2],add[N<<2];
void build(int l,int r,int p)
{
    if(l==r)
    {
        t[p]=add[p]=b[l];
        return;
    }
    int mid=(l+r)>>1;
    build(l,mid,p1);
    build(mid+1,r,p2);
    t[p]=t[p1]+t[p2];
}
void update(int l,int r,int x,int y,int z,int p)
{
    if(x>y) return;
    if(add[p]==z) return;
    if(x<=l&&r<=y)
    {
        add[p]=z;
        t[p]=z*(r-l+1);
        return;
    }
    int mid=(l+r)>>1;
    if(add[p]!=-1)
    {
        update(l,mid,l,mid,add[p],p1);
        update(mid+1,r,mid+1,r,add[p],p2);
        add[p]=-1;
    }
    if(x<=mid) update(l,mid,x,y,z,p1);
    if(y>mid) update(mid+1,r,x,y,z,p2);
    t[p]=t[p1]+t[p2];
    //cout<<l<<' '<<r<<' '<<t[p]<<endl;
}
int solve(int l,int r,int x,int y,int p)
{
    if(x<=l&&r<=y) return t[p];
    int ans=0,mid=(l+r)>>1;
    if(add[p]!=-1)
    {
        update(l,mid,l,mid,add[p],p1);
        update(mid+1,r,mid+1,r,add[p],p2);
        add[p]=-1;
    }
    if(x<=mid) ans+=solve(l,mid,x,y,p1);
    if(y>mid) ans+=solve(mid+1,r,x,y,p2);
    return ans;
}
bool check(int x)
{
    int i;
    for(i=1;i<=n;i++)
        if(a[i]>=x) b[i]=1;else b[i]=0;
    for(i=1;i<=n<<2;i++)
    {
        t[i]=0;
        add[i]=-1;
    }
    build(1,n,1);
    for(i=1;i<=m;i++)
    {
        int sum=solve(1,n,p[i].a,p[i].b,1);
        if(p[i].c==0)
        {
            sum=p[i].b-p[i].a+1-sum;
            update(1,n,p[i].a,p[i].a+sum-1,0,1);
            update(1,n,p[i].a+sum,p[i].b,1,1);
        } else
        {
            update(1,n,p[i].a,p[i].a+sum-1,1,1);
            update(1,n,p[i].a+sum,p[i].b,0,1);
        }
    }
    if(solve(1,n,k,k,1)==0) return 0;else return 1;
}
int erfen(int l,int r)
{
    if(l>r) return r;
    int mid=(l+r)>>1;
    if(check(mid)==0) return erfen(l,mid-1);else
    return erfen(mid+1,r);
}
int main()
{
    scanf("%d%d",&n,&m);
    for(i=1;i<=n;i++)
        scanf("%d",&a[i]);
    for(i=1;i<=m;i++)
        scanf("%d%d%d",&p[i].c,&p[i].a,&p[i].b);
    scanf("%d",&k);
    cout<<erfen(1,n);
    return 0;
}

 

转载于:https://www.cnblogs.com/lwq12138/p/5517560.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值