#define _USE_MATH_DEFINES #include 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 inline bool chmax(T& a, U b) { return a < b ? (a = b, true) : false; } template 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 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 Real to_real() const { return static_cast(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 Rational abs(Rational x) { if (x.num < 0) x.num = -x.num; return x; } template Rational max(const Rational& a, const Rational& b) { return a < b ? b : a; } template Rational min(const Rational& a, const Rational& b) { return a < b ? a : b; } template struct numeric_limits> { static constexpr Rational max() { return std::numeric_limits::max(); } static constexpr Rational lowest() { return std::numeric_limits::lowest(); } }; } // namespace std int main() { using rational = Rational<>; constexpr int N = 1000000, C = 641; vector 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 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; // } }