这三个题写了一天半,第一个题写了大概一整天。出题人劝我从后往前写,我不听,结果T1想了+调了一天QWQWQ

Task 1:序列

【问题描述】

 定义一个“好的序列”为一个长度为M的不下降序列,且序列中的元素均为1-N的正整数。现在我们随机生成一个“好的序列”,每个不同的“好的序列”出现的概率相同,求这个序列中众数的出现次数的期望。

【输入】

 输入文件含有多组数据

 每组数据读入一行两个整数M、N。

【输出】

 对于每组测试数据,输出一行一个实数表示答案,精确到小数点后4位。

【输入输出样例】

 1 5

 2 9

 3 3

 1.0000

 1.2000

 2.2000

【样例解释】

 令x表示序列中众数的出现次数。

 当M = 3,N = 3时可能的序列有:

  • 1 1 1(x=3) 1 1 2(x=2) 1 1 3(x=2)

  • 1 2 2(x=2) 1 2 3(x=1) 1 3 3(x=2)

  • 2 2 2(x=3) 2 2 3(x=2) 2 3 3(x=2)

  • 3 3 3(x=3)

 因此x的期望为2.2。

【数据范围】

 1 ≤ M≤ 250,1 ≤ N ≤ 10^9

 每个点的总数据组数不会超过15组。

【提示】

  • 题目难度与题目顺序无关。

  • 请注意精度问题。

 关键难点在于如何构建状态。

 这里我们设sum[ k ][ i ][ j ]为众数个数为[ 1 , k ],数列长度为i,选用元素为j的情况数量之和。

 这样设定的好处,看下面很快就会明白:

  • 对于数列本身,我们有这样的划分:

    • 对于某个特定长度的数列,我们考虑把它离散化处理。

      成多块,每块由相同元素组成,一共ki[ 1 , m ]块。

      • 在n个元素中选出k个填入这些块中,离散化的结果和数列本身效果一致,根据乘法原理,直接把可以选的个数C(n,k)与最终答案对应相乘即可。

        • 为了处理答案,我们维护一个前缀和:众数个数<=k的情况下序列总个数的和
        • 为什么维护前缀和而不是单独求本身某一个?因为求和制约条件更宽松,容易推导和判断。

更新时,累加前面两种情况:

  • 上一个选择的元素是最新的元素

  • 上一个选择的元素不是最新元素

需要判重的情况:

  • 众数个数超过k个的话。

  • 如果同一个元素在这里累积了k+1次,就应该去掉。

 所以转移方程就可以得到了:

sum[ k ][ i ][ j ]=sum[ k ][ i - 1 ][ j - 1 ]+sum[ k ][ i - 1 ][ j ];
if( k < i )sum[ k ][ i ][ j ]-=sum[ k ][ i - k - 1 ][ j - 1 ];

 为了保证答案精度,这里使用long double。

 组合数的计算:这里我选择线性递推法.

 由于数据有15组之多,我们选择先进行不同长度情况的预处理,对于不同的可选元素个数只需要搞个组合数上去就可以了。

 复杂度:O(n^3)

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std; int T,m,n;
long double C[255],sum[255][255][255],f[255][255][255],g[255]; int main(){
freopen("sequence.in","r",stdin);
freopen("sequence.out","w",stdout);
C[0]=true;
for(register int k=1;k<=251;++k){
sum[k][0][0]=1;
for(register int i=1;i<=251;++i){
for(register int j=1;j<=i;++j){
//i->序列长度 j->总共使用的元素个数 k->众数的个数
//sum存储的就是在当前的众数个数的情况下,可以构造出来的序列个数。
//由于是一步一步递推得到,可以直接适用于所有m∈[1,251]的情况。
sum[k][i][j]=sum[k][i-1][j-1]+sum[k][i-1][j];
//这里维护的是一个前缀和:求的是众数个数为k情况下序列总个数的和
//更新时,累加前面两种情况:
//->上一个选择的元素是最新的元素
//->上一个选择的元素不是最新元素
if(k<i){
sum[k][i][j]-=sum[k][i-k-1][j-1];
//需要判重的情况:众数个数超过k个的话。
//等效于->如果同一个元素在这里累积了k+1次,就应该去掉。
}
}
}
}
while(cin>>m>>n){
memset(g,0,sizeof(g));
for(register int i=1;i<=m;++i){
C[i]=C[i-1]*(n-i+1)/i;
}//O(m)递推求组合数
for(register int k=1;k<=m;++k){
for(register int j=1;j<=m;++j){
g[k]+=sum[k][m][j]*C[j];
}//计数求和。
}
long double ans=0;
//g表示的是总共可能出现的情况种类数,是累加
for(register int k=1;k<=m;++k){
ans+=k*(g[k]-g[k-1])/g[m];
} //给出答案
printf("%.4Lf\n",(double)ans);
}
return 0;
}

 T1相当有难度,但T2T3都是大水题。


Task 2:游戏

【问题描述】

 小H喜欢玩游戏。她会首先选择一个数

