题目链接:Click here~~
题意:
RT。
解题思路:
http://www.cnblogs.com/jbelial/archive/2011/08/08/2131165.html
求重心时需要注意的地方:
1、总面积为0的时候,因为总面积会作除数。(本题没有出现)
2、面积要用有向面积。(考虑凹多边形的情况,如下图)
#include <math.h>
#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
#define N 1000005
struct Point
{
int x,y;
}P[N];
const double eps = 1e-6;
int sgn(double x)
{
return fabs(x)<eps ? 0 : (x > eps ? 1 : -1);
}
double Dis(const Point& p1,const Point& p2)
{
return sqrt( (p1.x-p2.x)*(p1.x-p2.x) + (p1.y-p2.y)*(p1.y-p2.y)*1.0);
}
double Cross(const Point& p1,const Point& p2,const Point& p3,const Point& p4)
{
return (p2.x-p1.x)*(p4.y-p3.y) - (p2.y-p1.y)*(p4.x-p3.x);
}
double Area(const Point& p1,const Point& p2,const Point& p3)
{
return Cross(p1,p2,p1,p3);
}
double CGx(const Point& p1,const Point& p2,const Point& p3)
{
return (p1.x+p2.x+p3.x)/3.0;
}
double CGy(const Point& p1,const Point& p2,const Point& p3)
{
return (p1.y+p2.y+p3.y)/3.0;
}
int main()
{
int z,n;
double ansx,ansy,sum;
scanf("%d",&z);
while(z--)
{
scanf("%d",&n);
for(int i=0;i<n;i++)
scanf("%d%d",&P[i].x,&P[i].y);
ansx = ansy = sum = 0;
for(int i=1;i<n-1;i++)
{
double area = Area(P[0],P[i],P[i+1]);
ansx += CGx(P[0],P[i],P[i+1]) * area;
ansy += CGy(P[0],P[i],P[i+1]) * area;
sum += area;
}
printf("%.2f %.2f\n",ansx/sum,ansy/sum);
}
return 0;
}
分享到:
相关推荐
HDU的1250,主要是利用高精度加法,但是代码有点繁琐,效率不是很高
杭电ACMhdu1163
HDU1059的代码
hdu1001解题报告
hdu 1574 passed sorce
HDU的一题........HDU DP动态规
The least common multiple (LCM) of a set of positive integers is the smallest positive integer which is divisible by all the numbers in the set. For example, the LCM of 5, 7 and 15 is 105. Input Input...
hdu2101AC代码
hdu acm 教案 搜索入门 hdu acm 教案 搜索入门
搜索 dfs 解题代码 hdu1241
hdu 5007 Post Robot 字符串枚举。 暴力一下就可以了。
hdu acm 教案 动态规划(1) hdu acm 教案 动态规划(1)
ACM HDU题目分类,我自己总结的大概只有十来个吧
hdu 1166线段树代码
HDU最全ac代码
自己做的HDU ACM已经AC的题目
hdu动态规划算法集锦
hdu题目分类
HDU图论题目分类
hdu-acm源代码(上百题)hdu-acm源代码、hdu-acm源代码hdu-acm源代码