E.New Year and Castle Construction
题意
给定n个点,对于每个点(p),求出4-point 子集(该子集有四个点,并且围成的圈包含(p))的个数
数据给的点中没有三点共线的情况
这个题的题意并不是很好理解,补题过后还发现自己理解的题解是有问题的。
分析
-
四个点围成一个圈,第5个点在其中。所以总方案数为(n*C_{n-1}^4),枚举中间那个点,然后从剩余的(n-1)个点中选剩余的四个点。
-
考虑(n^2)的做法,正向枚举所有方案中可行的方案并不是很合适,因为要考虑更大的子集的情况,所以反向枚举。
-
不合法的方案:其他的四个点并不能包含点(p)。枚举这四个点中的其中一个点(q),计算取剩余3个点的方案数。考虑(p)与(q)连成的直线,可以从这条直线某一侧取出三个点,可以发现这样的组合一定不会使得四个点构成的圈包含(p)。
-
如何不重不漏的删除所有不合法的方案?对于每个(p),(q)作为四个点中,与(p)连线斜率最小的那个点,然后在这样的基础上,只能从(pq)连线一侧拿另外三个点。可以想到这样的方案对于(p)来讲是没有重复的。而(p)的枚举又是独立的,所以可以补充不漏的删除所有方案。
const int N = 3005;
long double x[N],y[N];
int main()
{
int n;scanf("%d",&n);
for(int i=1;i<=n;i++){
cin >> x[i] >> y[i];
}
ll res = 1ll * n * (n-1) * (n-2) * (n-3) * (n-4) / 24;
long double pi = acos(-1.0L);
for(int i=1;i<=n;i++){
vector<long double> v;
for(int j=1;j<=n;j++){
if(i == j)continue;
v.push_back(atan2(y[j]-y[i],x[j]-x[i]));
}
sort(v.begin(),v.end());
int m = n-1,index = 0;
for(int j=0;j<m;j++){
while(index < j + m){
long double ang = v[index%m] - v[j];
if(ang < 0) ang += 2 * pi;
if(ang < pi)index ++;
else break;
}
ll cnt = index - j - 1;
res -= 1ll * cnt * (cnt - 1) * (cnt - 2) / 6;
}
}
cout << res << endl;
return 0;
}
参考链接
B站up主:https://www.bilibili.com/video/av82161298?p=411
ps:题解有些许啰嗦,up主视频中讲解的很清楚,然后我对这个解法的正确性进行了一些额外的思考,表达能力欠佳,如果读者发现描述中有不准确或者不理解的地方,烦请在下面留言,谢谢!在最后非常感谢B站up主qscqesze的讲解。