#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #ifndef KK2_TEMPLATE_PROCON_HPP #define KK2_TEMPLATE_PROCON_HPP 1 #ifndef KK2_TEMPLATE_CONSTANT_HPP #define KK2_TEMPLATE_CONSTANT_HPP 1 #ifndef KK2_TEMPLATE_TYPE_ALIAS_HPP #define KK2_TEMPLATE_TYPE_ALIAS_HPP 1 using u32 = unsigned int; using i64 = long long; using u64 = unsigned long long; using i128 = __int128_t; using u128 = __uint128_t; using pi = std::pair; using pl = std::pair; using pil = std::pair; using pli = std::pair; template using vc = std::vector; template using vvc = std::vector>; template using vvvc = std::vector>; template using vvvvc = std::vector>; template using pq = std::priority_queue; template using pqi = std::priority_queue, std::greater>; #endif // KK2_TEMPLATE_TYPE_ALIAS_HPP template constexpr T infty = 0; template <> constexpr int infty = (1 << 30) - 123; template <> constexpr i64 infty = (1ll << 62) - (1ll << 31); template <> constexpr i128 infty = (i128(1) << 126) - (i128(1) << 63); template <> constexpr u32 infty = infty; template <> constexpr u64 infty = infty; template <> constexpr u128 infty = infty; template <> constexpr double infty = infty; template <> constexpr long double infty = infty; constexpr int mod = 998244353; constexpr int modu = 1e9 + 7; constexpr long double PI = 3.14159265358979323846; #endif // KK2_TEMPLATE_CONSTANT_HPP #ifndef KK2_TEMPLATE_FUNCTION_UTIL_HPP #define KK2_TEMPLATE_FUNCTION_UTIL_HPP 1 #ifndef KK2_MATH_MONOID_MAX_HPP #define KK2_MATH_MONOID_MAX_HPP 1 #ifndef KK2_TYPE_TRAITS_IO_HPP #define KK2_TYPE_TRAITS_IO_HPP 1 namespace kk2 { namespace type_traits { struct istream_tag {}; struct ostream_tag {}; } // namespace type_traits template using is_standard_istream = typename std::conditional::value || std::is_same::value, std::true_type, std::false_type>::type; template using is_standard_ostream = typename std::conditional::value || std::is_same::value, std::true_type, std::false_type>::type; template using is_user_defined_istream = std::is_base_of; template using is_user_defined_ostream = std::is_base_of; template using is_istream = typename std::conditional::value || is_user_defined_istream::value, std::true_type, std::false_type>::type; template using is_ostream = typename std::conditional::value || is_user_defined_ostream::value, std::true_type, std::false_type>::type; template using is_istream_t = std::enable_if_t::value>; template using is_ostream_t = std::enable_if_t::value>; } // namespace kk2 #endif // KK2_TYPE_TRAITS_IO_HPP namespace kk2 { namespace monoid { template > struct Max { static constexpr bool commutative = true; using M = Max; S a; bool is_unit; Max() : a(S()), is_unit(true) {} Max(S a_) : a(a_), is_unit(false) {} operator S() const { return a; } inline static M op(M l, M r) { if (l.is_unit or r.is_unit) return l.is_unit ? r : l; return Compare{}(l.a, r.a) ? r : l; } inline static M unit() { return M(); } bool operator==(const M &rhs) const { return is_unit == rhs.is_unit and (is_unit or a == rhs.a); } bool operator!=(const M &rhs) const { return is_unit != rhs.is_unit or (!is_unit and a != rhs.a); } template * = nullptr> friend OStream &operator<<(OStream &os, const M &x) { if (x.is_unit) os << "-inf"; else os << x.a; return os; } template * = nullptr> friend IStream &operator>>(IStream &is, M &x) { is >> x.a; x.is_unit = false; return is; } }; } // namespace monoid } // namespace kk2 #endif // MATH_MONOID_MAX_HPP #ifndef KK2_MATH_MONOID_MIN_HPP #define KK2_MATH_MONOID_MIN_HPP 1 namespace kk2 { namespace monoid { template > struct Min { static constexpr bool commutative = true; using M = Min; S a; bool is_unit; Min() : a(S()), is_unit(true) {} Min(S a_) : a(a_), is_unit(false) {} operator S() const { return a; } inline static M op(M l, M r) { if (l.is_unit or r.is_unit) return l.is_unit ? r : l; return Compare{}(l.a, r.a) ? l : r; } inline static M unit() { return M(); } bool operator==(const M &rhs) const { return is_unit == rhs.is_unit and (is_unit or a == rhs.a); } bool operator!=(const M &rhs) const { return is_unit != rhs.is_unit or (!is_unit and a != rhs.a); } template * = nullptr> friend OStream &operator<<(OStream &os, const M &x) { if (x.is_unit) os << "inf"; else os << x.a; return os; } template * = nullptr> friend IStream &operator>>(IStream &is, M &x) { is >> x.a; x.is_unit = false; return is; } }; } // namespace monoid } // namespace kk2 #endif // KK2_MATH_MONOID_MIN_HPP #ifndef KK2_TYPE_TRAITS_CONTAINER_TRAITS_HPP #define KK2_TYPE_TRAITS_CONTAINER_TRAITS_HPP 1 namespace kk2 { template struct is_vector : std::false_type {}; template struct is_vector> : std::true_type {}; template struct is_container : std::false_type {}; template struct is_container> : std::true_type { }; template struct is_container> : std::true_type {}; template struct is_container> : std::true_type {}; template struct is_container> : std::true_type {}; template struct is_container> : std::true_type {}; template using is_container_t = typename std::enable_if_t::value>; } // namespace kk2 #endif // KK2_TYPE_TRAITS_CONTAINER_TRAITS_HPP namespace kk2 { template auto make_vector(int first, Sizes... sizes) { if constexpr (sizeof...(sizes) == 0) { return std::vector(first); } else { return std::vector(sizes...))>(first, make_vector(sizes...)); } } template void fill_all(std::vector &v, const U &x) { if constexpr (is_vector::value) { for (auto &u : v) fill_all(u, x); } else { std::fill(v.begin(), v.end(), T(x)); } } template int iota_all(std::vector &v, U x, int offset = 0) { if constexpr (is_vector::value) { for (auto &u : v) offset += iota_all(u, x + offset); } else { for (auto &u : v) u = x++, ++offset; } return offset; } template int mysize(const C &c) { return size(c); } // T: commutative monoid, F: (U, T) -> U template U all_monoid_prod(const std::vector &v, U unit, const F &f) { U res = unit; if constexpr (is_vector::value) { for (const auto &x : v) res = f(res, all_monoid_prod(x, unit, f)); } else { for (const auto &x : v) res = f(res, x); } return res; } template U all_sum(const std::vector &v, U unit = U()) { return all_monoid_prod(v, unit, [](U a, U b) { return a + b; }); } template U all_prod(const std::vector &v, U unit = U(1)) { return all_monoid_prod(v, unit, [](U a, U b) { return a * b; }); } template U all_xor(const std::vector &v, U unit = U()) { return all_monoid_prod(v, unit, [](U a, U b) { return a ^ b; }); } template U all_and(const std::vector &v, U unit = U(-1)) { return all_monoid_prod(v, unit, [](U a, U b) { return a & b; }); } template U all_or(const std::vector &v, U unit = U()) { return all_monoid_prod(v, unit, [](U a, U b) { return a | b; }); } template U all_min(const std::vector &v) { return all_monoid_prod, T>(v, monoid::Min::unit(), monoid::Min::op); } template U all_max(const std::vector &v) { return all_monoid_prod, T>(v, monoid::Max::unit(), monoid::Max::op); } template U all_gcd(const std::vector &v, U unit = U()) { return all_monoid_prod(v, unit, [](U a, U b) { return std::gcd(a, b); }); } template U all_lcm(const std::vector &v, U unit = U(1)) { return all_monoid_prod(v, unit, [](U a, U b) { return std::lcm(a, b); }); } template int all_count(const std::vector &v, U x) { return all_monoid_prod(v, 0, [x](int a, U y) { return a + int(x == y); }); } } // namespace kk2 #endif // KK2_TEMPLATE_FUNCTION_UTIL_HPP #ifndef KK2_TEMPLATE_IO_UTIL_HPP #define KK2_TEMPLATE_IO_UTIL_HPP 1 // なんかoj verifyはプロトタイプ宣言が落ちる namespace impl { struct read { template inline static void all_read(IStream &is, T &x) { is >> x; } template inline static void all_read(IStream &is, std::pair &p) { all_read(is, p.first); all_read(is, p.second); } template inline static void all_read(IStream &is, std::vector &v) { for (T &x : v) all_read(is, x); } template inline static void all_read(IStream &is, std::array &a) { for (T &x : a) all_read(is, x); } }; struct write { template inline static void all_write(OStream &os, const T &x) { os << x; } template inline static void all_write(OStream &os, const std::pair &p) { all_write(os, p.first); all_write(os, ' '); all_write(os, p.second); } template inline static void all_write(OStream &os, const std::vector &v) { for (int i = 0; i < (int)v.size(); ++i) { if (i) all_write(os, ' '); all_write(os, v[i]); } } template inline static void all_write(OStream &os, const std::array &a) { for (int i = 0; i < (int)F; ++i) { if (i) all_write(os, ' '); all_write(os, a[i]); } } }; } // namespace impl template * = nullptr> IStream &operator>>(IStream &is, std::pair &p) { impl::read::all_read(is, p); return is; } template * = nullptr> IStream &operator>>(IStream &is, std::vector &v) { impl::read::all_read(is, v); return is; } template * = nullptr> IStream &operator>>(IStream &is, std::array &a) { impl::read::all_read(is, a); return is; } template * = nullptr> OStream &operator<<(OStream &os, const std::pair &p) { impl::write::all_write(os, p); return os; } template * = nullptr> OStream &operator<<(OStream &os, const std::vector &v) { impl::write::all_write(os, v); return os; } template * = nullptr> OStream &operator<<(OStream &os, const std::array &a) { impl::write::all_write(os, a); return os; } #endif // KK2_TEMPLATE_IO_UTIL_HPP #ifndef KK2_TEMPLATE_MACROS_HPP #define KK2_TEMPLATE_MACROS_HPP 1 #define rep1(a) for (long long _ = 0; _ < (long long)(a); ++_) #define rep2(i, a) for (long long i = 0; i < (long long)(a); ++i) #define rep3(i, a, b) for (long long i = (a); i < (long long)(b); ++i) #define repi2(i, a) for (long long i = (a) - 1; i >= 0; --i) #define repi3(i, a, b) for (long long i = (a) - 1; i >= (long long)(b); --i) #define overload3(a, b, c, d, ...) d #define rep(...) overload3(__VA_ARGS__, rep3, rep2, rep1)(__VA_ARGS__) #define repi(...) overload3(__VA_ARGS__, repi3, repi2, rep1)(__VA_ARGS__) #define fi first #define se second #define all(p) begin(p), end(p) #endif // KK2_TEMPLATE_MACROS_HPP struct FastIOSetUp { FastIOSetUp() { std::ios::sync_with_stdio(false); std::cin.tie(nullptr); } } fast_io_set_up; auto &kin = std::cin; auto &kout = std::cout; auto (*kendl)(std::ostream &) = std::endl>; void Yes(bool b = 1) { kout << (b ? "Yes\n" : "No\n"); } void No(bool b = 1) { kout << (b ? "No\n" : "Yes\n"); } void YES(bool b = 1) { kout << (b ? "YES\n" : "NO\n"); } void NO(bool b = 1) { kout << (b ? "NO\n" : "YES\n"); } void yes(bool b = 1) { kout << (b ? "yes\n" : "no\n"); } void no(bool b = 1) { kout << (b ? "no\n" : "yes\n"); } template inline bool chmax(T &a, const S &b) { return (a < b ? a = b, 1 : 0); } template inline bool chmin(T &a, const S &b) { return (a > b ? a = b, 1 : 0); } std::istream &operator>>(std::istream &is, u128 &x) { std::string s; is >> s; x = 0; for (char c : s) { assert('0' <= c && c <= '9'); x = x * 10 + c - '0'; } return is; } std::istream &operator>>(std::istream &is, i128 &x) { std::string s; is >> s; bool neg = s[0] == '-'; x = 0; for (int i = neg; i < (int)s.size(); i++) { assert('0' <= s[i] && s[i] <= '9'); x = x * 10 + s[i] - '0'; } if (neg) x = -x; return is; } std::ostream &operator<<(std::ostream &os, u128 x) { if (x == 0) return os << '0'; std::string s; while (x) { s.push_back('0' + x % 10); x /= 10; } std::reverse(s.begin(), s.end()); return os << s; } std::ostream &operator<<(std::ostream &os, i128 x) { if (x == 0) return os << '0'; if (x < 0) { os << '-'; x = -x; } std::string s; while (x) { s.push_back('0' + x % 10); x /= 10; } std::reverse(s.begin(), s.end()); return os << s; } #endif // KK2_TEMPLATE_PROCON_HPP #ifndef KK2_TEMPLATE_DEBUG_HPP #define KK2_TEMPLATE_DEBUG_HPP 1 #ifndef KK2_TYPE_TRAITS_MEMBER_HPP #define KK2_TYPE_TRAITS_MEMBER_HPP 1 namespace kk2 { #define HAS_MEMBER_FUNC(member) \ template struct has_member_func_##member##_impl { \ template \ static std::true_type check(decltype(std::declval().member(std::declval()...)) *); \ template static std::false_type check(...); \ using type = decltype(check(nullptr)); \ }; \ template struct has_member_func_##member \ : has_member_func_##member##_impl::type {}; \ template using has_member_func_##member##_t = \ std::enable_if_t::value>; \ template using not_has_member_func_##member##_t = \ std::enable_if_t::value>; #define HAS_MEMBER_VAR(member) \ template struct has_member_var_##member##_impl { \ template static std::true_type check(decltype(std::declval().member) *); \ template static std::false_type check(...); \ using type = decltype(check(nullptr)); \ }; \ template struct has_member_var_##member \ : has_member_var_##member##_impl::type {}; \ template using has_member_var_##member##_t = \ std::enable_if_t::value>; \ template using not_has_member_var_##member##_t = \ std::enable_if_t::value>; HAS_MEMBER_FUNC(debug_output) HAS_MEMBER_FUNC(val) #undef HAS_MEMBER_FUNC #undef HAS_MEMBER_VAR } // namespace kk2 #endif // KK2_TYPE_TRAITS_MEMBER_HPP namespace kk2 { namespace debug { #ifdef KK2 template *> void output(OStream &os); template *> void output(OStream &os, const T &t); template *> void output(OStream &os, const T &t, const Args &...args); template *> void output(OStream &os, const std::vector &v); template *> void output(OStream &os, const std::vector> &d); template *> void output(OStream &os, const std::array &a); template *> void output(OStream &os, const std::pair &p); template *> void output(OStream &os, const std::tuple &t); template *> void output(OStream &os, const std::queue &q); template *> void output(OStream &os, const std::priority_queue &q); template *> void output(OStream &os, const std::deque &d); template *> void output(OStream &os, const std::stack &s); template *> void output(OStream &os, const std::set &s); template *> void output(OStream &os, const std::multiset &s); template *> void output(OStream &os, const std::unordered_set &s); template *> void output(OStream &os, const std::unordered_multiset &s); template *> void output(OStream &os, const std::map &m); template *> void output(OStream &os, const std::unordered_map &m); template * = nullptr> void output(OStream &) {} template * = nullptr> void output(OStream &os, const T &t) { if constexpr (has_member_func_debug_output::value) { t.debug_output(os); } else { os << t; } } template * = nullptr> void output(OStream &os, const T &t, const Args &...args) { output(os, t); os << ", "; output(os, args...); } template * = nullptr> void output(OStream &os, const std::vector &v) { os << "["; for (int i = 0; i < (int)v.size(); i++) { output(os, v[i]); if (i + 1 != (int)v.size()) os << ", "; } os << "]"; } template * = nullptr> void output(OStream &os, const std::vector> &d) { os << "[\n"; for (int i = 0; i < (int)d.size(); i++) { output(os, d[i]); output(os, "\n"); } os << "]"; } template * = nullptr> void output(OStream &os, const std::array &a) { os << "["; for (int i = 0; i < (int)F; i++) { output(os, a[i]); if (i + 1 != (int)F) os << ", "; } os << "]"; } template * = nullptr> void output(OStream &os, const std::pair &p) { os << "("; output(os, p.first); os << ", "; output(os, p.second); os << ")"; } template * = nullptr> void output(OStream &os, const std::tuple &t) { os << "("; std::apply([&](const Args &...args) { output(os, args...); }, t); os << ")"; } template * = nullptr> void output(OStream &os, const std::queue &q) { os << "["; std::queue tmp = q; while (!tmp.empty()) { output(os, tmp.front()); tmp.pop(); if (!tmp.empty()) os << ", "; } os << "]"; } template * = nullptr> void output(OStream &os, const std::priority_queue &q) { os << "["; std::priority_queue tmp = q; while (!tmp.empty()) { output(os, tmp.top()); tmp.pop(); if (!tmp.empty()) os << ", "; } os << "]"; } template * = nullptr> void output(OStream &os, const std::deque &d) { os << "["; std::deque tmp = d; while (!tmp.empty()) { output(os, tmp.front()); tmp.pop_front(); if (!tmp.empty()) os << ", "; } os << "]"; } template * = nullptr> void output(OStream &os, const std::stack &s) { os << "["; std::stack tmp = s; std::vector v; while (!tmp.empty()) { v.push_back(tmp.top()); tmp.pop(); } for (int i = (int)v.size() - 1; i >= 0; i--) { output(os, v[i]); if (i != 0) os << ", "; } os << "]"; } template * = nullptr> void output(OStream &os, const std::set &s) { os << "{"; std::set tmp = s; for (auto it = tmp.begin(); it != tmp.end(); ++it) { output(os, *it); if (std::next(it) != tmp.end()) os << ", "; } os << "}"; } template * = nullptr> void output(OStream &os, const std::multiset &s) { os << "{"; std::multiset tmp = s; for (auto it = tmp.begin(); it != tmp.end(); ++it) { output(os, *it); if (std::next(it) != tmp.end()) os << ", "; } os << "}"; } template * = nullptr> void output(OStream &os, const std::unordered_set &s) { os << "{"; std::unordered_set tmp = s; for (auto it = tmp.begin(); it != tmp.end(); ++it) { output(os, *it); if (std::next(it) != tmp.end()) os << ", "; } os << "}"; } template * = nullptr> void output(OStream &os, const std::unordered_multiset &s) { os << "{"; std::unordered_multiset tmp = s; for (auto it = tmp.begin(); it != tmp.end(); ++it) { output(os, *it); if (std::next(it) != tmp.end()) os << ", "; } os << "}"; } template * = nullptr> void output(OStream &os, const std::map &m) { os << "{"; std::map tmp = m; for (auto it = tmp.begin(); it != tmp.end(); ++it) { output(os, it->first); os << ": "; output(os, it->second); if (std::next(it) != tmp.end()) os << ", "; } os << "}"; } template * = nullptr> void output(OStream &os, const std::unordered_map &m) { os << "{"; std::unordered_map tmp = m; for (auto it = tmp.begin(); it != tmp.end(); ++it) { output(os, it->first); os << ": "; output(os, it->second); if (std::next(it) != tmp.end()) os << ", "; } os << "}"; } template * = nullptr> void outputln(OStream &os) { os << '\n'; os.flush(); } template * = nullptr> void outputln(OStream &os, const T &t, const Args &...args) { output(os, t, args...); os << '\n'; os.flush(); } std::vector sep(const char *s) { std::vector res; std::string now; int dep = 0; while (true) { if (*s == '\0') { res.emplace_back(now); break; } if (*s == '(' or *s == '[' or *s == '{') dep++; if (*s == ')' or *s == ']' or *s == '}') dep--; if (dep == 0 and *s == ',') { res.emplace_back(now); now.clear(); } else if (!isspace(*s)) { now += *s; } s++; } return res; } void show_vars(const std::vector &, int) {} template void show_vars(const std::vector &name, int pos, const T &t, const Args &...args) { output(std::cerr, name[pos++] + ":", t); if (sizeof...(args) > 0) output(std::cerr, ", "); show_vars(name, pos, args...); } #define kdebug(...) \ { \ std::cerr << "line:" << __LINE__ << ' '; \ kk2::debug::show_vars(kk2::debug::sep(#__VA_ARGS__), 0, __VA_ARGS__); \ std::cerr << std::endl; \ } #define kput(s) \ { \ std::cerr << "line:" << __LINE__ << ' '; \ kk2::debug::outputln(std::cerr, s); \ } #else template * = nullptr> void output(OStream &, const Args &...) {} template * = nullptr> void outputln(OStream &, const Args &...) {} template void fix_warn(const Args &...) {} #define kdebug(...) kk2::debug::fix_warn(__VA_ARGS__); #define kput(s) kk2::debug::fix_warn(s) #endif // KK2 } // namespace debug } // namespace kk2 #endif // KK2_TEMPLATE_DEBUG_HPP using namespace std; void solve() { /* E: fixed 1/E * N^1/E * log(N) Eは40まで試せばいい sum_{i=L, ..., R} i^Eの計算は面倒くさい!!!! https://judge.yosupo.jp/problem/sum_of_exponential_times_polynomial E = 1はなぞ.因数分解とか find (l, r) r * (r + 1) / 2 - l * (l - 1) / 2 = n すみません面倒なので,AIにやらせます prompt: l + l + 1 + ... + r == nをみたすようなl, r (0 < l <= r) をすべて見つけたいです.これは,左辺を因数分解し,nを素因数分解すれば解けます. E = 2, 3は頑張る # >= 4は愚直でも二分探索する必要がない */ i64 n; kin >> n; vc> res; // E = 1: l + (l+1) + ... + r = n // (r-l+1)(l+r)/2 = n => (r-l+1)(l+r) = 2n // len = r-l+1, sum = l+r とすると len * sum = 2n // l = (sum - len + 1)/2, r = (sum + len - 1)/2 { i64 two_n = 2 * n; // two_n の約数を全列挙 vc divs; for (i64 d = 1; d * d <= two_n; d++) { if (two_n % d == 0) { divs.push_back(d); if (d != two_n / d) divs.push_back(two_n / d); } } for (auto len : divs) { i64 sum = two_n / len; // l = (sum - len + 1) / 2, r = (sum + len - 1) / 2 // 条件: sum - len + 1 が偶数(つまり sum と len の偶奇が同じ) if ((sum - len + 1) % 2 != 0) continue; i64 l = (sum - len + 1) / 2; i64 r = (sum + len - 1) / 2; if (l > 0 && l <= r) { res.emplace_back(1, l, r); } } } auto fast = [&](int e, auto calc) { // auto calc = [&](i64 l, i64 r) { // return i128(r) * (r + 1) * r * (r + 1) / 4 - i128(l - 1) * l * (l - 1) * l / 4; // }; auto is_ok = [&](i64 l, i64 r) -> bool { return n >= calc(l, r); }; kdebug(calc(1, 2)); kdebug(calc(1, 3)); rep (r, 1, 1e9) { if (calc(r, r) > n) break; int ng = 0, ok = r; while (ok - ng > 1) { int mid = (ng + ok) / 2; // kdebug(mid, r, calc(mid, r)); if (is_ok(mid, r)) ok = mid; else ng = mid; } // kdebug(r, ok); if (calc(ok, r) == n) { res.emplace_back(e, ok, r); } } }; fast(2, [](i64 l, i64 r) { return i64(r) * (r + 1) * (2 * r + 1) / 6 - i64(l - 1) * l * (2 * l - 1) / 6; }); fast(3, [](i64 l, i64 r) { return i128(r) * (r + 1) * r * (r + 1) / 4 - i128(l - 1) * l * (l - 1) * l / 4; }); rep (e, 4, 61) { vc pre(1e3 + 1, -1); pre[0] = 0; auto safe_pow = [&](int x) -> i64 { i64 pw = 1; rep (e) { if (pw > n / x) return -1; pw *= x; } return pw; }; rep (i, 1, 1e3 + 1) { auto tmp = safe_pow(i); if (tmp == -1) break; pre[i] = tmp + pre[i - 1]; } rep (r, 1, 1e3 + 1) rep (l, 1, r + 1) if (pre[r] != -1) { if (pre[r] - pre[l - 1] == n) res.emplace_back(e, l, r); } } sort(all(res)); kout << res.size() << "\n"; for (auto [e, l, r] : res) kout << e << " " << l << " " << r << "\n"; } int main() { #ifdef KK2 int t = 2; #else int t = 1; #endif // kin >> t; rep (t) solve(); return 0; } // Author: kk2 // converted by https://github.com/kk2a/cpp-bundle // 2025-11-02 22:32:39