标签:cos 1.0 int sum 解题 Incenters ca AGC039D sin
AGC039D Incenters 解题报告:
题意
给定单位圆上 \(n\) 个点,随机选三个点求三角形内心期望。
\(3\leqslant n\leqslant 3000\)。
分析
不会推式子/dk。
令 \(C=\frac{n(n-1)(n-2)}{6}\) 为总方案数。
\[\text{E}(P_x)\times C=\sum_{i=1}^n(((n+2i)\sum_{j=i+1}^n\cos a_j-\sum_{j=i+1}^n 2j\cos a_j)\cos a_i-((n+2i)\sum_{j=i+1}^n \sin a_j-\sum_{j=i+1}^n 2j\sin a_j)\sin a_i) \]同理可以推出:
\[\text{E}(P_y)\times C=\sum_{i=1}^n(((n+2i)\sum_{j=i+1}^n\cos a_j-\sum_{j=i+1}^n 2j\cos a_j)\sin a_i-((n+2i)\sum_{j=i+1}^n \sin a_j-\sum_{j=i+1}^n 2j\sin a_j)\cos a_i) \]时间复杂度 \(O(n)\)。
代码
#include<stdio.h>
#include<math.h>
const int maxn=3005;
const double pi=acos(-1.0);
int n,L;
int a[maxn];
double ansx,ansy,C,sc,sci,ss,ssi;
int main(){
scanf("%d%d",&n,&L);
for(int i=1;i<=n;i++)
scanf("%d",&a[i]);
for(int i=n;i>=1;i--){
double A=1.0*a[i]*pi/L,ca=cos(A),sa=sin(A),vc=(1.0*(n+i+i)*sc-2*sci),vs=(1.0*(n+i+i)*ss-2*ssi);
ansx+=vc*ca-vs*sa,ansy+=vc*sa+vs*ca,sc+=ca,sci+=ca*i,ss+=sa,ssi+=sa*i;
}
double C=1.0*n*(n-1)*(n-2)/6;
printf("%.15lf %.15lf\n",ansx/C,ansy/C);
return 0;
}
标签:cos,1.0,int,sum,解题,Incenters,ca,AGC039D,sin 来源: https://www.cnblogs.com/xiaoziyao/p/15940149.html
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。