[计蒜客T2238]礼物_线段树_归并排序_概率期望
礼物
题目大意:
数据范围:
题解:
这题有意思啊($md$卡常
直接做怎么做?
随便上个什么东西,维护一下矩阵乘和插入,比如说常数还算小的$KD-Tree$(反正我是没见人过过
我们漏掉了一个条件,就是所有二元组都是随机的。
这个条件很好,它几乎就保证了,任选一个区间的话,优秀二元组只有$log$个。
这是为什么呢?
其实区间内,优秀二元组的个数,就相当于把区间按照$x$排序后,$y$值是前缀最大值的期望个数。
因为二元组是随机的,所以$x$排序后,$y$仍然是随机的。
就是给定一个随机数列,求前缀最大值的期望个数。
这是调和级数的。
所以,我们就开一棵线段树,线段树上每个节点维护一个数组,存这个节点管辖区间内的优秀二元组。
合并用归并,复杂度是$O(log)$的。
所以每次查询的复杂度是$O(log^2n)$的。
总复杂度是$O(nlog^2n)$的,有点小卡常,加了输出优化才过(读入优化是必备。
代码:
#include <bits/stdc++.h> #define ls p << 1 #define rs p << 1 | 1 #define N 200010 using namespace std; const int mod = 1000000007 ; typedef long long ll; char *p1, *p2, buf[100000]; #define nc() (p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 100000, stdin), p1 == p2) ? EOF : *p1 ++ ) int rd() {
int x = 0;
char c = nc();
while (c < 48) {
c = nc();
}
while (c > 47) {
x = (((x << 2) + x) << 1) + (c ^ 48), c = nc();
}
return x;
} char pbuf[100000],*pp=pbuf;
void push(const char c) {
if(pp-pbuf==100000) fwrite(pbuf,1,100000,stdout),pp=pbuf;
*pp++=c;
}
void write(int x) {
static int sta[35];
int top=0;
do{sta[top++]=x%10,x/=10;}while(x);
while(top) push(sta[--top]+'0');
push('\n');
} struct Node {
int x, y;
friend bool operator < (const Node &a, const Node &b) {
return a.x == b.x ? a.y < b.y : a.x < b.x;
}
}num[N]; Node q[60]; struct Mode {
Node v[30];
int len;
int sum;
Mode() { len = 0, sum = 1; }
}a[N << 2]; inline Mode operator + (const Mode &a, const Mode &b) {
Mode re;
int cnt = 0;
int j = 1;
for (int i = 1; i <= a.len; i ++ ) {
while ((j <= b.len) && (a.v[i] < b.v[j])) {
q[ ++ cnt] = b.v[j];
j ++ ;
}
q[ ++ cnt] = a.v[i];
}
while (j <= b.len) {
q[ ++ cnt] = b.v[j];
j ++ ;
}
// int i = 1, j = 1;
// while (i <= a.len && j <= b.len) {
// if (a.v[i] < b.v[j]) {
// q[ ++ cnt] = a.v[i];
// i ++ ;
// }
// else {
// q[ ++ cnt] = b.v[j];
// j ++ ;
// }
// }
// while (i <= a.len) {
// q[ ++ cnt] = a.v[i];
// i ++ ;
// }
// while (j <= b.len) {
// q[ ++ cnt] = b.v[j];
// j ++ ;
// }
// for (int i = 1; i <= a.len; i ++ ) {
// q[ ++ cnt] = a.v[i];
// }
// for (int i = 1; i <= b.len; i ++ ) {
// q[ ++ cnt] = b.v[i];
// }
// sort(q + 1, q + cnt + 1);
int mx = 0;
for (int i = 1; i <= cnt; i ++ ) {
if (q[i].y > mx) {
re.v[ ++ re.len] = q[i];
re.sum = (ll)re.sum * (q[i].x ^ q[i].y % mod) % mod;
mx = q[i].y;
}
}
// reverse(re.v + 1, re.v + re.len + 1);
return re;
} void build(int l, int r, int p) {
if (l == r) {
a[p].v[ ++ a[p].len] = num[l];
a[p].sum = (num[l].x ^ num[l].y) % mod;
return;
}
int mid = (l + r) >> 1;
build(l, mid, ls);
build(mid + 1, r, rs);
a[p] = a[ls] + a[rs];
} Mode query(int x, int y, int l, int r, int p) {
if (x <= l && r <= y) {
return a[p];
}
int mid = (l + r) >> 1;
if (mid < x)
return query(x, y, mid + 1, r, rs);
else if(y <= mid)
return query(x, y, l, mid, ls);
else
return query(x, y, l, mid, ls) + query(x, y, mid + 1, r, rs);
} int main() {
int n = rd();
for (int i = 1; i <= n; i ++ ) {
num[i].x = rd();
num[i].y = rd();
}
build(1, n, 1);
int q = rd();
while(q -- ) {
int x = rd(), y = rd();
Mode now = query(x, y, 1, n, 1);
write(now.sum);
}
fwrite(pbuf,1,pp-pbuf,stdout);
return 0;
}
小结:这种期望的题还是要自己证才行,不然结论根本记不过来。
[计蒜客T2238]礼物_线段树_归并排序_概率期望的更多相关文章
- 【原创】tyvj1038 忠诚 & 计蒜客 管家的忠诚 & 线段树(单点更新,区间查询)
最简单的线段树之一,中文题目,不翻译.... 注释讲的比较少,这已经是最简单的线段树,如果看不懂真的说明最基础的理论没明白 推荐一篇文章http://www.cnblogs.com/liwenchi/ ...
- [计蒜客] 矿石采集【记搜、Tarjan缩点+期望Dp】
Online Judge:计蒜客信息学3月提高组模拟赛 Label:记搜,TarJan缩点,树状数组,期望Dp 题解 整个题目由毫无关联的两个问题组合成: part1 问题:对于每个询问的起点终点,求 ...
- 计蒜客 28315.Excellent Engineers-线段树(单点更新、区间最值) (Benelux Algorithm Programming Contest 2014 Final ACM-ICPC Asia Training League 暑假第一阶段第二场 E)
先写这几道题,比赛的时候有事就只签了个到. 题目传送门 E. Excellent Engineers 传送门 这个题的意思就是如果一个人的r1,r2,r3中的某一个比已存在的人中的小,就把这个人添加到 ...
- 计蒜客 38228. Max answer-线段树维护单调栈(The Preliminary Contest for ICPC China Nanchang National Invitational I. Max answer 南昌邀请赛网络赛) 2019ICPC南昌邀请赛网络赛
Max answer Alice has a magic array. She suggests that the value of a interval is equal to the sum of ...
- 计蒜客 41391.query-二维偏序+树状数组(预处理出来满足情况的gcd) (The Preliminary Contest for ICPC Asia Xuzhou 2019 I.) 2019年徐州网络赛)
query Given a permutation pp of length nn, you are asked to answer mm queries, each query can be rep ...
- 计蒜客 Prefix Free Code(字典树+树状数组)
Consider n initial strings of lower case letters, where no initial string is a prefix of any other i ...
- 计蒜客 NOIP 提高组模拟竞赛第一试 补记
计蒜客 NOIP 提高组模拟竞赛第一试 补记 A. 广场车神 题目大意: 一个\(n\times m(n,m\le2000)\)的网格,初始时位于左下角的\((1,1)\)处,终点在右上角的\((n, ...
- 计蒜客 28449.算个欧拉函数给大家助助兴-大数的因子个数 (HDU5649.DZY Loves Sorting) ( ACM训练联盟周赛 G)
ACM训练联盟周赛 这一场有几个数据结构的题,但是自己太菜,不会树套树,带插入的区间第K小-替罪羊套函数式线段树, 先立个flag,BZOJ3065: 带插入区间K小值 计蒜客 Zeratul与Xor ...
- 计蒜客 A1607 UVALive 8512 [ACM-ICPC 2017 Asia Xi'an]XOR
ICPC官网题面假的,要下载PDF,点了提交还找不到结果在哪看(我没找到),用VJ交还直接return 0;也能AC 计蒜客题面 这个好 Time limit 3000 ms OS Linux 题目来 ...
随机推荐
- Maven手动将jar导入本地仓库
1.使用cmd进入maven安装目录下的bin 2.运行mvn install:install-file -Dfile=jar包的路径 -DgroupId=gruopId中的内容 -Dartifact ...
- Neo4j (1)创建节点
节点模式的构成:(Variable:Lable1:Lable2{Key1:Value1,Key2,Value2}),实际上,每个节点都有一个整数ID,在创建新的节点时,Neo4j自动为节点设置ID值, ...
- load加载层
//eg1 var index = layer.load(); //eg2 var index = layer.load(1); //换了种风格 //eg3 var index = layer.loa ...
- linux下通过进程名查看其占用端口
linux下通过进程名查看其占用端口: 1.先查看进程pid ps -ef | grep 进程名 2.通过pid查看占用端口 netstat -nap | grep 进程pid 例:通过nginx进程 ...
- 实现多列等高布局_flex布局
详情参见此篇博客 http://www.w3cplus.com/css/creaet-equal-height-columns 建议掌握方法四.五 其实,利用最新的flex布局 http://www. ...
- leetcode题目142.环形链表Ⅱ(中等)
题目描述: 给定一个链表,返回链表开始入环的第一个节点. 如果链表无环,则返回 null. 为了表示给定链表中的环,我们使用整数 pos 来表示链表尾连接到链表中的位置(索引从 0 开始). 如果 p ...
- 同源策略和Ajax跨域访问
1. 什么是同源策略 理解跨域首先必须要了解同源策略.同源策略是浏览器上为安全性考虑实施的非常重要的安全策略. 何谓同源: URL由协议.域名.端口和路径组成,如果两个URL的协议 ...
- 代码bug管理工具bugfree与禅道
禅道官网有一键安装包,一键安装即可用 bugfree 搭建lnmp环境 详情参考zabbix 解压bugfree包 到nginx的html里 unzip bugfree3.0.4.zip mv b ...
- C#txt文件创建并写入信息
public static void LogWrite(string str) { //项目根目录 string path = HttpContext.Current.Server.MapPath(& ...
- 3.创建一个pod应用
创建一个应用:k8s增删查改: pod创建:kubectl run nginx-deploy --image=nginx:1.14-alpine --port=80 --replicas=1 [roo ...