CF 975E - 求质心+角度(数学)

解题思路：

质心：

一般的质心公式是：

X = (x1+x2+x3+..+xn)/n,Y = (y1+y2+y3+...+yn)/n,

但是这得是平均的多边形精度才是正确的.

求任意多边形质心步骤

1.将多边形划分为n-2个三角形根据上面的公式求得每个三角形的质心

2.根据叉积求每个三角形的面积

3.质心：X = (x1*m1+x2*m2+x3*m3+...+xn-2*mn-2)/(m1+m2+m3+...mn-2)

Y = (y1*m1+y2*m2+y3*m3+...+yn-2*mn-2)/(m1+m2+...+mn-2)，(mi表示第i个三角形的面积,xi,yi表示第i个三角形的质心坐标)

因为每次去掉一个钉子之后回复的平衡转态剩下那个钉子肯定就变成了最高点，那么根据物理学知识，质心受到的重力肯定与最高点的钉子的拉起是相反力，所以质心此时一定在最高点钉子的正下方。那么我们需做的就是预处理好每个顶点与质心的距离，其次就是求出两个顶点与质心的向量的夹角，这样就可以在知道最高点坐标的情况下知道任意顶点的坐标了。因为垂直关系我们可以根据夹角用三角函数求得顶点坐标。

 

#include<bits/stdc++.h>
#define pai 3.1415926536
using namespace std;
typedef long long ll;
const int mx = 2e5 + 10;
int n,m,ma=-1e9,pos;
double xi[mx],yi[mx],angle[mx],X,Y;
double mi,dis[mx];
double cross(int p1,int p2,int c)
{
    return (xi[p1]-xi[c])*(yi[p2]-yi[c]) - (xi[p2]-xi[c])*(yi[p1]-yi[c]);
}
int main()
{
    scanf("%d%d",&n,&m);
    int up = 0,a,b,c,down;
    for(int i=1;i<=n;i++) scanf("%lf%lf",xi+i,yi+i);
    double sum1 = 0,sum2 = 0,ret = 0;
    for(int i=2;i<n;i++){
        mi = fabs(cross(1,i,i+1));
        ret += mi;
        sum1 += (xi[1]+xi[i]+xi[i+1])/3*mi;
        sum2 += (yi[1]+yi[i]+yi[i+1])/3*mi;
    }
    X = sum1/ret;Y = sum2/ret;
    for(int i=1;i<=n;i++){
        angle[i] = atan2(yi[i]-Y,xi[i]-X);//求斜率
        if(angle[i]<0) angle[i] = 2*pai + angle[i];
    }
    for(int i=1;i<=n;i++) dis[i] = hypot(xi[i]-X,yi[i]-Y);//求直角三角形最长边
    while(m--){
        scanf("%d",&a);
        if(a==1){
            scanf("%d%d",&b,&c);
            swap(down,c);
            if(!up){
                up = b==1?2:1;
                X = xi[up],Y = yi[up] - dis[up];
            }else if(b==up){
                double v = angle[b]-angle[c];
                if(v<0) v = 2*pai+v;
                xi[c] = X + dis[c]*sin(v);
                yi[c] = Y + dis[c]*cos(v);
                X = xi[c],Y = yi[c] - dis[c];
                up = c;
            }
        }else{
            scanf("%d",&c);
            if(!up){
                printf("%.10f %.10f\n",xi[c],yi[c]);
                continue;
            }
            double v = angle[up]-angle[c];
            if(v<0) v = 2*pai+v;
            xi[c] = X + dis[c]*sin(v);
            yi[c] = Y + dis[c]*cos(v);
            printf("%.10f %.10f\n",xi[c],yi[c]);
        }
    }
	return 0;
}