二维偏序问题应用(二维数点)

已于 2022-05-12 16:48:08 修改
阅读量415 收藏
点赞数 1
分类专栏： 数学 文章标签： c++ 算法 数据结构
于 2022-05-12 16:07:58 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_52093121/article/details/124730910
版权

二维偏序问题的数学解答
将二维数点问题转移为二维偏序问题来做，CDQ分治主席树也可以做，但是本人太菜，只会树状数组做法：
题目集

无需离散化写法：(注意询问数组开四倍，树状数组0没有意义+1)
例题:园丁的烦恼

#include <bits/stdc++.h>
using namespace std;
#define x first
#define y second
# define rep(i,be,en) for(int i=be;i<=en;i++)
# define pre(i,be,en) for(int i=be;i>=en;i--)
#define ll long long
#define endl "\n"
#define LOCAL
#define pb push_back
typedef pair<ll, ll> PII;
#define eb emplace_back
#define sp(i) setprecision(i)
const int N = 5e5 + 10, INF = 0x3f3f3f3f;

int n;
PII a[N];
int tr[10000005];
int m;
ll ans[N];

struct ques
{
	int x, y, id, type;
	ques(int _x = 0, int _y = 0, int _id = 0, int _type = 0)
	{
		x = _x;
		y = _y;
		id = _id;
		type = _type;
	}
	bool operator<(const ques &t)const
	{
		if (x == t.x)
			{
				return y < t.y;
			}
		else
			{
				return x < t.x;
			}
	}
} q[N << 2];
int cnt;



inline int lowbit(int x)
{
	return x & -x;
}

void modify(int x, int t)
{
	for (int i = x; i < 10000005; i += lowbit(i))
		{
			tr[i] += t;
		}
}


int  query(int x)
{
	ll sum = 0;
	for (int i = x; i; i -= lowbit(i))
		{
			sum += tr[i];
		}
	return sum;
}

void solve()
{
	cnt = 0;
	n = read();
	m = read();
	for (int i = 1; i <= n; i++)
		{
			a[i].x = read() + 1;
			a[i].y = read() + 1;
		}

	for (int i = 1; i <= m; i++)
		{
			int x1, x2, y1, y2;
			x1 = read() + 1;
			y1 = read() + 1;
			x2 = read() + 1;
			y2 = read() + 1;
			q[++cnt] = ques(x2, y2, i, 1);
			q[++cnt] = ques(x2, y1 - 1, i, -1);
			q[++cnt] = ques(x1 - 1, y2, i, -1);
			q[++cnt] = ques(x1 - 1, y1 - 1, i, 1);
		}
	sort(a + 1, a + 1 + n);
	sort(q + 1, q + 1 + cnt);
	int now = 1;
	for (int i = 1; i <= cnt; i++)
		{
			while (a[now].x <= q[i].x && now <= n)
				{
					modify(a[now].y, 1);
					now++;
				}
			ans[q[i].id] += 1LL * q[i].type * query(q[i].y);
		}
	for (int i = 1; i <= m; i++)
		{
			write(ans[i]);
			puts("");
		}
}




signed main()
{
	//#ifdef LOCAL
	//freopen("data.in.txt","r",stdin);
	//freopen("data.out.txt","w",stdout);
	//#endif
	int __ = 1;
	//cin >> __;
	while (__--)
		{
			solve();
		}
	return 0;
}

需要离散化写法：
代码源oj 数数
这时要将两个数组放到一起
不需要加一，离散化会配好的：

/*********************************************************************
    程序名:
    版权:
    作者:	Joecai
    日期: 2022-05-12 14:44
    说明:
*********************************************************************/

#include <bits/stdc++.h>
using namespace std;
#define x first
#define y second
# define rep(i,be,en) for(int i=be;i<=en;i++)
# define pre(i,be,en) for(int i=be;i>=en;i--)
#define ll long long
#define endl "\n"
#define LOCAL
#define pb push_back
#define int	long long
typedef pair<ll, ll> PII;
#define eb emplace_back
#define sp(i) setprecision(i)
const int N = 1e6 + 10, INF = 0x3f3f3f3f;
int n, m;
int tr[N];
int t[N];
int ans[N];

inline int read()
{
	int X = 0;
	bool flag = 1;
	char ch = getchar();
	while (ch < '0' || ch > '9')
		{
			if (ch == '-')
				flag = 0;
			ch = getchar();
		}
	while (ch >= '0' && ch <= '9')
		{
			X = (X << 1) + (X << 3) + ch - '0';
			ch = getchar();
		}
	if (flag)
		return X;
	return ~(X - 1);
}

