传送门
解析:
读优没有写负数又被卡了半个小时。。。
这里采用
J
a
r
r
i
s
Jarris
Jarris步进法求凸包。。主要讲一讲怎么求多边形面积。
思路:
满足题意的显然是这些点的凸包,而我们要做的就是求出凸包面积。
那么怎么求多边形面积?
考虑三角剖分,我们将多边形由一个顶点发射出的所有对角线将该多边形剖分成若干个三角形。我们求每一个三角形的面积就是求多边形面积了。
好的问题就变成了怎么求三角形的面积。。。
我们还是不求三角形面积,改为求平行四边形面积。
求平行四边形面积就直接向量叉积就好了啊,证明什么的我才不会写
要证明的在评论区留言,我就抽空补了。
代码:
#include<cstdio>
#include<algorithm>
#include<cctype>
#include<cmath>
#include<iostream>
using namespace std;
#define ll long long
#define re register
#define gc getchar
#define pc putchar
#define cs const
inline
int getint(){
re int num;
re char c;
re bool f=0;
while(!isdigit(c=gc()))f^=c=='-';num=c^48;
while(isdigit(c=gc()))num=(num<<1)+(num<<3)+(c^48);
return f?-num:num;
}
cs int N=10002;
struct Point{
int x,y;
Point(int _x=0,int _y=0):x(_x),y(_y){}
Point operator+(cs Point &b)cs{return Point(x+b.x,y+b.y);}
Point operator-(cs Point &b)cs{return Point(x-b.x,y-b.y);}
Point operator*(cs int &b)cs{return Point(x*b,y*b);}
int operator*(cs Point &b)cs{return x*b.y-y*b.x;}
int dot()cs{return x*x+y*y;}
friend int dot(cs Point &a,cs Point &b){return a.x*b.x+a.y*b.y;}
friend int cross(cs Point &a,cs Point &b){return a.x*b.y-a.y*b.x;}
}O,q[N];
inline
int dist(cs Point &a){
return sqrt(a.dot());
}
inline
int dist(cs Point &a,cs Point &b){
return sqrt((a-b).dot());
}
inline
bool cmp(cs Point &a,cs Point &b){
return a.x==b.x?a.y<b.y:a.x<b.x;
}
struct Polygon{
Point p[N];
int n;
Polygon convex_hull(Point *q,int m){
sort(q+1,q+m+1,cmp);
n=0;
for(int re i=1;i<=m;++i){
while(n>=2&&cross(q[i]-p[n-1],p[n]-p[n-1])>=0)--n;
p[++n]=q[i];
}
int k=n;
for(int re i=m;i>=1;--i){
while(n>k&&cross(q[i]-p[n-1],p[n]-p[n-1])>=0)--n;
p[++n]=q[i];
}
if(m>1)--n;
return *this;
}
int area(){
int res=0;
for(int re i=2;i<=n;++i)res+=(p[i]-p[1])*(p[i-1]-p[1]);
return abs(res)/2;
}
}poly;
int m;
signed main(){
m=getint();
for(int re i=1;i<=m;++i)q[i].x=getint(),q[i].y=getint();
poly.convex_hull(q,m);
cout<<poly.area()/50;
return 0;
}