一个有趣的问题(随缘更新)

已于 2022-12-23 19:04:01 修改
阅读量225 收藏
点赞数 7
文章标签： c语言 算法
于 2022-11-13 07:50:40 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/only_me_canthis/article/details/127827988
版权

736aaaceca604d1f880a218b9d9dd100.jpg

 最近我发现了一个有趣的问题，问题如上。问题简化一下，大概就是一圆内有四个随机生成的点，求四个点在同一半圆的概率。

问题分析:

  1. 怎么判断四点是否在同一半圆内？
  2. 怎么判断一个点是否在一个凸多边形内？
  3. 怎么通过四个点生成一个凸多边形？
  4. 怎么定量判断一点是否在一个三角形内？

先解决第一个问题，那我们都知道，任意一天圆的直径都能将一个圆平分，而任意一条直径都过圆心，所以只要圆心不在四点所围成的凸多边形，那为什么是凸多边形呢？后面会有解答。

解决第一个问题后，随之衍生出第二个问题，先明确凸多边形是啥意思？凸多边形（Convex Polygon）指如果把一个多边形的所有边中，任意一条边向两方无限延长成为一直线时，其他各边都在此直线的同旁，那么这个多边形就叫做凸多边形。简单的说，就是内角都得小于180才叫凸多边形。

8355de56f0cb48eba432ca1d0b299e49.png

 再看一什么是非凸多边形

d5108568467e40bbaaf2016b3bcbc446.jpg

 可以看到，图中圆内四点所形成的就是一个非凸多边形，此时圆心确实也在该多边形的外面，但是却不能找到一条直径让这四点分布在同一侧，所以我们强调要凸多边形。

解决了上述问题，那就要思考一下给四个点怎么能生成一个凸多边形？我就想到，任一多边形都能由几个三角形拼接而成。所以，任给四个点，我们可以先生成C43也就是4个三角形，然后把这四个三角形围成的区域并起来就可以形成一个凸多边形，当然后面也证实这确实是我们所需要的凸多边形。类似下面这样。

b067851988334eda99a7de5592443ab5.jpg

 上图中，圆内四点ABCD，圆心O，然后是有ABC，ABD，ACD，BCD四个三角形，可以看到这四个三角形区域并起来围成的区域是三角形ABD，圆心位于该多边形内(三角形也是凸多边形)，这时就得到了正确的结论。

到现在问题已经转化到如何判断一点是否在一个三角形内了。因为，圆内生成的四点可以三三生成四个三角形，那么我们只需要保证圆心不在任意一个生成的三角形内就可以保证该四点可被一条直径划分到同一侧，即全部分布在一个半圆内。当然相反地，圆心只要分布在至少一个三角形内就说明没有一条直径能够将这四点划分在同一侧，即非全部分布在一个半圆内。那么到底该怎么判断勒？看下图

61fea6c5963f48b8ad76d794ac274d20.jpg

 上图呢，是一个任意的三角形ABC和内外两点P1，P2，将P1，P2都与三角形的三个顶点相连又可以得到三个三角形，观察得到，内点P1形成的三个三角形面积之和就等于ABC的面积，而外点P2形成的三个三角形面积之和比ABC的面积大一个P2BC的面积。所以得出结论就是只需要计算点与三角形顶点所形成的三个三角形面积之和再与原三角形面积比较即可得到该点是否在三角形内。

写到这里，如果不用概率论与数理统计的知识去解决，而单纯的用计算机模拟的话，问题已经变得很简单了。继问题四往后推，也就是需要实现一个求三角形面积的函数，实现如下
 

double Striangle(double *A,double *B,double *C)

{

    //输入三点ABC坐标(二维数组)

    //输出所围三角形面积

    double deta_x1,deta_x2,deta_y1,deta_y2;

    double Ax,Ay,Bx,By,Cx,Cy;

    double len_a,len_b,sineC;

    double S=-1;

    Ax=A[0];

    Ay=A[1];

    Bx=B[0];

    By=B[1];

    Cx=C[0];

    Cy=C[1];

    deta_x1=(Ax-Bx);

    deta_x2=(Ax-Cx);

    deta_y1=(Ay-By);

    deta_y2=(Ay-Cy);

    if(deta_x1*deta_x2!=0 && deta_y1/deta_x1==deta_y2/deta_x2) {

        return 0;

    } else if(deta_x1==0 && deta_y1==0) {

        return 0;

    } else if(deta_x2==0 && deta_y2==0) {

        return 0;

    } else {

        len_a=length(B,C);

        len_b=length(A,C);

        sineC=sine_angle(A,B,C);

        S=0.5*len_a*len_b*sineC;

        return S;

    }

}

