BZOJ 1027 合金

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

---

【分析】
计算几何
参考 http://blog.csdn.net/popoqqq/article/details/40539273
关键在于叉积判断点在直线的哪一侧与floyd求最小环
主函数里有个很迷的地方，我加了if(i!=j) 就错掉了，不加就对了…这精度卡的有点丧啊…

---

【代码】

//BZOJ 1027 合金
#include<iostream>
#include<cstring>
#include<cstdio>
#define eps 1e-6
#define inf 1e9+7
#define M(a) memset(a,0,sizeof a)
#define fo(i,j,k) for(i=j;i<=k;i++)
using namespace std;
const int mxn=505;
int n,m,ans=inf;
int f[mxn][mxn];
double x[mxn],y[mxn],a[mxn],b[mxn],l,r;
inline int floyd()
{
    int i,j,k;
    fo(k,1,m)
      fo(i,1,m)
        if(f[i][k]<inf)
          fo(j,1,m)
            f[i][j]=min(f[i][j],f[i][k]+f[k][j]);
    fo(i,1,m) ans=min(ans,f[i][i]);
    if(ans==inf) return -1;
    return ans;
}
int main()
{
    int i,j,k;
    scanf("%d%d",&m,&n);
    fo(i,1,m) scanf("%lf%lf%lf",&x[i],&y[i],&l);
    fo(i,1,n) scanf("%lf%lf%lf",&a[i],&b[i],&r);
    fo(i,1,m) fo(j,1,m) f[i][j]=inf;
    fo(i,1,m)
      fo(j,1,m)
      {
          bool flag=1;
          fo(k,1,n)
          {
              double x1=a[k]-x[i],y1=b[k]-y[i];
              double x2=a[k]-x[j],y2=b[k]-y[j];
              double tmp=x1*y2-x2*y1; //叉积 
              if(tmp>=eps) {flag=0;break;}
              if(tmp>-eps && tmp<eps && x1*x2+y1*y2>=eps) {flag=0;break;} 
          }
          if(flag) f[i][j]=1;
      }
    printf("%d\n",floyd());
    return 0;
}