基本算法如下，但不知main函数如何写
请各位高手速速给出解答，谢谢！
#include <iostream.h>
#include <math.h>
#include <iomanip.h>

struct point
{
double x;
double y;
};

//求多边形的重心算法
//说明:
//求多边形重心并不是简单的把求三角形的重心公式推广就行了
//我的算法是在平面上取一点(一般取原点, 这样可以减少很多计算, 而且使思路更清晰^_^)
//这样就得到了N个三角形OP[i]P[i+1](其中点的顺序要为逆时针的),
//分别求出这N个三角形的重心Ci和面积Ai(注意此处面积是又向面积, 就是用叉乘求面积时保留其正负号)
//在求出A = A1+A2+...+AN(同样保留正负号的代数相加)
//最终重心C = sigma(Ai+Ci)/A;
point gravity(point *p, int n)
{
double area = 0;
point center;
center.x = 0;
center.y = 0;

for (int i = 0; i < n-1; i++)
{
   area += (p[i].x*p[i+1].y - p[i+1].x*p[i].y)/2;
   center.x += (p[i].x*p[i+1].y - p[i+1].x*p[i].y) * (p[i].x + p[i+1].x);
   center.y += (p[i].x*p[i+1].y - p[i+1].x*p[i].y) * (p[i].y + p[i+1].y);
}

area += (p[n-1].x*p[0].y - p[0].x*p[n-1].y)/2;
center.x += (p[n-1].x*p[0].y - p[0].x*p[n-1].y) * (p[n-1].x + p[0].x);
center.y += (p[n-1].x*p[0].y - p[0].x*p[n-1].y) * (p[n-1].y + p[0].y);

center.x /= 6*area;
center.y /= 6*area;

return center;
}
   
各点如下：  
1    13.28    1.32
2    13.76    1.48
3    15.50    0.63
4    4.71    1.75
5    13.68    1.35
7    4.50    2.49
8    4.15    12.12
10    0.63    8.41
11    4.74    3.12
12    8.25    1.98
13    11.96    1.03

汗 子程序都写了 主程序写不成了？晕呢

我感觉这个算法似乎不对耶？御坂御坂挠着头

拜托啦，急求

有偿帮你，QQ:1247142440

你物理专业啊，果然牛逼，想法就是不一样~~

解释下，为什么求重心不能直接把各点相加求平均值啊？？？
主函数可以这样啊~
#include<iostream.h>
void main()
{
    const n=10;
    int g[n];
    int i;
    for(i=0;i<n;i++)
        cin>>g[n];
    point ans=gravity(g,n);
    cout<<"重心是:"<<'('<<ans.x<<','<<ans.y<<')'<<endl;
}

[ 本帖最后由 promising 于 2010-8-20 09:10 编辑 ]

你的main函数貌似有错误，根据你的我改了一下，能运行了，但输出有点问题

我改后的代码为
void main()
{
    const n=1;
    point g[n];
    int i;
    for(i=0;i<n;i++)
    {
        cin>>g[i].x;
        cin>>g[i].y;
    }
    point ans=gravity(g,n);
    cout<<"重心是:"<<'('<<ans.x<<','<<ans.y<<')'<<endl;
}

ok,好了，调试成功，6楼，感谢你的帮助

解释下：为什么求重心不能直接把各点坐标相加求平均值啊？？？从小到大我都是这么求的啊~~