可以看到在运行这个函数之前，还需要再实现两个函数求两点距离和三角形某个角正弦值，实现如下

double length(double *A,double *B)

{

    //输入两点AB坐标(二维数组)

    //输出两点距离(模长)

    double Ax,Ay,Bx,By,R;

    Ax=A[0];

    Ay=A[1];

    Bx=B[0];

    By=B[1];

    R=sqrt(pow(Ax-Bx,2) + pow(Ay-By,2));

    return R;

}

double sine_angle(double *A,double *B,double *C)

{

    //输入三点ABC坐标(二维数组)

    //输出角C的正弦值

    double len_a,len_b,len_c;

    double cosine2ab,cosine,sine;

    len_a=length(B,C);

    len_b=length(A,C);

    len_c=length(A,B);

    cosine2ab=pow(len_a,2)+pow(len_b,2)-pow(len_c,2);

    if(len_a*len_b!=0) {

        cosine=cosine2ab / (2*len_a*len_b);

    } else {

        cosine=1;

    }

    sine=sqrt(1-pow(cosine,2));

    return sine;

}

到这里就没啥难度了，下面再说一下怎么随机生成四点坐标，实现如下

void rand_circle_coord(*CC)

{

    //输入一个空二维数组

    //将随机数赋值给传入的指针

    //srand(time(0));

    CC[0]=rand()%(2*eps+1) - eps;

    CC[1]=rand()%(2*eps+1) - eps;

}



void rand_circle_coords(double *A,double *B,double *C,double *D)

{

    //输入四个二维数组

    //将随机坐标复制给传入的四个指针

    int intA[2],intB[2],intC[2],intD[2];

    double O[2]= {0,0};



    //srand(time(0));

    rand_circle_coord(intA);

    A[0]=intA[0] / (double)eps;

    A[1]=intA[1] / (double)eps;

    while(length(A,O)>=1.00) {

        rand_circle_coord(intA);

        A[0]=intA[0] / (double)eps;

        A[1]=intA[1] / (double)eps;

    }

    //printf("A是(%lf,%lf)\n",A[0],A[1]);



    rand_circle_coord(intB);

    B[0]=intB[0] / (double)eps;

    B[1]=intB[1] / (double)eps;

    while(length(A,O)>=1.00) {

        rand_circle_coord(intB);

        B[0]=intB[0] / (double)eps;

        B[1]=intB[1] / (double)eps;

    }

    //printf("B是(%lf,%lf)\n",B[0],B[1]);



    rand_circle_coord(intC);

    C[0]=intC[0] / (double)eps;

    C[1]=intC[1] / (double)eps;

    while(length(C,O)>=1.00) {

        rand_circle_coord(intC);

        C[0]=intC[0] / (double)eps;

        C[1]=intC[1] / (double)eps;

    }

    //printf("C是(%lf,%lf)\n",C[0],C[1]);



    rand_circle_coord(intD);

    D[0]=intD[0] / (double)eps;

    D[1]=intD[1] / (double)eps;

    while(length(D,O)>=1.00) {

        rand_circle_coord(intD);

        D[0]=intD[0] / (double)eps;

        D[1]=intD[1] / (double)eps;

    }

    //printf("D是(%lf,%lf)\n",D[0],D[1]);

}

到这里基本的一些函数就全部实现好了，剩下的就是将这些函数联合起来用，达到模拟的效果，整个代码如下

#include<stdio.h>
#include<math.h>
#include<stdlib.h>
#include<time.h>

#define eps 1000000

double length(double *A,double *B);
double sine_angle(double *A,double *B,double *C);
double Striangle(double *A,double *B,double *C);
int judgeP_in_triangle(double *A,double *B,double *C,double *P);
int judgeO_in_triangles(double *A,double *B,double *C,double *D);
void rand_circle_coord(int *CC);
void rand_circle_coords(double *A,double *B,double *C,double *D);

int main()
{
    double A[2];
    double B[2];
    double C[2];
    double D[2];
    int i,j,P=0,R[100],epochs=1000;
    srand(time(0));
    for(j=0; j<100; j++) {
        R[j]=0;
        for(i=0; i<epochs; i++) {
            rand_circle_coords(A,B,C,D);
            if(judgeO_in_triangles(A,B,C,D)) {
                R[j]+=1;
            }
        }
    }
    for(j=0;j<100;j++){
        P+=R[j];
    }

    printf("\n结果是%lf",(double)P / (epochs*100));
    return 0;
}

