bz1027 [JSOI2007]合金

由于本人傻×，给bzoj上贡献了满满一屏的wa，真是对不住了......

这个题其实挺水的，我就简单一说吧。可以看出选的原材料组成的凸包一定要包住需求的材料组成的凸包，那么我们枚举任意两个材料的点，如果他们的向量在所有材料组成的点的同侧的话，就给这两个点连一条1的边，否则连一条正无穷的边，那么用floyd跑一边最小环 就行了。注意特殊情况，如果只有1个点就能满足，那么这个点必须和所有需求的点重合，如果只有2个点就能满足，那么所有需求的点都应在这两个点的连线段上。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#define maxn 1000
#define inf 100000000
#define eps 1e-9
using namespace std;
struct point 
{
    double x,y;
    friend point operator -(const point &a,const point &b)
    {
        point tmp;
        tmp.x=a.x-b.x; tmp.y=a.y-b.y;
        return tmp;
    }
    friend double operator *(const point &a,const point &b)
    {
        return a.x*b.y-a.y*b.x;
    }
}a[maxn],b[maxn];
double z;
int n,m,ans;
int f[maxn][maxn],d[maxn][maxn];

double fabs(double a)
{
    return a>0?a:-a;
}

bool judge(int i,int j)
{
    for (int k=1;k<=n;k++)
        if (((a[j]-a[i])*(b[k]-a[j]))<-eps) return 0;
    return 1;
}

bool single(int i)
{ 
    for (int k=1;k<=n;k++)
        if (a[i].x!=b[k].x||a[i].y!=b[k].y) return 0;
    return 1;
}

bool check(int i,int j)
{
    for (int k=1;k<=n;k++)
    {
        if (b[k].x>max(a[i].x,a[j].x)||b[k].x<min(a[i].x,a[j].x)) return 0;
        if (b[k].y>max(a[i].y,a[j].y)||b[k].y<min(a[i].y,a[j].y)) return 0;
    }
    return 1;
}


int main()
{
    //freopen("metal.in","r",stdin);
    scanf("%d%d",&m,&n);
    for (int i=1;i<=m;i++) scanf("%lf%lf%lf",&a[i].x,&a[i].y,&z);
    for (int i=1;i<=n;i++) scanf("%lf%lf%lf",&b[i].x,&b[i].y,&z);
    for (int i=1;i<=m;i++)
        for (int j=1;j<=m;j++)
        {
            if (i==j) continue;
            if (judge(i,j)) d[i][j]=f[i][j]=1; else d[i][j]=f[i][j]=inf;
        }
    for (int i=1;i<=m;i++)
        if (single(i))
        {
            printf("1\n");
            return 0;
        }
    for (int i=1;i<=m;i++)
        for (int j=1;j<i;j++)
            if (f[i][j]==1&&f[j][i]==1&&check(i,j))
            {
                printf("2\n");
                return 0;
            }
    ans=inf;
    for (int k=1;k<=m;k++)
    {
        for (int i=1;i<k;i++)
            for (int j=1;j<i;j++)
                ans=min(ans,f[i][j]+d[j][k]+d[k][i]);    
        for (int i=1;i<=m;i++)
            for (int j=1;j<=m;j++)
                f[i][j]=min(f[i][j],f[i][k]+f[k][j]);
    }
    if (ans==inf) ans=-1;
    printf("%d\n",ans);
    return 0;
}