計算幾何總結

來源:互聯網
上載者:User

設平面上有N 個離散資料點( xi , yi ) ( i = 1, 2, ., n) , 其
多邊形重心G( . x1, . y1) 為:



  成累加和:




#include <iostream>#include <algorithm>#include <cstring>#include <cstdio>using namespace std;struct Point{double x,y;};//叉積求面積double Area(const Point p0,const Point p1,const Point p2){    // 另外在求解的過程中,不需要考慮點的輸入順序是順時針還是逆時針,相除後就抵消了。     return (p0.x*p1.y+p1.x*p2.y+p2.x*p0.y-p1.x*p0.y-p2.x*p1.y-p0.x*p2.y)/2.0;}int main(){    int T,n;    scanf("%d",&T);    while(T--)    {        Point p0,p1,p2;        double area = 0,suma = 0,sumx = 0,sumy = 0;        scanf("%d",&n);        scanf("%lf%lf",&p0.x,&p0.y);        scanf("%lf%lf",&p1.x,&p1.y);        for(int i = 2;i < n;++i)        {            scanf("%lf%lf",&p2.x,&p2.y);            area = Area(p0,p1,p2);            suma += area;            sumx += (p0.x+p1.x+p2.x)*area;            sumy += (p0.y+p1.y+p2.y)*area;            p1 = p2;        }           printf("%.2lf %.2lf\n",sumx/suma/3,sumy/suma/3);    }    return 0;}








相關文章

聯繫我們

該頁面正文內容均來源於網絡整理,並不代表阿里雲官方的觀點,該頁面所提到的產品和服務也與阿里云無關,如果該頁面內容對您造成了困擾,歡迎寫郵件給我們,收到郵件我們將在5個工作日內處理。

如果您發現本社區中有涉嫌抄襲的內容,歡迎發送郵件至: info-contact@alibabacloud.com 進行舉報並提供相關證據,工作人員會在 5 個工作天內聯絡您,一經查實,本站將立刻刪除涉嫌侵權內容。

A Free Trial That Lets You Build Big!

Start building with 50+ products and up to 12 months usage for Elastic Compute Service

  • Sales Support

    1 on 1 presale consultation

  • After-Sales Support

    24/7 Technical Support 6 Free Tickets per Quarter Faster Response

  • Alibaba Cloud offers highly flexible support services tailored to meet your exact needs.