1
0
Fork 0

Add number/fractional.cc

Signed-off-by: subcrip <contact@subc.rip>
This commit is contained in:
subcrip 2024-11-16 07:48:23 +00:00
parent 367c0bef37
commit cb03157987
1 changed files with 31 additions and 0 deletions

31
number/fractional.cc Normal file
View File

@ -0,0 +1,31 @@
template <typename T>
inline T mygcd(T a, T b) { return b == 0 ? a : mygcd(b, a % b); }
template<typename T>
struct fractional {
T p, q;
inline void reduce(void) {
if (q < 0) p = -p, q = -q;
if (p == 0) q = 1; else { T g = mygcd(p, q); p /= g; q /= g; }
}
fractional(void) : p(0), q(1) {}
template <typename U>
fractional(const U& p) : p(p), q(1) { reduce(); }
fractional(const T& p, const T& q) : p(p), q(q) { reduce(); }
inline fractional operator+(void) const { return *this; }
inline fractional operator-(void) const { return { -p, q }; }
inline fractional operator+(const fractional& rhs) const { return { p * rhs.q + q * rhs.p, q * rhs.q }; }
inline fractional operator-(const fractional& rhs) const { return *this + (-rhs); }
inline fractional operator*(const fractional& rhs) const { return { p * rhs.p, q * rhs.q }; }
inline fractional operator/(const fractional& rhs) const { return *this * fractional(rhs.q, rhs.p); }
inline fractional& operator+=(const fractional& rhs) { return *this = *this + rhs; }
inline fractional& operator-=(const fractional& rhs) { return *this = *this - rhs; }
inline fractional& operator*=(const fractional& rhs) { return *this = *this * rhs; }
inline fractional& operator/=(const fractional& rhs) { return *this = *this / rhs; }
inline bool operator==(const fractional& rhs) const { return p == rhs.p and q == rhs.q; }
inline bool operator!=(const fractional& rhs) const { return not (*this == rhs); }
inline bool operator<(const fractional& rhs) const { return (*this - rhs).p < 0; }
inline bool operator>=(const fractional& rhs) const { return not (*this < rhs); }
inline bool operator>(const fractional& rhs) const { return *this >= rhs and *this != rhs; }
inline bool operator<=(const fractional& rhs) const { return *this < rhs or *this == rhs; }
};