結果

問題 No.5009 Draw A Convex Polygon
コンテスト
ユーザー emthrm
提出日時 2022-12-02 18:21:38
言語 C++17(gcc12)
(gcc 12.3.0 + boost 1.87.0)
結果
AC  
実行時間 891 ms / 2,600 ms
コード長 5,230 bytes
コンパイル時間 2,098 ms
実行使用メモリ 24,448 KB
スコア 1,000,000
平均クエリ数 1000001.00
最終ジャッジ日時 2022-12-02 18:21:43
合計ジャッジ時間 4,880 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge12
純コード判定しない問題か言語
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 1
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#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 = 998244353;
// constexpr int MOD = 1000000007;
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 num * x.den < x.num * den; }
  // 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() {
  using rational = Rational<>;
  constexpr int N = 1000000, C = 641;
  vector<rational> v{0};
  v.reserve(C * C + 1);
  for (int num = 1; num <= C; ++num) {
    for (int den = 1; den <= C; ++den) {
      v.emplace_back(num, den);
    }
  }
  sort(ALL(v));
  v.erase(unique(ALL(v)), v.end());
  v.resize(N / 4);

  cout << N << '\n';
  int x = 0, y = -80159743;
  for (const rational& r : v) {
    x += r.den;
    y += r.num;
    cout << x << ' ' << y << '\n';
  }
  for (const rational& r : v) {
    x -= r.num;
    y += r.den;
    cout << x << ' ' << y << '\n';
  }
  for (const rational& r : v) {
    x -= r.den;
    y -= r.num;
    cout << x << ' ' << y << '\n';
  }
  for (const rational& r : v) {
    x += r.num;
    y -= r.den;
    cout << x << ' ' << y << '\n';
  }
  return 0;

  // int lb = llround(sqrt(N / 4)), ub = 10000;
  // while (ub - lb > 1) {
  //   const int mid = (lb + ub) / 2;
  //   vector<rational> v;
  //   v.reserve(mid * mid);
  //   for (int num = 1; num <= mid; ++num) {
  //     for (int den = 1; den <= mid; ++den) {
  //       v.emplace_back(num, den);
  //     }
  //   }
  //   sort(ALL(v));
  //   v.erase(unique(ALL(v)), v.end());
  //   (v.size() >= N / 4 - 2 ? ub : lb) = mid;
  // }
}
0