double length(double *A,double *B)
{
    //输入两点AB坐标(二维数组)
    //输出两点距离(模长)
    double Ax,Ay,Bx,By,R;
    Ax=A[0];
    Ay=A[1];
    Bx=B[0];
    By=B[1];
    R=sqrt(pow(Ax-Bx,2) + pow(Ay-By,2));
    return R;
}

double sine_angle(double *A,double *B,double *C)
{
    //输入三点ABC坐标(二维数组)
    //输出角C的正弦值
    double len_a,len_b,len_c;
    double cosine2ab,cosine,sine;
    len_a=length(B,C);
    len_b=length(A,C);
    len_c=length(A,B);
    cosine2ab=pow(len_a,2)+pow(len_b,2)-pow(len_c,2);
    if(len_a*len_b!=0) {
        cosine=cosine2ab / (2*len_a*len_b);
    } else {
        cosine=1;
    }
    sine=sqrt(1-pow(cosine,2));
    return sine;
}

double Striangle(double *A,double *B,double *C)
{
    //输入三点ABC坐标(二维数组)
    //输出所围三角形面积
    double deta_x1,deta_x2,deta_y1,deta_y2;
    double Ax,Ay,Bx,By,Cx,Cy;
    double len_a,len_b,sineC;
    double S=-1;
    Ax=A[0];
    Ay=A[1];
    Bx=B[0];
    By=B[1];
    Cx=C[0];
    Cy=C[1];
    deta_x1=(Ax-Bx);
    deta_x2=(Ax-Cx);
    deta_y1=(Ay-By);
    deta_y2=(Ay-Cy);
    if(deta_x1*deta_x2!=0 && deta_y1/deta_x1==deta_y2/deta_x2) {
        return 0;
    } else if(deta_x1==0 && deta_y1==0) {
        return 0;
    } else if(deta_x2==0 && deta_y2==0) {
        return 0;
    } else {
        len_a=length(B,C);
        len_b=length(A,C);
        sineC=sine_angle(A,B,C);
        S=0.5*len_a*len_b*sineC;

        return S;
    }
}

int judgeP_in_triangle(double *A,double *B,double *C,double *P)
{
    int flag;
    double Spab,Spac,Spbc,Sabc,deta_S;
    Spab=Striangle(A,B,P);
    Spac=Striangle(A,C,P);
    Spbc=Striangle(B,C,P);
    Sabc=Striangle(A,B,C);
    deta_S=Spab+Spac+Spbc-Sabc;
    if(deta_S<=1e-06) {
        flag=1;
    } else {
        flag=0;
    }

    return flag;
}

int judgeO_in_triangles(double *A,double *B,double *C,double *D)
{
    double O[2]= {0,0};
    int flag1,flag2,flag3,flag4;
    flag1=!judgeP_in_triangle(A,B,C,O);
    flag2=!judgeP_in_triangle(A,B,D,O);
    flag3=!judgeP_in_triangle(A,C,D,O);
    flag4=!judgeP_in_triangle(B,C,D,O);
    if(flag1 && flag2 && flag3 && flag4) {
        return 1;
    } else {
        return 0;
    }
}

void rand_circle_coord(int *CC)
{
    //输入一个空二维数组
    //将随机数赋值给传入的指针
    //srand(time(0));
    CC[0]=rand()%(2*eps+1) - eps;
    CC[1]=rand()%(2*eps+1) - eps;
}

