#DTOJ 5134 小h的几何

我是小 h，我是世界级几何大师。

我希望你了解一些定义。我可以证明三角形的三条高交于一点，我把这点称作垂心；我还可以证明三角形三边的中点、三条高的垂足、三角形三个顶点与垂心连线的中点、共九个点在一个圆上。我把这个圆叫九点圆，并记 $△ABC$ 的九点圆圆心为 $\Omega (A,B,C)$。

现在我在单位圆（即 { ( x , y ) ∣ x ∈ R , y ∈ R , x 2 + y 2 = 1 } \{(x, y)|x \in \mathbb R, y \in \mathbb R, x^2 + y^2 = 1 \} {(x,y)∣x∈R,y∈R,x2+y2=1}）上放置 $n$ 个两两不同的点 $P_1,P_2,\cdots ,P_n$，其中 $P_i(\cos {\theta }_i,\sin {\theta }_i)$。我还在这 $n$ 个点中等概率随机三个不同点，你能求出它们九点圆圆心的期望吗？

形式化的说，你需要输出$\frac{6}{n(n-1)(n-2)}\sum _{1\le i<j<k\le n}\Omega (P_i,P_j,P_k)$。

这里我定义两个点 $(a,b),(c,d)$ 的和为 $(a+c,b+d)$。

数据范围

对所有数据，保证 $3\le n\le 5\times 10^5,1\le q_i\le 2\times 10^9,q_i\ne q_j(\forall i\ne j)$。

子任务编号	$n\le$	特殊性质	分值
$1$	$20$		$20$
$2$	$200$		$30$
$3$	$2\times 10^3$		$20$
$4$	$5\times 10^5$	${\theta }_k=\frac{2k\pi }{n}$	$20$
$5$	$5\times 10^5$		$10$

题解

一道结论题：显然，其圆心在垂心和外心的中点上。
所以再猜个结论：圆心是其到三个顶点的向量之和。
于是这道题就写完了；

代码

#include<bits/stdc++.h>
#define ld long double
using namespace std;
const ld Pi=acos(-1);
const int N=5e5+10;
int n,q[N];
ld a[N],X=0,Y=0;
void puts(ld x,ld y){printf("%.15Lf %.15Lf\n",x,y);}
int main(){
//	freopen("geometry.in","r",stdin);
//	freopen("geometry.out","w",stdout);
//	ios::sync_with_stdio(false);
	cin>>n;
	for(int i=1;i<=n;++i)cin>>q[i],a[i]=(q[i]*Pi)*1.0/1e9,X+=cos(a[i])*3.0/n,Y+=sin(a[i])*3.0/n;
	puts(X/2,Y/2);
	return 0;
}