BZOJ1027[JSOI2007] 合金

bzoj链接：http://www.lydsy.com/JudgeOnline/problem.php?id=1027
洛谷链接：https://www.luogu.org/problemnew/show/P4049

合金

Description

　　某公司加工一种由铁、铝、锡组成的合金。他们的工作很简单。首先进口一些铁铝锡合金原材料，不同种类的
原材料中铁铝锡的比重不同。然后，将每种原材料取出一定量，经过融解、混合，得到新的合金。新的合金的铁铝
锡比重为用户所需要的比重。 现在，用户给出了n种他们需要的合金，以及每种合金中铁铝锡的比重。公司希望能
够订购最少种类的原材料，并且使用这些原材料可以加工出用户需要的所有种类的合金。

Input

　　第一行两个整数m和n（m, n ≤ 500），分别表示原材料种数和用户需要的合金种数。第2到m + 1行，每行三
个实数a, b, c（a, b, c ≥ 0 且 a + b + c = 1），分别表示铁铝锡在一种原材料中所占的比重。第m + 2到m +
n + 1行，每行三个实数a, b, c（a, b, c ≥ 0 且 a + b + c = 1），分别表示铁铝锡在一种用户需要的合金中
所占的比重。

Output

　　一个整数，表示最少需要的原材料种数。若无解，则输出–1。

Sample Input

10 10
0.1 0.2 0.7
0.2 0.3 0.5
0.3 0.4 0.3
0.4 0.5 0.1
0.5 0.1 0.4
0.6 0.2 0.2
0.7 0.3 0
0.8 0.1 0.1
0.9 0.1 0
1 0 0
0.1 0.2 0.7
0.2 0.3 0.5
0.3 0.4 0.3
0.4 0.5 0.1
0.5 0.1 0.4
0.6 0.2 0.2
0.7 0.3 0
0.8 0.1 0.1
0.9 0.1 0
1 0 0

Sample Output

5

题解

首先，很妙的一点是，由于三种元素比例相加等于一，所以第三种元素比例可以通过前两种算出来，因此第三维可以直接扔掉。
然后，如果把前两维作为坐标将每块合金映射到二维平面上，可以观察到非常显然的性质：
两块合金能炼成的合金在两点的连线上，三块即以上的合金能炼出的合金在这些点的凸包内。
因此，我们需要求点最少的能够包住客户需求的凸包。
考虑到数据量较小，可以用$O(n^2m)$求解。
这样我们思路就很清晰了，$O(m^2)$枚举每一条边，$O(n)$检验，如果客户需求都在它的一侧，那么就连一条边。最后，跑一遍floyed求最小环即可AC。

细节

虽然$O(2n^2m)$看起来可以稳过这道题，然而博主洛谷上T了四个点，加了各种常数优化才A，bzoj上无压力过呀。。。
另外，需要特判所有点共线或共点的情况，因为超时而且没有这样的数据，所以博主水着写的SPJ，居然A了。。。

代码

#include<bits/stdc++.h>
#define db double
using namespace std;
const int M=505;
const db eps=1e-8;
struct pt{db x,y;};
db operator * (pt a,pt b){return a.x*b.y-a.y*b.x;}
pt operator - (pt a,pt b){return (pt){a.x-b.x,a.y-b.y};}
db dot(pt a,pt b){return a.x*b.x+a.y*b.y;}
pt mat[M],req[M];
int dis[M][M],m,n;
int sig(db x){return (x>eps)-(x<-eps);}
void in()
{
    db kkk;
    scanf("%d%d",&m,&n);
    for(int i=1;i<=m;++i)
    scanf("%lf%lf%lf",&mat[i].x,&mat[i].y,&kkk);
    for(int i=1;i<=n;++i)
    scanf("%lf%lf%lf",&req[i].x,&req[i].y,&kkk);
}
void spj1()
{
    if(m==1){printf("1");exit(0);}
}
void build()
{
    db x;
    bool flag;
    memset(dis,127/3,sizeof(dis));
    for(int i=1;i<=m;++i)
    for(int j=1;j<=m;++j)
    {
        if(i==j) continue;
        flag=1;
        for(int k=1;k<=n;++k)
        {
            x=(mat[i]-req[k])*(mat[j]-req[k]);
            if(sig(x)>0){flag=0;break;}
            if(sig(x)==0&&sig(dot(mat[i]-req[k],mat[j]-req[k]))>0)
            {flag=0;break;}
        }
        if(flag) dis[i][j]=1;
    }
}
void ac()
{
    for(int k=1;k<=m;++k)
    for(int i=1;i<=m;++i)
    for(int j=1;j<=m;++j)
    dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);
    int ans=1e9;
    for(int i=1;i<=m;++i)
    ans=min(dis[i][i],ans);
    if(ans>=1e3) printf("-1");
    else printf("%d",ans);
}
int main()
{
    in();spj1();
    build();ac();
    return 0;
}