inline void write(int X)
{
	if (X < 0)
		{
			X = ~(X - 1);
			putchar('-');
		}
	if (X > 9)
		write(X / 10);
	putchar(X % 10 + '0');
}

struct node
{
	int x;
	int y;
	int id;
	int type;
	node(int _x = 0, int _y = 0, int _id = 0, int _type = 0)
	{
		x = _x;
		y = _y;
		id = _id;
		type = _type;
	}

	bool operator<(const node &w)const
	{
		if (x == w.x)
			{
				if (y == w.y)
					{
						return id < w.id;
					}
				else
					return y < w.y;
			}
		else
			return x < w.x;
	}
} q[N];

inline int lowbit(int x)
{
	return x & (-x);
}

void up(int x, int s)
{
	for (int i = x; i < N; i += lowbit(i))
		{
			tr[i] += s;
		}
}

int query(int x)
{
	int sum = 0;
	for (int i = x; i; i -= lowbit(i))
		{
			sum += tr[i];
		}
	return sum;
}


void solve()
{
	int cnt = 0;
	memset(tr, 0, sizeof(tr));
	memset(ans, 0, sizeof(ans));
	cin >> n >> m;
	for (int i = 1; i <= n; i++)
		{
			int x;
			x = read();
			q[++cnt].x = i;
			q[cnt].y = x;
			q[cnt].type = 0;
			q[cnt].id = 0;
		}

	for (int i = 1; i <= m; i++)
		{
			int x, y, z;
			x = read();
			y = read();
			z = read();
			q[++cnt] = node(y, z, i, 1);
			q[++cnt] = node(x - 1, -1, i, 1);
			q[++cnt] = node(x - 1, z, i, -1);
			q[++cnt] = node(y, -1, i, -1);
		}
	for (int i = 1; i <= cnt; i++)
		t[i] = q[i].y;
	sort(q + 1, q + 1 + cnt);
	sort(t + 1, t + 1 + cnt);
	int te = unique(t + 1, t + 1 + cnt) - t - 1;
	for (int i = 1; i <= cnt; i++)
		{
			q[i].y = lower_bound(t + 1, t + 1 + te, q[i].y) - t;
		}
	for (int i = 1; i <= cnt; i++)
		{
			if (q[i].id) //查询的点
				{
					ans[q[i].id] += q[i].type * query(q[i].y);
				}
			else
				{
					up(q[i].y, 1);
				}
		}
	for (int i = 1; i <= m; i++)
		{
			write(ans[i]);
			printf(" ");
		}
	cout << endl;
}




signed main()
{

	//#ifdef LOCAL
	//freopen("data.in.txt","r",stdin);
	//freopen("data.out.txt","w",stdout);
	//#endif
	int __ = 1;
	cin >> __;
	while (__--)
		{
			solve();
		}
	return 0;
}

1

/*********************************************************************
    程序名:
    版权:
    作者:	Joecai
    日期: 2022-05-12 16:20
    说明:
*********************************************************************/

#include <bits/stdc++.h>
using namespace std;
#define x first
#define y second
# define rep(i,be,en) for(int i=be;i<=en;i++)
# define pre(i,be,en) for(int i=be;i>=en;i--)
#define ll long long
#define endl "\n"
#define LOCAL
#define pb push_back
typedef pair<ll, ll> PII;
#define eb emplace_back
#define sp(i) setprecision(i)
const int N = 1e6 + 10, INF = 0x3f3f3f3f;


int n;
int tr[N];

inline int lowbit(int x)
{
	return x & (-x);
}

void up(int x, int s)
{
	for (int i = x; i < N; i += lowbit(i))
		{
			tr[i] += s;
		}
}

int query(int x)
{
	int sum = 0;
	for (int i = x; i; i -= lowbit(i))
		{
			sum += tr[i];
		}
	return sum;
}


struct node
{
	int x;
	int y;
	int id;
	int type;
	node(int _x = 0, int _y = 0, int _id = 0, int _type = 0)
	{
		x = _x;
		y = _y;
		id = _id;
		type = _type;
	}
	bool operator<(const node &w)const
	{
		if (x == w.x)
			{
				return y < w.y;
			}
		else
			return x < w.x;
	}
} q[N];

int t[N];
int ans[N];

