csp-s模拟测试101的T3代码+注释
因为题目过于大神所以单独拿出来说。而且既然下发std了颓代码貌似也不算可耻233
很难讲啊,所以还是写在代码注释里面吧
因为比较认真的写了不少注释,所以建议缩放到80%观看,或者拿到gedit上
1 //不要在意我写的是“气”,输入法找不到那个字。。。
2 #include<cstdio>
3 #include<algorithm>
4 long double C[205][205];int n,g,t,a[205],N;
5 struct Ex{//用于存储期望值
6 long double cnt,tot;//cnt表示到达这个状态的总方案数,tot表示这些方案所获得的神的总和
7 void add(long double e,long double c){tot+=e*c;cnt+=c;}//这里是正推的,所以概率加和的结果不为1,所以期望不能直接相加
8 //期望其实类似于加权平均数。现在你有cnt个方案,和为tot。要加入c个方案,这些方案的平均数是e,所以这c个方案的和为e*c
9 //那么总的方案数就是cnt+c。总的和就是tot+e*c
10 long double E(){return cnt?tot/cnt:0;}//期望值就是总量/方案数啊,可以重载成+=,但是因为参数有2个还得写pair什么的所以就没有重载
11 }w[205][205],f[205][205];//全局变量初值为0,所以刚开始的方案数和期望值都是0
12 long double cal(int l,int r,int ord){return ord>n?0:(l+r+1)/2.0;}//计算平均数,表示用了一份l的气所得到的神的期望。
13 //如果用的是原来不存在的气,就是你加的那些编号为[n+1,n+t]的气,这些气是不存在的,所以其实你放弃了这一份精,所以并不会得到神,返回0
14 int main(){
15 freopen("surmount.in","r",stdin);freopen("surmount.out","w",stdout);
16 scanf("%d%d%d",&n,&g,&t);//方便查阅:n是气的份数,g是单次最大的精量,t是轮数
17 for(int i=1;i<=n;++i)scanf("%d",&a[i]);
18 for(int i=1;i<=n;++i)if(a[i]>g)a[i]=g;//大于g的气是没有用的,直接重置为g
19 std::sort(a+1,a+1+n);N=n+t;//排序方便处理连续的区间
20 for(int i=n+1;i<=N;++i)a[i]=g;//添加一些足够大的气,这样在下面处理的时候就能保证每一份精都能被压制
21 for(int i=0;i<=t;++i)C[i][0]=1;
22 for(int i=1;i<=t;++i)for(int j=1;j<=i;++j)C[i][j]=C[i-1][j-1]+C[i-1][j];//处理组合数
23 //w[i][j]表示对于[i,j)这段区间的所有气,如果使用了它们所获得的神的期望值
24 for(int i=0;i<=N;++i)w[i][i].add(0,1);//赋初值:空区间的方案数是1,对神的贡献为0
25 for(int i=N;~i;--i)for(int j=i+1;j<=N&&j-i<=t;++j)for(int k=i;k<j;++k)
26 w[i][j].add(w[i][k].E()+w[k+1][j].E()+cal(a[k+1],a[i],k+1),w[i][k].cnt*w[k+1][j].cnt*(a[k+1]-a[i])*C[j-i-1][k-i]);
27 //现在的操作是要把[i,k)与[k+1,j)合并成[i,j)。假设最后一个使用的精是k。注意区间开闭。
28 //add的两个参数:期望值就是[i,k)和[k+1,j)这两段的期望值的和,然而还有使用k所得到的神,你使用它是用来压制[a[k+1],a[i]]的精,得到对应的神,详见cal函数
29 //总方案数就是两边的方案相乘(乘法计数原理)。还要乘上组合数,因为要考虑这些气的顺序可以颠倒。因为左右两边内部已经有序了,所以只需要把两个混合起来。
30 //就是类似于归并排序两个内部有序的数组,总方案数就是C[j-i-1][k-i],表示在所有j-i-1个位置里选出左边所用的k-i个元素的位置,注意第k个被你强制定为最后一个了所以是j-i还要-1
31 for(int i=0;i<=t;++i)f[i][i]=w[0][i];
32 //初值,表示前i-1份气没有用。超过t没有意义。
33 for(int i=0;i<=N;++i)for(int k=0;k<=t;++k)if(f[i][k].cnt)for(int j=i+1;j<=N&&k+j-i-1<=t;++j)
34 f[j][k+j-i-1].add(f[i][k].E()+w[i+1][j].E(),f[i][k].cnt*w[i+1][j].cnt*C[k+j-i-1][k]);
35 //i表示转移之前已经用了几份气(含),j表示转移之后用的气,k表示转移了几轮,所以k+j-i-1就是转移之后已经经过了的轮数
36 //f[a][b]表示用了前a份气,已经经过了b轮之后得到的神的期望,所以最后答案是f[N][t]
37 //转移就是,先看期望,就是已经做完的f[i][k]再加上新的一轮你用了[i+1,j)这一段的气,这样所获得的神
38 //总方案数其实也一样,和w的转移比较类似,也是排序了两个内部有序的过程
39 printf("%.10Lf\n",f[N][t].E());
40 }//By DeepinC究级压行+垃圾注释(硬核25行)
累。。。
csp-s模拟测试101的T3代码+注释的更多相关文章
- NOIP模拟测试18(T3待更新)
T1: 直接模拟,详见代码注释. 复杂度$O(NM)$. Code: #include<iostream> #include<cstdio> #include<vecto ...
- [考试反思]1105csp-s模拟测试101: 临别
先不改题,这次主要不在T3上. 这次有必要粘文件得分了. 临考前总解锁新锅我也不知道这是什么个事啊... T1宏定义写挂.因为原来在OJ上没事所以一直没注意.在Lemon评测下直接全部RE. GG在主 ...
- 利用Python中的mock库对Python代码进行模拟测试
这篇文章主要介绍了利用Python中的mock库对Python代码进行模拟测试,mock库自从Python3.3依赖成为了Python的内置库,本文也等于介绍了该库的用法,需要的朋友可以参考下 ...
- 【转】利用Python中的mock库对Python代码进行模拟测试
出处 https://www.toptal.com/python/an-introduction-to-mocking-in-python http://www.oschina.net/transla ...
- noi2019模拟测试赛(四十七)
noi2019模拟测试赛(四十七) T1与运算(and) 题意: 给你一个序列\(a_i\),定义\(f_i=a_1\&a_2\&\cdots\&a_i\),求这个序列的所 ...
- [考试反思]0729NOIP模拟测试10
安度因:哇哦. 安度因:谢谢你. 第三个rank1不知为什么就来了.迷之二连?也不知道哪里来的rp 连续两次考试数学都占了比较大的比重,所以我非常幸运的得以发挥我的优势(也许是优势吧,反正数学里基本没 ...
- csp-s模拟测试95
csp-s模拟测试95 去世场祭. $T1$:这不裸的除法分块吗. $T2$:这不裸的数据结构优化$Dp$吗. $T3$:这不裸的我什么都不会搜索骗$30$分吗. 几分钟后. 这除法分块太劲了..(你 ...
- 0823NOIP模拟测试赛后总结
考了两场感觉虚了... NOIP模拟测试30 分着考的. 就只有T2的美妙的暴力拿分了,60分rank10,挂了. T1是一道sb题,爆零了十分遗憾. 许多人都掉进了输出格式的坑里,C没大写.少个空格 ...
- 2019.8.14 NOIP模拟测试21 反思总结
模拟测试20的还没改完先咕着 各种细节问题=错失190pts T1大约三分钟搞出了式子,迅速码完,T2写了一半的时候怕最后被卡评测滚去交了,然后右端点没有初始化为n…但是这样还有80pts,而我后来还 ...
随机推荐
- 在 Cocos Creator 中使用 Protobufjs(一)
一. 环境准备 我一直在探索Cocos H5正确的开发姿势,目前做javascript项目已经离不开 nodejs.npm或grunt等脚手架工具了. 1.初始化package.json文件 npm ...
- 从零开始入门 K8s | 可观测性:你的应用健康吗?
作者 | 莫源 阿里巴巴技术专家 一.需求来源 首先来看一下,整个需求的来源:当把应用迁移到 Kubernetes 之后,要如何去保障应用的健康与稳定呢?其实很简单,可以从两个方面来进行增强: 首先是 ...
- Redis 的底层数据结构(跳跃表)
字典相对于数组,链表来说,是一种较高层次的数据结构,像我们的汉语字典一样,可以通过拼音或偏旁唯一确定一个汉字,在程序里我们管每一个映射关系叫做一个键值对,很多个键值对放在一起就构成了我们的字典结构. ...
- 10个值得深思的_PHP_面试问题
Q1 第一个问题关于弱类型 $str1 = 'yabadabadoo'; $str2 = 'yaba'; if (strpos($str1,$str2)) { echo "\"&q ...
- docker在Mac上的下载安装
在Mac上下载安装docker,下载链接:Stable 安装成功后启动终端,检查安装后的docker版本: yanguobindeMacBook-Pro:~ yanguobin$ docker --v ...
- 如何正确遍历删除List中的元素(普通for循环、增强for循环、迭代器iterator、removeIf+方法引用)
遍历删除List中符合条件的元素主要有以下几种方法: 普通for循环 增强for循环 foreach 迭代器iterator removeIf 和 方法引用 其中使用普通for循环容易造成遗漏元素的问 ...
- 【TencentOS tiny】 超详细的TencentOS tiny移植到STM32F103全教程
移植前的准备工作 1. 获取STM32的裸机工程模板 STM32的裸机工程模板直接使用野火STM32开发板配套的固件库例程即可.可以从我github上获取https://github.com/jiej ...
- 超详细!! sql server 同步数据库 发布 订阅 跨网段 无公网ip 常见问题
问题描述 主机1:发布端 阿里云服务器--有公网ip 主机2:订阅端 笔记本--无公网ip 数据量很小,主要是熟悉发布订阅的操作流程. 主机2仅仅作为主机1的本地备份,要求修改云服务器上数据后,能通过 ...
- 将自定义功能添加到Spring Data Repository
Spring Data非常方便,可以加快开发速度,避免使用样板代码. 但是,在某些情况下,注释查询不足,而无法达到您可能希望实现的自定义功能. 因此,Spring Data允许我们向Spring Da ...
- PMBOK(第六版) PMP笔记——《十一》第十一章(项目风险管理)
PMBOK(第六版) PMP笔记——<十一>第十一章(项目风险管理) 第十一章 风险管理: 项目的独特性导致项目充满风险,项目风险是一种不确定的事件或条件,可能发生.将 要发生,也可能不发 ...