码迷,mamicode.com
首页 > Web开发 > 详细

[bzoj1027][JSOI2007]合金

时间:2017-04-03 21:59:49      阅读:216      评论:0      收藏:0      [点我收藏+]

标签:题目   getchar   for   博客   zoj   允许   print   复杂   jsoi2007   

来自FallDream的博客,未经允许,请勿转载,谢谢。


 

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

n,m<=500

看到这道神题根本没头绪,只好看看题解。

因为三个参数加起来等于1,所以我们可以忽略第三个,用前两维表示一个平面上的点,两个材料能混合出来的材料是这两个点之间的线段,所以题目转换为求尽量少的点,使得这些点的凸包可以包住需要的材料的凸包。具体实现,如果所有需要的材料的点都在一条线段左边,那么这两个点之间连边,然后floyd求最小环。复杂度n^3

注意各种特判。如果全部都重点输出1,如果一条直线覆盖了所有的点判断一下是否都在线段上,是就输出2。最后要判断答案必须大于2.

#include<iostream>
#include<cstring>
#include<cstdio>
#define INF 1000000000
#define MN 500
#define eps 1e-8
using namespace std;
inline int read()
{
    int x = 0 , f = 1; 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;
}
 
int f[MN+5][MN+5];
int n,m;double Xmn=INF,Xmx=-INF,Ymn=INF,Ymx=-INF;
struct P{double x,y;
    P(double x=0,double y=0):x(x),y(y){}
    double operator^(P b){return x*b.y-y*b.x;}
    P operator-(P b){return P(x-b.x,y-b.y);}
}p[MN+5],q[MN+5];
 
bool check(double x1,double x2,double x3,double x4)
{
    if(x1>x2)swap(x1,x2);
    return x1<=x3&&x4<=x2;
}
 
int solve(int x,int y)
{
    int num1=0,num2=0;
    for(int i=1;i<=m;i++)
    {
        double crs=(p[y]-p[x])^(q[i]-p[x]);
        if(crs<-eps)num2++;
        if(crs>eps) num1++;
    }
    if(!num1&&!num2&&check(p[x].x,p[y].x,Xmn,Xmx)&&check(p[x].y,p[y].y,Ymn,Ymx))return -1;
    if(num1*num2)return 0;
    if(num1)return 2;
    if(num2)return 1;
    return 3;
}
 
bool special()
{
    for(int i=2;i<=n;i++)
        if(p[i].x!=p[1].x||p[i].y!=p[1].y)return 0;
    for(int i=1;i<=m;i++)
        if(q[i].x!=p[1].x||q[i].y!=p[1].y)return 0;
    return 1;
}
 
int main()
{
    n=read();m=read();double d;
    for(int i=1;i<=n;i++)
        scanf("%lf%lf%lf",&p[i].x,&p[i].y,&d);
    for(int i=1;i<=m;i++)
    {
        scanf("%lf%lf%lf",&q[i].x,&q[i].y,&d);
        Xmn=min(Xmn,q[i].x);Xmx=max(Xmx,q[i].x);
        Ymn=min(Ymn,q[i].y);Ymx=max(Ymx,q[i].y);
    }
    if(special())return 0*puts("1");
    memset(f,63,sizeof(f));
    for(int i=1;i<n;i++)
        for(int j=i+1;j<=n;j++)
        {
            int res=solve(i,j);
            if(res==-1)return 0*puts("2");
            if(res&1)f[i][j]=1;
            if(res&2)f[j][i]=1;
        }
    for(int k=1;k<=n;k++)
        for(int i=1;i<=n;i++)
             if(f[i][k]<=INF)
                for(int j=1;j<=n;j++)
                    f[i][j]=min(f[i][j],f[i][k]+f[k][j]);
    int minn=INF;
    for(int i=1;i<=n;i++) minn=min(minn,f[i][i]);
    if(minn>2&&minn<INF)printf("%d\n",minn);else puts("-1");
    return 0;
}

 

[bzoj1027][JSOI2007]合金

标签:题目   getchar   for   博客   zoj   允许   print   复杂   jsoi2007   

原文地址:http://www.cnblogs.com/FallDream/p/bzoj1027.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!