BZOJ 1670: [Usaco2006 Oct]Building the Moat护城河的挖掘

Description

求凸包周长.

Sol

凸包+计算几何.

这好像叫什么 Graham Scan 算法...

这个可以求凸包的周长,直径,面积.

选择一个基点,然后按极角排序,最后用一个栈一直维护方向单调.

极角排序就是先按与基点的向量和 \(x\) 轴的夹角排序,就是点积变一变.

维护方向的时候就是用叉积判断顺逆关系...

Code

/**************************************************************
    Problem: 1670
    User: BeiYu
    Language: C++
    Result: Accepted
    Time:36 ms
    Memory:1352 kb
****************************************************************/
 
#include
#include
#include
#include
#include
using namespace std;
 
#define mpr make_pair
#define sqr(x) ((x)*(x))
#define x first
#define y second
typedef pair< int,int > pr;
const int N = 5005;
int n,b;
pr g[N];
int stk[N],top;
 
inline int in(int x=0,char ch=getchar()){ while(ch>'9' || ch<'0') ch=getchar();
    while(ch>='0' && ch<='9') x=(x<<3)+(x<<1)+ch-'0',ch=getchar();return x; }
pr operator - (const pr &a,const pr &b){ return mpr(a.x-b.x,a.y-b.y); }
int operator * (const pr &a,const pr &b){ return a.x*b.y-b.x*a.y; }
double Length(const pr &a){ return sqrt(sqr(1.0*a.x)+sqr(1.0*a.y)); }
int cmp(const pr &a,const pr &b){
    pr v1=a-g[1],v2=b-g[1];double l1=Length(v1),l2=Length(v2);
    if(v1.x*l22 && (g[i]-g[stk[top]])*(g[stk[top]]-g[stk[top-1]])<0) --top;
        stk[++top]=i;
    }
//  cout<<"-----------"<

　　

 

转载于:https://www.cnblogs.com/beiyuoi/p/6052392.html