汕头市队赛 SRM 06 C 秀恩爱
C 秀恩爱 SRM 06
背景&&描述
KPM坐在直升机上俯瞰小渔村景象。
渔村可看作二维平面,密密麻麻地到处都是单身狗,KPM当前所在坐标为(sx,sy)。
KPM的后宫团们自发地聚集在一起为他送行,从空中看,后宫团形成了一个多边形。
当然了KPM是不在那个多边形内的。
直升机突然开始原地转圈,后宫团们因为想看着KPM的正脸,所以也跟着以KPM所在坐标为中心旋转。
后宫团所经之处单身狗尸横遍野。赶来救治伤员的医护人员想知道,多边形扫过的面积是多少。
注意,本题不保证横坐标互不相同、纵坐标互不相同什么的。
请注意运算过程中可能出现的/0等问题。
输入格式
第一行三个整数,n,sx,sy。n表示多边形的顶点数。
接下来n行每行俩整数,分别表示多边形一个顶点的横纵坐标。
(顶点是按照顺时针或者逆时针顺序给出的,并且所有点的坐标绝对值<=,保证不存在共线的三个顶点)
输出格式
一个整数,表示面积四舍五入为整数的结果。
样例输入
3 0 0 0 1 -1 2 1 2
样例输出
13
数据范围与约定
- 对于100%的数据:
样例解释
这道题主要是求最小半径以及最大半径
易证明最大半径mx一定在顶点位置
而最小半径也可能在两点的直线上
明白这一点后代码就很好实现了 不过要注意精度问题 今天似乎大爷们都被精度卡了 2333
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
#include<cstdio> #include<cstring> #include<algorithm> #include<cmath> #define LL long long using namespace std; const int M=155; const double P=3.141592653589793238462643383279502884; LL read(){ LL ans=0,f=1,c=getchar(); while(c<'0'||c>'9'){if(c=='-') f=-1; c=getchar();} while(c>='0'&&c<='9'){ans=ans*10+(c-'0'); c=getchar();} return ans*f; } double mx,mn=1e15; double k1,k2,b,L,mxx,mnx,nx; LL x[M],y[M],n,sx,sy; int main() { n=read(); sx=read(); sy=read(); for(int i=1;i<=n;i++){ x[i]=read()-sx; y[i]=read()-sy; L=1.0*(fabs(x[i])*fabs(x[i])+fabs(y[i])*fabs(y[i])); mx=max(mx,L); mn=min(mn,L); } for(int i=1;i<=n;i++){ int j=i+1<=n?i+1:1; if(x[i]==x[j]){ if((y[i]>0&&y[j]<0)||(y[i]<0&&y[j]>0)) mn=min(mn,(double)x[i]*x[j]); continue; } mxx=max(x[i],x[j]); mnx=min(x[i],x[j]); if(y[i]==y[j]){ if(mxx>0&&mnx<0) mn=min(mn,(double)y[i]*y[j]); continue; } k1=(double)(y[i]-y[j])/(x[i]-x[j]); b=(double) y[i]-k1*x[i]; k2=-1.0/k1; nx=b/(k2-k1); if(nx>mxx||nx<mnx) continue; L=b/sqrt(k1*k1+1); mn=min(mn,L*L); } printf("%lld",(long long)(P*(mx-mn)+0.5)); return 0; }