bzoj1027【JSOI2007】合金
题目描述
某公司加工一种由铁、铝、锡组成的合金。他们的工作很简单。首先进口一些铁铝锡合金原材料,不同种类的
原材料中铁铝锡的比重不同。然后,将每种原材料取出一定量,经过融解、混合,得到新的合金。新的合金的铁铝
锡比重为用户所需要的比重。 现在,用户给出了n种他们需要的合金,以及每种合金中铁铝锡的比重。公司希望能
够订购最少种类的原材料,并且使用这些原材料可以加工出用户需要的所有种类的合金。
输入格式
第一行两个整数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),分别表示铁铝锡在一种用户需要的合金中
所占的比重。
输出格式
一个整数,表示最少需要的原材料种数。若无解,则输出–1。
题解
- $a+b+c=1$所以忽略$c$,根据定比分点原则两个点可以表示的点组成它们连线段,推广一下;
- 所以对$n$个点求出凸包,问题即用$m$个点去圈住凸包求最小点数;
- 对两个材料$i,j$,如果凸包上的点都在连线$\vec{ij}$的左边就$i$向$j$连边,$floyd$求所有$dis[i][i]$即可;
- 特判:
- 当凸包是一个点时特判;
- 当凸包是一条直线时,如果叉积为$0$还需要盖住所有点才可以加边,用点积特判;
#include<bits/stdc++.h>
#define db double
#define inf 0x3f3f3f3f
#define il inline
#define eps 1e-7
using namespace std;
const int N=;
int n,m,top,dis[N][N];
il int dcmp(db x){return fabs(x)<eps?:x<?-:;}
struct poi{
db x,y;
poi(db _x=,db _y=):x(_x),y(_y){};
il poi operator -(const poi&A)const{return poi(x-A.x,y-A.y);}
il bool operator <(const poi&A)const{return x==A.x?y<A.y:x<A.x;}
}p1[N],p2[N],q[N];
il db crs(poi A,poi B){return A.x*B.y-A.y*B.x;}
il db dot(poi A,poi B){return A.x*B.x+A.y*B.y;}
il bool spj(){
int fg=;
for(int i=;i<=n;++i)if(dcmp(p1[].x-p1[i].x)||dcmp(p1[].y-p1[i].y)){fg=;break;}
for(int i=;i<=m;++i)if(dcmp(p1[].x-p2[i].x)||dcmp(p1[].y-p2[i].y)){fg=;break;}
if(!fg){puts("");return true;}
else return false;
}
il void convex(){
sort(p2+,p2+m+);
q[top=]=p2[];
if(m==)return;
for(int i=;i<=m;++i){
while(top>&&dcmp(crs(q[top]-q[top-],p2[i]-q[top]))<=)top--;
q[++top]=p2[i];
}
int now=top;
for(int i=m-;i;--i){
while(top>now&&dcmp(crs(q[top]-q[top-],p2[i]-q[top]))<=)top--;
q[++top]=p2[i];
}
top--;
}
il bool judge(int a,int b){
int i; poi p = p1[b]-p1[a];
for(i=;i<=top;++i){
int c = dcmp(crs(p, q[i]-p1[a]));
if(c>)break;
if(!c&&dcmp(dot(q[i]-p1[b], q[i]-p1[a]))>)break;
}
return i==top+?true:false;
}
int main(){
#ifndef ONLINE_JUDGE
freopen("bzoj1027.in","r",stdin);
freopen("bzoj1027.out","w",stdout);
#endif
db tmp;
scanf("%d%d",&n,&m);
for(int i=;i<=n;++i)scanf("%lf%lf%lf",&p1[i].x,&p1[i].y,&tmp);
for(int i=;i<=m;++i)scanf("%lf%lf%lf",&p2[i].x,&p2[i].y,&tmp);
if(spj())return ;
convex();
memset(dis,0x3f,sizeof(dis));
for(int i=;i<=n;++i)
for(int j=;j<=n;++j)if(i!=j){
if(judge(i,j))dis[i][j]=;
}
for(int k=;k<=n;++k)
for(int i=;i<=n;++i)
for(int j=;j<=n;++j){
if(dis[i][j]>dis[i][k]+dis[k][j])
dis[i][j] = dis[i][k]+dis[k][j];
}
int ans=inf;
for(int i=;i<=n;++i)ans = min(ans, dis[i][i]);
if(ans==inf)puts("-1"); else printf("%d\n",ans);
return ;
}bzoj1027
bzoj1027【JSOI2007】合金的更多相关文章
- bzoj1027 [JSOI2007]合金
1027: [JSOI2007]合金 Time Limit: 4 Sec Memory Limit: 162 MBSubmit: 2671 Solved: 703[Submit][Status][ ...
- bzoj千题计划123:bzoj1027: [JSOI2007]合金
http://www.lydsy.com/JudgeOnline/problem.php?id=1027 因为x+y+z=1,所以z=1-x-y 第三维可以忽略 将x,y 看做 平面上的点 简化问题: ...
- BZOJ1027 [JSOI2007]合金 【计算几何 + floyd】
题目 某公司加工一种由铁.铝.锡组成的合金.他们的工作很简单.首先进口一些铁铝锡合金原材料,不同种类的 原材料中铁铝锡的比重不同.然后,将每种原材料取出一定量,经过融解.混合,得到新的合金.新的合金的 ...
- 【bzoj1027】合金
[bzoj1027]合金 分析 数形结合+计算几何+Floyd最小环. http://blog.csdn.net/popoqqq/article/details/40539273 虽然这样占大家的很不 ...
- bzoj 1027 [JSOI2007]合金(计算几何+floyd最小环)
1027: [JSOI2007]合金 Time Limit: 4 Sec Memory Limit: 162 MBSubmit: 2970 Solved: 787[Submit][Status][ ...
- BZOJ 1027 [JSOI2007]合金
1027: [JSOI2007]合金 Time Limit: 4 Sec Memory Limit: 162 MBSubmit: 2605 Solved: 692[Submit][Status][ ...
- [bzoj 1027][JSOI2007]合金(解析几何+最小环)
题目:http://www.lydsy.com:808/JudgeOnline/problem.php?id=1027 分析: 首先因为一个合金的和为1,所以考虑2个材料合金能否合成一个需求合金的时候 ...
- 1027: [JSOI2007]合金 - BZOJ
Description 某公司加工一种由铁.铝.锡组成的合金.他们的工作很简单.首先进口一些铁铝锡合金原材料,不同种类的原材料中铁铝锡的比重不同.然后,将每种原材料取出一定量,经过融解.混合,得到新的 ...
- [JSOI2007]合金
Description 某公司加工一种由铁.铝.锡组成的合金.他们的工作很简单.首先进口一些铁铝锡合金原材料,不同种类的 原材料中铁铝锡的比重不同.然后,将每种原材料取出一定量,经过融解.混合,得到新 ...
- P4049 [JSOI2007]合金
传送门 我数学可能白学了-- 因为三个数加起来等于\(1\),那么只要用前两个数就能表示,那么就能把每一种金属看成一个二维向量.考虑只有两个向量的时候,设这两个向量为\(a,b\),那么一个向量\(c ...
随机推荐
- openstack系列文章(一)
学习openstack的系列文章-虚拟化 虚拟化 KVM CPU 虚拟化 KVM 内存虚拟化 全虚拟化 I/O 设备 半虚拟化 I/O 设备 I/O PCI PCIe 设备直接分配 SR-IOV 在 ...
- BOM 头是什么,怎么除去
WINDOWS自带的记事本,在保存一个以 UTF-8 编码的文件时,会在文件开始的地方插入三个不可见的字符(0xEF 0xBB 0xBF,即BOM).它是一串隐藏的字符,用于让记事本等编辑器识别这个文 ...
- nodejs 搭建自己的简易缓存cache管理模块
http://www.infoq.com/cn/articles/built-cache-management-module-in-nodejs/ 为什么要搭建自己的缓存管理模块? 这个问题其实也是在 ...
- ES6的新特性(23)——ArrayBuffer
ArrayBuffer ArrayBuffer对象.TypedArray视图和DataView视图是 JavaScript 操作二进制数据的一个接口.这些对象早就存在,属于独立的规格(2011 年 2 ...
- 浏览器播放rtmp流
我是利用flash插件实现的,需要以下几个文件: flowplayer-3.2.8.min.js flowplayer-3.2.18.swf flowplayer.rtmp-3.2.8.swf flo ...
- TCP系列52—拥塞控制—15、前向重传与RACK重传拥塞控制处理对比
一.概述 这里主要简单分析一个丢包重传并恢复的场景,通过不同的设置让这个相同的场景分别触发RACK重传和前向重传,通过对比说明以下问题: Forward Retransmit可以产生只有重传标记的数据 ...
- poi excel导入 数字自动加小数点
问题:导入excel表,若表格中为整数数字,不管单元格设置成数字格式还是文本格式,导入时都会出现小数点和0. 我遇到的问题是:一个名称,做测试数据的时候做了纯整形数字,发现了这个问题. 解决办法:在代 ...
- loadrunner在win10破解提示:Cannot save the license information because acceses to the registry is denied的解决办法
方法1 下图1-1中提示就是说明了破解的时候权限不足导致,解决办法就是使用管理员权限打开loadrunner破解就好了,但是右键“以管理员身份运行”选项打开loadrunner又是会提示1-2中的问题 ...
- springmvc+json 前后台数据交互
1. 配置(1) 文件配置参考这里(2) 导入jackson相关包:jackson-annotations-2.9.4.jar,jackson-core-2.9.4.jar,jackson-datab ...
- 随机场(Random field)
一.随机场定义 http://zh.wikipedia.org/zh-cn/随机场 随机场(Random field)定义如下: 在概率论中, 由样本空间Ω = {0, 1, …, G − 1}n取样 ...