2020ICPC·小米网络选拔赛热身赛K-Random Point in Triangle

最新推荐文章于 2021-09-27 19:55:57 发布

狙击美佐

最新推荐文章于 2021-09-27 19:55:57 发布

阅读量464

点赞数 1

分类专栏：计算几何文章标签： 1024程序员节

本文链接：https://blog.csdn.net/fztsilly/article/details/109260850

版权

计算几何专栏收录该内容

2 篇文章 0 订阅

订阅专栏

该博客讨论了一道计算几何题目，涉及在三角形内部随机选取一点，计算其与三角形三边形成的三角形最大面积期望。首先通过分析等边三角形的情况，得出期望值，然后推广到任意三角形，尽管未给出证明，但提出了大胆的猜想。博主给出了等边三角形的解析解，并提供了C++代码实现。

摘要由CSDN通过智能技术生成

2020ICPC·小米网络选拔赛热身赛K-Random Point in Triangle 计算几何

- 题意

#传送门： https://ac.nowcoder.com/acm/contest/8409/K

题意

$给3个点(x1,y1),(x2,y2),(x3,y3)，在三角形内部中取一点p，设E=max(S_{PAB},S_{PAC},S_{PBC}),求E的概率期望。$

$对于任何三角形，先看看特殊的等边三角形。然后在推广到任意三角形（ A C M 需要猜！）$
在这里插入图片描述

$O 点为三角形 A B C 的内心。所以对于 p 的位置的选取，只需要在三角形的一半，然后一半即可。即在三角形 O S B 中选取。$

$然后构造下图为：$
在这里插入图片描述
$设该等边三角形的边长为 2 a ， O N 为 y ， N P 为 x ，则根据相似可得：$

$PH=\frac{2}{3}a+x+\frac{\sqrt 3}{3}y$
$PG=\frac{\sqrt 3}{2}PH=\frac{\sqrt3}{3}a+\frac{\sqrt 3}{2}x+\frac{1}{2}y$

$则E=S_{PAC}=\frac{1}{2}*AH*PG=\frac{\sqrt3}{3}a^2+\frac{\sqrt 3}{2}ax+\frac{1}{2}ay$

$根据上图，P点的y可在[0，\frac{\sqrt a}{3}a]中选取，所以x可以在[0，\sqrt 3y]中选取。$

$对面积积分如下：$

$\int_{0}^{\frac{\sqrt a}{3}}dy\int_{0}^{\sqrt 3y}(\frac{\sqrt3}{3}a^2+\frac{\sqrt 3}{2}ax+\frac{1}{2}ay)dx=\frac{11}{36}a^4$

$而P点在该区域的选取面积为\frac{\sqrt 3}{6}a^2，期望为\frac{11}{18}\sqrt 3 a^2，而等边三角形的面积为\sqrt 3a^2$
$即期望为\frac{11}{18}S_{ABC}，在乘以36得22S_{ABC}。$

$把等边三角形推广到任意三角形，虽然我不会证明，大胆猜测总没错！$

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef long double ld;
typedef pair<int, int> pdd;

#define INF 0x3f3f3f3f
#define lowbit(x) x & (-x)
#define mem(a, b) memset(a , b , sizeof(a))
#define FOR(i, x, n) for(int i = x;i <= n; i++)

// const ll mod = 998244353;
 const ll mod = 1e9 + 7;
// const double eps = 1e-6;
// const double PI = acos(-1);
// const double R = 0.57721566490153286060651209;

void solve() {
    ll x1, x2, x3, y1, y2, y3;
    while(cin >> x1 >> y1 >> x2 >> y2 >> x3 >> y3) {
        cout <<  abs(x1 * y2 - x2 * y1 + x3 * y1 - x1 * y3 + x2 * y3 - x3 * y2) * 11 << endl;
    }
}

signed main() {
    ios_base::sync_with_stdio(false);
    //cin.tie(nullptr);
    //cout.tie(nullptr);
#ifdef FZT_ACM_LOCAL
    freopen("in.txt", "r", stdin);
    freopen("out.txt", "w", stdout);
    signed test_index_for_debug = 1;
    char acm_local_for_debug = 0;
    do {
        if (acm_local_for_debug == '$') exit(0);
        if (test_index_for_debug > 20)
            throw runtime_error("Check the stdin!!!");
        auto start_clock_for_debug = clock();
        solve();
        auto end_clock_for_debug = clock();
        cout << "Test " << test_index_for_debug << " successful" << endl;
        cerr << "Test " << test_index_for_debug++ << " Run Time: "
             << double(end_clock_for_debug - start_clock_for_debug) / CLOCKS_PER_SEC << "s" << endl;
        cout << "--------------------------------------------------" << endl;
    } while (cin >> acm_local_for_debug && cin.putback(acm_local_for_debug));
#else
    solve();
#endif
    return 0;
}