CDQ分治 陌上花开(三维偏序)
CDQ分治或树套树可以切掉
CDQ框架:
- 先分
- 计算左边对右边的贡献
- 再和
所以这个题可以一维排序,二维CDQ,三维树状数组统计
CDQ代码
# include <stdio.h>
# include <stdlib.h>
# include <iostream>
# include <algorithm>
# include <string.h>
# define IL inline
# define RG register
# define Fill(a, b) memset(a, b, sizeof(a))
using namespace std;
typedef long long ll;
IL ll Read(){
RG char c = getchar(); RG ll x = 0, z = 1;
for(; c < '0' || c > '9'; c = getchar()) z = c == '-' ? -1 : 1;
for(; c >= '0' && c <= '9'; c = getchar()) x = (x << 1) + (x << 3) + c - '0';
return x * z;
}
const int MAXN(100010), MAXK(200010);
int n, k, t[MAXK], m;
struct Point{
int a, b, c, cnt, id, ans;
IL bool operator ==(RG Point B) const{ return a == B.a && b == B.b && c == B.c; }
} p[MAXN], q[MAXN];
IL bool Cmp1(RG Point x, RG Point y){ return x.a < y.a || (x.a == y.a && (x.b < y.b || (x.b == y.b && x.c < y.c))); }
IL bool Cmp2(RG Point x, RG Point y){ return x.b < y.b || (x.b == y.b && (x.c < y.c || (x.c == y.c && x.a < y.a))); }
IL void Add(RG int x, RG int d){ for(; x <= k; x += x & -x) t[x] += d; }
IL int Query(RG int x){ RG int cnt = 0; for(; x; x -= x & -x) cnt += t[x]; return cnt; }
IL void CDQ(RG int l, RG int r){
if(l == r) return;
RG int mid = (l + r) >> 1; CDQ(l, mid); CDQ(mid + 1, r);
sort(p + l, p + r + 1, Cmp2);//可以归并排序
for(RG int i = l; i <= r; i++)
if(p[i].id <= mid) Add(p[i].c, p[i].cnt);
else p[i].ans += Query(p[i].c);
for(RG int i = l; i <= r; i++)
if(p[i].id <= mid) Add(p[i].c, -p[i].cnt);
}
int main(RG int argc, RG char* argv[]){
n = Read(); k = Read();
for(RG int i = 1; i <= n; i++) q[i] = (Point){Read(), Read(), Read()};
sort(q + 1, q + n + 1, Cmp1);
for(RG int i = 1; i <= n; i++)
if(q[i] == p[m]) p[m].cnt++;
else p[++m] = q[i], p[m].cnt = 1, p[m].id = m;
CDQ(1, m);
sort(p + 1, p + m + 1, Cmp1);
for(RG int i = 1; i <= m; i++) t[p[i].ans + p[i].cnt - 1] += p[i].cnt;
for(RG int i = 0; i < n; i++) printf("%d\n", t[i]);
return 0;
}
附上树套树
# include <stdio.h>
# include <stdlib.h>
# include <iostream>
# include <string.h>
# include <algorithm>
# define IL inline
# define RG register
# define Fill(a, b) memset(a, b, sizeof(a))
using namespace std;
typedef long long ll;
IL ll Read(){
RG char c = getchar(); RG ll x = 0, z = 1;
for(; c > '9' || c < '0'; c = getchar()) z = c == '-' ? -1 : 1;;
for(; c >= '0' && c <= '9'; c = getchar()) x = (x << 1) + (x << 3) + c - '0';
return x * z;
}
const int MAXN(100010);
int n, k, ans[MAXN], NOW;
struct YYB{
int a, b, c;
IL bool operator <(RG YYB B) const{
if(a != B.a) return a < B.a;
if(b != B.b) return b < B.b;
return c < B.c;
}
IL bool operator ==(RG YYB B) const{
return a == B.a && b == B.b && c == B.c;
}
} yyb[MAXN];
struct SBT{
SBT *ch[2], *fa;
int size, val, cnt;
IL SBT(RG SBT *f, RG int v, RG int d){ cnt = size = d; val = v; fa = f; }
} *rt[MAXN << 1];
IL int Son(RG SBT *x){ return x -> fa -> ch[1] == x; }
IL int Size(RG SBT *x){ return x == NULL ? 0 : x -> size; }
IL void Updata(RG SBT *x){ x -> size = Size(x -> ch[0]) + Size(x -> ch[1]) + x -> cnt; }
IL void Rot(RG SBT *x){
RG int c = !Son(x); RG SBT *y = x -> fa;
y -> ch[!c] = x -> ch[c];
if(x -> ch[c] != NULL) x -> ch[c] -> fa = y;
x -> fa = y -> fa;
if(y -> fa != NULL) y -> fa -> ch[Son(y)] = x;
x -> ch[c] = y; y -> fa = x;
Updata(y);
}
IL void Splay(RG SBT *x, RG SBT *f){
while(x -> fa != f){
RG SBT *y = x -> fa;
if(y -> fa != f)
Son(x) ^ Son(y) ? Rot(x) : Rot(y);
Rot(x);
}
Updata(x); rt[NOW] = x;
}
IL void Insert(RG SBT *x, RG int v, RG int d){
if(rt[NOW] == NULL){ rt[NOW] = new SBT(NULL, v, d); return; }
while(2333){
RG int id = v > x -> val;
if(v == x -> val){ x -> size += d; x -> cnt += d; break; }
if(x -> ch[id] == NULL){
x -> ch[id] = new SBT(x, v, d);
x = x -> ch[id];
break;
}
x = x -> ch[id];
}
Splay(x, NULL);
}
IL int Calc(RG SBT *x, RG int v){
RG int cnt = 0;
while(x != NULL){
if(v >= x -> val) cnt += Size(x -> ch[0]) + x -> cnt;
if(v == x -> val) break;
x = x -> ch[v > x -> val];
}
return cnt;
}
IL void Modify(RG int id, RG int d){
for(NOW = yyb[id].b; NOW <= k; NOW += NOW & -NOW) Insert(rt[NOW], yyb[id].c, d);
}
IL int Query(RG int id){
RG int cnt = 0;
for(RG int i = yyb[id].b; i; i -= i & -i) cnt += Calc(rt[i], yyb[id].c);
return cnt;
}
int main(){
n = Read(); k = Read();
for(RG int i = 1; i <= n; i++)
yyb[i] = (YYB){Read(), Read(), Read()};
sort(yyb + 1, yyb + n + 1);
for(RG int i = 1, sum = 0; i <= n; i++){
sum++;
if(yyb[i] == yyb[i + 1]) continue;
ans[Query(i) + sum - 1] += sum;
Modify(i, sum);
sum = 0;
}
for(RG int i = 0; i < n; i++) printf("%d\n", ans[i]);
return 0;
}
CDQ分治 陌上花开(三维偏序)的更多相关文章
- 并不对劲的cdq分治解三维偏序
为了反驳隔壁很对劲的太刀流,并不对劲的片手流决定与之针锋相对,先一步发表cdq分治解三维偏序. 很对劲的太刀流在这里-> 参照一.二维偏序的方法,会发现一位偏序就是直接排序,可以看成通过排序使 ...
- cdq分治解决三维偏序
问题背景 在三维坐标系中有n个点,坐标为(xi,yi,zi). 定义一个点A比一个点B小,当且仅当xA<=xB,yA<=yB,zA<=zB.问对于每个点,有多少个点比它小.(n< ...
- 【算法学习】【洛谷】cdq分治 & P3810 三维偏序
cdq是何许人也?请参看这篇:https://wenku.baidu.com/view/3b913556fd0a79563d1e7245.html. 在这篇论文中,cdq提出了对修改/询问型问题(Mo ...
- hdu5618(cdq分治求三维偏序)
题意:给n(n<=100000)个点,坐标为(xi,yi,zi)(1<=xi,yi,zi<=100000),定义一个点A比一个点B小,当且仅当xA<=xB,yA<=yB, ...
- BZOJ 3262: 陌上花开 (cdq分治,三维偏序)
#include <iostream> #include <stdio.h> #include <algorithm> using namespace std; c ...
- SPOJ:Another Longest Increasing Subsequence Problem(CDQ分治求三维偏序)
Given a sequence of N pairs of integers, find the length of the longest increasing subsequence of it ...
- Luogu 3810 & BZOJ 3262 陌上花开/三维偏序 | CDQ分治
Luogu 3810 & BZOJ 3263 陌上花开/三维偏序 | CDQ分治 题面 \(n\)个元素,每个元素有三个值:\(a_i\), \(b_i\) 和 \(c_i\).定义一个元素的 ...
- CDQ分治-陌上花开(附典型错误及原因)
CDQ分治-陌上花开 题目大意 对于给遗传给定的序列: \[ (x,y,z)_1, (x,y,z)_2, (x,y,z)_3, \cdots, (x,y,z)_n \] 求: \[ \sum_{x_i ...
- HDU4742 CDQ分治,三维LIS
HDU4742 CDQ分治,三维LIS 传送门:http://acm.hdu.edu.cn/showproblem.php?pid=4742 题意: 每个球都有三个属性值x,y,z,要求最长的lis的 ...
- BZOJ3262 陌上花开 —— 三维偏序 CDQ分治
题目链接:https://vjudge.net/problem/HYSBZ-3262 3262: 陌上花开 Time Limit: 20 Sec Memory Limit: 256 MBSubmit ...
随机推荐
- 关于 IO的同步异步间要描述
IO在计算机中指Input/Output,也就是输入和输出.由于程序和运行时数据是在内存中驻留,由CPU这个超快的计算核心来执行,涉及到数据交换的地方,通常是磁盘.网络等,就需要IO接口. 比如你打开 ...
- Servlet3.0+springmvc5+log4j2正确的开启姿势(WebLookUp)
前言 java社区占据市场份额比较大的日志组件由log4j 1.×,到logback,再到整合后的升级版 log4j 2.×,有网友测试后据说log4j2的性能最NB.于是开始往自己的springmv ...
- ActiveMq笔记2-消息持久化
为了避免意外宕机以后丢失信息,需要做到重启后可以恢复消息队列,消息系统一般都会采用持久化机制. ActiveMQ的消息持久化机制有JDBC,AMQ,KahaDB和LevelDB, 无论使用哪种持久化方 ...
- HttpURLConnection 411错误解决
1.问题描述 Java HttpURLConnection类发送Http请求链接外网返回ResponseCode为411,对方为.NET服务器,网上查明原因可能来自以下: 1.可能是没有传http中的 ...
- 基于全志H3芯片的ARM开发环境搭建
基于全志H3芯片的ARM开发环境搭建 最近买了个友善之臂的NanoPi M1板子,又在网上申请了个NanoPi NEO板子,这两个都是基于全志H3芯片的Crotex-A7四核ARM开发板,两个板子可以 ...
- 一文解决python模块导入
python 模块导入 原理 查找是按照 sys.path 中的路径挨个扫描.若都不存在则提示error. sys.path路径第一个是当前运行脚本所在的目录,其后是PYTHONPATH(一般若步专门 ...
- MinGW下简单编译FFmpeg
2009.03.21补充:ffmpeg-0.5正式发布,地址为:[url]http://www.ffmpeg.org/releases/ffmpeg-0.5.tar.bz2[/url].修改了第7步, ...
- Anaroid WebView API详解
android提供了webView控件专门用来浏览网页 使用WebView时,在程序中装载WebView控件,可以把WebView写到布局中,也可以动态New出WebView实例,可以设置属性,比如: ...
- Sublime codeIntel 配置支持php自动提示
Sublime codeIntel 配置支持php自动提示 下载地址:https://github.com/SublimeCodeIntel/SublimeCodeIntel 安装方法:下载后放到su ...
- linux下面的fd限制
如果不考虑内存大小的限制,在linux下面,fd (即file descriptor)的数量来自2个限制(阈值).其一:是操作系统的限制.这个限制主要是在linux内核中,我们知道,用户程序的fope ...