NYOJ3 3 多边形重心问题

#include<stdio.h>
#include<math.h>
#include<string.h>
#include<stdlib.h>
struct Point{
    double x;
    double y;
}point[10010];
int i,n,m;
int main()
{
    scanf("%d",&n);
    while(n--){
        scanf("%d",&m);
        double tmp,area=0,x=0,y=0;
        for(i=0;i<m;i++)
            scanf("%lf%lf",&point[i].x,&point[i].y);
        for(i=1;i<=m;i++){
            tmp=(point[i%m].x*point[i-1].y-point[i-1].x*point[i%m].y)/2.0;
            area+=tmp;
            x+=tmp*(point[i%m].x+point[i-1].x)/3.0;
            y+=tmp*(point[i%m].y+point[i-1].y)/3.0;
        }
        if(fabs(area-0)<0.0000001)
            puts("0.000 0.000");  
        else
            printf("%.3lf %.3lf\n",fabs(area),(x+y)/area);
    }
    //system("pause");
    return 0;
}

精度要求很高，把变量定义成float都不行，贡献了好几个WA！