http://www.lydsy.com/JudgeOnline/problem.php?id=3295

题意:简单明了。

思路:终于好像有点明白CDQ分治处理三维偏序了。把删除操作看作是插入操作,那么可以按照插入的时间顺序看作是一维x,插入的数在原本序列的下标是一维y,插入的数本身是一维z。那么问题可以转化成每插入一个数(xx,yy,zz),求有多少个数(x,y,z)使得 x < xx,y < yy,z > zz 。一开始先对 x 进行排序,然后进行CDQ分治。这样可以干掉一维,保证随着时间递增。在分治的时候,通过标记判断那一个点属于左半区间还是右半区间,然后对 y 进行排序。如果在左半区间,那么它的 x 必定是小于 右半区间的,它所修改的结果会影响右半区间的查询,因此要去更新左半区间的元素。因为 y 是升序的,那么正着查询大于该点的 z 值的个数,就是查询可以满足 y < yy, z > zz 的条件的个数了。反着查询小于该点的 z 值的个数,即满足 y > yy, z < zz 的条件的个数。这样就可以找全插入一个数对整个数组产生的逆序对的个数了。

 #include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long LL;
#define N 100010
struct node {
int x, y, z, f;
node () {}
node (int x, int y, int z) : x(x), y(y), z(z) {}
} p[N], s[N]; int bit[N], gap, num[N], Hash[N];
LL ans[N]; bool cmpx(const node &a, const node &b) {
return a.x < b.x;
}
bool cmpy(const node &a, const node &b) {
if(a.y == b.y) return a.z < b.z;
return a.y < b.y;
} int lowbit(int x) { return x & (-x); } LL query(int x) {
LL ans = ;
while(x) { ans += bit[x]; x -= lowbit(x); }
return ans;
} void update(int x, int w) {
while(x <= gap) { bit[x] += w; x += lowbit(x); }
} void CDQ(int l, int r) {
if(l == r) return ;
int m = (l + r) >> , cnt = ;
CDQ(l, m); CDQ(m + , r);
for(int i = l; i <= m; i++) s[++cnt] = p[i], s[cnt].f = ; // 在左半部分
for(int i = m + ; i <= r; i++) s[++cnt] = p[i], s[cnt].f = ; // 在右半部分
sort(s + , s + + cnt, cmpy); // 根据y排序
for(int i = ; i <= cnt; i++) { // 正着扫
if(!s[i].f) update(s[i].z, ); // 左半部分对右半部分的查询有影响因此更新
else ans[s[i].x] += query(gap) - query(s[i].z); // 在[m,r]区间查询大于它的z的数量
}
for(int i = ; i <= cnt; i++) if(!s[i].f) update(s[i].z, -);
for(int i = cnt; i >= ; i--) { // 逆着扫
if(!s[i].f) update(s[i].z, );
else ans[s[i].x] += query(s[i].z); // 在[m,r]区间查询小于它的z的数量
}
for(int i = ; i <= cnt; i++) if(!s[i].f) update(s[i].z, -);
} int main() {
int n, m;
while(~scanf("%d%d", &n, &m)) {
int a, cnt = ;
gap = ;
for(int i = ; i <= n; i++) {
scanf("%d", &num[i]);
Hash[num[i]] = i;
p[i] = node(, i, num[i]);
if(num[i] > gap) gap = num[i];
}
for(int i = ; i <= m; i++) {
scanf("%d", &a);
p[Hash[a]].x = n - i + ;
}
for(int i = ; i <= n; i++)
if(p[i].x == ) p[i].x = ++cnt;
sort(p + , p + + n, cmpx);
memset(bit, , sizeof(bit));
CDQ(, n);
LL res = ;
for(int i = ; i <= n; i++) res += ans[i];
for(int i = n; i > n - m; i--) {
printf("%lld\n", res);
res -= ans[i];
}
}
return ;
}

