昨晚乱入学弟的训练赛,想了一下这个题。推导的过程中,加深了对公理化的概率论理解。$\newcommand{\d}{\mathop{}\!\mathrm{d}}$


解法一

考虑 $ d < L$ 的情形。
\begin{equation*}
P(X = 1) = \frac{d}{L}
\end{equation*}
\begin{align*}
P(X = 2) &= \int_0^{L - d} \frac{\d x}{L} \frac{d}{L - x} \\
&= \frac{d}{L}\ln\frac{L}{d}
\end{align*}
\begin{align*}
P(X = 3) &= \int_{0}^{L-d}\frac{\d x}{L}\int_{0}^{L-x-d}\frac{\d y}{L-x}\frac{d}{L-x-y} \\
&= \int_{0}^{L-d} \frac{\d x}{L} \frac{d}{L - x} \ln \frac{L - x}{d} \\
&= \frac{d}{L}\frac{1}{2}\ln^2\frac{L}{d}
\end{align*}
\begin{align}
P(X = 4) &= \int_0^{L - d}\frac{\d x}{L} \int_0^{L - x - d} \frac{\d y}{L-x}\int_0^{L-x-y-d}\frac{\d z}{L - x -y}\frac{d}{L - x - y - z} \notag\\
&= \int_0^{L - d}\frac{\d x}{L} \int_0^{L - x - d} \frac{\d y}{L-x} \frac{d}{L - x -y} \ln \frac{L - x -y}{d} \notag\\
&= \int_0^{L - d}\frac{\d x}{L} \frac{d}{L-x}\frac{1}{2}\ln^2\frac{L-x}{d} \label{Int:1}
\end{align}
令 $u = \frac{L-x}{d}$ ,则 $\d x = -d\d u$ ,有
\begin{align*}
\eqref{Int:1} &= \int_1^\frac Ld\frac dL\frac{\d u}{u}\frac 12\ln^2u \\
&= \int_1^\frac Ld\frac dL\frac 16\d\ln^3u \\
&= \frac 16\frac dL\ln^3\frac Ld
\end{align*}
不难推出
\begin{equation*}
P(X = n) = \frac dL\frac1{(n-1)!}\ln^{n-1}\frac Ld
\end{equation*}
所以
\begin{align*}
E(X) &= \sum_{ n \ge 1 } n P(X=n) \\
&= \frac dL \sum_{n \ge 1} \frac n{(n-1)!}\ln^{n-1}\frac Ld \\
&= \frac dL \sum_{n \ge 0} \frac{n+1}{n!} \ln^n\frac Ld \\
&= \frac dL (\ln\frac Ld + 1) \mathrm{e}^{\ln\frac Ld} \\
&= \ln\frac Ld + 1
\end{align*}
上式中的求和用到了 $(x+1)\mathrm{e}^x$ 的 Maclaurin 展开:
\begin{equation*}
(x+1)\mathrm{e}^x = \sum_{n\ge 0} \frac{n + 1}{n!} x^n
\end{equation*}

解法二

用 $f(x)$ 表示绳长为 $x$ 时切割次数的期望,则有
$$
f(x) =
\begin{cases}
0, &&\text{if $x\le d$;} \\
1 + \int_0^x \frac{dy}{x}f(y), && \text{otherwise.}
\end{cases}
$$
考虑 $x>d$ 的情形,此时有
\begin{align}
f(x) &= 1 + \int_0^x \frac{\d y}{x}f(y) \notag\\
&= 1 + \int_0^d \frac{\d y}{x}f(y) + \int_d^x \frac{\d y}{x}f(y) \notag\\
&= 1 + \int_d^x \frac{\d y}{x}f(y) \label{Int:2}
\end{align}
对 \eqref{Int:2} 式两边求导,得
\begin{align*}
f'(x) &= \frac{f(x)}x - \frac1{x^2}\int_d^x\d yf(y) \\
&= \frac{f(x)}x - \frac1x(f(x) -1) \\
&= \frac1x
\end{align*}
又 $\lim\limits_{x\to d^+} f(x) = 1 $,得 $$ f(x) = \ln x + 1 - \ln d $$


解法二来自 Huo Chen

