「JSOI2014」矩形并

传送门

我们首先考虑怎么算这个期望比较好。

我们不难发现每一个矩形要和 \(n - 1\) 个矩形去交,而总共又有 \(n\) 个矩形,所以我们把矩形两两之间的交全部加起来再除以 \(n(n - 1)\) 就是答案。

至于算矩形之间的交我们可以考虑把每个矩形都视为在这个矩形范围内区间加上 \(1\) ,那么我们只需要查询一个矩形内的和 - 该矩形自身的贡献就可以算出一个矩形与其他矩形的交。

所以现在相当于我们只需要实现二维的区间加/查询。

但是数据范围很大我们不可能用二维树状数组搞 不然这题就被爆艹了 。

所以我们考虑扫描线 + 一维树状数组来搞。

我们把纵坐标用树状数组维护起来,横坐标用扫描线代替。

对于每个矩形弄四条扫描线,两条修改两条查询。

其中一条修改在左边界加,另一条在右边界右边一个位置减。

其中一个查询在左边界左边一个位置算负的贡献,另一个在右边界算正的贡献。

那么有一个问题,为什么我们要对于左边界左边搞一条扫描线算负贡献呢?

我们考虑差分时的效果,我们对于一条加的修改扫描线,我们是把这条扫描线以左的一整块矩形区域都加了 \(1\)

那么我们如果在加之前,把一段区间的贡献减掉,实际上是减掉了若干个两个不同矩形左边界之间的值,那么我们这时如果在右边碰到一条查询正贡献的扫描线,我们会发现这两条扫描线的贡献一结合刚好就得到了一块矩形交的贡献,也就是说我们这样算就可以很方便地算矩形之间的交了。

具体细节可以手画一下图,那么我们接下来就用二维树状数组的实现技巧来搞。

还有一些小问题:运算过程中可能会爆 long long 所以我们用 long double 来存;然后就是树状数组不能以零为下标,所以我们再输入时把所有矩形的横纵坐标都往各正方向平移 \(2\) 个单位。

参考代码:

#include <algorithm>
#include <cstdio>
#define rg register
#define file(x) freopen(x".in", "r", stdin), freopen(x".out", "w", stdout)
using namespace std;
template < class T > inline void read(T& s) {
s = 0; int f = 0; char c = getchar();
while ('0' > c || c > '9') f |= c == '-', c = getchar();
while ('0' <= c && c <= '9') s = s * 10 + c - 48, c = getchar();
s = f ? -s : s;
} typedef long double LL;
const int _ = 2e5 + 5, __ = 2e6 + 5; int n, num = 0; LL ans = 0;
struct node { int opt, x, l, r, v; } t[_ << 2];
inline bool cmp(const node& a, const node& b)
{ return a.x != b.x ? a.x < b.x : (a.opt != b.opt ? a.opt < b.opt : a.v < b.v); } struct BIT {
LL tr[__];
inline void update(int x, LL v) { for (rg int i = x; i < __; i += i & -i) tr[i] += v; }
inline LL query(int x) { LL res = 0; for (rg int i = x; i >= 1; i -= i & -i) res += tr[i]; return res; }
} tr1, tr2, tr3, tr4; inline void Update(int x, int y, int v) {
tr1.update(y, v), tr2.update(y, x * v), tr3.update(y, y * v), tr4.update(y, (LL) x * y * v);
} inline LL Query(int x, int y) {
return tr1.query(y) * (x + 1) * (y + 1) - tr2.query(y) * (y + 1) - tr3.query(y) * (x + 1) + tr4.query(y);
} int main() {
#ifndef ONLINE_JUDGE
file("cpp");
#endif
read(n);
for (rg int x, y, a, b, i = 1; i <= n; ++i) {
read(x), read(y), read(a), read(b);
x += 2, y += 2;
t[++num] = (node) { 0, x, y, y + b - 1, 1 };
t[++num] = (node) { 0, x + a, y, y + b - 1, -1 };
t[++num] = (node) { 1, x - 1, y, y + b - 1, -1 };
t[++num] = (node) { 1, x + a - 1, y, y + b - 1, 1 };
ans -= 1ll * a * b;
}
sort(t + 1, t + num + 1, cmp);
for (rg int i = 1; i <= num; ++i) {
if (t[i].opt == 0)
Update(t[i].x, t[i].l, t[i].v), Update(t[i].x, t[i].r + 1, -t[i].v);
else
ans += t[i].v * (Query(t[i].x, t[i].r) - Query(t[i].x, t[i].l - 1));
}
printf("%.9Lf\n", 1.0 * ans / n / (n - 1));
return 0;
}

