bzoj1027 [JSOI2007]合金（计算几何+Floyd）

因为a+b+c=1,所以可以只看前两维，忽略第三维。
把每种合金看做一个二维平面上的点，我们可以证明两种合金的连线上的合金都可以由这两种合成出来。进而我们可以发现，选择的原材料合金围成的凸包中的点我们都可以合成出来。

因而问题就转化为：在点集A中选择一个子集S，使得点集B中的所有点都处于子集S围成的凸包中，最小化子集S的大小。

我们可以利用floyd求最小环来解决：枚举点集A中的任意两点i,j(i可以等于j)，如果点集B中的所有点都处于i->j的同一侧（$\vec{ki}*\vec{kj}>0$，即都是逆时针）或者处于线段ij上（$\vec{ki}*\vec{kj}=0\land\vec{ki}·\vec{kj}<0$）就连有向边i->j，边权为1.最后求一个最小环就是答案了。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 510
#define eps 1e-7
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,m,mp[N][N];
struct P{
    double x,y;
    P(){}
    P(double _x,double _y){x=_x;y=_y;}
    friend double dot(P a,P b){return a.x*b.x+a.y*b.y;}
    friend double det(P a,P b){return a.x*b.y-a.y*b.x;}
    friend P operator-(P a,P b){return P(a.x-b.x,a.y-b.y);}
}a[N],b[N];
inline bool jud(int i,int j){
    for(int k=1;k<=m;++k){
        double res=det(b[k]-a[i],b[k]-a[j]);
        if(res>eps||fabs(res)<eps&&dot(b[k]-a[i],b[k]-a[j])<eps) continue;return 0;
    }return 1;
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();memset(mp,inf,sizeof(mp));
    for(int i=1;i<=n;++i) scanf("%lf%lf%*lf",&a[i].x,&a[i].y);
    for(int i=1;i<=m;++i) scanf("%lf%lf%*lf",&b[i].x,&b[i].y);
    for(int i=1;i<=n;++i)
        for(int j=1;j<=n;++j)
            if(jud(i,j)) mp[i][j]=1;
    for(int k=1;k<=n;++k)
        for(int i=1;i<=n;++i)
            for(int j=1;j<=n;++j)
                mp[i][j]=min(mp[i][j],mp[i][k]+mp[k][j]);
    int ans=inf;for(int i=1;i<=n;++i) ans=min(ans,mp[i][i]);
    if(ans==inf) puts("-1");
    else printf("%d\n",ans);
    return 0;
}