原文链接https://www.cnblogs.com/zhouzhendong/p/NowCoder-2018-Summer-Round4-E.html

题目传送门 - https://www.nowcoder.com/acm/contest/142/E

题意

  给定二维平面上的 $n$ 个点,第 $i$ 个点的坐标是 $(x_i,y_i)$ ,第 $i$ 个点出现的概率是 $a_i\times b_i^{-1}$ 。

  现在让你求  [满足 “存在 $i$ 使得 $x\leq x_i,y\leq y_i$ ”的点 $(x,y)$  的数量]  的期望。答案对于 $10^9+7$ 取模。

  多组数据,$n\leq 10^5,\sum n\leq 10^6, \rm{Time\ Limit = 10s}$

题解

  出题人大概是着急有事吧。他讲的太快了这题不知道说了什么。只记得我去存了一下队友比赛时候的 AC 代码然后回来发现他讲了好几题。

  下面是我们队的做法。

  首先,我们考虑补集转化: 期望点数 = 可能的点 - 不满足条件的点数的期望

  其中可能的点我们设定为 左下角为 (1,1) ,右上角为 $ (max_x,max_y)$ ,其中右上角的两个坐标值分别为输入点的两个坐标值的最大值。

  那么我们考虑一个点对于“不满足条件的点数的期望”的贡献。这个贡献显然是 $1\times $ “在这个点左上方向的所有点都不出现的概率” 。

  我们只需要用一个算法把所有的点的贡献加起来就可以了。

  首先,我们离散化一下。然后线段树扫描线,从上往下一行一行操作。一个点的"操作"就是对 $x$ 坐标在其左侧的区间乘上 它不出现的概率。这个显然可以线段树区间乘法。每操作完一行,统计一下这一行与下一行之间区域的答案,并更新总答案。

  时间复杂度 $O (n\log n)$ 。

  不进行补集转化应该也可以做吧,但是会难写一些吧。比如我一开始用树状数组写到弃疗……

代码

%:pragma GCC optimize("Ofast")
%:pragma GCC optimize("inline")
#include <bits/stdc++.h>
using namespace std;
const int N=100005,mod=1e9+7;
namespace HashTable{
int hs,Ha[N];
void clear(){hs=0;memset(Ha,0,sizeof Ha);}
void push(int x){Ha[++hs]=x;}
void HASH(){
sort(Ha+1,Ha+hs+1);
int _hs=1;
for (int i=2;i<=hs;i++)
if (Ha[i]!=Ha[i-1])
Ha[++_hs]=Ha[i];
hs=_hs;
}
int find(int x){return lower_bound(Ha+1,Ha+hs+1,x)-Ha;}
}
using namespace HashTable;
inline int read(){
int x=0;
char ch=getchar();
while (!isdigit(ch))
ch=getchar();
while (isdigit(ch))
x=(x<<1)+(x<<3)+ch-48,ch=getchar();
return x;
}
inline int Pow(int x,int y){
int ans=1;
for (;y;y>>=1,x=1LL*x*x%mod)
if (y&1)
ans=1LL*ans*x%mod;
return ans;
}
struct Point{
int x,y,p;
}a[N];
bool cmp(Point a,Point b){
if (a.y!=b.y)
return a.y>b.y;
return a.x>b.x;
}
struct Seg{
int v,add;
}t[N<<2];
inline void build(int rt,int L,int R){
t[rt].add=1;
if (L==R){
t[rt].v=Ha[L]-Ha[L-1];
return;
}
int mid=(L+R)>>1,ls=rt<<1,rs=ls|1;
build(ls,L,mid);
build(rs,mid+1,R);
t[rt].v=(t[ls].v+t[rs].v)%mod;
}
inline void pushdown(int rt){
int &v=t[rt].add;
if (v==1)
return;
int ls=rt<<1,rs=ls|1;
t[ls].v=1LL*t[ls].v*v%mod,t[ls].add=1LL*t[ls].add*v%mod;
t[rs].v=1LL*t[rs].v*v%mod,t[rs].add=1LL*t[rs].add*v%mod;
v=1;
}
inline void update(int rt,int L,int R,int x,int d){
if (R<=x){
t[rt].v=1LL*t[rt].v*d%mod;
t[rt].add=1LL*t[rt].add*d%mod;
return;
}
pushdown(rt);
int mid=(L+R)>>1,ls=rt<<1,rs=ls|1;
update(ls,L,mid,x,d);
if (x>mid)
update(rs,mid+1,R,x,d);
t[rt].v=(t[ls].v+t[rs].v)%mod;
}
inline void solve(){
int n=read();
clear();
for (int i=1;i<=n;i++){
a[i].x=read(),a[i].y=read();
push(a[i].x);
int A=read(),B=read();
a[i].p=1LL*A*Pow(B,mod-2)%mod;
}
HASH();
for (int i=1;i<=n;i++)
a[i].x=find(a[i].x);
sort(a+1,a+n+1,cmp);
build(1,1,hs);
int ans=0;
for (int i=1,j;i<=n;i=j+1){
for (j=i;j<n&&a[j+1].y==a[i].y;j++);
for (int k=i;k<=j;k++)
update(1,1,hs,a[k].x,(mod+1-a[k].p)%mod);
ans=(1LL*(a[i].y-a[j+1].y)*(Ha[hs]-t[1].v+mod)+ans)%mod;
}
printf("%d\n",ans);
}
int main(){
for (int T=read();T;T--)
solve();
return 0;
}

  

