初学CDQ分治-NEU1702

关于CDQ分治,首先需要明白分治的复杂度。

T(n) = 2T(n/2)+O(kn), T(n) = O(knlogn)

T(n) = 2T(n/2)+O(knlogn), T(n) = O(knlog^2n)

T(n) = 2T(n/2)+O(k), T(n) = O(kn)

那么我们要处理[l, r]内的询问,我们可以分别处理[l, m]和[m+1, r]的询问,然后以较小的复杂度计算出[l, m]对[m+1, r]的贡献。

最简单的cdq就是三维偏序问题。

两点(x1, y1, z1)和(x2, y2, z2),同时满足x1 < x2, y1 < y2, z1 < z2,则前面的点小于后面的点。

首先按第一维x排序。

则处理的问题变成对于排在前面的点,统计多少个点满足y维与z维同时小于该点。

CDQ分治。

假设已处理出[l, m]与[m+1, r]。对于[m+1, r]内的所有点,我们还要统计[l, m]内有多少个点相比它更小。

对[l, r]按y维排序,对z维用树状数组统计。

扫描一遍排序后的[l, r]。

若该点在排序前属于[l, m],树状数组单点修改;否则该点在排序前属于[m+1, r],统计一次。

复杂度为O(nlognlogn)

CDQ分治算法的核心就在于:去掉时间的限制,将所有查询要求发生的时刻同化,化动态修改为静态查询 
(其实对于有些问题来说可以把某一维的限制通过排序看作时间限制然后运用CDQ分治) 
这类分治的特殊性在于分治的左右两部分的合并,作用两部分在合并的时候作用是不同的,比如,通过左半部分的影响来更新右半部分,所以分治开始前都要按照某一个关键字排序,然后利用这个顺序,考虑一个区间[l, r]的两部分间的影响。

框架为

void cdq(int l, int r){
    if(l == r) return ;
    int m = (l+r)/2;
    cdq(l, m);
    cdq(m+1, r);

    //统计[l, m]对[m+1, r]的贡献。整体排序后统计。
    sort(pp+l, pp+r+1, yzx);
    for(int i = l; i <= r; i++)
        if(pp[i].x <= m)
            add(pp[i].z, 1);
        else
            ans[ pp[i].n ] += sum(pp[i].z);
    for(int i = l; i <= r; i++)
        if(pp[i].x <= m)
            add(pp[i].z, -1);
}

NEU1702

题意:一个人的魅力值是相对于周围人来说的,如果他的颜值,内涵和智慧值同时不低于另外一个人,那么他的魅力值就会加1,给你一些人的颜值,内涵,和智慧值,请输出这些人的魅力值。

 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 typedef long long ll;
 4 const int maxn = 1e5+10;
 5 struct p{
 6     int n, x, y, z;
 7     p(){}
 8     p(int n, int x, int y, int z): n(n), x(x), y(y), z(z){}
 9     bool operator <(const p& A) const{
10         if(x != A.x) return x < A.x;
11         return y != A.y? y < A.y: z < A.z;
12     }
13     bool operator ==(const p& A) const{
14         return x == A.x&&y == A.y&&z == A.z;
15     }
16 };
17 bool yzx(p A, p B){//y z x
18     if(A.y != B.y) return A.y < B.y;
19     return A.z != B.z? A.z < B.z : A.x < B.x;
20 }
21 bool cmpn(p A, p B){//n
22     return A.n < B.n;
23 }
24 p pp[maxn];
25
26 int c[maxn], Maxn;
27 int lowbit(int x){ return x&-x;}
28 int add(int x, int d){
29     for(int i = x; i <= Maxn; i += lowbit(i))
30         c[i] += d;
31 }
32 int sum(int x){
33     int ret = 0;
34     for(int i = x; i; i -= lowbit(i))
35         ret += c[i];
36     return ret;
37 }
38
39 int ans[maxn];
40
41 void cdq(int l, int r){
42     if(l == r) return ;
43     int m = (l+r)/2;
44     cdq(l, m);
45     cdq(m+1, r);
46     sort(pp+l, pp+r+1, yzx);
47     for(int i = l; i <= r; i++)
48         if(pp[i].x <= m)
49             add(pp[i].z, 1);
50         else
51             ans[ pp[i].n ] += sum(pp[i].z);
52     for(int i = l; i <= r; i++)
53         if(pp[i].x <= m)
54             add(pp[i].z, -1);
55 }
56
57 int same[maxn];// smae[i] 表示 下标为i的ans 与 下标为same[i]相同
58
59 int main(){
60     int T; scanf("%d", &T);
61     while(T--){
62         int n;scanf("%d", &n);
63         for(int i = 0; i < n ; i++){
64             pp[i].n = i;
65             scanf("%d%d%d", &pp[i].x, &pp[i].y, &pp[i].z);
66             Maxn = max(pp[i].z, Maxn);
67         }
68         sort(pp, pp+n);//x y z
69
70         for(int i = 0; i < n; ){
71             int j = i+1;
72             while(j < n&&pp[i] == pp[j]) j++;
73             while(i < j)
74                 same[ pp[i++].n ] = pp[j-1].n;
75         }
76         for(int i = 0; i < n; i++)
77             pp[i].x = i;
78
79         memset(ans, 0, sizeof(int)*(n+5) );
80         cdq(0, n-1);
81
82         sort(pp, pp+n, cmpn);
83         for(int i = 0; i < n; i++)
84             printf("%d\n", ans[ same[ pp[i].n ] ]);
85     }
86     return 0;
87
88 }

