poj3348 Cows (Graham Scan算法学习)

题目

　　http://poj.org/problem?id=3348

题解

　　题目其实就是让你求一个凸包，并求出它的面积。
　　我们使用Graham Scan算法

求解凸包问题的Graham Scan算法

　　我主要是看了哈工大出版社的“计算几何及其应用”学的，虽然最后书上的标程没看懂，但是我还是照自己的理解把程序写出来了。
　　感觉叉积是种很神奇的东西，假设有向量$\vec a,\vec b$，那么当$\vec a\times \vec b > 0$时，$\vec a$ $在 \vec b$的顺时针方向，如果 <0 <script type="math/tex" id="MathJax-Element-982"><0</script>则在逆时针方向，如果$=0$则向量共线；$S_{\triangle ABC}=\frac{|\vec{AB}\times\vec{AC}|}{2}$。
　　我们可以把向量坐标化，用$(x,y)$表示一个向量，那么$\vec a\times \vec b=x_a\times y_b-x_b\times y_a$。向量$\vec c = \vec a-\vec b=(x_b-x_a,y_b-y_a)$
　　显然，所有点中最左下角的点(即按先纵坐标、后横坐标排序得到的第一个点)必在凸包上。那么选取这个点，把它作为新的坐标系原点，重新建系。在新的坐标系中，把除原点之外的所有点按照其与原点所成向量和$x$轴的夹角排序(这夹角就是极角，取值范围$[0,2\pi)$，即从x轴逆时针旋转到向量的正方向所形成的角)，若极角相同，则按离原点距离的远近排序，近的优先。
　　排序好之后，先把原点放进栈，然后每次加入一个新的点，如果它和栈顶点所成的向量$\vec{v_1}$，与栈顶和次栈顶所成的向量$\vec{v_2}$成顺时针夹角(用叉积判断)，则弹出栈顶，用一个while直到上述条件不满足，最后将当前点加入。可以证明，排序后的最后一个点一定也在凸包上。
　　面积的计算：从原点引出若干条线段，把凸包划分成若干个三角形，每次用叉积计算面积，加到答案中去就好了。
　　相关名词：向量vector，叉积cross product

代码

//Graham Scan求凸包 
#include <cstdio>
#include <algorithm>
#define maxn 15000
using namespace std;
struct vec{double x, y;}p[maxn], s[maxn];
int N, top;
double S;
double abs(double x){return x>0?x:-x;}
double cp(vec p1, vec p2)
{return p1.x*p2.y-p2.x*p1.y;}
bool cmp(vec p1, vec p2)
{
    double x=p[1].x, y=p[1].y;
    p1.x-=x, p1.y-=y;
    p2.x-=x, p2.y-=y;
    return cp(p1,p2)==0?(p1.x==p2.x?p1.y<p2.y:p1.x<p2.x):cp(p1,p2)>0;
}
void Graham_Scan()
{
    int i, j;
    double x, y;
    vec v1, v2, pt;
    j=1,x=p[1].x,y=p[1].y;
    for(i=2;i<=N;i++)if(p[i].y<y or p[i].y==y and p[i].x==x)x=p[i].x,y=p[i].y,j=i;
    pt=p[j],p[j]=p[1],p[1]=pt;
    s[++top]=pt;
    sort(p+2,p+N+1,cmp);
    for(i=2;i<=N;i++)
    {
        v1=(vec){s[top].x-s[top-1].x,s[top].y-s[top-1].y};
        v2=(vec){p[i].x-s[top].x,p[i].y-s[top].y};
        while(top>1 and cp(v1,v2)<0)
        {
            top--;
            v1=(vec){s[top].x-s[top-1].x,s[top].y-s[top-1].y};
            v2=(vec){p[i].x-s[top].x,p[i].y-s[top].y};
        }
        s[++top]=p[i];
    }
}
void CalcS()
{
    int i;
    double x=s[1].x, y=s[1].y;
    for(i=3;i<=top;i++)
        S+=0.5*abs(cp((vec){s[i].x-x,s[i].y-y},(vec){s[i-1].x-x,s[i-1].y-y}));
}
int main()
{
    int i, j;
    scanf("%d",&N);
    for(i=1;i<=N;i++)scanf("%lf%lf",&p[i].x,&p[i].y);
    Graham_Scan();
    CalcS();
    printf("%d\n",int(S/50));
    return 0;
}