-
Notifications
You must be signed in to change notification settings - Fork 1
/
base_ld.hpp
54 lines (40 loc) · 1.08 KB
/
base_ld.hpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
#pragma once
#include "../template/template.hpp"
namespace lib {
using vec = complex<ld>;
const ld PI = acos(-1);
void ldout(int len = 20) {
cout << fixed << setprecision(len);
}
int sgn(ld a, const ld eps = 1e-7) {
return (a < -eps) ? -1 : (a > eps) ? 1 : 0;
}
bool same_vec(vec a, vec b) {
a -= b;
return sgn(a.real()) == 0 && sgn(a.imag()) == 0;
}
ld dot(const vec &a, const vec &b) {
return (conj(a) * b).real();
}
ld cross(const vec &a, const vec &b) {
return (conj(a) * b).imag();
}
int isp(const vec &a, const vec &b, const vec &c) {
int cross_sgn = sgn(cross(b - a, c - a));
if (cross_sgn == 0) {
if (sgn(dot(b - a, c - a)) < 0) return -2;
if (sgn(dot(a - b, c - b)) < 0) return 2;
}
return cross_sgn;
}
vec rot90(const vec &a) {
return {-a.imag(), a.real()};
}
vec rot(const vec &a, ld rad) {
return a * vec(cosl(rad), sinl(rad));
}
bool comp_for_argument_sort(const vec &lhs, const vec &rhs) {
// if (abs(arg(lhs)-arg(rhs)) < eps) return false; // need ?
return arg(lhs) < arg(rhs);
}
} // namespace lib