HDU 5839 Special Tetrahedron (2016CCPC网络赛08) (暴力+剪枝)

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=5839

在一个三维坐标，给你n个点，问你有多少个四面体(4个点，6条边) 且满足至少四边相等 其余两边不相邻。

暴力4重循环，但是在第3重循环的时候需要判断是否是等腰三角形，这便是一个剪枝。在第4重循环的时候判断4点是否共面 (叉乘)， 5或者6边相等就+1，4边相等就判断另外两边是否相交就行了。

赛后过的，觉得自己还是太菜了。

 //#pragma comment(linker, "/STACK:102400000, 102400000")8
 #include <algorithm>
 #include <iostream>
 #include <cstdlib>
 #include <cstring>
 #include <cstdio>
 #include <vector>
 #include <cmath>
 #include <ctime>
 #include <list>
 #include <set>
 #include <map>
 using namespace std;
 typedef long long LL;
 typedef pair <int, int> P;
 const int N = 2e2 + ;
 int cost[N][N]; //两点的边长的平方
 struct edge {
     int x, y, z;
 }a[N*N];

 int main()
 {
     int t;
     scanf("%d", &t);
     for(int ca = ; ca <= t; ++ca) {
         int n;
         scanf("%d", &n);
         for(int i = ; i <= n; ++i) {
             scanf("%d %d %d", &a[i].x, &a[i].y, &a[i].z);
             for(int j = ; j < i; ++j) {
                 cost[i][j] = cost[j][i] = (a[i].x - a[j].x)*(a[i].x - a[j].x) + (a[i].y - a[j].y)*(a[i].y - a[j].y) + (a[i].z - a[j].z)*(a[i].z - a[j].z);
             }
         }
         edge s1, s2, s3;
         int ans, res = , x1, x2, y1, y2, len1, len2; //x1 y1是一个三角面不同边的两端 len1为等腰三角形的腰
         for(int i1 = ; i1 < n; ++i1) {
             for(int i2 = i1 + ; i2 < n; ++i2) {
                 for(int i3 = i2 + ; i3 < n; ++i3) {
                     //判断3点是否形成等腰三角形
                     if(cost[i1][i2] != cost[i1][i3] && cost[i2][i3] != cost[i1][i2] && cost[i2][i3] != cost[i1][i3])
                         continue;
                     else if(cost[i1][i2] == cost[i1][i3] && cost[i1][i2] == cost[i2][i3]) {
                         len1 = cost[i1][i2], x1 = y1 = ;
                     }
                     else if(cost[i1][i2] == cost[i2][i3]) {
                         len1 = cost[i1][i2], x1 = i1, y1 = i3;
                     }
                     else if(cost[i1][i2] == cost[i1][i3]) {
                         len1 = cost[i1][i2], x1 = i2, y1 = i3;
                     }
                     else {
                         len1 = cost[i2][i3], x1 = i1, y1 = i2;
                     }
                     for(int i4 = i3 + ; i4 <= n; ++i4) {
                         if(cost[i1][i4] != cost[i2][i4] && cost[i1][i4] != cost[i3][i4] && cost[i3][i4] != cost[i2][i4])
                             continue;
                         else if(cost[i1][i4] == cost[i2][i4] && cost[i1][i4] == cost[i3][i4]) {
                             len2 = cost[i1][i4], x2 = y2 = ;
                         }
                         else if(cost[i1][i4] == cost[i2][i4]) {
                             len2 = cost[i1][i4], x2 = i3, y2 = i4;
                         }
                         else if(cost[i1][i4] == cost[i3][i4]) {
                             len2 = cost[i1][i4], x2 = i2, y2 = i4;
                         }
                         else {
                             len2 = cost[i2][i4], x2 = i1, y2 = i4;
                         }
                         s1.x=a[i2].x-a[i1].x;s1.y=a[i2].y-a[i1].y;s1.z=a[i2].z-a[i1].z;
                         s2.x=a[i3].x-a[i1].x;s2.y=a[i3].y-a[i1].y;s2.z=a[i3].z-a[i1].z;
                         s3.x=a[i4].x-a[i1].x;s3.y=a[i4].y-a[i1].y;s3.z=a[i4].z-a[i1].z;
                         ans=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(ans == ) //4点是否共面
                             continue;
                         if(!x1 || !x2) { //6边或5边相等
                             if(len1 == len2)
                                 res++;
                         }
                         else {
                             if(len1 == len2 && (x1 != x2 && x1 != y2 && y1 != x2 && y1 != y2)) //四边相等
                                 res++;
                         }
                     }
                 }
             }
         }
         printf("Case #%d: %d\n", ca, res);
     }
     return ;
 }