数学一本通——K. 点与三角形关系(叉积求面积)

简介: 数学一本通——K. 点与三角形关系(叉积求面积)

Description

已知一个平面坐标系中三角形三个点A,B,C的坐标,判断另外一个点D是否在三角形内(点在三角形边上也认为在三角形内)


Input

输入共四行,每行两个数,前三行表示A,B,C的坐标,第四行为D的坐标。


Output

输出一个字符串,“in”表示点D在三角形ABC内,“out”表示点D在三角形ABC外。


Samples

Input 复制

1.0 2.0

7.0 1.0

7.0 5.0

5.0 3.0

Output

in


思路:

运用三角形的面积关系来判断点与三角形的关系。

假设点为P,三角形的三个顶点分别为A、B、C。

分别计算三角形ABC,ABP,ACP,BCP的面积:

如果:

S(ABC)=S(ABP)+S(ACP)+S(BCP)

则说明点在三角形内或在三角形的边上。

所以问题就转化成了已知三角形的三个点坐标求三角形的面积。

在网上百度了几种方法,都有精度误差:

double cul(node a,node b,node c){
    double x1=a.x,y1=a.y;
    double x2=b.x,y2=b.y;
    double x3=c.x,y3=c.y;
    return x1*y2-x1*y3+x2*y3-x2*y1+x3*y1-x2*y2;
}
double cul(node a,node b,node c)
{
    double x1=a.x,y1=a.y;
    double x2=b.x,y2=b.y;
    double x3=c.x,y3=c.y;
    double aa=sqrt( (x1-x2)*(x1-x2) + (y1-y2)*(y1-y2) );
    double bb=sqrt( (x3-x2)*(x3-x2) + (y3-y2)*(y3-y2) );
    double cc=sqrt( (x1-x3)*(x1-x3) + (y1-y3)*(y1-y3) );
    double pp=(aa+bb+cc)/2;
    double S=sqrt( pp*(pp-aa)*(pp-bb)*(pp-cc));
    return S;
}

然后张巨说他要是能用三角形的面积过了就“打死我”。

张巨的思路是用叉积来求三角形的面积。

根据高中知识,假设三角形的两个边长分别为a,b,两边的夹角为c,那么三角形的面积为a ∗ b ∗ s i n c ∗ 1 2 a*b*sinc*\frac{1}{2}a∗b∗sinc∗∗1/2

而叉积的定义是两个矢量模的乘积再乘夹角正弦,转为坐标计算就可以了。

代码:

#pragma GCC optimize(2)
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<ll,ll>PLL;
typedef pair<int,int>PII;
typedef pair<double,double>PDD;
#define I_int ll
inline ll read()
{
    ll x=0,f=1;
    char ch=getchar();
    while(ch<'0'||ch>'9')
    {
        if(ch=='-')f=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
    {
        x=x*10+ch-'0';
        ch=getchar();
    }
    return x*f;
}
#define read read()
#define closeSync ios::sync_with_stdio(0);cin.tie(0);cout.tie(0)
#define multiCase int T;cin>>T;for(int t=1;t<=T;t++)
#define rep(i,a,b) for(int i=(a);i<=(b);i++)
#define repp(i,a,b) for(int i=(a);i<(b);i++)
#define per(i,a,b) for(int i=(a);i>=(b);i--)
#define perr(i,a,b) for(int i=(a);i>(b);i--)
ll ksm(ll a,ll b,ll p)
{
    ll res=1;
    while(b)
    {
        if(b&1)res=res*a%p;
        a=a*a%p;
        b>>=1;
    }
    return res;
}
#define PI acos(-1)
const double eps=1e-10;
struct node
{
    double x,y;
    node operator - (node &s){
        return (node){x-s.x,y-s.y};
    }
};
double operator*(node a,node b){
    return a.x*b.y-b.x*a.y;
}
double cul(node a,node b,node c){
    return fabs((b-a)*(c-a)/2);
}
int main()
{
    node a,b,c,d;
    cin>>a.x>>a.y>>b.x>>b.y>>c.x>>c.y>>d.x>>d.y;
    double abc=cul(a,b,c);
    double abd=cul(a,b,d);
    double acd=cul(a,c,d);
    double bcd=cul(b,c,d);
    if(abd+acd+bcd==abc) puts("in");
    else puts("out");
    return 0;
}
目录
相关文章
|
1月前
|
Java C++ Python
试题 基础练习 圆的面积
试题 基础练习 圆的面积
29 1
|
1月前
|
Java C++ Python
试题 基础练习 矩形面积交
试题 基础练习 矩形面积交
30 0
|
8月前
|
Serverless C语言 C++
【数学建模】利用C语言来实现 太阳赤纬 太阳高度角 太阳方位角 计算和求解分析 树木树冠阴影面积与种植间距的编程计算分析研究
【数学建模】利用C语言来实现 太阳赤纬 太阳高度角 太阳方位角 计算和求解分析 树木树冠阴影面积与种植间距的编程计算分析研究
150 1
|
1月前
|
算法 测试技术 C#
【数学】【计算几何】1453. 圆形靶内的最大飞镖数量
【数学】【计算几何】1453. 圆形靶内的最大飞镖数量
|
7月前
|
算法 搜索推荐 程序员
C语言第四练——利用海伦公式编程求三角形的面积
C语言第四练——利用海伦公式编程求三角形的面积
141 0
|
11月前
|
机器学习/深度学习 Python
【每周一坑】输出三角形
如果输出固定长度对你来说太简单了,可以增加一个输入 n(n为正整数且 n>3),作为输出三角形第一行星号的数量。
|
Go 开发者
经典案例打印空心金字塔 | 学习笔记
简介:快速学习经典案例打印空心金字塔
119 0
经典案例打印空心金字塔 | 学习笔记
宇宙无敌加法器 (20 分)(数组模拟)
宇宙无敌加法器 (20 分)(数组模拟)
64 0
矩形面积交[蓝桥杯]
平面上有两个矩形,它们的边平行于直角坐标系的X轴或Y轴。对于每个矩形,我们给出它的一对相对顶点的坐标,请你编程算出两个矩形的交的面积。
矩形面积交[蓝桥杯]
|
C语言
浙大版《C语言程序设计(第3版)》题目集 - 习题9-3 平面向量加法(15 分)
浙大版《C语言程序设计(第3版)》题目集 - 习题9-3 平面向量加法(15 分)
111 0