bzoj1091: [SCOI2003]切割多边形

Description

　　有一个凸p边形(p<=8)，我们希望通过切割得到它。一开始的时候，你有一个n*m的矩形，即它的四角的坐标分
别为(0,0), (0,m), (n,0), (n,m)。每次你可以选择一条直线把当前图形切割成两部分，保留其中一个部分（另一
部分扔掉）切割线的长度为此直线在多边形内部的部分的长度。求出最短的切割线总长度。下面是一个例子。我们
需要得到中间的多边形。

　　分别沿着直线1，2，3，4进行切割即可，得到中间的四边形。

Input

　　第一行有两个整数n, m(0 < n,m < 500)，第二行为一个整数p(3<=p<=8)。以下p行每行为两个整数x, y(0 < x
 < n, 0 < y < m)，为按顺时针给出的各顶点坐标。数据保证多边形的是凸的，无三点共线。输入数据无错误。

Output

　　仅一行，为最短切割线的总长度，四舍五入到小数点后3位。允许有0.001的误差。

枚举切割顺序，然后模拟加入直线，暴力计算与其余直线的交点

#include<cstdio>
#include<cmath>
#include<algorithm>
int n,id[],lp=;
double v1,v2,ans=1e10;
struct pos{
    double x,y;
    void init(){scanf("%lf%lf",&x,&y);}
    pos operator+(pos a){return (pos){x+a.x,y+a.y};}
    pos operator-(pos a){return (pos){x-a.x,y-a.y};}
    pos operator*(double a){return (pos){x*a,y*a};}
    double operator*(pos a){return x*a.y-y*a.x;}
    double dot(pos a){return x*a.x+y*a.y;}
    double abs(){return sqrt(x*x+y*y);}
}ps[];
double mn,mx;
struct line{
    pos a,b;
    void chk(line w){
        double c=w.b*b;
        if(c==)return;
        c=(a*w.b+w.b*w.a)/c;
        if(c>0.5)c<mx&&(mx=c);
        else c>mn&&(mn=c);
    }
}ls[],l0[];
int main(){
    scanf("%lf%lf%d",&v1,&v2,&n);
    for(int i=;i<=n;++i)ps[i].init(),id[i]=i;
    ps[n+]=ps[];
    pos p1=(pos){,},p2=(pos){v1,},p3=(pos){v1,v2},p4=(pos){,v2};
    ls[lp++]=(line){p1,p2-p1};
    ls[lp++]=(line){p2,p3-p2};
    ls[lp++]=(line){p3,p4-p3};
    ls[lp++]=(line){p4,p1-p4};
    for(int i=;i<=n;++i)l0[i]=(line){ps[i],ps[i+]-ps[i]};
    do{
        lp=;
        double s=;
        for(int i=;i<=n;++i){
            int w=id[i];
            mn=-1e10,mx=1e10;
            for(int j=;j<lp;++j)l0[w].chk(ls[j]);
            ls[lp++]=l0[w];
            s+=(mx-mn)*l0[w].b.abs();
        }
        if(s<ans)ans=s;
    }while(std::next_permutation(id+,id+n+));
    printf("%.3f",ans);
    return ;
}