Description
某公司加工一种由铁、铝、锡组成的合金。他们的工作很简单。首先进口一些铁铝锡合金原材料,不同种类的原材料中铁铝锡的比重不同。然后,将每种原材料取出一定量,经过融解、混合,得到新的合金。新的合金的铁铝锡比重为用户所需要的比重。 现在,用户给出了n种他们需要的合金,以及每种合金中铁铝锡的比重。公司希望能够订购最少种类的原材料,并且使用这些原材料可以加工出用户需要的所有种类的合金。
题解:
好题啊!可以发现,c是没用的,因为c=1-a-b。然后对于一个原材料我们可以看做是平面上的一个点(a,b),对于两个原材料,他们能加工出的点一定在他们连的线段上,而多个的话就是在他们围成的凸包中。对于两个原材料,若所有需要的点都在这条线段的某一边(可以左可以右,但一定要统一),或在线段上,那么连一条边权为1的有向边,最后用floyd求最小环,就是答案。
代码:
#include<bits/stdc++.h>
using namespace std;
const int Maxn=510;
const double eps=1e-8;
struct Point{double x,y;}p[Maxn],P[Maxn];
double multi(Point p1,Point p2,Point p0)
{
double x1=p1.x-p0.x,y1=p1.y-p0.y;
double x2=p2.x-p0.x,y2=p2.y-p0.y;
return x1*y2-x2*y1;
}
int m,n,f[Maxn][Maxn];
bool check(Point p1,Point p2,Point p3)
{
if(abs(multi(p1,p2,p3))<=eps&&
min(p1.x,p2.x)<=p3.x&&p3.x<=max(p1.x,p2.x)&&
min(p1.y,p2.y)<=p3.y&&p3.y<=max(p1.y,p2.y))return true;
return multi(p1,p2,p3)>=eps;
}
int main()
{
memset(f,63,sizeof(f));
scanf("%d%d",&m,&n);
for(int i=1;i<=m;i++)
{
double c;
scanf("%lf%lf%lf",&p[i].x,&p[i].y,&c);
}
for(int i=1;i<=n;i++)
{
double c;
scanf("%lf%lf%lf",&P[i].x,&P[i].y,&c);
}
for(int i=1;i<=m;i++)
for(int j=1;j<=m;j++)
{
bool tf=false;
for(int k=1;k<=n;k++)
if(!check(p[i],p[j],P[k])){tf=true;break;}
if(!tf)f[i][j]=1;
}
for(int k=1;k<=m;k++)
for(int i=1;i<=m;i++)
for(int j=1;j<=m;j++)
if(f[i][k]+f[k][j]<f[i][j])f[i][j]=f[i][k]+f[k][j];
int ans=1061109567;
for(int i=1;i<=m;i++)ans=min(ans,f[i][i]);
if(ans==1061109567)puts("-1");
else printf("%d",ans);
}