FZU 2035 Axial symmetry（判断多边形是否对称）

题目链接：http://acm.fzu.edu.cn/problem.php?pid=2035

这个题目没发现平行x轴线和平行y轴条件的作用，直接按照随便的多边形思想来，没想到精度控高了，一直没能AC，后来把

1e-10改成1e-6就AC了！

解题思路，首先多边形的对称轴一定经过重心，所以先求出边形的重心

然后开始从第枚举对称点，从第一个枚举（为了防止枚举的正好是对称轴，特意枚举了第二个点）

枚举的时候直接判断两点的中点与重心的连线是否与这两点连成的直线垂直，再判断夹在这两点中间

的和两边的是不是对称，因为找到一对称点之后，夹在之间的就一定对称（还有两点之外的，程序中

枚举1的时候使用个小技巧，将最后位置设置为第一个点，这样可以减少很多代码量）！

#include <string.h>
#include <algorithm>
#include <stdio.h>
#include <cmath>
#include <iostream>
using namespace std;
#define maxn 600
#define eps 1e-6
bool zero(double a) {return (fabs(a)<eps?1:0);}
struct point{
    double x,y;
}po[maxn],center;
int n;
struct line{
    point a,b;
};
double cross1(point &a,point &b,point &c){
    return (b.x-a.x)*(c.y-a.y)-(c.x-a.x)*(b.y-a.y);
}
double cross(point p1,point p2,point p0){
    return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}
point intersection(line u,line v){
    point ret=u.a;
    double t=((u.a.x-v.a.x)*(v.a.y-v.b.y)-(u.a.y-v.a.y)*(v.a.x-v.b.x))/((u.a.x-u.b.x)*(v.a.y-v.b.y)-(u.a.y-u.b.y)*(v.a.x-v.b.x));
    ret.x+=(u.b.x-u.a.x)*t;
    ret.y+=(u.b.y-u.a.y)*t;
    return ret;
}
point barycenter(point a,point b,point c){
    line u,v;
    u.a.x=(a.x+b.x)/2;
    u.a.y=(a.y+b.y)/2;
    u.b=c;
    v.a.x=(a.x+c.x)/2;
    v.a.y=(a.y+c.y)/2;
    v.b=b;
    return intersection(u,v);
}
point barycenter(int n,point *p){
    point ret,t;
    double t1=0,t2;
    int i;
    ret.x=ret.y=0;
    for(i=1;i<n-1;i++){
        if(fabs(t2=cross(p[0],p[i],p[i+1]))>eps){
            t=barycenter(p[0],p[i],p[i+1]);
            ret.x+=t.x*t2;
            ret.y+=t.y*t2;
            t1+=t2;
        }
    }
    if(fabs(t1)>eps)
        ret.x/=t1,ret.y/=t1;
    return ret;
}
double dmult(point &a,point &b){
    return a.x*b.x+a.y*b.y;
}
double dmult(double a,double b,double c,double d){
    return a*c+b*d;
}
bool is_ok(int l,int r,int a,int b){
    if(r<l)return true;
    point temp,t;
    while(l<r){
        temp.x=(po[l].x+po[r].x)/2;
        temp.y=(po[l].y+po[r].y)/2;
         if(!zero(dmult(po[l].x-po[r].x,po[l].y-po[r].y,temp.x-center.x,temp.y-center.y)))
            return false;
         l++,r--;
    }
    temp.x=(po[a].x+po[b].x)/2;
    temp.y=(po[a].y+po[b].y)/2;
    if(l==r)
    return zero(cross(temp,center,po[l]));
    return true;
}
bool sovle(){
    int i,k;
    point temp;
    for(i=1;i<n;i++){
        temp.x=(po[0].x+po[i].x)/2;
        temp.y=(po[0].y+po[i].y)/2;
        if(zero(dmult(po[i].x-po[0].x,po[i].y-po[0].y,temp.x-center.x,temp.y-center.y))){
                if(is_ok(1,i-1,0,i)&& is_ok(i+1,n-1,0,i))
                return true;
        }
    }
     for(i=2;i<n;i++){
        temp.x=(po[1].x+po[i].x)/2;
        temp.y=(po[1].y+po[i].y)/2;
        if(zero(dmult(po[i].x-po[1].x,po[i].y-po[1].y,temp.x-center.x,temp.y-center.y))){
            //printf("HHHH%lf %lf\n",po[i].x,po[i].y);
                if(is_ok(2,i-1,1,i)&& is_ok(i+1,n,1,i))
                return true;
        }
    }
    return false;
}
int main(){
    int i,j,k,t,cas=0;
    scanf("%d",&t);
    while(t--){

        scanf("%d",&n);
        for(i=0;i<n;i++)
        scanf("%lf%lf",&po[i].x,&po[i].y);
        k=1;
        for(i=2;i<n;i++){
            if(zero(cross(po[k-1],po[k],po[i])))
            po[k]=po[i];
            else po[++k]=po[i];
        }
        k++;
        if(zero(cross(po[0],po[k-1],po[k-2])))
        k--;
        if(zero(cross(po[k-1],po[0],po[1]))){
                for(i=1;i<k;i++)
                 po[i-1]=po[i];
                 k--;
        }
        n=k;
        printf("Case %d: ",++cas);
        if(k<=2){
            printf("YES\n");
            continue;
        }
        po[n]=po[0];
        center=barycenter(n,po);
        if(sovle())
            printf("YES\n");
        else
            printf("NO\n");
    }
    return 0;
}