BZOJ1913: [Apio2010]signaling 信号覆盖

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

计算几何

我们考虑用总贡献/总方案数 先把圆上的3个点提出来

考虑剩下点的贡献， 对于一个凹四边形 贡献只有1（凹进去的在圆内）

对于一个凸四边形 贡献为2（无四点共圆）

枚举一个点，极角排序， 枚举起点，用双指针卡出一个范围，组合数统计答案即可

P.S：本题无SPJ，输出6位小数

#include <bits/stdc++.h>

using namespace std;

const long long MAXN = 1510;
const double pi = acos( -1 );

struct Point
{
	double x, y;
}p[MAXN];

long long n, top;

long long ans, tot, ret;

double st[MAXN << 1];

inline long long C(long long x, long long y)
{
	if( y == 1 ) return x;
	if( y == 2 ) return x * ( x - 1 ) / 2;
	if( y == 3 ) return x * ( x - 1 ) * ( x - 2 ) / 6;
	if( y == 4 ) return x * ( x - 1 ) * ( x - 2 ) * ( x - 3 ) / 24;
}

inline long long cal(long long x)
{
	top = ret = 0;
	for( long long i = 1 ; i <= n ; i++ ) if( i ^ x ) st[ ++top ] = atan2( p[ i ].y - p[ x ].y, p[ i ].x - p[ x ].x );
	sort( st + 1, st + top + 1 );
	long long t = top;
	for( long long i = 1 ; i <= t ; i++  ) st[ ++top ] = st[ i ] + pi * 2; 
	long long p = 1;
	for( long long i = 1 ; i <= t ; i++ )
	{
		p = max( p, i + 1 );
		while( p <= top && st[ p ] < st[ i ] + pi ) p++;
		if( p - i - 1 >= 2 ) ret += C( p - i - 1, 2 );
	}
	return C( n - 1, 3 ) - ret;
}

int main()
{
	scanf( "%d", &n );
	for( long long i = 1 ; i <= n ; i++ ) scanf( "%lf%lf", &p[ i ].x, &p[ i ].y );
	if( n <= 3 ) return puts( "0" ), 0;
	tot = C( n, 3 );
	for( long long i = 1 ; i <= n ; i++ )
		ans += cal( i );
	long long ans2 = C( n, 4 ) - ans;
	printf( "%.6lf\n", ( double )( ans + ans2 * 2 ) / tot + 3.0 );
}