【bzoj1027】合金
【bzoj1027】合金
分析
数形结合+计算几何+Floyd最小环。
http://blog.csdn.net/popoqqq/article/details/40539273
虽然这样占大家的很不好,但是我的确需要这么做。
小结
(1)最小环的写法
(2)凸包的一种简易算法
RT。
(3)关于计算几何的练习
通过做这道题,计算几何的模板也熟悉了一些。
计算几何如何调试?
这个么...好像只能瞪眼法了。
代码
#include <cstdio>
#include <climits>
#include <cmath>
#include <algorithm>
using namespace std;
#define rep(i,a,b) for (int i=(a);i<=(b);i++)
const int M=512;
const int N=512;
const int MAX=INT_MAX>>1;
const double EPS=1e-15;
int cmp(double x);
int m,n;
struct point
{
double x,y;
point(double _x=0,double _y=0) {
x=_x,y=_y;
}
friend point operator - (point a,point b) {
return point(a.x-b.x,a.y-b.y);
}
friend point operator * (point a,point b) {
return point(a.x*b.x,a.y*b.y);
}
friend int operator == (point a,point b) {
return !cmp(a.x-b.x)&&!cmp(a.y-b.y);
}
friend int operator != (point a,point b) {
return !(a==b);
}
}mer[M],prd[N];
int f[M][M];
int res;
void Init(void) {
scanf("%d%d",&m,&n);
rep(i,1,m) {
double x,y; scanf("%lf%lf%*lf",&x,&y);
mer[i]=point(x,y);
}
rep(i,1,n) {
double x,y; scanf("%lf%lf%*lf",&x,&y);
prd[i]=point(x,y);
}
}
int cmp(double x)
{
if(fabs(x)<EPS) return 0;
return x>0?1:-1;
}
double det(point a,point b) {
return a.x*b.y-a.y*b.x;
}
int Left(point k,point a,point b) {
double t=det(b-a,k-a);
return t>0;
}
int On(point k,point a,point b) {
if (cmp(det(k-a,b-a))!=0) return 0;
if (!(cmp(min(a.x,b.x)-k.x)<=0&&cmp(k.x-max(a.x,b.x))<=0)) return 0;
if (!(cmp(min(a.y,b.y)-k.y)<=0&&cmp(k.y-max(a.y,b.y))<=0)) return 0;
return 1;
}
/*
bool On(const point &a,const point &b,const point &c) {
return !cmp(det(a-b,c-b))&&cmp((a.x-b.x)*(a.x-c.x))<=0&&cmp((a.y-b.y)*(a.y-c.y))<=0;
}
*/
int Check(point a,point b) {
rep(i,1,n) {
int t1=Left(prd[i],a,b);
int t2=On(prd[i],a,b);
if (!t1&&!t2) return 0;
}
return 1;
}
void MakeGraph(void) {
rep(i,1,m) rep(j,1,m) f[i][j]=MAX;
rep(i,1,m) rep(j,1,m)
if (Check(mer[i],mer[j]))
f[i][j]=1;
}
int Floyd(void) {
rep(i,1,m) {
int mrk=1;
rep(j,1,n) if (mer[i]!=prd[j]) {
mrk=0;
break;
}
if (mrk) return 1;
}
rep(i,1,m) rep(j,1,m) {
int mrk=1;
rep(k,1,n)
if (!On(prd[k],mer[i],mer[j])) {
mrk=0;
break;
}
if (mrk) return 2;
}
/*
rep(i,1,m) {
rep(j,1,m)
if (f[i][j]!=MAX)
printf("%d ",f[i][j]);
else printf("0 ");
printf("\n");
}
*/
int ans=MAX;
rep(k,1,m) rep(i,1,m) rep(j,1,m)
f[i][j]=min(f[i][j],f[i][k]+f[k][j]);
rep(i,1,m) ans=min(ans,f[i][i]);
return ans!=MAX?ans:-1;
}
int main(void) {
#ifndef ONLINE_JUDGE
freopen("bzoj1027.in","r",stdin);
freopen("bzoj1027.out","w",stdout);
#endif
Init();
MakeGraph();
res=Floyd();
printf("%d\n",res);
return 0;
}
【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】
题目 某公司加工一种由铁.铝.锡组成的合金.他们的工作很简单.首先进口一些铁铝锡合金原材料,不同种类的 原材料中铁铝锡的比重不同.然后,将每种原材料取出一定量,经过融解.混合,得到新的合金.新的合金的 ...
- [bzoj 1027][JSOI2007]合金(解析几何+最小环)
题目:http://www.lydsy.com:808/JudgeOnline/problem.php?id=1027 分析: 首先因为一个合金的和为1,所以考虑2个材料合金能否合成一个需求合金的时候 ...
- 合金装备V 幻痛 制作技术特辑
合金装备V:幻痛 制作特辑 资料原文出自日版CGWORLD2015年10月号 在[合金装备4(Metal Gear Solid IV)]7年后,序章作品[合金装备5 :原爆点 (Metal Gea ...
- 【BZOJ】【1027】【JSOI2007】合金
计算几何/凸包/Floyd Orz rausen大爷太强辣 计算几何题目果然不会做>_> 这个题……虽然他给了3个坐标,但实际上是个二维的计算几何题= =因为第三维坐标可以直接用前两维坐标 ...
- 1027: [JSOI2007]合金 - BZOJ
Description 某公司加工一种由铁.铝.锡组成的合金.他们的工作很简单.首先进口一些铁铝锡合金原材料,不同种类的原材料中铁铝锡的比重不同.然后,将每种原材料取出一定量,经过融解.混合,得到新的 ...
- BZOJ_1027_[JSOI2007]_合金_(计算几何+Floyd求最小环)
描述 http://www.lydsy.com/JudgeOnline/problem.php?id=1027 共三种金属,\(m\)种材料,给出每种材料中三种金属的占比. 给出\(n\)种合金的三种 ...
随机推荐
- winform 发邮件
System.Net.Mail.SmtpClient client = new SmtpClient("smtp.163.com"); client.UseDefaultCre ...
- 抓取Js动态生成数据且以滚动页面方式分页的网页
代码也可以从我的开源项目HtmlExtractor中获取. 当我们在进行数据抓取的时候,如果目标网站是以Js的方式动态生成数据且以滚动页面的方式进行分页,那么我们该如何抓取呢? 如类似今日头条这样的网 ...
- SqlSever基础 lower函数 返回字符串的小写形式
镇场诗:---大梦谁觉,水月中建博客.百千磨难,才知世事无常.---今持佛语,技术无量愿学.愿尽所学,铸一良心博客.------------------------------------------ ...
- activiti学习总结
Activiti界面元素的使用总结 一.图形设计中元素的使用 1.SequenceFlow:连接线,可以连接两个任务,来管理流程实例的流向 -----General -----id:流程的id,用与程 ...
- BASE64,MD5,SHA,HMAC加密與解密算法(java)
package com.ice.webos.util.security; import java.io.UnsupportedEncodingException; import java.math.B ...
- 我的android学习经历15
利用Intent实现有返回结果的页面跳转 主要用的方法: (1)Intent的构造方法:intent(当前界面对象,要跳转的界面.class); (2)接受结果的方法onActivityResult( ...
- BZOJ 3514 Codechef MARCH14 GERALD07加强版
题目链接:http://www.lydsy.com:808/JudgeOnline/problem.php?id=3514 题意:给出一个图m条边.每次询问只加入编号在区间[L,R]之内的边有多少连通 ...
- Using Pre-Form Trigger In Oracle Forms
Pre-Form trigger in Oracle Forms fires during the form start-up, before forms navigates to the first ...
- Hibernate 查询语句
Hibernate查询语言(HQL)是一种面向对象的查询语言,类似于SQL,但不是对表和列操作,HQL适用于持久对象和它们的属性. HQL查询由Hibernate转换成传统的SQL查询,这在圈上的数据 ...
- 07.常用的SQL语句
常用SQL语句,仿照着写(只改动字段和表名,库名) 查看数据库引擎 SHOW ENGINES;(默认InnoDB) 用数据库生成UUID selec ...