void rand_circle_coords(double *A,double *B,double *C,double *D)
{
    //输入四个二维数组
    //将随机坐标复制给传入的四个指针
    int intA[2],intB[2],intC[2],intD[2];
    double O[2]= {0,0};

    //srand(time(0));
    rand_circle_coord(intA);
    A[0]=intA[0] / (double)eps;
    A[1]=intA[1] / (double)eps;
    while(length(A,O)>=1.00) {
        rand_circle_coord(intA);
        A[0]=intA[0] / (double)eps;
        A[1]=intA[1] / (double)eps;
    }
    //printf("A是(%lf,%lf)\n",A[0],A[1]);

    rand_circle_coord(intB);
    B[0]=intB[0] / (double)eps;
    B[1]=intB[1] / (double)eps;
    while(length(A,O)>=1.00) {
        rand_circle_coord(intB);
        B[0]=intB[0] / (double)eps;
        B[1]=intB[1] / (double)eps;
    }
    //printf("B是(%lf,%lf)\n",B[0],B[1]);

    rand_circle_coord(intC);
    C[0]=intC[0] / (double)eps;
    C[1]=intC[1] / (double)eps;
    while(length(C,O)>=1.00) {
        rand_circle_coord(intC);
        C[0]=intC[0] / (double)eps;
        C[1]=intC[1] / (double)eps;
    }
    //printf("C是(%lf,%lf)\n",C[0],C[1]);

    rand_circle_coord(intD);
    D[0]=intD[0] / (double)eps;
    D[1]=intD[1] / (double)eps;
    while(length(D,O)>=1.00) {
        rand_circle_coord(intD);
        D[0]=intD[0] / (double)eps;
        D[1]=intD[1] / (double)eps;
    }
    //printf("D是(%lf,%lf)\n",D[0],D[1]);
}

这里我们模拟1000次，重复实验100次然后取平均结果如下

56738624b56841b5a3a6771406075ab6.jpg

 大概是0.5的概率，一半一半吧。

好了。整个问题我的想法就是这样了。如有错误欢迎指正，一起进步，加油💪💪💪

