【bzoj1913】 Apio2010—signaling 信号覆盖
2016-12-25 10:10
417 查看
http://www.lydsy.com/JudgeOnline/problem.php?id=1913 (题目链接)
题意
给出一个平面上n个点,求任选3个点画一个圆所包含的点的期望值。Solution
这个问题可以转化为凹凸多边形的问题求解(当然我是没想到的)。。左转题解:http://blog.csdn.net/regina8023/article/details/45556321细节
注意存放极角的数组要开成2倍。代码
// bzoj1913 #include<algorithm> #include<iostream> #include<cstdlib> #include<cstring> #include<cstdio> #include<cmath> #include<queue> #define LL long long #define inf 2147483640 #define eps 1e-6 #define Pi acos(-1.0) #define free(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout); using namespace std; const int maxn=2000; struct point {int x,y;}p[maxn]; LL C[maxn][maxn]; double a[maxn<<1]; int n; int main() { scanf("%d",&n); for (int i=1;i<=n;i++) scanf("%d%d",&p[i].x,&p[i].y); for (int i=0;i<=n;i++) C[i][0]=1; for (int i=1;i<=n;i++) for (int j=1;j<=i;j++) C[i][j]=C[i-1][j]+C[i-1][j-1]; LL t2,t1=0; for (int i=1;i<=n;i++) { int cnt=0; for (int j=1;j<=n;j++) if (i!=j) { a[++cnt]=atan2(p[j].y-p[i].y,p[j].x-p[i].x); if (a[cnt]<0) a[cnt]+=2*Pi; } sort(a+1,a+cnt+1); for (int j=1;j<=cnt;j++) a[j+cnt]=a[j]+2*Pi; for (int j=1,k=1;j<=cnt;j++) { while (a[k+1]-a[j]<Pi && k<2*cnt) k++; t1-=C[k-j][2]; //凸多边形个数 } t1+=C[n-1][3]; } t2=C [4]-t1; printf("%.6lf",(double)(t1+t2*2)/C [3]+3); return 0; }
相关文章推荐
- bzoj1913: [Apio2010]signaling 信号覆盖
- 【BZOJ1913】[Apio2010]signaling 信号覆盖
- BZOJ 1913: [Apio2010]signaling 信号覆盖
- BZOJ 1913: [Apio2010]signaling 信号覆盖
- BZOJ1913: [Apio2010]signaling 信号覆盖
- [BZOJ1913][Apio2010]signaling 信号覆盖
- [BZOJ 1913][Apio2010]signaling 信号覆盖:计算几何
- [BZOJ1913][Apio2010]signaling 信号覆盖(计算几何+组合数学)
- bzoj 1913: [Apio2010]signaling 信号覆盖
- 【BZOJ 1913】 [Apio2010]signaling 信号覆盖
- bzoj 1913: [Apio2010]signaling 信号覆盖【旋转卡壳(?)】
- bzoj1913[Apio2010]signaling 信号覆盖 计算几何
- [BZOJ1913][APIO2010]信号覆盖(计算几何+计数)
- 1913: [Apio2010]signaling 信号覆盖
- BZOJ 1913 [Apio]signaling 信号覆盖 方法简析
- [BZOJ 1913][APIO 2011]信号覆盖(计算几何)
- 【BZOJ-1913】signaling信号覆盖 极角排序 + 组合
- BZOJ 1913 signaling 信号覆盖
- bzoj 1913 signaling 信号覆盖 (极角排序 组合数)
- bzoj1913: [Apio2010]signaling 信号覆盖