BJ模拟 帐篷【凸包】

题目描述：

n,m<=400

解题思路：

当n,m很小时，直接暴力枚举起点，用三角剖分的方式，状态记录当前点和上一个点dp，复杂度为$O(n^4)$
当n,m较大时，h已经远超1e9z，直接求周长最小凸包即可，注意周长最小时自动满足凸性，所以可以只枚举起点，状态记录当前点dp即可，也可以用floyed求最小环的方式dp，时间复杂度皆为$O(n^3)$

#include<bits/stdc++.h>
using namespace std;
int getint()
{
    int i=0,f=1;char c;
    for(c=getchar();c!='-'&&(c<'0'||c>'9');c=getchar());
    if(c=='-')f=-1,c=getchar();
    for(;c>='0'&&c<='9';c=getchar())i=(i<<3)+(i<<1)+c-'0';
    return i*f;
}

const int N=1005;
const double INF=1e50;
struct point
{
    double x,y;int id;
    point(){}
    point(double _x,double _y):x(_x),y(_y){}
    inline friend point operator - (const point &a,const point &b){return point(a.x-b.x,a.y-b.y);}
    inline friend double operator * (const point &a,const point &b){return a.x*b.y-a.y*b.x;}
    inline double dis(){return sqrt(x*x+y*y);}
}a[N],b[N],p[N];
int n,m;
bool ok[N][N];
double f[N][N],g[N][N],fib[N],h;
inline bool cmp(const point &a,const point &b)
{
    double det=(a-p[1])*(b-p[1]);
    if(det)return det>0;
    else return (a-p[1]).dis()<(b-p[1]).dis();
}
double cs(int i,int j)
{
    double d=(p[j]*p[i])/(p[j]-p[i]).dis();
    return d*d;
}
double calc(int o)
{
    int cnt=1;p[1]=a[o],p[1].id=o;
    for(int i=1;i<=n;i++)if(a[i].x>a[o].x||a[i].x==a[o].x&&a[i].y>a[o].y)
        p[++cnt]=a[i],p[cnt].id=i;
    sort(p+2,p+cnt+1,cmp);
    for(int i=1;i<=cnt;i++)
        for(int j=1;j<=i;j++)
            f[i][j]=INF;
    f[1][1]=0;
    for(int i=1;i<=cnt;i++)
        for(int j=1;j<i;j++)if(ok[p[j].id][p[i].id])
            for(int k=1;k<=j;k++)
                if((p[j]-p[k])*(p[i]-p[k])>=0)
                    f[i][j]=min(f[i][j],f[j][k]+(p[i]-p[j]).dis()*sqrt(h*h+cs(i,j))/h); 
    double res=INF; 
    for(int i=2;i<=cnt;i++) if(ok[p[i].id][p[1].id])
            for(int j=1;j<i;++j)
                if((p[i]-p[j])*(p[1]-p[j])>=0)
                    res=min(res,f[i][j]+(p[1]-p[i]).dis()*sqrt(h*h+cs(1,i))/h);
    return res;
}
void solve1()
{
    fib[0]=fib[1]=1;
    for(int i=2;i<=n+m;i++)fib[i]=fib[i-1]+fib[i-2];
    double ans=INF;h=fib[n+m];
    for(int i=1;i<=n;i++)ans=min(ans,calc(i));
    if(ans==INF)puts("NO SOLUTION");
    else printf("%0.3f\n",ans/2);
}
void solve2()
{
    double ans=INF;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
        {
            f[i][j]=g[i][j]=INF;
            if(ok[i][j])f[i][j]=g[i][j]=(a[i]-a[j]).dis();
        }
    for(int k=1;k<=n;k++)
    {
        for(int i=1;i<k;i++)
            for(int j=1;j<i;j++)
                ans=min(ans,min(f[i][j]+g[j][k]+g[k][i],f[j][i]+g[i][k]+g[k][j]));
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                f[i][j]=min(f[i][j],f[i][k]+f[k][j]);
    }
    if(ans==INF)puts("NO SOLUTION");
    else printf("%0.3f\n",ans/2);
}
int main()
{
    //freopen("lx.in","r",stdin);
    //freopen("lx.out","w",stdout);
    while(scanf("%d %d",&n,&m)!=EOF)
    {
        for(int i=1;i<=n;i++)a[i].x=getint(),a[i].y=getint();
        for(int i=1;i<=m;i++)b[i].x=getint(),b[i].y=getint();
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)if(i!=j)
            {
                ok[i][j]=1;
                for(int k=1;k<=m&&ok[i][j];k++)
                    if((a[j]-a[i])*(b[k]-a[i])<=0)ok[i][j]=0;
            }
        if(n+m<=60)solve1();
        else solve2();  
    }
    return 0;
}