結果
| 問題 |
No.1936 Rational Approximation
|
| コンテスト | |
| ユーザー |
emthrm
|
| 提出日時 | 2022-05-13 22:03:07 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 15 ms / 2,000 ms |
| コード長 | 4,652 bytes |
| コンパイル時間 | 1,801 ms |
| コンパイル使用メモリ | 192,792 KB |
| 最終ジャッジ日時 | 2025-01-29 07:00:55 |
|
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 1 |
| other | AC * 14 |
ソースコード
#define _USE_MATH_DEFINES
#include <bits/stdc++.h>
using namespace std;
#define FOR(i,m,n) for(int i=(m);i<(n);++i)
#define REP(i,n) FOR(i,0,n)
#define ALL(v) (v).begin(),(v).end()
using ll = long long;
constexpr int INF = 0x3f3f3f3f;
constexpr long long LINF = 0x3f3f3f3f3f3f3f3fLL;
constexpr double EPS = 1e-8;
constexpr int MOD = 1000000007;
// constexpr int MOD = 998244353;
constexpr int DY4[]{1, 0, -1, 0}, DX4[]{0, -1, 0, 1};
constexpr int DY8[]{1, 1, 0, -1, -1, -1, 0, 1};
constexpr int DX8[]{0, -1, -1, -1, 0, 1, 1, 1};
template <typename T, typename U>
inline bool chmax(T& a, U b) { return a < b ? (a = b, true) : false; }
template <typename T, typename U>
inline bool chmin(T& a, U b) { return a > b ? (a = b, true) : false; }
struct IOSetup {
IOSetup() {
std::cin.tie(nullptr);
std::ios_base::sync_with_stdio(false);
std::cout << fixed << setprecision(20);
}
} iosetup;
template <typename T = long long>
struct Rational {
T num, den;
Rational() : num(0), den(1) {}
Rational(const T num, const T den = 1) : num(num), den(den) {
// assert(den != 0);
reduce();
}
template <typename Real = long double>
Real to_real() const { return static_cast<Real>(num) / den; }
Rational& operator+=(const Rational& x) {
const T g = std::__gcd(den, x.den);
num = num * (x.den / g) + x.num * (den / g);
den *= x.den / g;
reduce();
return *this;
}
Rational& operator-=(const Rational& x) { return *this += -x; }
Rational& operator*=(const Rational& x) {
const T g1 = std::__gcd(num, x.den), g2 = std::__gcd(den, x.num);
num = (num / g1) * (x.num / g2);
den = (den / g2) * (x.den / g1);
reduce();
return *this;
}
Rational& operator/=(const Rational& x) {
return *this *= Rational(x.den, x.num);
}
bool operator==(const Rational& x) const {
return num == x.num && den == x.den;
}
bool operator!=(const Rational& x) const { return !(*this == x); }
bool operator<(const Rational& x) const { return (x - *this).num > 0; }
bool operator<=(const Rational& x) const { return !(x < *this); }
bool operator>(const Rational& x) const { return x < *this; }
bool operator>=(const Rational& x) const { return !(*this < x); }
Rational& operator++() {
if ((num += den) == 0) den = 1;
return *this;
}
Rational operator++(int) {
const Rational res = *this;
++*this;
return res;
}
Rational& operator--() {
if ((num -= den) == 0) den = 1;
return *this;
}
Rational operator--(int) {
const Rational res = *this;
--*this;
return res;
}
Rational operator+() const { return *this; }
Rational operator-() const { return Rational(-num, den); }
Rational operator+(const Rational& x) const { return Rational(*this) += x; }
Rational operator-(const Rational& x) const { return Rational(*this) -= x; }
Rational operator*(const Rational& x) const { return Rational(*this) *= x; }
Rational operator/(const Rational& x) const { return Rational(*this) /= x; }
friend std::ostream& operator<<(std::ostream& os, const Rational& x) {
if (x.den == 1) return os << x.num;
return os << x.num << '/' << x.den;
}
private:
void reduce() {
const T g = std::__gcd(num, den);
num /= g;
den /= g;
if (den < 0) {
num = -num;
den = -den;
}
}
};
namespace std {
template <typename T>
Rational<T> abs(Rational<T> x) {
if (x.num < 0) x.num = -x.num;
return x;
}
template <typename T>
Rational<T> max(const Rational<T>& a, const Rational<T>& b) {
return a < b ? b : a;
}
template <typename T>
Rational<T> min(const Rational<T>& a, const Rational<T>& b) {
return a < b ? a : b;
}
template <typename T> struct numeric_limits<Rational<T>> {
static constexpr Rational<T> max() { return std::numeric_limits<T>::max(); }
static constexpr Rational<T> lowest() {
return std::numeric_limits<T>::lowest();
}
};
} // namespace std
int main() {
int p, q; cin >> p >> q;
cout << p + q << '\n';
return 0;
// using rational = Rational<>;
// FOR(P, 2, M) FOR(Q, P + 1, M) {
// if (gcd(P, Q) > 1) continue;
// rational r(INF, 1);
// for (int q = Q - 1; q >= 1; --q) {
// const int p = (P * q + Q - 1) / Q;
// chmin(r, rational(p, q));
// // cout << p << '/' << q << '\n';
// }
// rational l(0);
// for (int q = Q - 1; q >= 1; --q) {
// const int p = P * q / Q - (P * q % Q == 0);
// chmax(l, rational(p, q));
// // cout << p << '/' << q << '\n';
// }
// cout << P << ' ' << Q << ": " << l.num + l.den + r.num + r.den << '\n';
// // cout << P << ' ' << Q << ": " << l << ' ' << r << '\n';
// }
}
emthrm