题解 \(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$的更多相关文章

  1. NOIP 模拟 $34\; \rm Equation$

    题解 \(by\;zj\varphi\) 发现每个点的权值都可以表示成 \(\rm k\pm x\). 那么对于新增的方程,\(\rm x_u+x_v=k\pm x/0\) 且 \(\rm x_u+x ...

  2. NOIP 模拟 $34\; \rm Merchant$

    题解 \(by\;zj\varphi\) 对于选的物品,总值一定有在前一段区间递减,后一段递增的性质,那么就可以二分. check()时只递归归并大的一段,用nth_element即可 Code #i ...

  3. noip模拟34[惨败]

    noip模拟34 solutions 我从来不为失败找借口,因为败了就是败了,没人听你诉说任何事情 今天很伤感,以来考试没考好,二来改题改半天也改不出来 这次算是炸出来了我经常范的一些错误,比如除以0 ...

  4. NOIP模拟 1

    NOIP模拟1,到现在时间已经比较长了.. 那天是6.14,今天7.18了 //然鹅我看着最前边缺失的模拟1,还是终于忍不住把它补上,为了保持顺序2345重新发布了一遍.. #   用  户  名   ...

  5. 2021.5.22 noip模拟1

    这场考试考得很烂 连暴力都没打好 只拿了25分,,,,,,,,好好总结 T1序列 A. 序列 题目描述 HZ每周一都要举行升旗仪式,国旗班会站成一整列整齐的向前行进. 郭神作为摄像师想要选取其中一段照 ...

  6. NOIP模拟3

    期望得分:30+90+100=220 实际得分:30+0+10=40 T1智障错误:n*m是n行m列,硬是做成了m行n列 T2智障错误:读入三个数写了两个%d T3智障错误:数值相同不代表是同一个数 ...

  7. 7.22 NOIP模拟7

    又是炸掉的一次考试 T1.方程的解 本次考试最容易骗分的一道题,但是由于T2花的时间太多,我竟然连a+b=c都没判..暴力掉了40分. 首先a+b=c,只有一组解. 然后是a=1,b=1,答案是c-1 ...

  8. 20190725 NOIP模拟8

    今天起来就是虚的一批,然后7.15开始考试,整个前半个小时异常的困,然后一看题,T1一眼就看出了是KMP,但是完了,自己KMP的打法忘的一干二净,然后开始打T2,T2肝了一个tarjan点双就扔上去了 ...

  9. 20190902+0903合集-NOIP模拟

    一直没时间写QwQ 于是补一下. Day 1 晚饭吃的有点恶心…… $1s\,2s\,5s$ 还开 -O2 ?? 有点恐怖. T1 猛的一想: 把外面设成一个点, 向入口连一条权为排队时间的边 从出口 ...

随机推荐

  1. Vue3 + TypeScript 开发实践总结

    前言 迟来的Vue3文章,其实早在今年3月份时就把Vue3过了一遍.在去年年末又把 TypeScript 重新学了一遍,为了上 Vue3 的车,更好的开车.在上家公司4月份时,上级领导分配了一个内部的 ...

  2. ARTS第六周

    第六周.后期补完,太忙了. 1.Algorithm:每周至少做一个 leetcode 的算法题2.Review:阅读并点评至少一篇英文技术文章3.Tip:学习至少一个技术技巧4.Share:分享一篇有 ...

  3. Django基础012--接口开发

    全局参数(get,post,put,delete) 接口:/api/parameter GET:获取全局参数的所有数据 POST:创建全局参数 PUT:更新全局参数 DELETE:删除全局参数 1.创 ...

  4. etcd的使用

    etcd的使用 什么是etcd etcd的特点 etcd的应用场景 服务注册与发现 消息发布和订阅 负载均衡 分布式通知与协调 分布式锁 分布式队列 集群监控与Leader竞选 参考 etcd的使用 ...

  5. vue2.x移动端ui框架选型

    前言 最近公司准备做移动端spa项目,需要选一个ui框架.优先考虑谷歌Material Design设计风格.针对市面上的框架进行了一次调研,简单总结如下. 选型原则:1. 优先考虑md风格. 2. ...

  6. python编程面试题

    # 实现需求为 注册.登录.查看昵称的功能 # def userN(): #     username = input("请输入账号: \n") #     password =  ...

  7. MySQL触发器笔记

    当操作了某张数据表时,希望同时触发一些动作或行为,就可以使用触发器完成. 当操作微博表时,同时生成一条日志记录 -- 插入时触发 create trigger tri_weiboAdd after i ...

  8. GeoServer Rest服务启动匿名认证的配置方法

    GeoServer Rest服务数据默认需要进行用户名.密码的认证,如不需进行该认证,则启动匿名认证即可,配置方式如下(针对war包发布的GeoServer应用): 在GeoServer war包的解 ...

  9. ESLint自用规则

    官方文档地址 rules: { // allow async-await 'generator-star-spacing': 'off', // allow debugger during devel ...

  10. oracle Group by 分组查询后,分页

    ------------恢复内容开始------------ 1.分页查询 select count(*) times,title from menulog group by title order ...