hdu 5649 DZY Loves Sorting 二分+线段树

题目链接

给一个序列, 两种操作, 一种是将[l, r]里所有数升序排列, 一种是降序排列。
所有操作完了之后, 问你a[k]等于多少。




真心是涨见识了这题..好厉害。
因为最后只询问一个位置, 所以我们二分这个位置的值。 将所有大于等于它的值赋为1, 小于的赋为0. 然后现在整个序列只有01, 更改什么的线段树就很好搞。
如果a[k]最后为1, 那么我们增加l, 否则减少r。


那么为什么可以这样呢。
因为二分mid的时候, 将等于mid的也赋为了1, 所以如果a[k]是0的话,代表mid比答案大, 所以减少r。
那么a[k]为1的时候有两种, 一种是mid刚好是答案, 另一种是mid比答案小, 所以我们增大l逼近答案。

#include <iostream>
#include <vector>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <complex>
#include <cmath>
#include <map>
#include <set>
#include <string>
#include <queue>
#include <stack>
#include <bitset>
using namespace std;
#define pb(x) push_back(x)
#define ll long long
#define mk(x, y) make_pair(x, y)
#define lson l, m, rt<<1
#define mem(a) memset(a, 0, sizeof(a))
#define rson m+1, r, rt<<1|1
#define mem1(a) memset(a, -1, sizeof(a))
#define mem2(a) memset(a, 0x3f, sizeof(a))
#define rep(i, n, a) for(int i = a; i<n; i++)
#define fi first
#define se second
typedef complex <double> cmx;
typedef pair<int, int> pll;
const double PI = acos(-1.0);
const double eps = 1e-8;
const int mod = 1e9+7;
const int inf = 1061109567;
const int dir[][2] = { {-1, 0}, {1, 0}, {0, -1}, {0, 1} };
const int maxn = 1e5+5;
int sum[maxn<<2], a[maxn], cover[maxn<<2];
pair<int, pll> q[maxn];
void push_down(int rt, int m) {
    if(~cover[rt]) {
        sum[rt<<1] = cover[rt]*(m-(m>>1));
        sum[rt<<1|1] = cover[rt]*(m>>1);
        cover[rt<<1] = cover[rt<<1|1] = cover[rt];
        cover[rt] = -1;
    }
}
void update(int L, int R, int l, int r, int rt, int val) {
    if(L>R)
        return ;
    if(L<=l&&R>=r) {
        cover[rt] = val;
        sum[rt] = val*(r-l+1);
        return ;
    }
    push_down(rt, r-l+1);
    int m = l+r>>1;
    if(L<=m)
        update(L, R, lson, val);
    if(R>m)
        update(L, R, rson, val);
    sum[rt] = sum[rt<<1] + sum[rt<<1|1];
}
int query(int L, int R, int l, int r, int rt) {
    if(L<=l&&R>=r) {
        return sum[rt];
    }
    push_down(rt, r-l+1);
    int m = l+r>>1, ret = 0;
    if(L<=m)
        ret += query(L, R, lson);
    if(R>m)
        ret += query(L, R, rson);
    return ret;
}
void check(int x, int n, int m) {
    mem(sum); mem1(cover);
    for(int i = 1; i <= n; i++)
        if(a[i]>=x)
            update(i, i, 1, n, 1, 1);
    for(int i = 0; i < m; i++) {
        int tmp = query(q[i].fi, q[i].se.fi, 1, n, 1);
        if(q[i].se.se) {
            update(q[i].fi, q[i].fi+tmp-1, 1, n, 1, 1);
            update(q[i].fi+tmp, q[i].se.fi, 1, n, 1, 0);
        } else {
            update(q[i].se.fi-tmp+1, q[i].se.fi, 1, n, 1, 1);
            update(q[i].fi, q[i].se.fi-tmp, 1, n, 1, 0);
        }

    }

}
void solve() {
    int n, m, k;
    cin>>n>>m;
    for(int i = 1; i <= n; i++) {
        scanf("%d", &a[i]);
    }
    for(int i = 0; i < m; i++) {
        scanf("%d%d%d", &q[i].se.se, &q[i].fi, &q[i].se.fi);
    }
    cin>>k;
    int l = 1, r = n, ans;
    while(l<=r) {
        int mid = l+r>>1;
        check(mid, n, m);
        if(query(k, k, 1, n, 1)) {
            ans = mid;
            l = mid+1;
        } else {
            r = mid-1;
        }
    }
    cout<<ans<<endl;
}
int main()
{
    int t;
    cin>>t;
    while(t--) {
        solve();
    }
    return 0;
}

posted on 2016-03-28 19:29  yohaha  阅读(153)  评论(0编辑  收藏  举报

导航