ICPC Northwestern European Regional Programming Contest B Bottle Flip

Paddi930

已于 2024-10-05 17:19:49 修改

阅读量730

点赞数 9

文章标签： c++

于 2024-10-05 17:18:39 首次发布

本文链接：https://blog.csdn.net/2302_80542709/article/details/142716588

版权

ICPC Northwestern European Regional Programming Contest B Bottle Flip

题意：

给你一个高度为 $H$ ,半径为 $r$ 的水杯，并告诉你空气密度 $\rho_{a}$ 和水的密度 $\rho_{w}$ ,求使得整个杯子装水后质心最低的水量

大致思路：

前置物理知识：如果处理多个部分的重心比较困难，但整体处理比较简单时，我们可以将求出这一部分物体的质心和质量，这两部分在处理整个物体的质心时时等价的。

即若把物体分成 $a, b$ 两个部分有：
$X=\frac{\sum_{i=0}^k x_{i}m_{i}}{\sum_{i=1}^k m_{i}}\\ X_{a}={{\sum_{i=1}^{k_a} x_i m_i}\over \sum_{i=1}^{k_a} m_i}\\ X_{b}={{\sum_{i=1}^{k_b} x_i m_i}\over \sum_{i=1}^{k_b} m_i}\\ M_a=\sum_{i=1}^{k_a} m_i\\ M_b=\sum_{i=1}^{k_b} m_i\\ X={X_aM_a+X_bM_b\over M_a+M_b}$
我们设水的高度为 $h$ ，那么重心高度则有
$X_w(h)= {h\over 2}\\ X_a(h)={H+h \over 2}\\ M_a(h)= S(H-h)\rho_a\\ M_w(h)=Sh\rho_w\\ X(h)=\frac{{1 \over 2} (Sh^2 \rho_w+(H^2-h^2) \rho_a)}{S(H-h)\rho_a+Sh\rho_w}= {1 \over 2}\frac{ (h^2 \rho_w+(H^2-h^2) \rho_a)}{(H-h)\rho_a+h\rho_w}$
数学较好的同学可以将 $X (h)$ 对 $h$ 求导。

另外，我们可以发现 $X (h)$ 是一个单峰函数，故三分可以通过本题

代码：

#include <bits/stdc++.h>
using namespace std;
#define int long long int
double H, y, Pa, Pw;
double eps = 1e-10;
double check(double h)
{
    double a1 = (h * h * Pw + (H * H - h * h) * Pa) * 0.5;
    double w = h * Pw + (H - h) * Pa;
    return a1 / w;
}
signed main()
{
    cin >> H >> y >> Pa >> Pw;
    double l = 0, r = H;
    while (r - l >= eps)
    {
        double midl, midr;
        midl = l + (r - l) / 3;
        midr = r - (r - l) / 3;
        if (check(midl) > check(midr))
            l = midl;
        else
            r = midr;
    }

    cout << fixed << setprecision(10) << l << endl;
    return 0;
}