bzoj1027 [JSOI2007]合金 floyd+计算几何

Description

---

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

第一行两个整数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），分别表示铁铝锡在一种用户需要的合金中
所占的比重。

Solution

---

首先只有两维有用
注意到两个东西融合只能造出它们连成线段上的点，三个不共线的则可以造出它们构成三角形内部的点，推广一下就是凸包的点集可以造出内部的任意点

考虑怎么求最小点集使得凸包包含所有关键点。我们枚举任意两点连一条边，所有关键点都在这条边的一侧说明它可能出现在凸包上，记边权为1，否则为INF
用建出来的图跑floyd，找到一个最小环就是答案

这题的特判很多，比如全部都在一个点，比如只能选出一条线

Code

---

#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <math.h>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define fill(x,t) memset(x,t,sizeof(x))

const double eps=1e-10;
const int N=2005;

struct pos {
	double x,y;
	pos operator -(pos b) {
		return (pos) {x-b.x,y-b.y};
	}
	double operator *(pos b) {
		return b.y*x-y*b.x;
	}
} a[N],b[N];

int f[N][N],n,m;

bool pd(pos x,pos y) {
	if (x.x>y.x) std:: swap(x,y);
	rep(i,1,m) {
		if (b[i].x<x.x||b[i].x>y.x) return 0;
	}
	if (x.y>y.y) std:: swap(x,y);
	rep(i,1,m) {
		if (b[i].y<x.y||b[i].y>y.y) return 0;
	}
	return 1;
}

int check(pos x,pos y) {
	int cnt1=0,cnt2=0;
	rep(i,1,m) {
		double tmp=(y-x)*(b[i]-x);
		if (tmp>eps) cnt1++;
		if (tmp<-eps) cnt2++;
		if (cnt1&&cnt2) return 0;
	}
	if (!cnt1&&!cnt2&&pd(x,y)) {
		puts("2"); return -1;
	}
	if (cnt1) return 1;
	if (cnt2) return 2;
	return 3;
}

int main(void) {
	freopen("data.in","r",stdin);
	scanf("%d%d",&n,&m); double wjp;
	bool flag=false;
	rep(i,1,n) {
		scanf("%lf%lf%lf",&a[i].x,&a[i].y,&wjp);
		if (fabs(a[i].x-a[1].x)>eps) flag=true;
		if (fabs(a[i].y-a[1].y)>eps) flag=true;
	}
	rep(i,1,m) {
		scanf("%lf%lf%lf",&b[i].x,&b[i].y,&wjp);
		if (fabs(b[i].x-a[1].x)>eps) flag=true;
		if (fabs(b[i].y-a[1].y)>eps) flag=true;
	}
	if (!flag) return 0&puts("1");
	fill(f,31);
	rep(i,1,n) rep(j,i+1,n) {
		int res=check(a[i],a[j]);
		if (res==-1) return 0;
		if (res==1) f[i][j]=1;
		if (res==2) f[j][i]=1;
		if (res==3) f[i][j]=f[j][i]=1;
	}
	rep(k,1,n) rep(i,1,n) rep(j,1,n) {
		f[i][j]=std:: min(f[i][j],f[i][k]+f[k][j]);
	}
	int ans=f[0][0];
	rep(i,1,n) ans=std:: min(ans,f[i][i]);
	if (ans==f[0][0]||ans<=2) ans=-1;
	printf("%d\n", ans);
	return 0;
}