NOIP 模拟 $34\; \rm Rectangle$
题解 \(by\;zj\varphi\)
对于没有在同一行或同一列的情况,直接枚举右边界,左边界从大到小,用树状数组维护上下边界即可。
而对于有多个在一列或一行的情况,这些点将左右分成了几个区间,枚举上边界在哪个区间,同时维护下边界。
Code
#include<bits/stdc++.h>
#define ri register signed
#define p(i) ++i
namespace IO{
char buf[1<<21],*p1=buf,*p2=buf;
#define gc() p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?(-1):*p1++
struct nanfeng_stream{
template<typename T>inline nanfeng_stream operator>>(T &x) {
ri f=0;x=0;register char ch=gc();
while(!isdigit(ch)) f|=ch=='-',ch=gc();
while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=gc();
return x=f?-x:x,*this;
}
}cin;
}
using IO::cin;
namespace nanfeng{
#define Node(x,y) (Node){x,y}
#define FI FILE *IN
#define FO FILE *OUT
template<typename T>inline T cmax(T x,T y) {return x>y?x:y;}
template<typename T>inline T cmin(T x,T y) {return x>y?y:x;}
typedef long long ll;
static const int N=1e5+7,MOD=1e9+7;
int cut[N],Y[N],vis[N],X[N],cnt,xm,tot,n;
ll ans;
struct node{int x,y;}pnt[N];
struct Node{ll x,y;};
inline bool operator<(const node &n1,const node &n2) {return n1.x==n2.x?n1.y<n2.y:n1.x>n2.x;}
struct BIT{
#define lowbit(x) ((x)&-(x))
int c[2507],b[2507];
inline void init() {
memset(c,0,sizeof(int)*(xm+7));
memset(b,0,sizeof(int)*(xm+7));
}
inline void update(int x,int k) {for (ri i(x);i<=xm;i+=lowbit(i)) c[i]+=k,p(b[i]);}
inline Node query(int x) {
register ll res(0),nm(0);
for (ri i(x);i;i-=lowbit(i)) res+=c[i],nm+=b[i];
return Node(res,nm);
}
}B;
inline int main() {
//FI=freopen("nanfeng.in","r",stdin);
//FO=freopen("nanfeng.out","w",stdout);
cin >> n;
for (ri i(1);i<=n;p(i)) cin >> pnt[i].x >> pnt[i].y,xm=cmax(xm,pnt[i].y);
std::sort(pnt+1,pnt+n+1);
cut[p(tot)]=1;
X[tot]=pnt[1].x;
for (ri i(1);i<=n;p(i)) {
if (pnt[i].x!=pnt[i-1].x&&i!=1) Y[p(cnt)]=xm+1,cut[p(tot)]=cnt+1,X[tot]=pnt[i].x;
Y[p(cnt)]=pnt[i].y;
}
Y[p(cnt)]=xm+1;
cut[p(tot)]=cnt+1;
for (ri i(1);i<tot;p(i)) {
register ll tmp=ans;
B.init();
memset(vis,0,sizeof(int)*(xm+7));
for (ri j(cut[i]);j<cut[i+1]-1;p(j))
if (!vis[Y[j]]) B.update(Y[j],Y[j]),vis[Y[j]]=1;
for (ri j(i+1);j<tot;p(j)) {
for (ri k(cut[j]);k<cut[j+1]-1;p(k))
if (!vis[Y[k]]) B.update(Y[k],Y[k]),vis[Y[k]]=1;
ri l1(cut[i]),l2(cut[j]);
ri mn(cmax(Y[l1],Y[l2]));
while(Y[l1+1]<=mn) p(l1);
while(Y[l2+1]<=mn) p(l2);
while(l1<cut[i+1]-1&&l2<cut[j+1]-1) {
int tmp=cmin(Y[l1+1],Y[l2+1]);
register Node tp1=B.query(tmp-1),tp2=B.query(mn-1),tp3=B.query(cmin(Y[l1],Y[l2]));
mn=tmp;
ans=(ans+(ll)(X[i]-X[j])*((tp1.x-tp2.x)*tp3.y-(tp1.y-tp2.y)*tp3.x))%MOD;
if (Y[l1+1]<=mn) p(l1);
if (Y[l2+1]<=mn) p(l2);
}
}
}
printf("%lld\n",ans);
return 0;
}
}
int main() {return nanfeng::main();}
NOIP 模拟 $34\; \rm Rectangle$的更多相关文章
- NOIP 模拟 $34\; \rm Equation$
题解 \(by\;zj\varphi\) 发现每个点的权值都可以表示成 \(\rm k\pm x\). 那么对于新增的方程,\(\rm x_u+x_v=k\pm x/0\) 且 \(\rm x_u+x ...
- NOIP 模拟 $34\; \rm Merchant$
题解 \(by\;zj\varphi\) 对于选的物品,总值一定有在前一段区间递减,后一段递增的性质,那么就可以二分. check()时只递归归并大的一段,用nth_element即可 Code #i ...
- noip模拟34[惨败]
noip模拟34 solutions 我从来不为失败找借口,因为败了就是败了,没人听你诉说任何事情 今天很伤感,以来考试没考好,二来改题改半天也改不出来 这次算是炸出来了我经常范的一些错误,比如除以0 ...
- NOIP模拟 1
NOIP模拟1,到现在时间已经比较长了.. 那天是6.14,今天7.18了 //然鹅我看着最前边缺失的模拟1,还是终于忍不住把它补上,为了保持顺序2345重新发布了一遍.. # 用 户 名 ...
- 2021.5.22 noip模拟1
这场考试考得很烂 连暴力都没打好 只拿了25分,,,,,,,,好好总结 T1序列 A. 序列 题目描述 HZ每周一都要举行升旗仪式,国旗班会站成一整列整齐的向前行进. 郭神作为摄像师想要选取其中一段照 ...
- NOIP模拟3
期望得分:30+90+100=220 实际得分:30+0+10=40 T1智障错误:n*m是n行m列,硬是做成了m行n列 T2智障错误:读入三个数写了两个%d T3智障错误:数值相同不代表是同一个数 ...
- 7.22 NOIP模拟7
又是炸掉的一次考试 T1.方程的解 本次考试最容易骗分的一道题,但是由于T2花的时间太多,我竟然连a+b=c都没判..暴力掉了40分. 首先a+b=c,只有一组解. 然后是a=1,b=1,答案是c-1 ...
- 20190725 NOIP模拟8
今天起来就是虚的一批,然后7.15开始考试,整个前半个小时异常的困,然后一看题,T1一眼就看出了是KMP,但是完了,自己KMP的打法忘的一干二净,然后开始打T2,T2肝了一个tarjan点双就扔上去了 ...
- 20190902+0903合集-NOIP模拟
一直没时间写QwQ 于是补一下. Day 1 晚饭吃的有点恶心…… $1s\,2s\,5s$ 还开 -O2 ?? 有点恐怖. T1 猛的一想: 把外面设成一个点, 向入口连一条权为排队时间的边 从出口 ...
随机推荐
- 全网唯一开源java开发的支持高扩展,高性能的Mqtt集群broker!
SMQTT是一款开源的MQTT消息代理Broker, SMQTT基于Netty开发,底层采用Reactor3反应堆模型,支持单机部署,支持容器化部署,具备低延迟,高吞吐量,支持百万TCP连接,同时支持 ...
- ROS2学习之旅(15)——编写简单的服务和客户节点(C++)
当节点使用服务进行通信时,发送数据请求的节点称为客户节点,响应请求的节点称为服务节点.请求和响应的结构由.srv文件决定. 本文的例子是一个简单的整数加法系统:一个节点请求两个整数的和,另一个节点响应 ...
- 使用 Java 和 Maven (JBake) 生成静态网站
使用 JBake("mvn generate-resources")构建您的静态网站或博客.使用布局.宏和数据文件. 我们迁移了整个www.optaplanner.org网站(13 ...
- JMeter之Throughput Controller吞吐量控制器
吞吐量控制器,它是用来控制该控制器下面元件的执行次数,与控制吞吐量的功能无关.(注:用Constant Throughput Timer可以控制吞吐量tps) 作用:控制其子节点的执行次数与负载比例分 ...
- React事件绑定的方式
一.是什么 在react应用中,事件名都是用小驼峰格式进行书写,例如onclick要改写成onClick 最简单的事件绑定如下: class ShowAlert extends React.Compo ...
- 【剑指offer】42.和为S的两个数字
42.和为S的两个数字 题目描述 输入一个递增排序的数组和一个数字S,在数组中查找两个数,使得他们的和正好是S,如果有多对数字的和等于S,输出两个数的乘积最小的. 示例: 输入:[1,2,4,7,11 ...
- PAT乙级:1090危险品装箱(25分)
PAT乙级:1090危险品装箱(25分) 题干 集装箱运输货物时,我们必须特别小心,不能把不相容的货物装在一只箱子里.比如氧化剂绝对不能跟易燃液体同箱,否则很容易造成爆炸. 本题给定一张不相容物品的清 ...
- 第五篇--Chorme浏览器主页被篡改
解决方法:关闭谷歌浏览器,右击桌面快捷方式,查看属性,然后将target后面的网址删掉.并且任务栏的google打开方式,最好也把流氓网址删掉.之后就正常了.
- 【原创】case、casez和casex谁是谁
在Verilog中case语句经常用于多分支表决的结构,case后的表达式会与各分支表达式"全等"那么对应的分支会被执行.其基本结构如下: case(expression) exp ...
- include 与 Widget
include :include 引入的页面动态加载数据时,需要在当前页面中写数据加载,而不是只在引入的页面中就可以了: