BZOJ 1027-合金（凸包+最小环）

1027: [JSOI2007]合金

Time Limit: 4 Sec   Memory Limit: 162 MB
Submit: 4064   Solved: 1191
[ Submit][ Status][ Discuss]

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

HINT

Source

刚开始写着到题想着二分答案随便搞搞，一搞就是几个小时，后来惊觉自己就是个sb。。。。

这题确实是好题，不看题解我绝不会往集合空间上想，更不会想到会和凸包有联系（虽然我当时不会凸包）

不过这道题并没有用到凸包的算法，但是我觉得学一下会更好：凸包的五种解法

既然前边提到要把这道题转化到空间中还是跟凸包有关，呢肯定是搞到二维平面上。

但是题目中明明是三个变量，不要忘了x+y+z==1，故，第三维没啥意义。。。。

为什么说会跟凸包有关呢：考虑一下这些二维平面上的点的意义……

假如我们现在手里有一种原材料 (x1,y1) 那么我们就只能搞出这一种合金=

那么如果有两种原材料，我们配一配，就可以配出 (x1,y1) 到 (x2,y2) 这条线段上的所有合金

如果有三种，我们能配出的就是那个三角形区域的所有种类的合金……

以此类推……我们有n种原材料，我们能搞出的合金就是这n个点的凸包围出来的区域。

floyd求最小环。。。具体就是如果目标点在线段a,b的一侧则dis[a][b]=1

否则dis[a][b]=inf，答案是min{dis[i][i]}

知识拓展：

                                                 判断 某一点在直线左右侧

          

      左右方向是相对前进方向的,只要指定了前进方向就可以知道左右(比如指定前进方向是从直线的起点到终点).

        

      判断点在直线的左侧还是右侧是计算几何里面的一个最基本算法.使用矢量来判断. 

      定义：平面上的三点P1(x1,y1),P2(x2,y2),P3(x3,y3)的面积量：

      S(P1,P2,P3)=|y1 y2 y3|= (x1-x3)*(y2-y3)-(y1-y3)*(x2-x3) 

     当P1P2P3逆时针时S为正的，当P1P2P3顺时针时S为负的。 

     令矢量的起点为A，终点为B，判断的点为C， 

      如果S（A，B，C）为正数，则C在矢量AB的左侧； 

      如果S（A，B，C）为负数，则C在矢量AB的右侧； 

      如果S（A，B，C）为0，则C在直线AB上。

     

      具体判断方法：点与直线的位置关系

#include<map>
#include<stack>
#include<queue>
#include<vector>
#include<math.h>
#include<stdio.h>
#include<iostream>
#include<string.h>
#include<stdlib.h>
#include<algorithm>
using namespace std;
typedef long long  ll;
#define inf 10000000
#define  maxn  505
#define  eps 1e-10
int d[maxn][maxn],n,m;
struct node
{
    double x,y;
}a[maxn],b[maxn];
bool tepan()//所有点重合要特判掉
{
    int i;
    for(i=1;i<=n;i++)
        if(fabs(a[i].x-a[1].x)>eps || fabs(a[i].y-a[1].y)>eps)
            return 0;
    for(i=1;i<=m;i++)
        if(fabs(b[i].x-a[1].x)>eps || fabs(b[i].y-a[1].y)>eps)
            return 0;
    printf("1\n");
    return 1;
}
double work(node a,node b, node c)
{
    return (a.x-c.x)*(b.y-c.y)-(a.y-c.y)*(b.x-c.x);
}
bool judge(node x,node y)//判断呢些点是不是在这两点构成的线段上
{
    if(x.x>y.x)
        swap(x,y);
    for(int i=1;i<=m;i++)
        if(b[i].x<x.x || b[i].x>y.x)
            return 0;
    if(x.y>y.y)
        swap(x,y);
    for(int i=1;i<=m;i++)
        if(b[i].y<x.y || b[i].y>y.y)
            return 0;
    printf("2\n");
    return 1;
}
int check(node x,node y)
{
    int i,t1=0,t2=0;
    for(i=1;i<=m;i++)
    {
        double tmp=work(x,y,b[i]);
        if(tmp>eps)//点在直线左侧
            t1++;
        if(tmp<-eps)//点在直线右侧
            t2++;
        if(t1 && t2)//该直线上这两点不是凸包上的点
            return 0;
    }
    if(t1)
        return 1;
    if(t2)
        return 2;
    if(judge(x,y))
        return -1;
    return 3;
}
void floyd()
{
    int i,j,k,ans=inf;
    for(k=1;k<=n;k++)
        for(i=1;i<=n;i++)
            for(j=1;j<=n;j++)
                d[i][j]=min(d[i][j],d[i][k]+d[k][j]);
    for(i=1;i<=n;i++)
        ans=min(ans,d[i][i]);
    if(ans>=inf || ans<=2)
        printf("-1\n");
    else
        printf("%d\n",ans);
}
void solve()
{
    int i,j;
    for(i=1;i<=n;i++)
        for(j=i+1;j<=n;j++)
        {
            int mark=check(a[i],a[j]);
            if(mark==-1)
                return;
            if(mark==1)
                d[i][j]=1;
            if(mark==2)
                d[j][i]=1;
            if(mark==3)
                d[i][j]=d[j][i]=1;
        }
    floyd();
}
int main(void)
{
    int i,j;
    double xx;
    scanf("%d%d",&n,&m);
    for(i=1;i<=n;i++)
        for(j=1;j<=n;j++)
            d[i][j]=inf;
    for(i=1;i<=n;i++)
        scanf("%lf%lf%lf",&a[i].x,&a[i].y,&xx);//xx可忽略
    for(i=1;i<=m;i++)
        scanf("%lf%lf%lf",&b[i].x,&b[i].y,&xx);
    if(tepan())
        return 0;
    solve();
    return 0;
}