1. 程式人生 > >51Nod 1265 四點共面(計算幾何)

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

題解:A、B、C、D四點共面則三向量滿足(AB×AC)*AD=0(AB×AC結果是該平面的法向量);

向量a×向量b= 

| i j k |

|a1 b1 c1|

|a2 b2 c2| 

=(b1c2-b2c1,c1a2-a1c2,a1b2-a2b1) 

#include<iostream>
#include<string.h>
#include<algorithm>
using namespace std;
struct node{
        int x,y,z;
}p[4];//向量
int main()
{
        int T;
        scanf("%d",&T);
        int a,b,c,x,y,z;
        while(T--)
        {
                scanf("%d%d%d",&a,&b,&c);
                for(int i=1;i<=3;i++){
                        scanf("%d%d%d",&x,&y,&z);
                        p[i].x=x-a;
                        p[i].y=y-b;
                        p[i].z=z-c;
                }
                x=p[1].y*p[2].z-p[1].z*p[2].y;
                y=p[1].z*p[2].x-p[1].x*p[2].z;
                z=p[1].x*p[2].y-p[1].y*p[2].x;//printf("%d*%d*%d\n",x,y,z);
                if((p[3].x*x+p[3].y*y+p[3].z*z)==0)
                        printf("Yes\n");
                else
                        printf("No\n");
        }
        return 0;
}