[Luogu P4049] [JSOI2007]合金

洛谷传送门

题目描述

某公司加工一种由铁、铝、锡组成的合金。他们的工作很简单。首先进口一些铁铝锡合金原材料，不同种类的原材料中铁铝锡的比重不同。然后，将每种原材料取出一定量，经过融解、混合，得到新的合金。新的合金的铁铝锡比重为用户所需要的比重。

现在，用户给出了n种他们需要的合金，以及每种合金中铁铝锡的比重。公司希望能够订购最少种类的原材料，并且使用这些原材料可以加工出用户需要的所有种类的合金。

输入输出格式

输入格式：

第一行两个整数m和n（$m,n\le 500$），分别表示原材料种数和用户需要的合金种数。

第2到m + 1行，每行三个实数a, b, c（$a,b,c\ge 0$ 且 $a+b+c=1$），分别表示铁铝锡在一种原材料中所占的比重。

第m + 2到m +n + 1行，每行三个实数a, b, c（$a, b, c ≥ 0 $且$ a + b + c = 1$），分别表示铁铝锡在一种用户需要的合金中

所占的比重。

输出格式：

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

输入输出样例

输入样例#1:

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

输出样例#1:

5

解题分析：

我们可以发现其实知道两种金属的比例即可确定金属， 所以我们可以将比例转化到二维图上分析。

如果我们以横坐标表示铁的比重， 用纵坐标表示铝的比重， 那么显然如果我们能造出某种比重的金属， 其表示的“产品”点一定落在这些点围成的封闭图形中。

于是我们的问题就变成了： 删去尽量多的“原材料”点使得所要求的点落在凸包中。这时我们可以枚举最大凸包上以及其中的“原材料”点， 看是否所有“产品点”都在它的一侧， 若成立则在这两个点间连一条边长为1 的边（注意我们必须规定方向， 否则无法分辨在凸包外还是凸包内）。 最后使用floyd跑一边求出最小环即可得到最小的“原材料点”个数。

#include <cstdio>
#include <cstring>
#include <cctype>
#include <cstdlib>
#include <cmath>
#include <algorithm>
#define R register
#define W while
#define IN inline
#define gc getchar()
namespace Geometry
{
    using std::min;
    #define EPS 1e-8
    #define MX 505
    #define db double
    struct pt
    {
        db x, y;
        IN bool operator < (const pt &a) const
        {
            return x==a.x ? y < a.y : x < a.x;
        }
    };
    IN pt operator + (const pt &x, const pt &y){return (pt){x.x + y.x, x.y + y.y};}
    IN pt operator - (const pt &x, const pt &y){return (pt){x.x - y.x, x.y - y.y};}
    IN db operator * (const pt &x, const pt &y){return x.x * y.y - x.y * y.x;}
    IN db operator ^ (const pt &x, const pt &y){return x.x * y.x + x.y * y.y;}
    IN int sig (const db &x){return (x > -EPS) - (x < EPS);}
    IN db dist (const pt &x){return sqrt(x.x * x.x + x.y * x.y);}
    pt mat[MX];
    pt req[MX];
    int mat_num, req_num;
    int mp[MX][MX], ans = 2e9;
    void floyd()
    {
        for (R int i = 1; i <= mat_num; ++i)
            for (R int j = 1; j <= mat_num; ++j)
                for (R int k = 1; k <= mat_num; ++k)
                    if (mp[j][k] > mp[j][i] + mp[i][k]) mp[j][k] = mp[j][i] + mp[i][k];
        for (R int i = 1; i <= mat_num; ++i)
            ans = min(ans, mp[i][i]);
    }
}
using namespace std;
using namespace Geometry;
int main()
{
    db cross;
    scanf("%d%d", &mat_num, &req_num);
    for (R int i = 1; i <= mat_num; ++i)
        scanf("%lf%lf%lf", &mat[i].x, &mat[i].y, &cross);
    for (R int i = 1; i <= req_num; ++i)
        scanf("%lf%lf%lf", &req[i].x, &req[i].y, &cross);
    memset(mp, 63, sizeof(mp));
    for (R int i = 1; i <= mat_num; ++i)
        for (R int j = 1; j <= mat_num; ++j)
        {
            for (R int k = 1; k <= req_num; ++k)
            {
                cross = (mat[i] - req[k]) * (mat[j] - req[k]);
                if(cross > EPS) goto e1;//方向反了
                if(!sig(cross) && ((mat[i] - req[k]) ^ (mat[j] - req[k])) > EPS) goto e1;//穿过“产品点”
            }
            mp[i][j] = 1;
            e1:
            ;
        }
        floyd();
        if(ans > 1e9) printf("-1");
        else printf("%d", ans);
        return 0;
}