洛谷4049 [jsoi2007]合金

建图后利用floyd原理求最小代价

Posted by yjjr's blog on October 21, 2018

标签:计算几何,floyd

题目

题目传送门

题目描述

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

现在,用户给出了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),分别表示铁铝锡在一种用户需要的合金中

所占的比重。

输出格式

一个整数,表示最少需要的原材料种数。若无解,则输出–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,所以可以将其缩成二维,看成平面上的一个点

两种原料能够配成的产品一定在这两点的连线/线段上,可以先预处理出两种原料的连线

于是本题转化为在m个点中选取最少的点,将n个目标点包围起来

然后可以用特殊的floyd求最小代价

需要特判下所有点共线/重合的情况

code

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<algorithm>
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define dep(i,a,b) for(int i=a;i>=b;i--)
#define ll long long
#define mem(x,num) memset(x,num,sizeof x)
#define reg(x) for(int i=last[x];i;i=e[i].next)
using namespace std;
inline ll read(){
	ll f=1,x=0;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;
}
//******head by yjjr******
#define inf 0x3f3f3f
#define eps 1e-9
const int maxn=1e3+6;
struct P{double x,y;}a[maxn],b[maxn];
int Map[maxn][maxn],n,m,ans;double k;
double operator * (P a,P b){return a.x*b.y-a.y*b.x;}
P operator - (P a,P b){P t;t.x=a.x-b.x;t.y=a.y-b.y;return t;}
inline bool col(P u,P v){
	if(u.x>v.x)swap(u,v);
	rep(i,1,m)if(b[i].x<u.x||b[i].x>v.x)return 0;
	if(u.y>v.y)swap(u,v);
	rep(i,1,m)if(b[i].y<u.y||b[i].y>v.y)return 0;
	return 1;
}
inline int judge(P x,P y){
	int c1=0,c2=0;
	rep(i,1,m){
		double t=(y-x)*(b[i]-x);
		if(t>eps)c1++;
		if(t<-eps)c2++;
		if(c1*c2)return 0;
	}
	if(!c1&&!c2&&col(x,y)){ans=2;return -1;}
	if(c1)return 1;if(c2)return 2;return 3;
}
void solve(){
	rep(i,1,n)
		rep(j,i+1,n){
			int flag=judge(a[i],a[j]);
			if(flag==-1)return;
			if(flag==1)Map[i][j]=1;
			else if(flag==2)Map[j][i]=1;
			else if(flag==3)Map[i][j]=Map[j][i]=1;
		}
}
void floyd(){
	rep(k,1,n)rep(i,1,n)rep(j,1,n)Map[i][j]=min(Map[i][j],Map[i][k]+Map[k][j]);
}
int main(){
	n=read(),m=read();
	rep(i,1,n)scanf("%lf%lf%lf",&a[i].x,&a[i].y,&k);
	rep(i,1,m)scanf("%lf%lf%lf",&b[i].x,&b[i].y,&k);
	int flag=1;mem(Map,inf);ans=1e9;
	rep(i,1,n)if(fabs(a[i].x-a[1].x)>eps||fabs(a[i].y-a[1].y)>eps)flag=0;
	rep(i,1,m)if(fabs(b[i].x-a[1].x)>eps||fabs(b[i].y-a[1].y)>eps)flag=0;
	if(flag)ans=1;
	solve();floyd();
	if(ans<1e9){cout<<ans<<endl;return 0;}
	rep(i,1,n)ans=min(ans,Map[i][i]);
	if(ans==1e9||ans<=2)puts("-1");else cout<<ans<<endl;
	return 0;
}
本文可以转载,但必须附上原文链接,否则你会终生找不到妹子!!!欢迎关注我的CSDN: ahyjjr