时间: 2024-11-02 01:18:52

初学CDQ分治-NEU1702的相关文章

ACdream1157 Segments(CDQ分治 + 线段树)

题目这么说的: 进行如下3种类型操作:1)D L R(1 <= L <= R <= 1000000000) 增加一条线段[L,R]2)C i (1-base) 删除第i条增加的线段,保证每条插入线段最多插入一次,且这次删除操作一定合法3) Q L R(1 <= L <= R <= 1000000000) 查询目前存在的线段中有多少条线段完全包含[L,R]这个线段,线段X被线段Y完全包含即LY <= LX <= RX <= RY) 初学CDQ分治是看了B

CDQ分治与整体二分小结

前言 这是一波强行总结. 下面是一波瞎比比. 这几天做了几道CDQ/整体二分,感觉自己做题速度好慢啊. 很多很显然的东西都看不出来 分治分不出来 打不出来 调不对 上午下午晚上的效率完全不一样啊. 完蛋.jpg 绝望.jpg. 关于CDQ分治 CDQ分治,求的是三维偏序问题都知道的. 求法呢,就是在分治外面先把一维变成有序 然后分治下去,左边(l,mid)关于右边(mid+1,r)就不存在某一维的逆序了,所以只有两维偏序了. 这个时候来一波"树状数组求逆序对"的操作搞一下二维偏序 就可

【BZOJ3963】[WF2011]MachineWorks cdq分治+斜率优化

[BZOJ3963][WF2011]MachineWorks Description 你是任意性复杂机器公司(Arbitrarily Complex Machines, ACM)的经理,公司使用更加先进的机械设备生产先进的机器.原来的那一台生产机器已经坏了,所以你要去为公司买一台新的生产机器.你的任务是在转型期内尽可能得到更大的收益.在这段时间内,你要买卖机器,并且当机器被ACM公司拥有的时候,操控这些机器以获取利润.因为空间的限制,ACM公司在任何时候都只能最多拥有一台机器. 在转型期内,有若

BZOJ 2225 [Spoj 2371]Another Longest Increasing(CDQ分治)

[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=2225 [题目大意] 给定N个数对(xi,yi),求最长上升子序列的长度. 上升序列定义为{(xi,yi)}满足对i<j有xi<xj且yi<yj. [题解] CDQ分治,将每个区间按照a排序,用区间左边的数据来更新右边的最长上升序列, 为排除a相等但是b上升情况的误统计,在排序时加入下标作为第二关键字, 使得a相等的情况下标小的后更新. [代码] #include <cs

BZOJ 3262: 陌上花开 [CDQ分治 三维偏序]

Description 有n朵花,每朵花有三个属性:花形(s).颜色(c).气味(m),又三个整数表示.现要对每朵花评级,一朵花的级别是它拥有的美丽能超过的花的数量.定义一朵花A比另一朵花B要美丽,当且仅当Sa>=Sb,Ca>=Cb,Ma>=Mb.显然,两朵花可能有同样的属性.需要统计出评出每个等级的花的数量. Input 第一行为N,K (1 <= N <= 100,000, 1 <= K <= 200,000 ), 分别表示花的数量和最大属性值. 以下N行,每

HDU 5618:Jam&#39;s problem again(CDQ分治+树状数组处理三维偏序)

http://acm.hdu.edu.cn/showproblem.php?pid=5618 题意:-- 思路:和NEUOJ那题一样的.重新写了遍理解了一下,算作处理三维偏序的模板了. 1 #include <cstdio> 2 #include <algorithm> 3 #include <iostream> 4 #include <cstring> 5 using namespace std; 6 #define INF 0x3f3f3f3f 7 #d

BZOJ 2726: [SDOI2012]任务安排( dp + cdq分治 )

考虑每批任务对后面任务都有贡献, dp(i) = min( dp(j) + F(i) * (T(i) - T(j) + S) ) (i < j <= N)  F, T均为后缀和. 与j有关的量只有t = dp(j) - F(i) * T(j) , 我们要最小化它. dp(j)->y, T(j)->x, 那么y = F(i) * x + t, 就是给一些点和一个斜率...然后最小化截距, 显然维护下凸包就可以了. 然后因为无比坑爹的出题人....时间可以为负数, 所以要用平衡树维护(

SPOJ LIS2 Another Longest Increasing Subsequence Problem 三维偏序最长链 CDQ分治

Another Longest Increasing Subsequence Problem Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://acm.hust.edu.cn/vjudge/problem/visitOriginUrl.action?id=19929 Description Given a sequence of N pairs of integers, find the length of the longest incre

CDQ分治与整体二分总结

Cdq分治和整体二分是两个很奇妙的东西.他们都是通过离线的思想来进行优化,从而更快的求出解. 整体二分通俗的讲就是二分答案,但是它了不起的地方是一下子把所有的答案都二分出来了,从而可以一下子得出所有查询. CDQ分治通俗的讲就是二分查询.通常的做法是把所有的查询分成两半,然后通过递归先计算出左边一半的所有的查询,然后通过这些已知的左半边的值来更新右半边的值.这里,最最重要的思想是通过左半边来更新右半边.更具体一点,就是用左半边的修改来更新右半边的查询. 重要的事情说话三遍: CDQ分治就是通过左