1913: [Apio2010]signaling 信号覆盖

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]

思路很不错，就是想不到。。。
先任意取四个点组成的四边形，
然后在这个四边形任取三个点构成三角形，做这个三角形的外接圆，考虑能否覆盖第四个点
若这个四边形是凹四边形，显然只有取非凹点外的三个点，这样的方案能覆盖凹点
若这个四边形是凸四边形，能覆盖的方案有且仅有两种
证明的话，考虑这张图
若一个点在一个三角形的外接圆内，必有$\alpha + \beta > \pi$
显然，对于一个凸多边形，满足对角和大于$\pi$的有且仅有一对(题目已告知不存在四点共圆)
假设所有任取四个点构成的四边形中，有$A$个凸四边形，$B$个凹四边形
那么显然有$Ans = \frac {2A + B}{C_n^3} + 3$
凸四边形的个数不太好统计，那么就统计出凹四边形的数量，然后容斥一下
枚举凹四边形的凹点，记为$o$，每次对剩下的点极角排序
枚举个点$x$，确定四边形的一条边为$o\rightarrow x$
找到点$y$使得$x，y$的极角差值恰好不超过$pi$
则在$x \rightarrow 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;
}