[CSP-S模拟测试]:Rectangle(模拟+树状数组)
题目描述
平面上有$n$个点,第$i$个点的坐标为$X_i,Y_i$。对于其中的一个非空点集$S$,定义$f(S)$为一个最小矩形,满足:
$\bullet$覆盖$S$中所有的点(在边界上也算覆盖);
$\bullet$边与坐标轴平行。
求所有不同的$f(S)$的面积和对$10^9+7$取模的结果。两个矩形被认为是不同的,当且仅当它们顶点坐标不同。
输入格式
从文件$rectangle.in$中读入数据。
第一行一个整数$n$。
接下来$n$行,每行两个整数$X_i,Y_i$。
输出格式
输出到文件$rectangle.out$中。
一行一个整数表示答案。
样例
样例输入:
4
1 2
3 1
4 4
5 1
样例输出:
45
数据范围与提示
样例解释:
有$8$个面积大于$0$的不同矩形,以下是它们左下角和右上角的坐标:
$(1,1),(3,2);(1,1),(4,4);(1,1),(5,2);(1,1),(5,4)$
$(1,2),(4,4);(3,1),(4,4);(3,1),(5,4);(4,1),(5,4)$
数据范围:
对于所有数据,满足$2\leqslant n\leqslant 10^4,1\leqslant X_i,Y_i\leqslant 2500$,没有重复的点。
$\bullet Subtask1(13\%)$,$n\leqslant 18$。
$\bullet Subtask2(9\%)$,$n\leqslant 50$。
$\bullet Subtask3(25\%)$,$n\leqslant 300$。
$\bullet Subtask4(21\%)$,$n\leqslant 2500,X_i\neq X_j,Y_i\neq Y_j$。
$\bullet Subtask5(19\%)$,$n\leqslant 2500$。
$\bullet Subtask6(13\%)$,没有特殊的约束。
题解
先来考虑$21\%$的$X_i\neq X_j,Y_i\neq Y_j$的情况。
我们可以$n^2$枚举左右边界,那么设边界上的点为$(L,y_1)$和$(R,y_2)$。
那么只有位于$(L,R)$且纵坐标$>\max(y_1,y_2)$和$<\min(y_1,y_2)$的点才能做贡献,我们可以考虑树状数组,存储$\sum y$即可(长度是变化的,但是高度不变)。
现在来考虑一般情况,每个$L$和$R$上可能有很多的点,我们依次枚举计数即可。
但是可能会出现如下图中的情况:
显然,我们在统计答案点$1,3$和点$2,3$的贡献的时候会将紫色矩阵算重,不用担心,我们只需要将纵坐标最靠下的统计就好了。
代码实现稍繁琐。
时间复杂度:$\Theta(nm\log m)$。
期望得分:$100$分。
实际得分:$100$分。
代码时刻
#include<bits/stdc++.h>
using namespace std;
const int mod=1000000007;
int n;
int Map[2501][2501];
int tr[2][2501][2501];
bool vis[2501][2501];
long long ans;
int lowbit(int x){return x&-x;}
void add(int id,int k,int x,int w)
{
for(int i=x;i<=2500;i+=lowbit(i))
tr[id][k][i]+=w;
}
int ask(int id,int k,int x)
{
int res=0;
for(int i=x;i;i-=lowbit(i))res+=tr[id][k][i];
return res;
}
int main()
{
scanf("%d",&n);
for(int i=1;i<=n;i++)
{
int x,y;
scanf("%d%d",&x,&y);
Map[x][++Map[x][0]]=y;
}
for(int i=1;i<=2500;i++)
{
sort(Map[i]+1,Map[i]+Map[i][0]+1);
Map[i][Map[i][0]+1]=2501;
}
for(int i=1;i<=2500;i++)
{
if(!Map[i][0])continue;
for(int j=1;j<=Map[i][0];j++)
if(!vis[i][Map[i][j]])
{
vis[i][Map[i][j]]=1;
add(1,i,Map[i][j],1);
add(0,i,Map[i][j],Map[i][j]);
}
for(int j=i-1;j;j--)
{
if(!Map[j][0])continue;
int l1=1,l2=1;
for(int k=1;k<=Map[j][0];k++)
if(!vis[i][Map[j][k]])
{
vis[i][Map[j][k]]=1;
add(1,i,Map[j][k],1);
add(0,i,Map[j][k],Map[j][k]);
}
int wzc=max(Map[i][1],Map[j][1]);
while(Map[i][l1+1]<=wzc)l1++;
while(Map[j][l2+1]<=wzc)l2++;
while(l1<=Map[i][0]&&l2<=Map[j][0])
{
int flag=min(Map[i][l1+1],Map[j][l2+1]);
ans=(ans+1LL*(i-j)*((ask(0,i,flag-1)-ask(0,i,wzc-1))*ask(1,i,min(Map[i][l1],Map[j][l2]))-((ask(1,i,flag-1)-ask(1,i,wzc-1))*ask(0,i,min(Map[i][l1],Map[j][l2])))))%mod;
wzc=flag;
if(Map[i][l1+1]<=wzc)l1++;
if(Map[j][l2+1]<=wzc)l2++;
}
}
}
printf("%lld\n",ans);
return 0;
}
rp++
[CSP-S模拟测试]:Rectangle(模拟+树状数组)的更多相关文章
- [CSP-S模拟测试]:柱状图(树状数组+二分+三分)
题目描述 $WTH$获得了一个柱状图,这个柱状图一共有$N$个柱子,最开始第$i$根柱子的高度为$x_i$,他现在要将这个柱状图排成一个屋顶的形状,屋顶的定义如下:$1.$屋顶存在一个最高的柱子,假设 ...
- [CSP-S模拟测试]:影魔(树状数组+线段树合并)
题目背景 影魔,奈文摩尔,据说有着一个诗人的灵魂.事实上,他吞噬的诗人灵魂早已成千上万.千百年来,他收集了各式各样的灵魂,包括诗人.牧师.帝王.乞丐.奴隶.罪人,当然,还有英雄.每一个灵魂,都有着自己 ...
- [CSP-S模拟测试]:统计(树状数组+乱搞)
题目传送门(内部题120) 输入格式 第一行,两个正整数$n,m$. 第二行,$n$个正整数$a_1,a_2,...,a_n$,保证$1\leqslant a_i\leqslant n$,可能存在相同 ...
- poj 2299 树状数组求逆序对数+离散化
Ultra-QuickSort Time Limit: 7000MS Memory Limit: 65536K Total Submissions: 54883 Accepted: 20184 ...
- The Stream of Corning 2( 权值线段树/(树状数组+二分) )
题意: 有两种操作:1.在[l,r]上插入一条值为val的线段 2.问p位置上值第k小的线段的值(是否存在) 特别的,询问的时候l和p合起来是一个递增序列 1<=l,r<=1e9:1< ...
- 【CSP模拟赛】奇怪的队列(树状数组 &二分&贪心)
题目描述 nodgd的粉丝太多了,每天都会有很多人排队要签名. 今天有n个人排队,每个人的身高都是一个整数,且互不相同.很不巧,nodgd今天去忙别的事情去了,就只好让这些粉丝们明天再来.同时nod ...
- [CSP-S模拟测试]:序列(二分答案+树状数组)
题目传送门(内部题98) 输入格式 第一行一个整数$n$,第二行$n$个整数$a_1\sim a_n$,第三行$n$个整数$b_1\sim b_n$. 输出格式 一行一个整数表示$\max(r-l+1 ...
- [CSP-S模拟测试]:小P的单调数列(树状数组+DP)
题目描述 小$P$最近喜欢上了单调数列,他觉得单调的数列具有非常多优美的性质.经过小$P$复杂的数学推导,他计算出了一个单调增数列的艺术价值等于该数列中所有书的总和.并且以这个为基础,小$P$还可以求 ...
- [CSP-S模拟测试]:降雷皇(DP+树状数组)
题目描述 降雷皇哈蒙很喜欢雷电,他想找到神奇的电光.哈蒙有$n$条导线排成一排,每条导线有一个电阻值,神奇的电光只能从一根导线传到电阻比它大的上面,而且必须从左边向右传导,当然导线不必是连续的.哈蒙想 ...
随机推荐
- java配置详解
JAVA_HOMED:\JavaTools\Java\jdk1.7.0_80\ D:\JavaEnvironment\Java\jdk1.7.0_71D:\JavaEnvironment\Java\j ...
- [Python3 填坑] 007 多才多艺的 len()
目录 1. print( 坑的信息 ) 2. 开始填坑 (1) 总的来说 (2) 举例说明 (3) 后记 1. print( 坑的信息 ) 挖坑时间:2019/01/10 明细 坑的编码 内容 Py0 ...
- [19/05/20-星期一] CSS_css的盒子模型
一.盒子模型 <!DOCTYPE html> <html> <!-- div 标签 块级标签 主要用来网页布局, 会将其中的子元素内容作为一个独立的整体存在 默认宽度是页 ...
- OOM排除与JVM调优
仅先记录,后续整理 1. 常用命令: jstat gcutil jmap 2. 打印GC执行情况: 通过执行jinfo -flag +PrintGCDetails <pid>直接动态开启, ...
- hdu-2389.rain on your parade(二分匹配HK算法)
Rain on your Parade Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 655350/165535 K (Java/Ot ...
- SCUT - 11 - 被钦定的选手 - 质因数分解
https://scut.online/p/11 T了好多次,还想用mutimap暴力分解每个数的质因数.后来记录每个数的最小质因子过了. #include <bits/stdc++.h> ...
- neo4j 的cql 语句,增、删、改、查(条件查询)(持续更新)
前言 因为做一个比赛的项目 ,需要用到 neo4j 数据库,所以要学习其语言cql,特来整理一下他的基本语言. 整片的语句是按照 了 Neo4j 数据库自带的示例 Movie Graph 来写的. 直 ...
- IOC详解
Ioc--控制反转详解(转载 http://www.cnblogs.com/qinqinmeiren/archive/2011/04/02/2151697.html) 本文转载与百度知道,简单例子让 ...
- Q开头的类找不到,无法加载插件:com.mysema.maven:apt-maven-plugin
http://www.jspxcms.com/documentation/297.html 如果出现无法加载com.mysema.maven:apt-maven-plugin插件的情况,通常是由于ma ...
- 10-基于TMS320C6678+XC7K325T的6U CPCI Full Camera Link图像处理平台
基于TMS320C6678+XC7K325T的6U CPCI Full Camera Link图像处理平台 1.板卡概述 板卡由我公司自主研发,基于6UCPCI架构,处理板包含一片TI DSPTMS3 ...