UVa 13099 Tobby and the Line Game

原创已于 2025-11-21 10:25:27 修改 · 777 阅读

12 ·

本内容遵循CC 4.0 BY-SA版权协议

GEO检测

标签

#算法 #c++ #青少年编程

于 2025-10-24 21:51:21 首次发布

组合数学专栏收录该内容

96 篇文章

订阅专栏

题目分析

问题描述

$Tobby\texttt{Tobby}$ 参加一个在矩形场地上进行的游戏。场地左下角坐标为 $x_L, y_L)$ ，右上角坐标为 $x_R, y_R)$ 。游戏规则是：每个参与者投掷两枚硬币，硬币落点均匀分布在矩形区域内（包括边界），且两次投掷相互独立。游戏的获胜者是两枚硬币落点之间线段长度最长的参与者。

$Tobby\texttt{Tobby}$ 想知道他获胜的机会有多大，因此需要计算两枚硬币落点之间距离平方的期望值：

$E[∣∣p1−p2∣∣2]=E[(x1−x2)2+(y1−y2)2]E[||\mathbf{p}_1 - \mathbf{p}_2||^2] = E[(x_1 - x_2)^2 + (y_1 - y_2)^2]$

关键观察

独立性： $x_1$ 与 $x_2$ 独立， $y_1$ 与 $y_2$ 独立，且 $x$ 坐标与 $y$ 坐标相互独立
均匀分布：所有坐标在各自区间内均匀分布
可分离性：期望运算可以分解为 $x$ 方向和 $y$ 方向分别计算

数学推导

$x$ 方向的计算

设 $X1,X2∼Uniform[a,b]X_1, X_2 \sim \text{Uniform}[a,b]$ ，其中 $a = x_L$ ， $b = x_R$ 。

均匀分布的均值和方差为：

均值： $\frac{a+b}{2}$
方差： $Var(X)=(b−a)212\text{Var}(X) = \frac{(b-a)^2}{12}$

对于两个独立随机变量 $X_1$ 和 $X_2$ ：
$E[(X1−X2)2]=Var(X1−X2)+[E(X1−X2)]2E[(X_1 - X_2)^2] = \text{Var}(X_1 - X_2) + [E(X_1 - X_2)]^2$

由于 $E[X_1 - X_2] = 0$ ，且方差具有可加性：
$Var(X1−X2)=Var(X1)+Var(X2)=(b−a)212+(b−a)212=(b−a)26\text{Var}(X_1 - X_2) = \text{Var}(X_1) + \text{Var}(X_2) = \frac{(b-a)^2}{12} + \frac{(b-a)^2}{12} = \frac{(b-a)^2}{6}$

因此：
$E[(x1−x2)2]=(xR−xL)26E[(x_1 - x_2)^2] = \frac{(x_R - x_L)^2}{6}$

$y$ 方向的计算

同理可得：
$E[(y1−y2)2]=(yR−yL)26E[(y_1 - y_2)^2] = \frac{(y_R - y_L)^2}{6}$

最终公式

将两个方向的结果相加：
$E[distance2]=(xR−xL)2+(yR−yL)26E[\text{distance}^2] = \frac{(x_R - x_L)^2 + (y_R - y_L)^2}{6}$

算法设计

基于上述推导，算法非常简单：

读取每个测试用例的 $x_L, y_L, x_R, y_R$
计算矩形宽度 $dx = x_R - x_L$ 和高度 $dy = y_R - y_L$
应用公式： $\frac{dx^2 + dy^2}{6}$
输出结果，精度满足 $10^{-5}$ 的要求

复杂度分析

时间复杂度： $O (1)$ 每个测试用例
空间复杂度： $O (1)$

样例验证

样例 1

输入：-100 -50 23 400
dx = 123, dy = 450
结果 = (123² + 450²) / 6 = (15129 + 202500) / 6 = 217629 / 6 = 36271.5

样例 2

输入：0 0 500 1
dx = 500, dy = 1
结果 = (500² + 1²) / 6 = (250000 + 1) / 6 = 250001 / 6 ≈ 41666.8333333333

代码实现

// Tobby and the Line Game
// UVa ID: 13099
// Verdict: Accepted
// Submission Date: 2025-10-24
// UVa Run Time: 0.000s
//
// 版权所有（C）2025，邱秋。metaphysis # yeah dot net

#include <iostream>
#include <cstdio>

using namespace std;

int main() {
    int xL, yL, xR, yR;
    while (scanf("%d%d%d%d", &xL, &yL, &xR, &yR) != EOF) {
        double dx = xR - xL;
        double dy = yR - yL;
        double result = (dx * dx + dy * dy) / 6.0;
        printf("%.15f\n", result);
    }
    return 0;
}