Central Europe Regional Contest 2016 C.Convex Contour

EM……凸包的题（暴力也可以写）

题目大意：

给你若干个正三角形，正方形，圆形，求包围他们一周所需要的最小长度。用凸包写完特殊判断下端点情况即可。

用的是Graham扫描算法

#include <iostream>
#include <bits/stdc++.h>
using namespace std;
const double PI = acos(-1.0);

double eps=1e-10;
double add(double a,double b)
{
    if(abs(a+b)<eps*(abs(a)+abs(b))) return 0;     /**考虑误差的加法运算**/
    return a+b;
}
struct P{
    double x,y;
    P(){}
    P(double x,double y):x(x),y(y){}
    P operator - (P p)                            /**重载-号**/
    {
        return P(add(x,-p.x),add(y,-p.y));
    }
    double det(P p)                              /**定义外积**/
    {
        return add(x*p.y,-y*p.x);
    }
};

bool cmp_x(const P& p,const P& q)
{
    if(p.x!=q.x) return p.x<q.x;
    return p.y<q.y;
}
vector<P> convex_hull(P *ps,int n)            /**抄模板**/
{
    sort(ps,ps+n,cmp_x);
    int k=0;
    vector<P> qs(n*2);
    for(int i=0;i<n;i++)
    {
        while(k>1 && (qs[k-1]-qs[k-2]).det(ps[i]-qs[k-1])<=0) k--;
        qs[k++]=ps[i];
    }
    for(int i=n-2,t=k;i>=0;i--)
    {
        while(k>t && (qs[k-1]-qs[k-2]).det(ps[i]-qs[k-1])<=0) k--;
        qs[k++]=ps[i];
    }
    qs.resize(k-1);
    return qs;
}
double dis(vector<P> v)                        /**计算凸包边长**/
{
    double ans=0;
    for(int i=0;i<v.size()-1;i++)
    {
        ans+=sqrt((v[i].x-v[i+1].x)*(v[i].x-v[i+1].x)+(v[i].y-v[i+1].y)*(v[i].y-v[i+1].y));
    }
    ans+=sqrt((v[0].x-v[v.size()-1].x)*(v[0].x-v[v.size()-1].x)+(v[0].y-v[v.size()-1].y)*(v[0].y-v[v.size()-1].y));
    return ans;
}
char s[25];
int N;
int n;
int findnotT(int dir)    // 1 left 2 right      
{                                       /**从左往右或从右往左找不是三角形的形状**/
    int ans=-1;
    if(dir==1)
    {
        for(int i=n-1;i>=0;i--)
        {
            if(s[i]!='T')
            {
                ans=i;
                break;
            }
        }
    }else{
        for(int i=0;i<n;i++)
        {
            if(s[i]!='T')
            {
                ans=i;
                break;
            }
        }
    }
    return ans;
}
P ps[100];
int main()
{
    while(~scanf("%d",&n))
    {
        scanf("%s",s);
        int cnt=0;
        for(int pos=0;pos<n;pos++)
        {
            if(s[pos]=='C')                      /**圆形的情况将最高点和最低点放入**/
            {
                ps[cnt].x=pos+0.5;
                ps[cnt++].y=0;
                ps[cnt].x=pos+0.5;
                ps[cnt++].y=1;
            }else if(s[pos]=='T')
            {
                ps[cnt].x=pos;
                ps[cnt++].y=0;
                ps[cnt].x=pos+0.5;
                ps[cnt++].y=sqrt(3)/2;
                ps[cnt].x=pos+1;
                ps[cnt++].y=0;
            }else{
                ps[cnt].x=pos;
                ps[cnt++].y=0;
                ps[cnt].x=pos;
                ps[cnt++].y=1;
                ps[cnt].x=pos+1;
                ps[cnt++].y=0;
                ps[cnt].x=pos+1;
                ps[cnt++].y=1;
            }
        }
        vector<P> pp=convex_hull(ps,cnt);
        double ans=dis(pp);
        if(s[0]=='C')                       /**因为上面计算的时候把圆当做线段处理所以这里修正一下**/
        {
            ans=ans-1+PI/2;
        }
        if(s[n-1]=='C')
        {
            ans=ans-1+PI/2;
        }
        if(s[0]=='T')
        {
            int index=findnotT(2);
            if(index!=-1)
            {
                if(s[index]=='C')            /**计算三角形和圆形所形成的长度**/
                {
                    ans=ans-sqrt(index*index+(1-sqrt(3)/2)*(1-sqrt(3)/2));          
                    double t0=sqrt(index*index+(0.5-sqrt(3)/2)*(0.5-sqrt(3)/2));
                    double af1=atan((sqrt(3)/2-0.5)/index);
                    double af2=acos(0.5/t0);
                    double af=PI/2-af1-af2;
                    ans+=sqrt(t0*t0-0.25);
                    ans+=af/2;
                }
            }
        }
        if(s[n-1]=='T')
        {
            int index=findnotT(1);
            if(index!=-1)
            {
                if(s[index]=='C')
                {
                    ans=ans-sqrt((n-index-1)*(n-index-1)+(1-sqrt(3)/2)*(1-sqrt(3)/2));
                    double t0=sqrt((n-index-1)*(n-index-1)+(0.5-sqrt(3)/2)*(0.5-sqrt(3)/2));
                    double af1=atan((sqrt(3)/2-0.5)/(n-index-1));
                    double af2=acos(0.5/t0);
                    double af=PI/2-af1-af2;
                    ans+=sqrt(t0*t0-0.25);
                    ans+=af/2;
                }
            }
        }
        printf("%.12f\n",ans);
    }
    return 0;
}

就这样~