src/geo2d/dis/pp/lib.hpp
| Line | Branch | Exec | Source |
|---|---|---|---|
| 1 | #pragma once | ||
| 2 | |||
| 3 | #include "../../ds/p/lib.hpp" | ||
| 4 | |||
| 5 | namespace tifa_libs { | ||
| 6 | |||
| 7 | // distance of two points | ||
| 8 | // @return $(|x_1-x_2|^p + |y_1-y_2|^p)^{1/p}$, p = 0 means $\infty$ | ||
| 9 | template <class FP, u32 p = 2> | ||
| 10 | 80908993 | CEXP FP dist_PP(point<FP> CR p1, point<FP> CR p2) NE { | |
| 11 | static_assert(p < 2 || std::floating_point<FP>); | ||
| 12 | 231 | if CEXP (p == 0) return max(abs(p1.x - p2.x), abs(p1.y - p2.y)); // Chebyshev | |
| 13 | 471 | else if CEXP (p == 1) return abs(p1.x - p2.x) + abs(p1.y - p2.y); // Manhattan | |
| 14 | 80906471 | else if CEXP (p == 2) return (p1 - p2).norm(); // Euclidian | |
| 15 | 1820 | else return std::pow(std::pow(abs(p1.x - p2.x), p) + std::pow(abs(p1.y - p2.y), p), FP{1} / p); | |
| 16 | } | ||
| 17 | |||
| 18 | } // namespace tifa_libs | ||
| 19 |