BZOJ 1027: [JSOI2007]合金 (计算几何+Floyd求最小环)
同样的方法…我200ms,他2000ms.
常数的幽怨…
CODE
#include <bits/stdc++.h>
using namespace std;
const int MAXN = 505;
const double eps = 1e-8;
struct Point {
double x, y;
Point(){}
Point(double x, double y):x(x), y(y){}
inline Point operator -(const Point &o)const { return Point(x-o.x, y-o.y); }
inline double operator *(const Point &o)const { return x*o.y - y*o.x; }
}a[MAXN], b[MAXN];
struct Line {
Point p, v;
Line(){}
Line(const Point &p, const Point &v):p(p), v(v){}
};
inline bool On_Right(const Line &l, const Point &p) {
return (p - l.p) * l.v > eps;
}
inline int dcmp(const double &x) {
return x < -eps ? -1 : x < eps ? 0 : 1;
}
int m, n, f[MAXN][MAXN];
int main () {
scanf("%d%d", &m, &n); double o;
for(int i = 1; i <= m; ++i) scanf("%lf%lf%lf", &a[i].x, &a[i].y, &o);
for(int i = 1; i <= n; ++i) scanf("%lf%lf%lf", &b[i].x, &b[i].y, &o);
for(int i = 1; i <= m; ++i) {
for(int j = 1; j <= m; ++j)
f[i][j] = m+1;
bool flg = 1;
for(int j = 1; j <= n && flg; ++j)
if(dcmp(b[j].x-a[i].x) || dcmp(b[j].y-a[i].y)) flg = 0;
if(flg) return puts("1"), 0;
}
for(int i = 1; i <= m; ++i)
for(int j = 1; j <= m; ++j)
if(i != j && !(a[i].x == a[j].x && a[i].y == a[j].y)) { //i->j
Line tmp = Line(a[i], a[j] - a[i]);
bool flg = 1;
for(int k = 1; k <= n && flg; ++k)
if(On_Right(tmp, b[k])) flg = 0;
if(flg) {
for(int k = 1; k <= n && flg; ++k)
if
(
(!dcmp((b[k] - tmp.p) * tmp.v))
&&
((b[k].x < a[i].x && b[k].x < a[j].x)
|| (b[k].x > a[i].x && b[k].x > a[j].x)
|| (b[k].y < a[i].y && b[k].y < a[j].y)
|| (b[k].y > a[i].y && b[k].y > a[j].y)
)
) flg = 0;
if(flg) f[i][j] = 1;
}
}
for(int k = 1; k <= m; ++k)
for(int i = 1; i <= m; ++i) if(f[i][k] <= m)
for(int j = 1; j <= m; ++j) if(f[k][j]+f[i][k] <= m)
f[i][j] = min(f[i][j], f[i][k] + f[k][j]);
int ans = m+1;
for(int i = 1; i <= m; ans = min(ans, f[i][i]), ++i);
printf("%d\n", ans > m ? -1 : ans);
}
BZOJ 1027: [JSOI2007]合金 (计算几何+Floyd求最小环)的更多相关文章
- BZOJ 1027 JSOI2007 合金 计算几何+Floyd
题目大意:给定一些合金,选择最少的合金,使这些合金能够按比例合成要求的合金 首先这题的想法特别奇异 看这题干怎么会想到计算几何 并且计算几何又怎么会跟Floyd挂边 好强大 首先因为a+b+c=1 所 ...
- 【BZOJ 1027】 (凸包+floyd求最小环)
[题意] 某公司加工一种由铁.铝.锡组成的合金.他们的工作很简单.首先进口一些铁铝锡合金原材料,不同种类的原材料中铁铝锡的比重不同.然后,将每种原材料取出一定量,经过融解.混合,得到新的合金.新的合金 ...
- BZOJ 1027 [JSOI2007]合金 ——计算几何
我们可以把每一种金属拆成一个二维向量,显然第三维可以计算出来,是无关的. 我们只需要考虑前两维的情况,显然可以构成点集所形成的凸包内. 然后我们枚举两两的情况,然后可以发现如果所有的点都在一侧是可以选 ...
- bzoj 1027 [JSOI2007]合金(计算几何+floyd最小环)
1027: [JSOI2007]合金 Time Limit: 4 Sec Memory Limit: 162 MBSubmit: 2970 Solved: 787[Submit][Status][ ...
- BZOJ_1027_[JSOI2007]_合金_(计算几何+Floyd求最小环)
描述 http://www.lydsy.com/JudgeOnline/problem.php?id=1027 共三种金属,\(m\)种材料,给出每种材料中三种金属的占比. 给出\(n\)种合金的三种 ...
- bzoj 1027: [JSOI2007]合金【凸包+Floyd】
参考:https://www.cnblogs.com/zhuohan123/p/3237246.html 因为一c可以由1-a-b得出,所以删掉c,把a,b抽象成二维平面上的点.首先考虑一个客户需求能 ...
- 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 某公司加工一种由铁.铝.锡组成的合金.他们的工作很简单.首先进口一些铁铝锡合金原材料,不同种类的原材料中铁铝锡的比重不同.然后,将每种原材料取出一定量,经过融解.混合,得到新的 ...
随机推荐
- [Cometoj#4 C]方块切割_质因数分解_贪心
方块切割 题目链接:https://cometoj.com/contest/39/problem/C?problem_id=1583 数据范围:略. 题解: 首先,如果我们知道了多少道在行上,多少刀在 ...
- spark 执行报错 java.io.EOFException: Premature EOF from inputStream
使用spark2.4跟spark2.3 做替代公司现有的hive选项. 跑个别任务spark有以下错误 java.io.EOFException: Premature EOF from inputSt ...
- Python 命名规范总结
Python推荐命名规范: 模块名和包名采用小写字母并且以下划线分隔单词的形式: 如:browser_driver 类名或异常名采用每个单词首字母大写的方式: 如:BasePage, Keyboard ...
- java中public protected friendly private作用域
1.public:public表明该数据成员.成员函数是对所有用户开放的,所有用户都可以直接进行调用 2.private:private表示私有,私有的意思就是除了class自己之外,任何人都不可以直 ...
- MySql设计表中的create_time和update_time字段
一般create_time和update_time字段的类型为datetime类型,长度为0
- [转]Life of a binary
发现一篇讲程序生命周期的文章,感觉蛮不错.
- 深入理解C++11 C3
继承构造函数 class A { public: A(int i):m_i(i) {} A(double d, int i):m_d(d),m_i(i){} private: int m_i{0}; ...
- 牛客 72C 小H和游戏 (动态点分治)
大意: 给定树, 每个点初始权值0, 每次询问给出$x$, $x$权值+1, 求距离$x$不超过2的权值和. 这题数据范围过大, 动态点分治卡不过去, 考虑其他做法 考虑每次只加范围$1$, c[0] ...
- IntelliJ IDEA 统一设置编码为utf-8编码
问题一: File->Settings->Editor->File Encodings 问题二: File->Other Settings->Default Settin ...
- el-table 单元格样式修改
<el-table :cell-style="set_cell_style"> set_cell_style({row, column, rowIndex, colum ...