51nod 1265 四点共面——计算几何

题目链接：http://www.51nod.com/Challenge/Problem.html#!#problemId=1265

以其中某一点向其它三点连向量，若四点共面，这三个向量定义的平行六面体的体积为零。

而这个体积等于这三个向量的混合积，所以只要验证混合积是否为零。

　　　　　　　　　　　　|ax ay az|

而混合积(a,b,c)=(a×b)·c=   |bx by bz|

　　　　　　　　　　　　|cx cy cz|

证明也不难：https://baike.baidu.com/item/%E6%B7%B7%E5%90%88%E7%A7%AF/10564182?fr=kg_general

所以只需验证ax*by*cz+ay*bz*cx+bx*cy*az-az*by*cx-ay*bx*cz-ax*bz*cy的值即可。

再提一句：空间向量a、b的叉积：a×b=(ax,ay,az)×(bx,by,bz)=(ay*bz-az*by , az*bx-ax*bz , ax*by-ay*bx)

 #include<cstdio>
 using namespace std;
 struct Point{ int x,y,z; }p[];
 int main(){
     int t;
     int x,y,z;
     int ax,ay,az,bx,by,bz,cx,cy,cz;
     scanf("%d", &t);
     for(int k = ; k <= t; ++k){
         scanf("%d%d%d",&x,&y,&z);//第一个点坐标
         for(int i = ; i < ; ++i)//后三个点坐标
             scanf("%d%d%d",&p[i].x,&p[i].y,&p[i].z);
         ///构建向量
         ax=p[].x-x; ay=p[].y-y; az=p[].z-z;
         bx=p[].x-x; by=p[].y-y; bz=p[].z-z;
         cx=p[].x-x; cy=p[].y-y; cz=p[].z-z;
         ///混合积
         if(ax*by*cz+ay*bz*cx+bx*cy*az-az*by*cx-ay*bx*cz-ax*bz*cy == )
             puts("Yes");
         else puts("No");
     }
     return ;
 }