[Swust OJ 249]--凸包面积

area+=(pos.x*tmp.y-pos.y*tmp.x)/2.0;(pos当前点，tmp下一个点，最后一个点和第一个点在来一次)

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<stack>
#include<cmath>
usingnamespacestd;
structnode{
doublex,y;
}point[110],pos,tmp;
doubledis(nodea,nodeb){
returnpow((a.x-b.y),2)+pow((a.y-b.y),2);
}
//按点集分布排序,利用向量平行关系判断(极角排序)，角度相同则距离小的在前面
boolcmp(nodea,nodeb){
doublepovit=(a.x-pos.x)*(b.y-pos.y)-(b.x-pos.x)*(a.y-pos.y);
if(povit>0||!povit&&(dis(a,pos)<dis(b,pos)))
returntrue;
returnfalse;
}
//当前点是否在点集左侧,利用叉乘比较3个点两条线的斜率关系
boolturn_left(nodep1,nodep2,nodep3){
return(p2.x*p1.y+p3.x*p2.y+p1.x*p3.y-p3.x*p1.y-p1.x*p2.y-p2.x*p3.y)>0?true:false;
}
intmain(){
inti,sign,n,t;
doublearea;
cin>>t;
while(t--){
cin>>n;
for(i=0;i<n;i++)
cin>>point[i].x>>point[i].y;
if(n<=2){
cout<<"0.0\n";
continue;
}
stack<node>Q;
sign=0;
pos=point[0];
for(i=1;i<n;i++){
if(pos.y==point[i].y&&pos.x>point[i].x||point[i].y<pos.y){
pos=point[i];
sign=i;
}
}
swap(point[0],point[sign]);
sort(point+1,point+n,cmp);
Q.push(point[0]),Q.push(point[1]),Q.push(point[2]);
for(i=3;i<n;i++){
while(!Q.empty()){
tmp=Q.top();
Q.pop();
if(turn_left(tmp,Q.top(),point[i])){
Q.push(tmp);
break;
}
}
Q.push(point[i]);
}
area=0;
tmp=Q.top(),Q.pop();
area+=(pos.x*tmp.y-pos.y*tmp.x)/2.0;
while(!Q.empty()){
area+=(tmp.x*Q.top().y-tmp.y*Q.top().x)/2.0;
tmp=Q.top();
Q.pop();
}
printf("%.1lf\n",fabs(area));
}
return0;
}