【清北学堂2018-刷题冲刺】Contest 5的更多相关文章

  1. 2017 清北济南考前刷题Day 7 afternoon

    期望得分:100+100+30=230 实际得分:100+100+30=230 1. 三向城 题目描述 三向城是一个巨大的城市,之所以叫这个名字,是因为城市中遍布着数不尽的三岔路口.(来自取名力为0的 ...

  2. 2017 清北济南考前刷题Day 1 afternoon

    期望得分:80+30+70=180 实际得分:10+30+70=110 T1 水题(water) Time Limit:1000ms   Memory Limit:128MB 题目描述 LYK出了道水 ...

  3. 2017 清北济南考前刷题Day 3 morning

    实际得分:100+0+0=100 T1 右上角是必败态,然后推下去 发现同行全是必胜态或全是必败态,不同行必胜必败交叉 列同行 所以n,m 只要有一个是偶数,先手必胜 #include<cstd ...

  4. 2017 清北济南考前刷题Day 3 afternoon

    期望得分:100+40+100=240 实际得分:100+40+100=240 将每个联通块的贡献乘起来就是答案 如果一个联通块的边数>点数 ,那么无解 如果边数=点数,那么贡献是 2 如果边数 ...

  5. 2017 清北济南考前刷题Day 4 afternoon

    期望得分:30+50+30=110 实际得分:40+0+0=40 并查集合并再次写炸... 模拟更相减损术的过程 更相减损术,差一定比被减数小,当被减数=减数时,停止 对于同一个减数来说,会被减 第1 ...

  6. 2017 清北济南考前刷题Day 7 morning

    期望得分:100+50+20=170 实际得分:10+50+20=80 1. 纸牌 题目描述 在桌面上放着n张纸牌,每张纸牌有两面,每面都写着一个非负整数.你的邪王真眼可以看到所有牌朝上的一面和朝下的 ...

  7. 2017 清北济南考前刷题Day 6 afternoon

    期望得分:100+100+30=230 实际得分: 正解: 枚举最高的位,这一位m是1但实际用了0 然后剩余的低位肯定是 正数就用1,负数用0 考场思路:数位DP #include<cstdio ...

  8. 2017 清北济南考前刷题Day 6 morning

    T1 贪心 10 元先找5元 20元 先找10+5,再找3张5 #include<cstdio> using namespace std; int m5,m10,m20; int main ...

  9. 2017 清北济南考前刷题Day 5 afternoon

    期望得分:100+100+30=230 实际得分:0+0+0=30 T1 直接模拟 #include<cstdio> #include<iostream> using name ...

  10. 2017 清北济南考前刷题Day 5 morning

    期望得分:100+100+0=200 实际得分: 坐标的每一位不是0就是1,所以答案就是 C(n,k) #include<cstdio> #include<iostream> ...

随机推荐

  1. 安装 BizTalk Server 2016

    在单台计算机上安装 BizTalk Server. 开始操作之前       系统管理员 – 安装 SQL Server 时,安装程序会自动向登录的帐户授予系统管理员权限. 由于安装 BizTalk ...

  2. Nginx 慢启动与拥塞窗口

    L:127

  3. ubuntu系统安装mysql(deb-bundle包)

    由于某些原因,又要在ubuntu系统中安装mysql了,之前曾经安装过好多次.都没记下来 以前一直动用源码包来安装,基于两个原因:1.一直用Python写代码.2.想使用文件来安装,而不是通过api ...

  4. nginx 正向代理上网

    配置文件: server { #resolver 21.202.152.10; #指定DNS服务器IP地址 |如果指定IP$scheme://22.2.65.214$request_uri 可以不指定 ...

  5. vue——script内容详解

    <script> export default { name: "Home", data() { return {}; }, methods: { // 组件的方法 } ...

  6. springboot 学习

    参考:http://www.cnblogs.com/sam-uncle/p/8796212.html spring boot 系列之一:spring boot 入门 注意:main启动类和contro ...

  7. 图灵机器人API接口

    调用图灵API接口实现人机交互 流程一: 注册 图灵机器人官网: http://www.tuling123.com/ 第一步: 先注册, 然后创建机器人, 拿到一个32位的key 编码方式 UTF-8 ...

  8. 数据库SQL SELECT查询的工作原理

    一般开发员只会应用SQL的四条经典语句:select,insert,delete,update.但是我从来没有研究过它们的工作原理,这篇我想说一说select在数据库中的工作原理. B/S架构中最经典 ...

  9. 洛谷3822 [NOI2017] 整数 【线段树】【位运算】

    题目分析: 首先这题的询问和位(bit)有关,不难想到是用线段树维护位运算. 现在我们压32位再来看这道题. 对于一个加法操作,它的添加位置可以得到,剩下的就是做不超过32的位移.这样根据压位的理论. ...

  10. Django+Vue打造购物网站(一)

    环境搭建 python == 3.6 Django == 2.0 创建工程 django-admin startproject MxShop 配置setting.py文件 # 数据库 DATABASE ...