51Nod 四点共面

1265 四点共面 

基准时间限制：1 秒 空间限制：131072 KB 分值: 0  难度：基础题

 收藏

 关注

给出三维空间上的四个点（点与点的位置均不相同），判断这4个点是否在同一个平面内（4点共线也算共面）。如果共面，输出"Yes"，否则输出"No"。

Input

第1行：一个数T，表示输入的测试数量(1 <= T <= 1000)
第2 - 4T + 1行：每行4行表示一组数据，每行3个数，x, y, z, 表示该点的位置坐标(-1000 <= x, y, z <= 1000)。

Output

输出共T行，如果共面输出"Yes"，否则输出"No"。

Input示例

1
1 2 0
2 3 0
4 0 0
0 0 0

Output示例

Yes

四个点构造三个向量，形成一个行列式 
若行列式的值为0，则共面；不为0，则不共面 

这一段是转自https://blog.csdn.net/rcy_zhu/article/details/74930995

这是方法一：代码如下：

#include<iostream>
#include<algorithm>
using namespace std;
struct Point {
	double x;
	double y;
	double z;
};
int main()
{
	int t;
	Point a[4];
	cin>>t;
	while(t--)
	{
		for(int i=0;i<4;i++)
		cin>>a[i].x>>a[i].y>>a[i].z;
		Point s1, s2, s3;
		s1.x=a[1].x-a[0].x;
		s1.y=a[1].y-a[0].y;
		s1.z=a[1].z-a[0].z;
		s2.x=a[2].x-a[1].x;
		s2.y=a[2].y-a[1].y;
		s2.z=a[2].z-a[1].z;
		s3.x=a[3].x-a[2].x;
		s3.y=a[3].y-a[2].y;
		s3.z=a[3].z-a[2].z;
		int num=(s1.x*s2.y*s3.z+s1.y*s2.z*s3.x+s1.z*s2.x*s3.y-s1.z*s2.y*s3.x-s1.x*s2.z*s3.y-s1.y*s2.x*s3.z);
		if(num==0)
		printf("Yes\n");
		else printf("No\n");
	}	
}

方法二：

先解出三个点的平面方程，再把第四个点带入，就可以了；

要先学一下向量叉乘的概念。

来源于百度；

首先，我们利用三个点，得到AB，AC，然后利用叉乘，计算出平面ABC的法向量。叉乘的公式在上面。得到法向量后，就可以设出平面ABC的方程了。最后再把点带进去。手算过程如下：

因为 AB=（-3,-3,3）,AC=（0,-2,3）,
因此平面的法向量为 n=AB×AC=（-3,9,6）,
所以,由平面方程的点法向式可得,所求平面方程为 -3(x-1)+9(y-1)+6(z+1)=0 ,
化简得 x-3y-2z=0 .

代码如下：

#include<iostream>
#include<algorithm>
using namespace std;
struct Point {
	double x;
	double y;
	double z;
};
Point xiangliangchacheng(Point a, Point b)
{
	Point c;
	c.x=a.y*b.z-a.z*b.y;
	c.y=a.z*b.x-a.x*b.z;
	c.z=a.x*b.y-a.y*b.x;
	return c;
}
int main()
{
	int t;
	Point a[4];
	cin>>t;
	while(t--)
	{
		for(int i=0;i<4;i++)
		cin>>a[i].x>>a[i].y>>a[i].z;
		Point s1, s2;
		s1.x=a[1].x-a[0].x;
		s1.y=a[1].y-a[0].y;
		s1.z=a[1].z-a[0].z;
		s2.x=a[2].x-a[1].x;
		s2.y=a[2].y-a[1].y;
		s2.z=a[2].z-a[1].z;
		Point s3;
		s3=xiangliangchacheng(s1,s2);
		int num=s3.x*(a[3].x-a[0].x)+s3.y*(a[3].y-a[0].y)+s3.z*(a[3].z-a[0].z);
		if(num==0)
		printf("Yes\n");
		else printf("No\n");
	}	
}