ZOJ 3871 Convex Hull(计算几何、凸包)

题意:给n个点,|x[i]|,|y[i]| <= 1e9。求在所有情况下的子集下(子集点数>=3),凸包的面积和。

这题主要有几个方面,一个是凸包的面积,可以直接用线段的有向面积和求得,这个当时没想到。还有就是,在180度以内的点数。

所以实际上是枚举2个点作为某个凸包的一条边,看有多少个这样的凸包。那个2^num - 1是保证除了2个点外至少还需1个点才能构成凸包。

时间复杂度O(n*n*logn)

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <cmath>
using namespace std; const double pi = acos(-1.0);
#define ll long long
#define maxn 1010
#define mod 998244353
int area2(int ax, int ay, int bx, int by){
return ((1ll*ax*by%mod-1ll*ay*bx%mod)%mod+mod)%mod;
}
struct Pnt{
int x,y;
double ang;
bool operator < (const Pnt &b) const{
return ang < b.ang;
}
}vec[maxn<<1];
int p2[maxn];
int main(){
p2[0]=1;
for(int i=1;i<=1000;++i) p2[i] = (p2[i-1]<<1)%mod;
int t;
scanf("%d",&t);
while(t--){
int n;
scanf("%d",&n);
int x[maxn],y[maxn];
for(int i=0;i<n;++i){
scanf("%d%d",x+i,y+i);
}
int ans = 0;
for(int i=0;i<n;++i){
for(int j=0;j<n;++j){
vec[j].x = x[j];
vec[j].y = y[j];
vec[j].ang = atan2(x[j]-x[i],y[j]-y[i]);
}
vec[i] = vec[n-1];
for(int j=0;j<n-1;++j){
vec[j+n-1] = vec[j];
vec[j+n-1].ang += pi*2;
}
int nn = n-1+n-1;
sort(vec,vec+nn);
int l=0,r=0;
while(l<n-1){
while(r<nn && vec[r].ang - vec[l].ang < pi) r++;
int area = area2(x[i],y[i],vec[l].x,vec[l].y);
int cnt = (p2[r-l-1]-1+mod)%mod;
ans = (ans+1ll*area*cnt%mod)%mod;
++l;
}
}
printf("%d\n",(mod-ans)%mod);
}
return 0;
}
上一篇:20165213 java学习第一周


下一篇:java子类继承关系