1913: [Apio2010]signaling 信号覆盖
Time Limit: 20 Sec Memory Limit: 64 MB
Submit: 1420 Solved: 573
[Submit][Status][Discuss]
Description
Input
输入第一行包含一个正整数 n, 表示房子的总数。接下来有 n 行,分别表示 每一个房子的位置。对于 i = 1, 2, .., n, 第i 个房子的坐标用一对整数 xi和yi来表 示,中间用空格隔开。
Output
输出文件包含一个实数,表示平均有多少个房子被信号所覆盖,需保证输出 结果与精确值的绝对误差不超过0.01。
Sample Input
4
0 2
4 4
0 0
2 0
Sample Output
3.500
HINT
3.5, 3.50, 3.500, … 中的任何一个输出均为正确。此外,3.49, 3.51,
3.499999,…等也都是可被接受的输出。
【数据范围】
100%的数据保证,对于 i = 1, 2, .., n, 第 i 个房子的坐标(xi, yi)为整数且
–1,000,000 ≤ xi, yi ≤ 1,000,000. 任何三个房子不在同一条直线上,任何四个房子不
在同一个圆上;
40%的数据,n ≤ 100;
70%的数据,n ≤ 500;
100%的数据,3 ≤ n ≤ 1,500。
Source
[Submit][Status][Discuss]
思路很不错,就是想不到。。。
先任意取四个点组成的四边形,
然后在这个四边形任取三个点构成三角形,做这个三角形的外接圆,考虑能否覆盖第四个点
若这个四边形是凹四边形,显然只有取非凹点外的三个点,这样的方案能覆盖凹点
若这个四边形是凸四边形,能覆盖的方案有且仅有两种
证明的话,考虑这张图
若一个点在一个三角形的外接圆内,必有 α+β>π
显然,对于一个凸多边形,满足对角和大于 π 的有且仅有一对(题目已告知不存在四点共圆)
假设所有任取四个点构成的四边形中,有 A 个凸四边形,B 个凹四边形
那么显然有 Ans=2A+BC3n+3
凸四边形的个数不太好统计,那么就统计出凹四边形的数量,然后容斥一下
枚举凹四边形的凹点,记为 o ,每次对剩下的点极角排序
枚举个点x ,确定四边形的一条边为 o→x
找到点 y 使得x,y 的极角差值恰好不超过 pi
则在 x→y 中的点(不包括 x )里任取两个,再与o,x 一起
这样构成的一定不是以 o <script type="math/tex" id="MathJax-Element-1399">o</script>点为凹点的凹四边形
容斥一下就知道凹四边形的数量了
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
const int maxn = 1505;
typedef double DB;
const DB PI = acos(-1.00);
const DB PI2 = PI * 2.00;
struct Point{
int x,y; Point(){}
Point(int x,int y): x(x),y(y){}
}p[maxn];
int n;
DB A,B,D[maxn*2];
inline DB C(int N,int M)
{
if (M > N) return 0; DB ret = 1;
for (int i = N; i > N - M; i--) ret *= (DB)(i);
for (int i = 1; i <= M; i++) ret /= (DB)(i); return ret;
}
inline void Calc(int k)
{
int tot = 0,tail = 1; DB sum = 0;
for (int i = 1; i < k; i++) D[++tot] = atan2(p[i].y - p[k].y,p[i].x - p[k].x);
for (int i = k + 1; i <= n; i++) D[++tot] = atan2(p[i].y - p[k].y,p[i].x - p[k].x);
sort(D + 1,D + tot + 1); D[2 * tot + 1] = 3.00001 * PI;
for (int i = 1; i <= tot; i++) D[i + tot] = D[i] + PI2;
for (int i = 1; i <= tot; i++)
{
while (D[tail + 1] - D[i] < PI) ++tail;
sum += C(tail - i,2);
}
B += C(tot,3) - sum;
}
int main()
{
#ifdef DMC
freopen("DMC.txt","r",stdin);
freopen("1.out","w",stdout);
#endif
cin >> n;
if (n == 3) {puts("3.000000"); return 0;}
for (int i = 1; i <= n; i++)
{
int x,y; scanf("%d%d",&x,&y);
p[i] = Point(x,y);
}
for (int i = 1; i <= n; i++) Calc(i);
A = C(n,4) - B; //cout << A << endl << B << endl;
printf("%.6lf\n",(2.00 * A + B) / C(n,3) + 3.00);
return 0;
}