加权法 多边形重心问题代码 求解释

听自己心跳的声音 2014-04-27 09:04:29
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;

#define INF 0.0000001
const int N = 10010;

struct point
{
double x, y;
point():x(0),y(0){}
}p[N];

int main()
{
int ncase;
scanf("%d", &ncase);
while(ncase--)
{
int n;
double result = 0;
point ans;
scanf("%d", &n);
for(int i = 0; i < n; ++i)
scanf("%lf%lf", &p[i].x, &p[i].y);
for(int i = 1; i <= n; ++i)
{
double temp = (p[i % n].x * p[i -1].y - p[i % n].y * p[i -1].x) / 2.0;
result += temp;
ans.x += temp * (p[i % n].x + p[i - 1].x) / 3.0;
ans.y += temp * (p[i % n].y + p[i - 1].y) / 3.0;
}
if(fabs(result - 0) < INF)
puts("0.000 0.000");
else
printf("%.3lf %.3lf\n", fabs(result), (ans.x + ans.y) / result); //result取绝对值
}
return 0;
}

其中这两行代码怎么理解?
ans.x += temp * (p[i % n].x + p[i - 1].x) / 3.0;
ans.y += temp * (p[i % n].y + p[i - 1].y) / 3.0;
重心与面积的权值,第一,可上边x或y相加的时候为什么只加了两个,是把另一个当做0了吗?temp所对应的面积与三角形重心是怎么对应的?比如有5条边,对应顶点分别为x0,y0 x1,y1 x2,y2 x3,y3 x4,y4 ,每一个叉积除以2所求的面积分别是三个顶点所在的三角形,可是这样对应下去加起来的面积貌似感觉就不是多边形的面积,而是所有的相邻三个顶点构成三角形的面积之和。
...全文
74 1 打赏 收藏 转发到动态 举报
写回复
用AI写文章
1 条回复
切换为时间正序
请发表友善的回复…
发表回复
zlm0913 2014-04-27
  • 打赏
  • 举报
回复
http://www.cnblogs.com/jbelial/archive/2011/08/08/2131165.html

64,644

社区成员

发帖
与我相关
我的任务
社区描述
C++ 语言相关问题讨论,技术干货分享,前沿动态等
c++ 技术论坛(原bbs)
社区管理员
  • C++ 语言社区
  • encoderlee
  • paschen
加入社区
  • 近7日
  • 近30日
  • 至今
社区公告
  1. 请不要发布与C++技术无关的贴子
  2. 请不要发布与技术无关的招聘、广告的帖子
  3. 请尽可能的描述清楚你的问题,如果涉及到代码请尽可能的格式化一下

试试用AI创作助手写篇文章吧