Tsinsen-A1490 osu! 【数学期望】
问题描述
一张osu!的地图是由若干个“音”排列而成的。在本题中,对于每个音我们只需要考虑成功点击和错过(miss)这两种情况。对于一张osu!地图,玩家的完成情况可以用一个01串表示(0代表miss,1代表成功)。在本题中,使用如下计分规则:将玩家完成一张地图的01串中所有的0删去,则这个串可能会断裂成若干段连续的1。对于一段长度为L的1(L≥1),你的总分会增加L^2+L+1。例如:一张地图有10个音,某玩家完成情况为1011101110,则删除所有0后得到的是“1”“111”和“111”。因此这个玩家的得分为(1+1+1)+(9+3+1)+(9+3+1)=29。
ACMonster要给Sandytea出一张osu!的地图。在一张图中,不同音对于Sandytea而言难度可能是不同的。我们定义一个音的难度系数为Sandytea成功完成这个音的概率,因此这个难度系数是介于0和1之间的。
现在ACMonster写下了一个包含N个音的序列,但是他不想直接把这个序列做成地图,而是选择其中的某个片段。设S(X,Y)代表序列上第X个音到第Y个音构成的序列,ACMonster想知道如果把S(X,Y)对应的序列做成地图,Sandytea的期望得分是多少。有时ACMonster会觉得某个音的难度系数不太合理,因此要进行修改。请你想办法处理ACMonster的修改,并回答他提出的问题。
第二行至第(N+1)行每行包含一个实数,第i行的实数表示第(i-1)个音的难度系数。
第(N+2)行至第(N+M+1)行每行包含三个数Type,X和Y。
如果Type=0,代表ACMonster询问S(X,Y)的期望得分,保证X和Y为整数,1≤X≤Y≤N。
如果Type=1,代表ACMonster要把第X个音的难度系数修改为Y,保证X为整数,1≤X≤N,0≤Y≤1。
0.5
0.5
0 1 2
1 1 0
0 1 2
1.50
对于第二次询问,00,01这两种情况出现的概率均为1/2,得分分别为0,3。因此期望得分为(0+3)/2=1.50。
60%的数据满足N,M≤50000;
100%的数据满足N,M≤500000。
![](/assets/blank.gif)
![](/assets/blank.gif)
1 #include <bits/stdc++.h> 2 #define rep(i, a, b) for (int i = a; i <= b; i++) 3 #define drep(i, a, b) for (int i = a; i >= b; i--) 4 #define REP(i, a, b) for (int i = a; i < b; i++) 5 #define mp make_pair 6 #define pb push_back 7 #define clr(x) memset(x, 0, sizeof(x)) 8 #define xx first 9 #define yy second 10 using namespace std; 11 typedef long long i64; 12 typedef pair<int, int> pii; 13 const int inf = ~0U >> 1; 14 const i64 INF = ~0ULL >> 1; 15 //******************************* 16 17 const int maxn = 500005; 18 19 struct node { 20 double p; 21 double a21, a22, a31, a32; 22 double mul1; 23 node() {} 24 node(double _p, double _a21, double _a22, double _a31, double _a32, double _mul1) : 25 p(_p), a21(_a21), a22(_a22), a31(_a31), a32(_a32), mul1(_mul1) {} 26 27 }; 28 29 double ans0, ans1, ans3[3]; 30 31 /* the matrix 32 {sum[i - 1], E[i - 1], 1} * {1 , 0 , 0} = {sum[i], E[i], 1} 33 {2pi, pi, 0} 34 {pi , pi, 1} 35 */ 36 37 double pi[maxn]; 38 39 struct Seg_Tree { 40 node T[maxn << 2]; 41 42 void Push_up(int o) { 43 T[o].p = T[o << 1].p + T[o << 1 | 1].p, T[o].mul1 = T[o << 1].mul1 + T[o << 1 | 1].mul1; 44 T[o].a21 = T[o << 1].a21 + T[o << 1].a22 * T[o << 1 | 1].a21; 45 T[o].a22 = T[o << 1].a22 * T[o << 1| 1].a22; 46 T[o].a31 = T[o << 1].a31 + T[o << 1].a32 * T[o << 1 | 1].a21 + T[o << 1 | 1].a31; 47 T[o].a32 = T[o << 1].a32 * T[o << 1 | 1].a22 + T[o << 1 | 1].a32; 48 } 49 50 void update(int o, int l, int r, int x, double p) { 51 if (l == r) { 52 T[o] = node(p, 2 * p, p, p, p, (1 - pi[l - 1]) * p); 53 return; 54 } 55 int mid = l + r >> 1; 56 if (x > mid) update(o << 1 | 1, mid + 1, r, x, p); 57 else update(o << 1, l, mid, x, p); 58 Push_up(o); 59 } 60 61 void query(int o, int l, int r, int ql, int qr) { 62 if (ql <= l && r <= qr) { 63 ans1 += T[o].p; 64 ans0 += T[o].mul1; 65 ans3[1] = ans3[1] + ans3[2] * T[o].a21 + T[o].a31; 66 ans3[2] = ans3[2] * T[o].a22 + T[o].a32; 67 return; 68 } 69 int mid = l + r >> 1; 70 if (ql <= mid) query(o << 1, l, mid, ql, qr); 71 if (qr > mid) query(o << 1 | 1, mid + 1, r, ql, qr); 72 } 73 } T; 74 75 int main() { 76 int n, m; 77 scanf("%d%d", &n, &m); 78 rep(i, 1, n) scanf("%lf", &pi[i]), T.update(1, 1, n, i, pi[i]); 79 while (m--) { 80 int op; 81 scanf("%d", &op); 82 if (!op) { 83 int x, y; scanf("%d%d", &x, &y); 84 ans1 = 0, ans0 = 0, ans3[1] = ans3[2] = 0; 85 T.query(1, 1, n, x, y); 86 printf("%.2lf\n", ans1 + ans0 + pi[x - 1] * pi[x] + ans3[1]); 87 } 88 else { 89 double w; int x; 90 scanf("%d%lf", &x, &w); 91 pi[x] = w, T.update(1, 1, n, x, w); 92 if (x != n) T.update(1, 1, n, x + 1, pi[x + 1]); 93 } 94 } 95 return 0; 96 }
View Code
转载于:https://www.cnblogs.com/y7070/p/5042928.html
Tsinsen-A1490 osu! 【数学期望】相关推荐
- 数论小白都能看懂的数学期望讲解
-1.灌水 这里阅读应该效果更佳 想了解更多关于数论的内容,可戳这里 感谢@command_block 大佬提出宝贵建议 也感谢洛谷及UVA的相关题目 如果有小瑕疵可以在评论区提出 内容可能有点多但很 ...
- 随机变量的数字特征(数学期望,方差,协方差与相关系数)
戳这里:概率论思维导图 !!! 数学期望 离散型随机变量的数学期望 (这里要求级数绝对收敛,若不绝对收敛,则E(X)不存在) 如果有绝对收敛,则有 ,其中 连续型随机变量的数学期望 (这里要求绝对收敛 ...
- 解题报告(一)F、(2018 ACM - ICPC shenyang I)Distance Between Sweethearts(数学期望 + 乘法原理 + FWT)(4.5)
繁凡出品的全新系列:解题报告系列 -- 超高质量算法题单,配套我写的超高质量题解和代码,题目难度不一定按照题号排序,我会在每道题后面加上题目难度指数(1∼51 \sim 51∼5),以模板题难度 11 ...
- R语言:求二维变量数学期望
想做一个二维变量数学期望实验, 查看若干资料终于找到方法 先看这篇文章熟悉一下R的函数 http://www.cyclismo.org/tutorial/R/tables.html 构造数据 通过下面 ...
- hdu 5419(数学期望)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5419 题解: 首先分母是C(m,3),考虑如何求出分子 考虑数学期望的独立性,我们首先可以用线性的时间 ...
- 概率论-2.2 随机变量的数学期望(重点:随机变量X的期望)
分布有关的特征数:均值,方差,分位数等 期望的定义: 设离散随机变量X的分布列为pi=p(xi)=P(X=xi),i=1,2,-,n 若Sum(| xi |*p(xi))收敛(等价于Sum( xi * ...
- jzoj3801-[NOIP2014模拟8.23]骰子【数学期望】
正题 题目链接:https://jzoj.net/senior/#main/show/3801 题目大意 mmm面的骰子是1∼m1\sim m1∼m,然后丢nnn次,求最大值的数学期望. 解题思路 若 ...
- HDU - 4586 数学期望
题意 有一个骰子有n个面,掷到每一个面的概率是相等的,每一个面上都有相应的钱数.其中当你掷到某些面 一共m个面之一时,你有多掷一次的机会.问最后所得钱数的期望. 分析 数学期望是什么 数学期望就是一种 ...
- 【BZOJ - 3450】Tyvj1952 Easy(数学期望,期望的线性性)
题干: 某一天WJMZBMR在打osu~~~但是他太弱逼了,有些地方完全靠运气:( 我们来简化一下这个游戏的规则 有n次点击要做,成功了就是o,失败了就是x,分数是按comb计算的,连续a个comb就 ...
- 几何分布的期望和方差公式推导_超几何分布的数学期望与方差推导
考虑 个外表相同的物品,其中有 个同类物品与另一类的 个物品:抽取 个物品,每个物品的抽取等概率随机. 上述便是一个超几何分布(Hypergeometric Distribution)的基本模型. 抽 ...
最新文章
- 多核时代 .NET Framework 4 中的并行编程6---并行LINQ
- 参赛队员对于第十六届全国大学智能车竞赛的建议-05月
- Android中的BuildConfig类怎么来的
- c语言数据结构-算法篇之冒泡排序
- chrome java上传文件_springboot+webuploader 实现大文件切片上传,兼容IE8+,chrome等浏览器,可运行...
- [SpringBoot2]自动配置
- 阿里达摩院python_阿里达摩院出品的735集的python教程
- Kafka集群全部断开,然后重启时报“The broker is trying to join the wrong cluster. Configured zookeeper.connect……”问题
- yum 崩溃的解决方法
- 什么是订单管理系统OMS
- 关于补码--看锆石fpga计算峰峰值时的疑问启示
- 加百列丨浅谈大水面水产养殖的几个误区及对策
- 高低频磨皮(商业磨皮一种)
- basICColor catch 5 Mac(光谱分析颜色测量软件) v5.0.7破解版
- python加减乘除运算代码_四则运算python版
- 【Paper】Learning Phrase Representations using RNN Encoder-Decoder for Statistical Machine Translation
- 数字时代,企业如何选择适合自己的转型平台?
- 软件测试--软件测试执行管理
- 【Java核心技术大会 PPT分享】陈阳:深入理解 Java 虚拟机编译原理
- 河南理工大学计算机学院保研,#研历言谈#优秀保研生介绍之测绘篇