Topcoder SRM 566 DIV2 1000 FencingPenguinsEasy

FencingPenguinsEasy

Task:
给定$n$个栅栏每个栅栏与原点的距离都为$r$，第一个栅栏坐标为$(r,0)$，其余栅栏均匀排布为一正多边形。
现有$p$个企鹅，你需要通过连线栅栏形成一个封闭图形，将这些企鹅全部包含在内。求此围栏的最小面积，无法完成输出-1。（$3<=n<=222,5<=r<=100000,1<=p<=50$）
Solution:
首先我们先考虑如何判断企鹅和多边形一条边的位置关系。
较为简单的一种方法是考虑直线的方向，如果直线、企鹅成逆时针，则说明在多边形内侧。如图所示：

判断逆时针可以采用向量的叉积（详见代码）。
而计算最终结果很显然就是一个区间dp了，每次引入的新点带来的面积增量就是（l,nl,r)这个三角形，用向量算算就好了。注意最终点还要检验能否连回第一个点。

class FencingPenguinsEasy {
public:
    vector<int>X,Y;
    double dp[250][250];
    bool judge_line(double x1,double y1,double x2,double y2){
        for(int i=0;i<X.size();i++){
            int x3=X[i],y3=Y[i];
            double t=(x1-x3)*(y2-y1)-(y1-y3)*(x2-x1);
            if(t<0)return false;
        }
        return true;
    }
    void Chk(double &a,double b){
        if(a>b)a=b;
    }
    double calcArea(double x1, double y1, double x2, double y2, double x3, double y3) {
        x1-=x3;y1-=y3;x2-=x3;y2-=y3;
        return(x1*y2-x2*y1)/2;
    }
    double calculateMinArea(int n, int r, vector <int> x, vector <int> y) {
        const double INF=1e50;
        double pi=acos(-1);
        for(int i=0;i<x.size();i++)X.push_back(x[i]);
        for(int i=0;i<y.size();i++)Y.push_back(y[i]);
        double Px[250],Py[250];
        double degree=2*pi/n;
        int sz=0;
        for(int i=0;i<n;i++){//各栅栏点坐标 
            sz++;
            Px[sz]=(r*cos(degree*i));
            Py[sz]=(r*sin(degree*i));
        }
        for(int i=2;i<=sz;i++)
            if(!judge_line(Px[i-1],Py[i-1],Px[i],Py[i]))return -1.0;
        if(!judge_line(Px[sz],Py[sz],Px[1],Py[1]))return -1.0;
        for(int i=1;i<=sz;i++)
            for(int j=1;j<=sz;j++)
                if(i==j)dp[i][j]=0;
                else dp[i][j]=INF;
        for(int r=2;r<=sz;r++){
            for(int l=r-1;l>=1;l--)
                for(int k=l+1;k<=r;k++)
                    if(judge_line(Px[l],Py[l],Px[k],Py[k]))Chk(dp[l][r],dp[k][r]+abs(calcArea(Px[k],Py[k],Px[l],Py[l],Px[r],Py[r])));
        }
        double ans=INF;
        for(int i=1;i<=sz;i++)
            for(int j=i+2;j<=sz;j++)
                if(judge_line(Px[j],Py[j],Px[i],Py[i]))Chk(ans,dp[i][j]);//j-i连线是否合法 
        return ans;
    }
};