BZOJ 3295:[Cqoi2011]动态逆序对(三维偏序 CDQ分治+树状数组)的更多相关文章

  1. BZOJ 3295 [CQOI2011]动态逆序对 (三维偏序CDQ+树状数组)

    题目大意: 题面传送门 还是一道三维偏序题 每次操作都可以看成这样一个三元组 $<x,w,t>$ ,操作的位置,权值,修改时间 一开始的序列看成n次插入操作 我们先求出不删除时的逆序对总数 ...

  2. BZOJ 3295: [Cqoi2011]动态逆序对

    3295: [Cqoi2011]动态逆序对 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 3865  Solved: 1298[Submit][Sta ...

  3. Bzoj 3295: [Cqoi2011]动态逆序对 分块,树状数组,逆序对

    3295: [Cqoi2011]动态逆序对 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 2886  Solved: 924[Submit][Stat ...

  4. 【刷题】BZOJ 3295 [Cqoi2011]动态逆序对

    Description 对于序列A,它的逆序对数定义为满足i<j,且Ai>Aj的数对(i,j)的个数.给1到n的一个排列,按照某种顺序依次删除m个元素,你的任务是在每次删除一个元素之前统计 ...

  5. BZOJ 3295: [Cqoi2011]动态逆序对 [CDQ分治]

    RT 传送门 首先可以看成倒着插入,求逆序对数 每个数分配时间(注意每个数都要一个时间)$t$,$x$位置,$y$数值 $CDQ(l,r)$时归并排序$x$ 然后用$[l,mid]$的加入更新$[mi ...

  6. bzoj 3295 [Cqoi2011]动态逆序对(cdq分治,BIT)

    [题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=3295 [题意] n个元素依次删除m个元素,求删除元素之前序列有多少个逆序对. [思路] ...

  7. bzoj 3295: [Cqoi2011]动态逆序对(树套树 or CDQ分治)

    Description 对于序列A,它的逆序对数定义为满足i<j,且Ai>Aj的数对(i,j)的个数.给1到n的一个排列,按照某种顺序依次删除m个元素,你的任务是在每次删除一个元素之前统计 ...

  8. BZOJ 3295 [Cqoi2011]动态逆序对 ——CDQ分治

    时间.位置.数字为三个属性. 排序时间,CDQ位置,树状数组处理数字即可. #include <cstdio> #include <cstring> #include < ...

  9. 2018.07.01 BZOJ3295: [Cqoi2011]动态逆序对(带修主席树)

    3295: [Cqoi2011]动态逆序对 **Time Limit: 10 Sec Memory Limit: 128 MB Description 对于序列A,它的逆序对数定义为满足i<j& ...

随机推荐

  1. PHP模拟单链表的数据结构

    <?php /*** * 单链表 */ //节点,下标,节点名称,下一个节点的地址 class Node { public $id; public $name; public $next; pu ...

  2. C#原子性运算 interlocked.compareExchanged

    缘起: 假设有一个类myClass, myclass里有一个count属性. 在多线程的环境下 每个线程中 直接使用count++,  如果两个线程并行执行时, 两个线程中的一个的结果会被覆掉, 非线 ...

  3. SQLServer 可更新订阅数据在线架构更改(增加字段)方案

    原文:SQLServer 可更新订阅数据在线架构更改(增加字段)方案 之前一直查找冲突发布和订阅数据不一致的原因,后来发现多少数据库升级引起,因为一直以来都是在发布数据库增加字段,订阅也会自动同步.在 ...

  4. C# Lambda表达式Contains方法 like

    原文:Lambda表达式Contains方法 like 1.使用Contains方法的必备条件: Contains等价于SQL中的like语句.不过Contains只针对于字符串(string)类型的 ...

  5. mysql三种修改密码的方式

    [root@MySQL ~]# mysqladmin -uroot -proot -S /data/3307/mysql.sock password '123'; 其中-p是现在的密码,passwor ...

  6. Android零基础入门第89节:Fragment回退栈及弹出方法

    在上一期分享的文章末尾留了一个课后作业,有去思考如何解决吗?如果已经会了那么恭喜你,如果还不会也没关系,本期一起来学习. 一.回退栈 在前面两期的示例中,当我们完成一些操作后,如果想要回到操作之前的状 ...

  7. Android零基础入门第85节:Fragment使用起来非常简单

    Fragment创建完成后并不能单独使用,还需要将Fragment加载到Activity中,在Activity中添加Fragment的方式有两种:静态加载和动态加载,接下来分别进行学习. 一.静态加载 ...

  8. Linux --- 程序后台运行的几种方法

    有时候我们运行一个程序,耗时比较长,所以在快下班的时候或是网络不稳定的时候就比较抓狂. 今天分享几个我在工作中用到的把程序放在后台运行的方法. nohup $ nohup --h Usage: noh ...

  9. Delphi 项目失败的总结

    随着项目的失败,这些天一直在总结失败的原因,到底是为什么? 一.技术层面        1.少用指针类型,多用类.            虽然指针类型能有效的节约内存和加快运行速度,但指针远没有类来得 ...

  10. 闰平年简介及计算过程描述 - Java代码实现

    import java.util.Scanner; /** * @author Shelwin Wei * 分析过程请参照<闰平年简介及计算过程描述>,网址 http://www.cnbl ...