2016 年 ACM/ICPC 青岛区域赛 Problem C Pocky的更多相关文章

  1. 2015 ACM / ICPC 亚洲区域赛总结(长春站&北京站)

    队名:Unlimited Code Works(无尽编码)  队员:Wu.Wang.Zhou 先说一下队伍:Wu是大三学长:Wang高中noip省一:我最渣,去年来大学开始学的a+b,参加今年区域赛之 ...

  2. ACM/ICPC2016 青岛区域赛

    A(hdu5982).(模拟) 题意:输入n对数,将每对数相乘并相加 分析:模拟 B(hdu5983).(模拟) 题意:给你一个二阶魔方,问能否通过一次旋转使得给定魔方的每个面颜色相同 分析:模拟 C ...

  3. Known Notation括号匹配类问题(2014年ACM/ICPC 亚洲区域赛牡丹江)

    题意: 给你数字或 * 的串,你可以交换一个*和数字.在最前面添1.在一个地方插入*,问你使串满足入栈出栈的(RNP)运算法则. 思路: 引用:https://blog.csdn.net/u01158 ...

  4. Digit sum (第 44 届 ACM/ICPC 亚洲区域赛(上海)网络赛)进制预处理水题

    131072K   A digit sum S_b(n)Sb​(n) is a sum of the base-bb digits of nn. Such as S_{10}(233) = 2 + 3 ...

  5. 2016年11月ACM/ICPC亚洲区北京赛赛后总结

    2016年11月12到11月13为期两天的比赛,这是我们这个对第一次去打亚洲区域赛,经过这次比赛,我认识到了自己与别人的差距,也许我们与别人的起点不同,但这不是理由. 这次的比赛12号的热身赛两点开始 ...

  6. 2019 ICPC 上海区域赛总结

    2019上海区域赛现场赛总结 补题情况(以下通过率为牛客提交): 题号 标题 已通过代码 通过率 我的状态 A Mr. Panda and Dominoes 点击查看 5/29 未通过 B Prefi ...

  7. 2014ACM/ICPC亚洲区域赛牡丹江站汇总

    球队内线我也总水平,这所学校得到了前所未有的8地方,因为只有两个少年队.因此,我们13并且可以被分配到的地方,因为13和非常大的数目.据领队谁oj在之上a谁去让更多的冠军.我和tyh,sxk,doub ...

  8. 2016青岛区域赛.Coding Contest(费用流 + 概率计算转换为加法计算)

    Coding Contest Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)To ...

  9. 2018 ACM-ICPC 亚洲区域赛青岛现场赛 —— Problem F. Tournament

    题面:http://acm.zju.edu.cn/contest-materials/qd2018/qd2018_problems.pdf 题意: n个骑士决斗K轮 要求是每个骑士只能跟另外一个骑士决 ...

随机推荐

  1. 如何处理错误消息Please install the Linux kernel header files

    Please install the Linux kernel "header" files matching the current kernel 当我启动minilkube时遇 ...

  2. Vector 容器简单介绍

    # Vector STL简要介绍 关于STL中的vector容器,以下做一些相关介绍. #### vector 简要概述 vector 称作向量类,属于容器类,实现了动态的数组,用于元素数量变化的对象 ...

  3. anaconda 安装各种库

    在anaconda prompt 添加清华源 https://mirrors.tuna.tsinghua.edu.cn/help/anaconda/ conda config --add channe ...

  4. Beta冲刺(周四)

    这个作业属于哪个课程 https://edu.cnblogs.com/campus/xnsy/SoftwareEngineeringClass1 这个作业要求在哪里 https://edu.cnblo ...

  5. CentOS配置主机名和主机映射

    1.修改本机主机名 vi /etc/sysconfig/network 修改hostname HOSTNAME=s0 2.配置主机映射 vi /etc/hosts 修改内容如下 192.168.32. ...

  6. SC || Chapter 5 复习向

    可复用性 ┉┉∞ ∞┉┉┉┉∞ ∞┉┉┉∞ ∞┉┉ 行为子结构 对于父子的继承关系的要求: ·子类可以增加方法,但不可以删 ·子类需实现抽象类型中未实现的方法 ·子类重写(override)的方法必须 ...

  7. $|^|\z|\Z|/a|/l

    #!/usr/bin/perl use strict; use warnings; foreach(<>) { if (/(\w*)/a){print "$1\n";} ...

  8. VS连接SQL Server 2008,并实现登录和注册功能

    --------------------- 作者:Cambridge 来源:CSDN 原文:https://blog.csdn.net/cambridgeacm/article/details/797 ...

  9. pycharm安装 suds模块报错:AttributeError: module 'pip' has no attribute 'main'

    需求:安装suds模块 遇到的问题: 一.报错信息:[file][Default Settint]---Project Interpreter 点击 搜索suds安装模块报错 解决:依据上图提示找到C ...

  10. 03等待多个线程返回WaitForMultipleObject

    二. WaitForMultipleObject 等待单个线程返回 1. 函数原型 DWORD WINAPI WaitForMultipleObjects( _In_ DWORD nCount, _I ...