CF 54E Vacuum Сleaner

对于一个凸包形状的吸尘器，问你这个吸尘器在清理矩形的角落时，遗留下的最小面积，可以旋转；

题解是计算从第i个点到第j个点所围成的面积，（从i到j 向量旋转不能超过90度），算出凸包的i与j连线在矩形中对应的长与宽，用长，宽，以及i，j连线构成的三角形减去前面的面积，就是够不到的角落，最后取最小值即可，注意顺时针逆时针都要来一次

#include<stdio.h>
#include<math.h>
#include<iostream>
#include<algorithm>
using namespace std;
struct P
{
    double x,y;
    P(double a=0,double b=0):x(a),y(b){}
    P operator +(const P &a)const{
        return P(x+a.x,y+a.y);
    }
    P operator -(const P &a)const{
        return P(x-a.x,y-a.y);
    }
    P operator *(const double a)const{
        return P(x*a,y*a);
    }
    P operator /(const double a)const{
        return P(x/a,y/a);
    }
    double operator *(const P &a)const{
        return x*a.x+y*a.y;
    }
    double operator ^(const P &a)const{
        return x*a.y-y*a.x;
    }
};
double dis(P a,P b)
{
    return sqrt((b-a)*(b-a));
}
P p[100100];
int n;
double sum[100100];
double solve()
{
    for(int i=n;i<n+n;i++) p[i]=p[i-n];
    sum[0]=0;
    for(int i=1;i<2*n;i++) {sum[i]=sum[i-1]+(p[i]^p[i-1])/2.0;}
    int j=0;
    double ans=1e18;
    for(int i=0;i<n;i++)
    {
        if(j<=i) j=i+1;
        while(j+1<2*n&&(((p[j+1]-p[i])*(p[i+1]-p[i]))>=((p[j]-p[i])*(p[i+1]-p[i])))) j++;
        //cout<<"j="<<j<<endl;
        if(j==i+1) ans=0;
        if(j>=2*n) break;
        double s=sum[j]-sum[i];

        s+=(p[i]^p[j])/2.0;

        double x=((p[j]-p[i])*(p[i+1]-p[i]))/dis(p[i],p[i+1]);
        double y=sqrt(((p[i]-p[j])*(p[i]-p[j]))-x*x);
        //cout<<"s="<<s<<" x="<<x;
        //printf(" y=%.10lf\n",y);
        ans=min(ans,x*y/2.0-fabs(s));
    }
    return ans;
}
int main()
{
    while(scanf("%d",&n)!=EOF)
    {
        for(int i=0;i<n;i++)
            scanf("%lf%lf",&p[i].x,&p[i].y);
        double ans=solve();
        for(int i=0;i<n/2;i++)
            swap(p[i],p[n-i-1]);
        ans=min(ans,solve());
        //printf("%.10lf\n",solve());
        printf("%.10lf\n",ans);
    }
    return 0;
}