【bzoj 1027】 合金 【JSOI2007】

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

这道题可以将每一种原料和需求的前两个数看作二维平面上两个点，那么任意三种原料所能制造的合金就在这三个点所围成的面上，这题就转化为从N个原料代表的点中选择尽量少的点围住所有原料代表的点，用floyd求最小环即可，下面是程序：

#include<stdio.h>
#include<string.h>
#include<iostream>
using namespace std;
const int N=1005,INF=20000;
const double eps=1e-8;
struct Point{
    double x,y;
    double operator *(const Point p)const{
        return x*p.y-y*p.x;
    }
    double operator ^(const Point p)const{
        return x*p.x+y*p.y;
    }
    Point operator -(const Point p)const{
        Point s;
        s.x=x-p.x;
        s.y=y-p.y;
        return s;
    }
}a[N],b[N];
int n,m,f[N][N];
double abs(double n){
	return n>0?n:-n;
}
bool check(int &x,int &y){
    int i;
    for(i=1;i<=m;i++){
        double t=(a[x]-b[i])*(a[y]-b[i]);
        if(t>eps){
            return 0;
        }
        if(abs(t)<eps&&((a[x]-b[i])^(a[y]-b[i]))>eps){
        	return 0;
		}
    }
    return 1;
}
bool SPJ(){
	int i;
	for(i=1;i<=n;i++){
		if(abs(a[i].x-a[1].x)>eps||abs(a[i].y-a[1].y)>eps){
			return 0;
		}
	}
	for(i=1;i<=m;i++){
		if(abs(b[i].x-a[1].x)>eps||abs(b[i].y-a[1].y)>eps){
			return 0;
		}
	}
	return 1;
}
int main(){
    int i,j,k;
    double tp;
    scanf("%d%d",&n,&m);
    for(i=1;i<=n;i++){
        scanf("%lf%lf%lf",&a[i].x,&a[i].y,&tp);
    }
    for(i=1;i<=m;i++){
        scanf("%lf%lf%lf",&b[i].x,&b[i].y,&tp);
    }
    if(!m){
        printf("0\n");
        return 0;
    }
    if(SPJ()){
    	printf("1\n");
    	return 0;
	}
    for(i=1;i<=n;i++){
        for(j=1;j<=n;j++){
            f[i][j]=(i!=j&&check(i,j))?1:INF;
        }
    }
    int ans=INF;
    for(k=1;k<=n;k++){
        for(i=1;i<=n;i++){
            for(j=1;j<=n;j++){
                f[i][j]=min(f[i][j],f[i][k]+f[k][j]);
            }
        }
    }
    for(i=1;i<=n;i++){
    	ans=min(ans,f[i][i]);
	}
    printf("%d\n",ans==INF?-1:ans);
    return ~~(0-0);
}