SwayDy_水滴
关注
  • 7
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
随缘网络留言板 v2.0 build 20120127
10-18
留言板自带后台，可批量管理留言，适合各中小企业将其整合植入网站或者做为一个单独留言系统交流使用。 随缘网络留言板v2.0 2012-1-27更新： 1、针对前一段时间网站明文密码安全问题对密码进行了md5加密。 2、...
记录ubuntu中遇到的问题（待续，随缘更新）
ping96的博客
03-30 97
1. ubuntu系统卡死后无法正常开机 2.谷歌浏览器的安装
【计算机网络】章节重点、习题问题（随缘更新）
人生若只如初见
03-26 998
方便自己，也方便他人
懒狗日记（随缘更新）
ch216003的博客
11-18 643
由于是新的服务器,为了更快的安装，需要将源换成国内源，具体操作可以参考【Linux教程】Ubuntu Linux 更换源教程老板要求点云数据，本来想使用python处理，但是没找到好的代码，反而c++的代码很多，因此被迫营业～由于是新买的服务器，里面的环境和配置基本上没有安装，所以有许多的安装步骤，希望对你的研究有所帮助。
[随缘一题]-LintCode-螺旋矩阵
呼延十
12-27 586
PS 本题代码来源于九章算法. 来源 lintcode-螺旋矩阵 描述 给出整数 n, 返回一个大小为 n * n 的螺旋矩阵 样例 给出 n = 3 则螺旋矩阵为: [ [1,2,3] [8,9,4] [7,6,5] ] 给出 n = 5 则螺旋矩阵为: [ [1,2,3,4,5] [16,17,18,19,6] [15,24,25,20,7] [14,23,22,21,8] [...
BugkuCTF 部分题解(随缘更新)
volcano的博客
12-30 7552
之前做的题在 佛系更新 等假期抽空做 bugku佛系更新MISC简单取证1南城旧梦成果狗成果狗 MISC 简单取证1 下载得到windows系统下一个目录，获取用户名和密码需要用SAM和system两个文件。 把SAM和SYSTEM文件放到Win32文件夹下，运行mimikatz，执行命令 所以flag{administrator-QQAAzz_forensics} 南城旧梦 mmz.bmp文件尾有一段DE@@=<6J:DB625K4，rot47解码后得到stoolkeyisqeadzc 意思是使
[随缘一题]-LintCode-平面列表
呼延十
12-26 767
来源 来源: lintcode-平面列表 描述 给定一个列表，该列表中的每个要素要么是个列表，要么是整数。将其变成一个只包含整数的简单列表。 样例 给定 [1,2,[1,2]]，返回 [1,2,1,2]。 给定 [4,[3,[2,[1]]]]，返回 [4,3,2,1]。 挑战 请用非递归方法尝试解答这道题。 解题思路 这道题一看就是用递归解决啦~,好,那我们就用递归. 啥玩意你不让用??? 那...
使用Idea过程中遇到的问题总结（随缘更新）
I'm Hexesdesu
06-07 202
坑
一份latex的学习总结（随缘持续更新…）
weixin_46390192的博客
02-17 865
一份Latex的学习笔记，随缘记得，随缘看吧，或许以后用到了新的东西还会再往里面写 写的感觉有的常识性的东西写的不是很清楚，如果看不明白评论一下吧（或者直接GPT更快）
开发感悟记录（后续随缘更新）
Sun_Raiser的博客
03-13 366
团队开发感悟
图论模板，随缘不定期更新
Ninght
01-06 1157
图论模板，随缘不定期更新网络流最大流dinic（更新于2021/1/6）hlpp（更新于2021/1/6）最小费用最大流（更新于2021/1/6）无源汇有上下界可行流（更新于2021/1/6）有源汇上下界最大流（更新于2021/1/6）最短路径dijkstra（更新于2021/1/6）spfa（更新于2021/1/6）最小生成树kruskal（更新于2021/1/6）prime（更新于2021/1/6）欧拉通路与欧拉回路并查集（更新于2021/1/6）判断存在欧拉回路（更新于2021/1/6）寻找欧拉路径（
随缘网络留言板 v2.0 build 20120127.rar
07-05
留言板自带后台，可批量管理留言，适合各中小企业将其整合植入网站或者做为一个单独留言系统交流使用。 2012-1-27更新： 1、针对前一段时间网站明文密码安全问题对密码进行了md5加密。 2、针对原验证码在留言提交...
随缘网络留言板 2.0 build 20120128.rar
05-24
留言板自带后台，可批量管理留言，适合各中小企业将其整合植入网站或者做为一个单独留言系统交流使用。 随缘网络留言板 2.0 build 20120128 更新内容: 1、针对前一段时间网站明文密码安全问题对密码进行了md5...
随缘文章网 v1.0源代码
03-23
因此我把一个很不错的新闻系统整理了一下，做了这个文章站，大家直接下载就可以用，不用在添加文章了，里面已经自带了进1000篇文章。不过这个程序有一个小小的不足，那就是不能添加分类。如果适合你的话就下载用吧。...
随缘网络php企业网站管理系统 v2.0.rar
07-16
4、后台采用最新版xhEditor编辑器，xhEditor是一个基于jQuery开发的跨平台开源迷你XHTML编辑器组件。 5、后台图片栏目可以对相关的幻灯，友情链接，logo，banner等进行修改设置。 5、后台可批量删除产品新闻。
力扣70 爬楼梯 C语言 动态规划 递归
XinxingZh的博客
05-07 323
爬 0 层和爬 1 层都只有一种情况， 但是爬两层有两种：一次爬一层一共爬两次、一次爬两层一共爬一次，爬三层有三种：一次爬一层一共爬三次、先爬一层再爬两层一共爬两次、先爬两层再爬一层一共爬两次。所以 f(0) = 1, f(1) = 1, f(2) = 2, f(3) = 3, f(4) = 5。规律是 f(n) = f(n-1) + f(n-2)，因为爬到第 n 阶有两种情况，分别是站在第 n-1 阶爬一层和站在第 n-2 阶爬两层，所以就是 f(n-1) 和 f(n-2)的和。有两种方法可以爬到楼顶。
c语言从入门到函数速成（2）
2302_80418176的博客
05-04 861
温馨提醒：本篇文章适合人群：刚学c又感觉那个地方不怎么懂的同学以及以及学了一些因为自身原因停学一段时间后又继续学​​​c的学好，正片开始！
【C语言】动态内存管理
2401_82759674的博客
05-07 929
本篇博客将C语言中的动态内存函数整体梳理一遍，能帮助大家快速掌握，并将其中容易犯的错误一一列举出来，使用题目加强理解
Java 7大排序
最新发布
m0_74270127的博客
05-10 705
插入排序、希尔排序、选择排序、堆排序、冒泡排序、快速排序、归并排序
若我们对k最短路问题添加一个约束:在最终得到的包含最短路、次短路、……、第k短路的集合A中，对于任意两条路径P,P∈A，我们要求Sim(P,P))≤t。优化目标为集合A中所有路径的总长度最短。该问题可称为考虑路径多样性的k最短路问题。请从贪心算法（(Greedyalgorithm)、分治算法（Divide and Conquer）以及动态规划（Dynamic Programming)中选取·种，设计相应算法求解该问题，不要求得到最优解。
02-06
1. 建立一个最短路径的堆，初始时堆中只包含起点到终点的最短路径。 2. 取出堆顶元素（即当前最短路径），并把它加入最终的路径集合A中。 3. 对于堆中剩余的每一条路径，如果它与当前最短路径的相似度小于等于t，...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值