POJ 1113

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

求凸包....

完全上模板的，先过了，然后自己在仔细看看Graham算法

#include <stdio.h>
#include <iostream>
#include <math.h>
#include <algorithm>
using namespace std;
const int maxn =1005;
const float pi = 3.14159265;
struct Point
{
    float x,y;
}p[maxn],ch[maxn];

double xmult(Point p1,Point p2,Point p0){
    return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}
double dis(Point p, Point q){
    return (p.x-q.x)*(p.x-q.x)+(p.y-q.y)*(p.y-q.y);
}
int cmp(Point a,Point b){
    float mu=xmult(a,b,p[0]);
    if(mu>0) return 1;
    if((mu==0)&&(dis(p[0],a)<dis(p[0],b))) return 1;
    return 0;
}
int Graham(int n){
    int i,top,k;
    Point tmp;
    k=0; top=2;
    if(n<=3){
        for(i=0;i<4;i++) ch[i]=p[i];
        return n;
    }
    for(i=1;i<n;i++)
        if((p[i].y<p[k].y)||((p[i].y==p[k].y)&&(p[i].x>p[k].x))) k=i;
    tmp=p[0]; p[0]=p[k]; p[k]=tmp;
    sort(p+1,p+n,cmp);
    ch[0]=p[0]; ch[1]=p[1]; ch[2]=p[2];
    for(i=3;i<n;i++){
        if((i+1<n)&&(xmult(p[i],p[i+1],p[0])==0)) 
            continue;
        while(xmult(p[i],ch[top],ch[top-1])>=0) top--;
        ch[++top]=p[i];
    }
    return ++top;
}

int main()
{
    int n,l,i,top;
    double num;
    scanf("%d%d",&n,&l);
    for(i=0;i<n;i++)
        scanf("%f%f",&p[i].x,&p[i].y);    
    top=Graham(n);
    num=0;
    for(i=0;i<top;i++)
    {
        num+=sqrt(dis(ch[i%top],ch[(i+1)%top]));
    }
    num+=2*pi*l;
    if(num-(int)num >= 0.5)
        printf("%d\n",(int)num+1);
    else
        printf("%d\n",(int)num);
    return 0;
}

本来一直错，因为原来的模板是qsort，被我改成了sort。

然后在sort()中写错了长度，还有自己对Graham中的一些步骤刚开始的时候有些不了解....

刚开始的时候以为当共线时只计算一次，且将共线的距离按照从大到小的排列。

实际上共线时，距离是从小到大，然后是不计算前面的，只有在最后一个不是共线的点的时候才将这个点放入点集中。