UVALive 5871 Arnooks's Defensive Line (cdq分治)

最新推荐文章于 2018-03-21 20:59:17 发布

yamiedie_

最新推荐文章于 2018-03-21 20:59:17 发布

阅读量621

点赞数

本文链接：https://blog.csdn.net/u013654696/article/details/46229135

版权

题意：

有两种操作，一个是插入一个[l,r]的区间，另一个是询问一个区间[l,r]，前面有多少个区间完全包含了这个区间。

思路：

首先离散化。。

很容易想到的一个搞法是线段树套平衡树或者树状数组套平衡树，也就是对于插入的区间，把[1,l]都插入一个r，然后对询问的区间就是看[1,l]里面有多少个数>=r，然而这样的空间复杂度是nlognlogn，n=50w明显要跪了。。

然后就是对序列进行分治。

对一个操作序列[a,b], 设m=(a+b)/2，先递归处理[a,m], [m+1,b]，然后考虑[a,m]里面的插入区间的操作对[m+1,b]里面询问的影响。也就是转化成静态的问题。可以先把区间按照左端点排序，然后从左向右扫，扫到一个插入区间用树状数组就在r处+1，询问就是答案加上sum(n)-sum(r-1)。

#include <iostream>
#include <cstring>
#include <string>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstdlib>
#include <stack>
#include <queue>
#include <map>
#include <vector>
#include <cstdlib>
using namespace std;

#define LL long long
#define N 500020
#define M 200200
#define eps 1e-8
#define MP make_pair
#define Pi acos(-1.0)
#pragma comment(linker, "/STACK:1024000000,1024000000")
#define inf 0x3f3f3f3f
#define ls (p[i].ch[0])
#define rs (p[i].ch[1])
#define Ls (i << 1)
#define Rs (Ls | 1)
#define md ((ll + rr) >> 1)
#define lson ll, md, Ls
#define rson md + 1, rr, Rs

struct node {
	int l, r, t, id;
	bool operator < (const node &b) const {
		return l < b.l;
	}
}p[N], b[N];
int s[N*2];

int san[N*2], cnt;

int ans[N];

int haxi(int x) {
    return lower_bound(san + 1, san + cnt + 1, x) - san;
}
void add(int x, int v) {
	while(x) {
		s[x] += v;
		x -= x & -x;
	}
}
int query(int x) {
	int ret = 0;
	while(x <= cnt) {
		ret += s[x];
		x += x & -x;
	}
	return ret;
}
void solve(int l, int r) {
	if(l >= r) return;
	int m = (l + r) / 2;
	solve(l, m);
	solve(m + 1, r);
	int c = 0;
	for(int i = l; i <= r; ++i)
		b[++c] = p[i];
	sort(b + 1, b + c + 1);

	for(int i = 1; i <= c; ++i) {
		int j = i;
		while(j <= c && b[j].l == b[i].l) ++j;
		--j;

		for(int k = i; k <= j; ++k) {
			if(b[k].t == 1 && b[k].id <= m) {
				add(b[k].r, 1);
			}
		}
		for(int k = i; k <= j; ++k) {
			if(b[k].t == 2 && b[k].id > m) {
				ans[b[k].id] += query(b[k].r);
			}
		}
		i = j;
	}
	for(int i = 1; i <= c; ++i) {
		if(b[i].t == 1 && b[i].id <= m)
			add(b[i].r, -1);
	}
}

int main() {
    //freopen("tt.txt", "r", stdin);
    int n;
    while(scanf("%d", &n) != EOF) {

        cnt = 0;
        for(int i = 1; i <= n; ++i){
            char op[4];
            int l, r;
            scanf("%s%d%d", op, &l, &r);
            if(op[0] == '+') p[i].t = 1;
            else p[i].t = 2;
			p[i].id = i;
            p[i].l = l, p[i].r = r;
            san[++cnt] = l;
            san[++cnt] = r;
        }
		memset(ans, 0, sizeof ans);
		memset(s, 0, sizeof s);
        sort(san + 1, san + cnt + 1);
        cnt = unique(san + 1, san + cnt + 1) - san - 1;
		for(int i = 1; i <= n; ++i)
			p[i].l = haxi(p[i].l), p[i].r = haxi(p[i].r);
		solve(1, n);
		for(int i = 1; i <= n; ++i)
			if(p[i].t == 2)
				printf("%d\n", ans[i]);
		return 0;
    }
    return 0;
}

搜索

复制