nyoj 1099 Lan Xiang's Square（判断正方形问题）

给四个点判断是否能构成正方形。

正方形判定定理：

1：对角线相等的菱形是正方形
2：对角线互相垂直的矩形是正方形,正方形是一种特殊的矩形
3：四边相等,有一个角是直角的四边形是正方形
4：一组邻边相等的矩形是正方形
5：一组邻边相等且有一个角是直角的平行四边形是正方形
6：四边均相等,对角线互相垂直平分且相等的平面四边形

没考略到sqrt()损失精度问题wrong了几次，四点的排序问题wrong了几次（不是严格按顺序来排的）。哎。。。。

#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
struct point
{
    double x,y;
} a[4];
bool cmp(point a,point b)
{
    if(a.x!=b.x)
        return a.x>b.x;
    return a.y<b.y;
}
int main()
{
    int t,k;
    double s1,s2,s3,s4;
    cin>>t;
    while(t--)
    {
        for(int i=0; i<4; i++)
            cin>>a[i].x>>a[i].y;
        sort(a,a+4,cmp);
        s1=sqrt((a[0].x-a[2].x)*(a[0].x-a[2].x)+(a[0].y-a[2].y)*(a[0].y-a[2].y));
        s2=sqrt((a[0].x-a[1].x)*(a[0].x-a[1].x)+(a[0].y-a[1].y)*(a[0].y-a[1].y));
        s3=sqrt((a[3].x-a[1].x)*(a[3].x-a[1].x)+(a[3].y-a[1].y)*(a[3].y-a[1].y));
        s4=sqrt((a[2].x-a[3].x)*(a[2].x-a[3].x)+(a[2].y-a[3].y)*(a[2].y-a[3].y));
        k=(a[0].x-a[2].x)*(a[0].x-a[1].x)+(a[0].y-a[1].y)*(a[0].y-a[2].y);
        if(s1==s2&&s3==s4&&s1==s2&&s1!=0&&fabs(k)<0.00001)
            cout<<"Yes"<<endl;
        else
            cout<<"No"<<endl;
    }
}