2018牛客网暑假ACM多校训练赛(第四场)E Skyline 线段树 扫描线的更多相关文章

  1. 2018牛客网暑假ACM多校训练赛(第二场)E tree 动态规划

    原文链接https://www.cnblogs.com/zhouzhendong/p/NowCoder-2018-Summer-Round2-E.html 题目传送门 - 2018牛客多校赛第二场 E ...

  2. 2018牛客网暑假ACM多校训练赛(第三场)I Expected Size of Random Convex Hull 计算几何,凸包,其他

    原文链接https://www.cnblogs.com/zhouzhendong/p/NowCoder-2018-Summer-Round3-I.html 题目传送门 - 2018牛客多校赛第三场 I ...

  3. 2018牛客网暑假ACM多校训练赛(第三场)G Coloring Tree 计数,bfs

    原文链接https://www.cnblogs.com/zhouzhendong/p/NowCoder-2018-Summer-Round3-G.html 题目传送门 - 2018牛客多校赛第三场 G ...

  4. 2018牛客网暑假ACM多校训练赛(第三场)D Encrypted String Matching 多项式 FFT

    原文链接https://www.cnblogs.com/zhouzhendong/p/NowCoder-2018-Summer-Round3-D.html 题目传送门 - 2018牛客多校赛第三场 D ...

  5. 2018牛客网暑假ACM多校训练赛(第七场)I Tree Subset Diameter 动态规划 长链剖分 线段树

    原文链接https://www.cnblogs.com/zhouzhendong/p/NowCoder-2018-Summer-Round7-I.html 题目传送门 -  https://www.n ...

  6. 2018牛客网暑假ACM多校训练赛(第六场)I Team Rocket 线段树

    原文链接https://www.cnblogs.com/zhouzhendong/p/NowCoder-2018-Summer-Round6-I.html 题目传送门 - https://www.no ...

  7. 2018牛客网暑假ACM多校训练赛(第五场)F take 树状数组,期望

    原文链接https://www.cnblogs.com/zhouzhendong/p/NowCoder-2018-Summer-Round5-F.html 题目传送门 - https://www.no ...

  8. 2018牛客网暑假ACM多校训练赛(第四场)B Interval Revisited 动态规划

    原文链接https://www.cnblogs.com/zhouzhendong/p/NowCoder-2018-Summer-Round4-B.html 题目传送门 - https://www.no ...

  9. 2018牛客网暑假ACM多校训练赛(第十场)H Rikka with Ants 类欧几里德算法

    原文链接https://www.cnblogs.com/zhouzhendong/p/NowCoder-2018-Summer-Round10-H.html 题目传送门 - https://www.n ...

随机推荐

  1. Springboot:开发时工程启动慢

    个人的解决办法是删除workspace中的histor 然后重启Eclipse,发现一切正常.

  2. 转载:UML学习(二)-----类图(silent)

    原文:http://www.cnblogs.com/huiy/p/8552607.html 1.什么是类图 类图(Class diagram)主要用于描述系统的结构化设计.类图也是最常用的UML图,用 ...

  3. Bootstrap -- 模态框实现拖拽移动

    ### 这里实现这个效果 需要引入 jquery-ui.min.js类库 jquery-ui.min.css样式 使用它提供的draggable()方法实现 ### 菜鸟教程 http://www.r ...

  4. Shiro配置URL过滤

      常用过滤器: anon     不需要认证 authc     需要认证 user     验证通过或RememberMe登录的都可以   URL说明: /admin?=authc      表示 ...

  5. Python-爬虫-租房Ziroom

    目标站点需求分析 涉及的库 import requestsimport timeimport pymongofrom lxml import etreefrom requests.exceptions ...

  6. python-异常处理、元类

    一.异常处理 1.异常处理介绍: 异常是错误发生的信号,一旦程序出错就会产生一个异常,如果该异常没有被应用程序处理,那么该异常就会被抛出来,程序执行随之停止 2.异常通常包含三个部分 1.traceb ...

  7. STM32L476应用开发之二:模拟量数据采集

    采集模拟量数据在一台一起中是必不可少的功能.在本次实验中我们要采集的模拟量值主要包括氧气传感器的输出以及压力变送器的输出. 1硬件设计 我们需要采集数据对精度有一定的要求,而STM32L476自带AD ...

  8. Confluence 6 管理协同编辑 - 修改编辑模式

    编辑模式确定了你站点所有用户使用协同编辑的体验,这个是你对协同编辑进行启用和关闭的地方. 希望修改编辑模式: 进入  > 基本配置(General Configuration) > 协同编 ...

  9. mysql数据库1

      desc 表名;  显示表结构

  10. LoadRunner监控window系统各项指标详解

    一.监控系统时,需要监控的项 System 系统 Processor 处理器 Memory 内存 PhysicalDisk 磁盘 Server 服务器 二.指标详解 (一). PhysicalDisk ...