[JSOI2007]合金，bzoj1027，非常巧妙的计算几何

导语

      出题人牛逼。

正题

      题目给出的就是原材料和目标合金，通过观察我们发现第三种金属是没有用的，因为其他两种金属得到之后，第三种金属自然然就出来了，所以就变成了两个金属。

      我们把这两个金属当成横纵坐标画在平面直角坐标系上面发现，两种原材料合成的合金好在这两种原材料的线段上面。

      那么几种原材料合成的新合金必定在这凸包里面（证明大概是就是可以知道边边上的可以合成，中间的可以通过边边上的两个点所构成的线段覆盖。

      问题就变成了，求覆盖n个点的最小凸包。

      考虑叉积+floyd，对于所有目标点都在其逆时针的线段，把它连边，这个东西可以通过叉积n^3得到。

      然后就是floyd跑最小环即可。

      注意只有两个条件就是k!=i && k!=j 。

      判断一些特殊的情况了。

      虽然我也A不了。

      下面这个是比较靠谱的std。。

 
#include<cstdio>
#include<cstring>
#include<cmath>
#define eps 1e-8
using namespace std;
struct node
{
    double x,y,z;
}a[510],b[510];
double multi(node p1,node p2,node p0)
{
    double x1=p1.x-p0.x;
    double y1=p1.y-p0.y;
    double x2=p2.x-p0.x;
    double y2=p2.y-p0.y;
    return x1*y2-x2*y1;
}
double dis(node p1,node p2)
{
    return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}
int d[510][510];
int main()
{
    int n,m;
    scanf("%d%d",&m,&n);
    for(int i=1;i<=m;i++) scanf("%lf%lf%lf",&a[i].x,&a[i].y,&a[i].z);
    for(int i=1;i<=n;i++) scanf("%lf%lf%lf",&b[i].x,&b[i].y,&b[i].z);
    int i,j,k;double t;
    memset(d,63,sizeof(d));
    for(i=1;i<=m;i++)
        for(j=1;j<=m;j++)
        {
            for(k=1;k<=n;k++)
            {
                t=multi(a[j],b[k],a[i]);
                if(t<-eps) break;
                if( fabs(t)<eps && dis(a[i],b[k])>dis(a[i],a[j]) )break;
            }
            if(k==n+1) d[i][j]=1;
        }
    for(k=1;k<=m;k++)
        for(i=1;i<=m;i++)if(i!=k)
            for(j=1;j<=m;j++)if(j!=k)
                if( d[i][j]>d[i][k]+d[k][j] )d[i][j]=d[i][k]+d[k][j];
    int  ans=m+1;
    for(i=1;i<=m;i++) if(ans>d[i][i])ans=d[i][i];
    if(ans==m+1)printf("-1\n");else printf("%d\n",ans);
    return 0;
}

下面是我的，有问题请麻烦一定一定告诉我。（拜托

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<cmath>
using namespace std;
 
int m,n;
struct op{
    double x,y;
    double operator^(op a)const{
        return x*a.y-a.x*y;
    }
}s[510],p[510];
int d[510][510];
double eps=1e-8;
 
double dist(op a){
    return sqrt(a.x*a.x+a.y*a.y);
}
 
int main(){
    scanf("%d %d",&m,&n);
    for(int i=1;i<=m;i++) scanf("%lf %lf %lf",&s[i].x,&s[i].y,&s[i+1].x);
    for(int i=1;i<=n;i++) scanf("%lf %lf %lf",&p[i].x,&p[i].y,&p[i+1].x);
    op l1,l2;
    memset(d,63,sizeof(d));
    for(int i=1;i<=m;i++)
        for(int j=1;j<=m;j++){
            bool tf=true;
            for(int k=1;k<=n;k++){
                l1=(op){s[j].x-s[i].x,s[j].y-s[i].y},l2=(op){p[k].x-s[i].x,p[k].y-s[i].y};
                if((l1^l2)<-eps || (fabs(l1^l2)<eps && dist(l1)<dist(l2))) {
                    tf=false;break;
                }
            }
            if(tf) d[i][j]=1;
        }
    for(int k=1;k<=m;k++)
        for(int i=1;i<=m;i++)if(i!=k)
            for(int j=1;j<=m;j++)if(j!=k)
                d[i][j]=min(d[i][j],d[i][k]+d[k][j]);
    int ans=m+1;
    for(int i=1;i<=m;i++)
        ans=min(ans,d[i][i]);
    if(ans!=m+1) printf("%d",ans);
    else printf("-1");
}