「JSOI2014」矩形并的更多相关文章

  1. 「JSOI2014」打兔子

    「JSOI2014」打兔子 传送门 首先要特判 \(k \ge \lceil \frac{n}{2} \rceil\) 的情况,因为此时显然可以消灭所有的兔子,也就是再环上隔一个点打一枪. 但是我们又 ...

  2. 「JSOI2014」电信网络

    「JSOI2014」电信网络 传送门 一个点选了就必须选若干个点,最大化点权之和,显然最大权闭合子图问题. 一个点向它范围内所有点连边,直接跑最大权闭合子图即可. 参考代码: #include < ...

  3. 「JSOI2014」学生选课

    「JSOI2014」学生选课 传送门 看到这题首先可以二分. 考虑对于当前的 \(mid\) 如何 \(\text{check}\) 我们用 \(f_{i,j}\) 来表示 \(i\) 对 \(j\) ...

  4. 「JSOI2014」歌剧表演

    「JSOI2014」歌剧表演 传送门 没想到吧我半夜切的 这道题应该算是 \(\text{JSOI2014}\) 里面比较简单的吧... 考虑用集合关系来表示分辨关系,具体地说就是我们把所有演员分成若 ...

  5. 「JSOI2014」支线剧情2

    「JSOI2014」支线剧情2 传送门 不难发现原图是一个以 \(1\) 为根的有根树,所以我们考虑树形 \(\text{DP}\). 设 \(f_i\) 表示暴力地走完以 \(i\) 为根的子树的最 ...

  6. 「JSOI2014」强连通图

    「JSOI2014」强连通图 传送门 第一问很显然就是最大的强连通分量的大小. 对于第二问,我们先把原图进行缩点,得到 \(\text{DAG}\) 后,统计出入度为零的点的个数和出度为零的点的个数, ...

  7. 「JSOI2014」序列维护

    「JSOI2014」序列维护 传送门 其实这题就是luogu的模板线段树2,之所以要发题解就是因为学到了一种比较NB的 \(\text{update}\) 的方式.(参见这题) 我们可以把修改操作统一 ...

  8. 「AHOI2014/JSOI2014」拼图

    「AHOI2014/JSOI2014」拼图 传送门 看到 \(n \times m \le 10^5\) ,考虑根号分治. 对于 \(n < m\) 的情况,我们可以枚举最终矩形的上下边界 \( ...

  9. 「AHOI2014/JSOI2014」宅男计划

    「AHOI2014/JSOI2014」宅男计划 传送门 我们首先要发现一个性质:存货天数随买食物的次数的变化类似于单峰函数. 具体证明不会啊,好像是二分加三分来证明?但是没有找到明确的严格证明. 感性 ...

随机推荐

  1. JS-对象常用方法整理

    查看对象的方法,继续控制台输出,如图: hasOwnProperty():返回一个布尔值,指示对象自身属性中是否具有指定的属性(也就是,是否有指定的键). let object1 = new Obje ...

  2. [Vue源码]一起来学Vue双向绑定原理-数据劫持和发布订阅

    有一段时间没有更新技术博文了,因为这段时间埋下头来看Vue源码了.本文我们一起通过学习双向绑定原理来分析Vue源码.预计接下来会围绕Vue源码来整理一些文章,如下. 一起来学Vue双向绑定原理-数据劫 ...

  3. IntelliJ IDEA 2017.3尚硅谷-----界面展示

  4. pycharm的一些快捷键[转]

    编辑类: Ctrl + Space 基本的代码完成(类.方法.属性)Ctrl + Alt + Space 类名完成Ctrl + Shift + Enter 语句完成Ctrl + P 参数信息(在方法中 ...

  5. JSON--WEB SERVICE

    Query ajax webservice:get 和 post 一.GET 方式 客户端 复制代码代码如下: var data = { classCode: "0001"}; / ...

  6. Sql 中获取年月日时分秒的函数

    getdate():获取系统当前时间 dateadd(datepart,number,date):计算在一个时间的基础上增加一个时间后的新时间值,比如:dateadd(yy,30,getdate()) ...

  7. Bugku-CTF之多次

    Day33   多次 http://123.206.87.240:9004 本题有2个flag flag均为小写 flag格式 flag{}  

  8. Codeforces Round #598 (Div. 3) C. Platforms Jumping

    There is a river of width nn. The left bank of the river is cell 00 and the right bank is cell n+1n+ ...

  9. ACM-ICPC实验室20.2.21测试-图论(二)

    A.患者的编号 给出一个有向图,要求你输出字典序最小的拓扑排序. 常规拓扑排序是做不了的,正解是反向建图,同时用大根堆的优先队列维护,保证每次优先访问编号大的结点,再反向输出~ #include< ...

  10. 003 CSS汇总

    字体属性:(font) 大小 {font-size: x-large;}(特大) xx-small;(极小) 一般中文用不到,只要用数值就可以,单位:PX.PD 样式 {font-style: obl ...