【题解】CF#229 E-Gifts
尽管是一道E题,但真心并不很难~不难发现,有一些物品是一定要被选择的,我们所需要决策的仅仅只有那几个有重复价值的物品。
而不同名字之间的概率并不互相影响,所以我们有 \(f[i][j]\) 表示名字为 \(i\) 的物品呼唤 \(j\) 次恰好获得前 \(j\) 大的价值的物品的概率。转移方程为:
\(f[i][j] = f[i][j - 1] * j * \frac{1}{a[i][0]−j+1}\)
为什么要\(*j\) 呢?因为这第 \(j\) 个物品的排列顺序并不是固定的。
要把这 \(n\) 个物品结合起来,我们可以再建立一个 dp 数组,\(g[i][[j]\) 表示前 \(i\) 个名字中,呼唤得到恰好 \(j\) 个有重复价值的物品。我们有转移方程:
\(g[i][j] = \sum g[i - 1][j - 1] * f[i][rec[j] +1]\)
与 \(g[i][j] = \sum g[i - 1][j] * f[i][rec[i]]\)
以上两个分别表示当前名字是否呼唤到一个重复价值的物品。
有没有感觉到有什么不对?没错,在计算的时候,我们的 \(f[i][k]\) 前面是没有带系数的,也就是我们并没有去统计以这样的方式去呼唤的概率是多少。但题目中明确说明当有几种可能呼唤到最高价值的物品时,我们会等概率的任选一种。所以我们可以考虑算出总的方案数 \(c[i][j]\) ,然后再除去这个方案数,即 \(ans =\frac{g[m][cnt]}{c[m][cnt]}\)。这个的转移很简单,可以看一下代码。表面 \(n ^{3}\) ,但第二维的枚举总数限定了范围,所以完全可以承受。
不过我也很好奇……为什么 \(c[i][j]\) 一定要开 double 类型呢?不开就WA了……求解释呀,有知道的还请回复我,私信也可以呀!感激不尽QAQ
#include <bits/stdc++.h>
using namespace std;
#define maxn 2500
#define db long double
int n, m, tot, cnt, rec[maxn];
int a[maxn][maxn], b[maxn];
db f[maxn][maxn], g[maxn][maxn], c[maxn][maxn]; int read()
{
int x = , k = ;
char c; c = getchar();
while(c < '' || c > '') { if(c == '-') k = -; c = getchar(); }
while(c >= '' && c <= '') x = x * + c - '', c = getchar();
return x * k;
} bool cmp(int x, int y) { return x > y; }
void Up(db &x, db y) { x = x + y; } int main()
{
n = read(), m = read();
for(int i = ; i <= m; i ++)
{
a[i][] = read();
for(int j = ; j <= a[i][]; j ++)
a[i][j] = read(), b[++ tot] = a[i][j];
sort(a[i] + , a[i] + + a[i][], cmp);
}
sort(b + , b + + tot, cmp);
for(int i = n; i; i --)
if(b[i] == b[i - ]) cnt ++;
else break;
cnt += ; int K = b[n]; c[][] = ;
for(int i = ; i <= m; i ++)
{
f[i][] = ;
for(int j = ; j <= a[i][]; j ++)
{
if(a[i][j] < K) break;
if(a[i][j] > K) rec[i] = j;
f[i][j] = (db) f[i][j - ] * (db) j * ((db) / (db) (a[i][] - j + ));
}
}
for(int i = , tem = , up = ; i <= m; i ++)
{
int r1 = ;
for(int j = ; j <= up; j ++) c[i][j] = c[i - ][j];
for(int j = rec[i] + ; j <= a[i][]; j ++)
{
if(a[i][j] < K) break;
int t = j - rec[i]; r1 ++;
for(int k = ; k <= up; k ++)
c[i][k + t] = (c[i][k + t] + c[i - ][k]);
}
up += r1;
}
g[][] = ;
for(int i = ; i <= m; i ++)
for(int j = ; j <= cnt; j ++)
{
if(j) Up(g[i][j], g[i - ][j - ] * f[i][rec[i] + ]);
Up(g[i][j], g[i - ][j] * f[i][rec[i]]);
}
cout << fixed << setprecision() << (g[m][cnt] / (db) c[m][cnt]) << endl;
return ;
}
【题解】CF#229 E-Gifts的更多相关文章
- 竞赛题解 - CF Round #524 Div.2
CF Round #524 Div.2 - 竞赛题解 不容易CF有一场下午的比赛,开心的和一个神犇一起报了名 被虐爆--前两题水过去,第三题卡了好久,第四题毫无头绪QwQ Codeforces 传送门 ...
- 题解——CF Manthan, Codefest 18 (rated, Div. 1 + Div. 2) T5(思维)
还是dfs? 好像自己写的有锅 过不去 看了题解修改了才过qwq #include <cstdio> #include <algorithm> #include <cst ...
- 竞赛题解 - [CF 1080D]Olya and magical square
Olya and magical square - 竞赛题解 借鉴了一下神犇tly的博客QwQ(还是打一下广告) 终于弄懂了 Codeforces 传送门 『题目』(直接上翻译了) 给一个边长为 \( ...
- [题解] [CF 1250J] The Parade
题面 题目大意: 给定一个 \(n\) , 所有军人的数量均在 \([1, n]\) 给定 \(a_i\) 代表高度为 \(i\) 的军人的个数 你要将这些军人分成 \(k\) 行, 满足下面两个条件 ...
- 题解 CF 1372 B
题目 传送门 题意 给出 \(n\),输出 \(a\) ,\(b\) (\(0 < a \leq b < n\)),使\(a+b=n\)且 \(\operatorname{lcm}(a,b ...
- 题解——CF Manthan, Codefest 18 (rated, Div. 1 + Div. 2) T4(模拟)
随便模拟下就过了qwq 然后忘了特判WA了QwQ #include <cstdio> #include <algorithm> #include <cstring> ...
- 题解——CF Manthan, Codefest 18 (rated, Div. 1 + Div. 2) T3(贪心)
是一道水题 虽然看起来像是DP,但其实是贪心 扫一遍就A了 QwQ #include <cstdio> #include <algorithm> #include <cs ...
- 题解——CF Manthan, Codefest 18 (rated, Div. 1 + Div. 2) T2(模拟)
题目要求很简单,做法很粗暴 直接扫一遍即可 注意结果会爆int #include <cstdio> #include <algorithm> #include <cstr ...
- 题解——CF Manthan, Codefest 18 (rated, Div. 1 + Div. 2) T1(找规律)
就是找一下规律 但是奈何昨天晚上脑子抽 推错了一项QwQ 然后重新一想 A掉了QwQ #include <cstdio> #include <algorithm> #inclu ...
随机推荐
- LeetCode:34. Search for a Range(Medium)
1. 原题链接 https://leetcode.com/problems/search-for-a-range/description/ 2. 题目要求 给定一个按升序排列的整型数组nums[ ]和 ...
- 二进制描述子 BRIEF(ORB), BRISK, FREAK
二进制描述子设计原则体现在三个部分: 采样pattern 方向orientation compensation 配对sampling pairs ORB基于BRIEF: BRISK是用于OKVIS的描 ...
- DirectX11与DirectX12在古墓丽影暗影中的表现
最近在关注这两个图形API,因为感兴趣,也算是初学者. 以下内容仅供参考. 使用古墓丽影暗影游戏,分别对这两个进行比较,得出的结论如下图(此笔记本散热很差,更改散热应该比下图结果好些): 首先看可以很 ...
- 【icon】 图标组件说明
小程序默认了几种类型图标,其组件原型如下: <icon type="[success | success_no_circle | info | warn | waiting | can ...
- 孤荷凌寒自学python第八十三天初次接触ocr配置tesseract环境
孤荷凌寒自学python第八十三天初次接触ocr配置tesseract环境 (完整学习过程屏幕记录视频地址在文末) 学习Python我肯定不会错过图片文字的识别,当然更重要的是简单的验证码识别了,今天 ...
- Open MPI集群运行
部署完之后,代码也能正确跑起来了,也确实集群分散了.跑一下各种各样的代码,发现了一个错误: $ ~/OpenMpi/bin/mpiexec -np ~/NetWorkTest My rank is M ...
- CDH/Hadoop 5.15 installation steps
I will talk the main steps to install CDH 5.15 on Linux(CENT OS 6.10). The installation method is M ...
- 2.安装hdfs yarn
下载hadoop压缩包设置hadoop环境变量设置hdfs环境变量设置yarn环境变量设置mapreduce环境变量修改hadoop配置设置core-site.xml设置hdfs-site.xml设置 ...
- Alpha 冲刺(3/10)
队名 火箭少男100 组长博客 林燊大哥 作业博客 Alpha 冲鸭鸭鸭! 成员冲刺阶段情况 林燊(组长) 过去两天完成了哪些任务 协调各成员之间的工作 协助后端界面的开发 搭建项目运行的服务器环境 ...
- 从零讲JAVA ,给你一条 清晰地学习道路!该学什么就学什么!!
1.计算机基础: 1.1数据机构基础: 主要学习:1.向量,链表,栈,队列和堆,词典.熟悉2.树,二叉搜索树.熟悉3.图,有向图,无向图,基本概念4.二叉搜索A,B,C类熟练,9大排序熟悉.5.树的前 ...