void solve()
{
	int n;
	cin >> n;
	int cnt = 0;
	for (int i = 1; i <= n; i++)
		{
			int x, y;
			cin >> x >> y;
			q[++cnt].x = x;
			q[cnt].y = y;
			q[++cnt] = node(y, y, i, 1);
			q[++cnt] = node(x - 1, x - 1, i, 1);
			q[++cnt] = node(x - 1, y, i, -1);
			q[++cnt] = node(y, x - 1, i, -1);
		}
	sort(q + 1, q + 1 + cnt);
	for (int i = 1; i <= cnt; i++)
		t[i] = q[i].y;
	sort(t + 1, t + 1 + cnt);
	int te = unique(t + 1, t + 1 + cnt) - t - 1;
	for (int i = 1; i <= cnt; i++)
		{
			q[i].y = lower_bound(t + 1, t + 1 + te, q[i].y) - t;
		}

	for (int i = 1; i <= cnt; i++)
		{
			if (q[i].id)
				{
					ans[q[i].id] += q[i].type * query(q[i].y);
				}
			else
				up(q[i].y, 1);
		}
	for (int i = 1; i <= n; i++)
		{
			cout << ans[i] - 1 << endl;
		}
}




signed main()
{
	std::ios::sync_with_stdio(false);
	std::cin.tie(nullptr);
	//#ifdef LOCAL
	//freopen("data.in.txt","r",stdin);
	//freopen("data.out.txt","w",stdout);
	//#endif
	int __ = 1;
	//cin>>__;
	while (__--)
		{
			solve();
		}
	return 0;
}
蔡珏
关注
  • 1
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
二维数点问题
weixin_30423977的博客
04-16 2382
二维数点问题: 给定平面上的$n$个点$(x_i,y_i)$, 权值$f(x_i,y_i)$, $m$次矩形查询$\sum\limits_{\substack{a \le i \le b \\ c \le j \le d}}f(i,j)$ 以下记$S(a,b)=\sum\limits_{\substack{i \le a \\ j \le b}}f(i,j)$ 一般解法: 先...
二维偏序入门
gml1999的博客
10-29 1494
一维sort； 多维cdq学不会； 那就学二维吧。 二维偏序：二维偏序就可以看作是一个坐标轴中的点（x,y）就是查找有多少个点x,y都比该点小。 一般都是排序一维，另一维树状数组更新维护 例题1：http://poj.org/problem?id=2352 POJ 2352 Stars 题意：给你N个坐标然后问你每个坐标的等级（等级的意思就是X,Y都小于当前（X,Y）的坐...
二维偏序专题
最新发布
Unlimitedz的博客
12-02 206
二维偏序专题
二维偏序/数点——学习笔记
Zenith_Habitant的博客
11-23 1234
二维偏序/数点 对于可离线的二维偏序问题一般都是先定一维的顺序，然后利用树状数组维护求解 应用 一. 题目链接：https://www.luogu.com.cn/problem/P2163 解题思路：二维数点的经典模板题，离线询问之后利用二维前缀和方法用树状数组维护即可 //#define LOCAL #include <bits/stdc++.h> using namespace std; #define ll long long #define mem(a, b) memset(
二维数点问题 (从线段树到CDQ分治)
蒟蒻的最后的倔强
08-21 2407
http://acm.hdu.edu.cn/showproblem.php?pid=1541 #include<bits/stdc++.h> using namespace std; typedef long long ll; const int maxn = 1e5 + 10; int n, cnt; struct node { int x, y; bool operator ...
【数据结构】二维数点/二维偏序
分享XCPC常用算法和相关技巧
03-29 5678
二维数点又称二维偏序，它是这样一类问题，给出一个二维平面內的若干个点，多次询问某个矩形区域內包含多少个点（边界也算）。又或者，给一个长为n的序列，多次询问区间[l,r]中值在[x,y]内的元素个数。能解决这一问题的数据结构较多，包括树状数组/线段树、K-D Tree、可持久化线段树等，运用CDQ分治可解决更高维的偏序问题。以下着重讲解扫描线思想+树状数组解决二维偏序问题的方法。
二维平面数点
di6499的博客
08-19 791
目录 CDQ分治 BZOJ1935 园丁的烦恼 BZOJ1176 Mokia Distance(2019年牛客多校第八场D题+CDQ+树状数组) 线段树 PUBG 1V3 CSUST2015 题意 思路 ...
二维偏序问题
z2656453168的博客
10-06 1220
二维偏序问题 偏序问题的概念 设RRR是集合AAA上的一个二元关系，若RRR满足： 自反性：对任意x∈Ax∈Ax∈A，有xRxxRxxRx 反对称性（即反对称关系）：对任意x,y∈Ax,y∈Ax,y∈A，若xRyxRyxRy，且yRxyRxyRx，则x=yx=yx=y 传递性：对任意x,y,z∈Ax, y,z∈Ax,y,z∈A，若xRyxRyxRy，且yRzyRzyRz，则xRzxRzxRz 则称RRR为AAA上的偏序关系，通常记作≼≼≼。注意这里的≼≼≼不必是指一般意义上的“小于或等于” 若然有x≼y
扫描线 二维数点
William_Wang_的博客
07-31 285
扫描线
牛客练习赛4 - A - Laptop（二维偏序）
sugarbliss
09-25 218
题目链接：https://ac.nowcoder.com/acm/contest/16/A 思路：二维偏序问题，按照第一维排序，再用树状数组处理第二维即可。 #include <bits/stdc++.h> using namespace std; typedef long long ll; typedef unsigned long long ul; typedef pair&...
AtCoder Beginner Contest 238 F - Two Exams「dp + 二维偏序」
Suryxin's blogs
03-16 334
F - Two Exams 题目描述： n个考生，参加了两次考试，第一次的排名是Pi，第二次是Qi,现在需要选m名考生去参加活动，必须保证不能存在未被挑中的人x的第一次的和第二次的排名比选中的人y的第一次和第二次的排名都靠前，也就是不可以出现如下的情况p[x]<p[y]且q[x]<q[y],x为被选中，y被选中p[x]<p[y]且q[x]<q[y],x为被选中，y被选中p[x]<p[y]且q[x]<q[y],x为被选中，y被选中，问有多少种挑选方法，模 99824435
Stars--HDU1541（二维偏序-树状数组）
孤风的博客
09-06 793
Stars Time Limit: 2000/1000 MS (Java/Others)Memory Limit: 65536/32768 K (Java/Others)题目链接http://acm.hdu.edu.cn/showproblem.php?pid=1541 Problem Description Astronomers often examine star maps w...
【无码专区8】三角形二维数点——计数有多少个给定点落在三角形区域内
11-09 6585
计数计数计数，怎么老是计数！我擦，我怎么感觉自己比题解还行。
[JOI Open 2014 day1]Fortune Telling 2(线段树+二维数点)
domy79360的博客
02-19 214
题目 NNN张牌，每张牌有AAA、BBB面。AAA面的值为A[i]，BBB面的值为B[i]。初始AAA面朝上。 有KKK次操...
二维偏序
qq_17080419的博客
10-24 1240
The Preliminary Contest for ICPC Asia Xuzhou 2019 I. query 题意：给定一个n的排列，给定询问区间[l,r]，问在区间内有多少对倍数关系 思路：可以用前缀和的思想，处理出[1,x]位置的所有倍数关系。答案就是[1,r]的倍数关系数减去[1,l]，同时还要减去其中一个点在l左边，一个点在l右边的情况 #include <iostream&...
20200607 练习：CQOI2017day2（二维数点）（网络流）（树形dp）
hongkongreporter
06-08 181
T1 P3755 [CQOI2017]老C的任务 思路： 二维数点 离线后线段树维护 代码： #include<bits/stdc++.h> using namespace std; #define in Read() #define int long long #define lch p<<1 #define rch p<<1|1 inline int in{ int s=0,f=1;char x; for(x=getchar();x<'0'||x>'
loj#2251. 「ZJOI2017」树状数组（二维数点，树套树维护概率）
juruo? juruo!
03-14 1111
先放代码，日后更。（*2） ===========================2018.3.21UPD=========================== 题面在这里 做法 容易发现只要将query的l−1,rl−1,rl-1,r改成r+1,lr+1,lr+1,l就是正确的了。所以最后答案之和l−1,rl−1,rl-1,r这两个点的奇偶性是否相同有关。问题转化为维护两个位置的奇偶性相...
主席树处理二维数点问题 园丁的烦恼
sophilex的博客
01-11 313
不妨先把所有点的坐标按照横坐标来排序，那么对于每次询问，我们可以很快获得点的序列范围，比如二分。然后考虑如何处理纵坐标上满足条件的点的数量。假设纵坐标的范围是[a,b],而我们确定的横坐标满足范围的点的序号范围是[c,d],那么我们按照顺序建立主席树，依次取出rt[c-1]和rt[d],每次获得一个区间和即可。主席树是满足可加性的，所以我们直接将两个结果减一下就好了。二维坐标里有n个点，q次询问，每次询问给定一个矩形的左下角和右上角，问在这个矩形范围内有多少个点。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值