結果
| 問題 | No.3522 冪乗乗 |
| コンテスト | |
| ユーザー |
 hitonanode
|
| 提出日時 | 2026-05-01 22:52:18 |
| 言語 | C++23 (gcc 15.2.0 + boost 1.89.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 12,174 bytes |
| 記録 | |
| コンパイル時間 | 6,202 ms |
| コンパイル使用メモリ | 508,852 KB |
| 実行使用メモリ | 30,320 KB |
| 平均クエリ数 | 23.57 |
| 最終ジャッジ日時 | 2026-05-01 22:52:37 |
| 合計ジャッジ時間 | 12,186 ms |
|
ジャッジサーバーID （参考情報） |
judge1_0 / judge3_0 |
（要ログイン）
| ファイルパターン | 結果 |
|---|---|
| other | AC * 51 WA * 12 |
ソースコード
#include <algorithm>
#include <array>
#include <bitset>
#include <cassert>
#include <chrono>
#include <cmath>
#include <complex>
#include <deque>
#include <forward_list>
#include <fstream>
#include <functional>
#include <iomanip>
#include <ios>
#include <iostream>
#include <limits>
#include <list>
#include <map>
#include <memory>
#include <numeric>
#include <optional>
#include <queue>
#include <random>
#include <set>
#include <sstream>
#include <stack>
#include <string>
#include <tuple>
#include <type_traits>
#include <unordered_map>
#include <unordered_set>
#include <utility>
#include <vector>
using namespace std;
using lint = long long;
using pint = pair<int, int>;
using plint = pair<lint, lint>;
struct fast_ios { fast_ios(){ cin.tie(nullptr), ios::sync_with_stdio(false), cout << fixed << setprecision(20); }; } fast_ios_;
#define ALL(x) (x).begin(), (x).end()
#define FOR(i, begin, end) for(int i=(begin),i##_end_=(end);i<i##_end_;i++)
#define IFOR(i, begin, end) for(int i=(end)-1,i##_begin_=(begin);i>=i##_begin_;i--)
#define REP(i, n) FOR(i,0,n)
#define IREP(i, n) IFOR(i,0,n)
template <typename T> bool chmax(T &m, const T q) { return m < q ? (m = q, true) : false; }
template <typename T> bool chmin(T &m, const T q) { return m > q ? (m = q, true) : false; }
const std::vector<std::pair<int, int>> grid_dxs{{1, 0}, {-1, 0}, {0, 1}, {0, -1}};
int floor_lg(long long x) { return x <= 0 ? -1 : 63 - __builtin_clzll(x); }
template <class T1, class T2> T1 floor_div(T1 num, T2 den) { return (num > 0 ? num / den : -((-num + den - 1) / den)); }
template <class T1, class T2> std::pair<T1, T2> operator+(const std::pair<T1, T2> &l, const std::pair<T1, T2> &r) { return std::make_pair(l.first + r.first, l.second + r.second); }
template <class T1, class T2> std::pair<T1, T2> operator-(const std::pair<T1, T2> &l, const std::pair<T1, T2> &r) { return std::make_pair(l.first - r.first, l.second - r.second); }
template <class T> std::vector<T> sort_unique(std::vector<T> vec) { sort(vec.begin(), vec.end()), vec.erase(unique(vec.begin(), vec.end()), vec.end()); return vec; }
template <class T> int arglb(const std::vector<T> &v, const T &x) { return std::distance(v.begin(), std::lower_bound(v.begin(), v.end(), x)); }
template <class T> int argub(const std::vector<T> &v, const T &x) { return std::distance(v.begin(), std::upper_bound(v.begin(), v.end(), x)); }
template <class IStream, class T> IStream &operator>>(IStream &is, std::vector<T> &vec) { for (auto &v : vec) is >> v; return is; }
template <class OStream, class T> OStream &operator<<(OStream &os, const std::vector<T> &vec);
template <class OStream, class T, size_t sz> OStream &operator<<(OStream &os, const std::array<T, sz> &arr);
template <class OStream, class T, class TH> OStream &operator<<(OStream &os, const std::unordered_set<T, TH> &vec);
template <class OStream, class T, class U> OStream &operator<<(OStream &os, const pair<T, U> &pa);
template <class OStream, class T> OStream &operator<<(OStream &os, const std::deque<T> &vec);
template <class OStream, class T> OStream &operator<<(OStream &os, const std::set<T> &vec);
template <class OStream, class T> OStream &operator<<(OStream &os, const std::multiset<T> &vec);
template <class OStream, class T> OStream &operator<<(OStream &os, const std::unordered_multiset<T> &vec);
template <class OStream, class T, class U> OStream &operator<<(OStream &os, const std::pair<T, U> &pa);
template <class OStream, class TK, class TV> OStream &operator<<(OStream &os, const std::map<TK, TV> &mp);
template <class OStream, class TK, class TV, class TH> OStream &operator<<(OStream &os, const std::unordered_map<TK, TV, TH> &mp);
template <class OStream, class... T> OStream &operator<<(OStream &os, const std::tuple<T...> &tpl);
template <class OStream, class T> OStream &operator<<(OStream &os, const std::vector<T> &vec) { os << '['; for (auto v : vec) os << v << ','; os << ']'; return os; }
template <class OStream, class T, size_t sz> OStream &operator<<(OStream &os, const std::array<T, sz> &arr) { os << '['; for (auto v : arr) os << v << ','; os << ']'; return os; }
template <class... T> std::istream &operator>>(std::istream &is, std::tuple<T...> &tpl) { std::apply([&is](auto &&... args) { ((is >> args), ...);}, tpl); return is; }
template <class OStream, class... T> OStream &operator<<(OStream &os, const std::tuple<T...> &tpl) { os << '('; std::apply([&os](auto &&... args) { ((os << args << ','), ...);}, tpl); return os << ')'; }
template <class OStream, class T, class TH> OStream &operator<<(OStream &os, const std::unordered_set<T, TH> &vec) { os << '{'; for (auto v : vec) os << v << ','; os << '}'; return os; }
template <class OStream, class T> OStream &operator<<(OStream &os, const std::deque<T> &vec) { os << "deq["; for (auto v : vec) os << v << ','; os << ']'; return os; }
template <class OStream, class T> OStream &operator<<(OStream &os, const std::set<T> &vec) { os << '{'; for (auto v : vec) os << v << ','; os << '}'; return os; }
template <class OStream, class T> OStream &operator<<(OStream &os, const std::multiset<T> &vec) { os << '{'; for (auto v : vec) os << v << ','; os << '}'; return os; }
template <class OStream, class T> OStream &operator<<(OStream &os, const std::unordered_multiset<T> &vec) { os << '{'; for (auto v : vec) os << v << ','; os << '}'; return os; }
template <class OStream, class T, class U> OStream &operator<<(OStream &os, const std::pair<T, U> &pa) { return os << '(' << pa.first << ',' << pa.second << ')'; }
template <class OStream, class TK, class TV> OStream &operator<<(OStream &os, const std::map<TK, TV> &mp) { os << '{'; for (auto v : mp) os << v.first << "=>" << v.second << ','; os << '}'; return os; }
template <class OStream, class TK, class TV, class TH> OStream &operator<<(OStream &os, const std::unordered_map<TK, TV, TH> &mp) { os << '{'; for (auto v : mp) os << v.first << "=>" << v.second << ','; os << '}'; return os; }
#ifdef HITONANODE_LOCAL
const string COLOR_RESET = "\033[0m", BRIGHT_CYAN = "\033[1;36m", NORMAL_FAINT = "\033[0;2m";
#define dbg(x) std::cerr << BRIGHT_CYAN << #x << COLOR_RESET << " = " << (x) << NORMAL_FAINT << " (L" << __LINE__ << ") " << __FILE__ << COLOR_RESET << std::endl
#define dbgif(cond, x) ((cond) ? std::cerr << BRIGHT_CYAN << #x << COLOR_RESET << " = " << (x) << NORMAL_FAINT << " (L" << __LINE__ << ") " << __FILE__ << COLOR_RESET << std::endl : std::cerr)
#else
#define dbg(x) ((void)0)
#define dbgif(cond, x) ((void)0)
#endif
// https://hitonanode.github.io/cplib-cpp/number/factorize.hpp
#include <algorithm>
#include <cassert>
#include <cstdint>
#include <numeric>
#include <vector>
namespace SPRP {
// http://miller-rabin.appspot.com/
const std::vector<std::vector<__int128>> bases{
{126401071349994536}, // < 291831
{336781006125, 9639812373923155}, // < 1050535501 (1e9)
{2, 2570940, 211991001, 3749873356}, // < 47636622961201 (4e13)
{2, 325, 9375, 28178, 450775, 9780504, 1795265022} // <= 2^64
};
inline int get_id(long long n) {
if (n < 291831) {
return 0;
} else if (n < 1050535501) {
return 1;
} else if (n < 47636622961201)
return 2;
else { return 3; }
}
} // namespace SPRP
// Miller-Rabin primality test
// https://ja.wikipedia.org/wiki/%E3%83%9F%E3%83%A9%E3%83%BC%E2%80%93%E3%83%A9%E3%83%93%E3%83%B3%E7%B4%A0%E6%95%B0%E5%88%A4%E5%AE%9A%E6%B3%95
// Complexity: O(lg n) per query
inline struct {
long long modpow(__int128 x, __int128 n, long long mod) noexcept {
__int128 ret = 1;
for (x %= mod; n; x = x * x % mod, n >>= 1) ret = (n & 1) ? ret * x % mod : ret;
return ret;
}
bool operator()(long long n) noexcept {
if (n < 2) return false;
if (n % 2 == 0) return n == 2;
int s = __builtin_ctzll(n - 1);
for (__int128 a : SPRP::bases[SPRP::get_id(n)]) {
if (a % n == 0) continue;
a = modpow(a, (n - 1) >> s, n);
bool may_composite = true;
if (a == 1) continue;
for (int r = s; r--; a = a * a % n) {
if (a == n - 1) may_composite = false;
}
if (may_composite) return false;
}
return true;
}
} is_prime;
// https://hitonanode.github.io/cplib-cpp/random/rand_nondeterministic.hpp
#include <chrono>
#include <random>
struct rand_int_ {
using lint = long long;
std::mt19937 mt;
rand_int_() : mt(42) {}
// rand_int_() : mt(std::chrono::steady_clock::now().time_since_epoch().count()) {}
lint operator()(lint x) { return this->operator()(0, x); } // [0, x)
lint operator()(lint l, lint r) {
std::uniform_int_distribution<lint> d(l, r - 1);
return d(mt);
}
} rnd;
#include <boost/multiprecision/cpp_int.hpp>
using mpint = boost::multiprecision::cpp_int;
lint power(lint x, mpint n, lint MOD)
{
lint ans = 1;
while (n>0)
{
if (n & 1) (ans *= x) %= MOD;
(x *= x) %= MOD;
n >>= 1;
}
return ans;
}
// https://hitonanode.github.io/cplib-cpp/number/euler_totient_phi.hpp
#include <numeric>
#include <vector>
// Euler's totient function
// Complexity: O(NlgN)
std::vector<int> euler_phi_table(int N) {
std::vector<int> ret(N + 1);
std::iota(ret.begin(), ret.end(), 0);
for (int p = 2; p <= N; p++) {
if (ret[p] == p) {
ret[p] = p - 1;
for (int i = p * 2; i <= N; i += p) ret[i] = ret[i] / p * (p - 1);
}
}
return ret;
}
// Euler's totient function
// Complexity: O(sqrtN)
// https://csacademy.com/contest/fii-code-2021-round-2/task/clown-fiesta/
// https://kopricky.github.io/code/Computation_Advanced/tetration.html
template <typename Int> Int euler_phi(Int n) {
Int ret = n;
if (n % 2 == 0) {
ret /= 2;
do { n /= 2; } while (n % 2 == 0);
}
if (n % 3 == 0) {
ret = ret / 3 * 2;
do { n /= 3; } while (n % 3 == 0);
}
for (Int i = 5; i * i <= n; i += 4) {
if (n % i == 0) {
ret = ret / i * (i - 1);
do { n /= i; } while (n % i == 0);
}
i += 2;
if (n % i == 0) {
ret = ret / i * (i - 1);
do { n /= i; } while (n % i == 0);
}
}
if (n != 1) ret = ret / n * (n - 1);
return ret;
}
lint B_tmp = 2;
lint Ask(lint i, int j) {
cout << "? " << i << ' ' << j << endl;
lint resp = 1 % B_tmp;
cin >> resp;
// REP(_, j) resp = resp * i % B_tmp;
return resp;
}
int main() {
auto Solve = [&](int N, int M, int L) -> lint {
if (L == 0) return Ask(N, 1);
if (L == 1) return Ask(N, M);
if (M == 0 and L > 0) return Ask(N, 0);
if (N == 0) return Ask(0, 1);
lint ng = 0;
lint ok = 1e9;
while (abs(ok - ng) > 1) {
const lint mid = (ok + ng) / 2;
auto a = Ask(mid, 1);
(a == mid ? ng : ok) = mid;
}
if (ok == 1) return 0;
// REP(_, 50) {
// int p = 4;
// while (!is_prime(p)) p = rnd(0.9e5, 1e5);
// // dbg(p);
// lint m = lint(p) * p - Ask(p, 2);
// gcd_b = gcd(gcd_b, m);
// }
// dbg(gcd_b);
mpint d = 1;
REP(_, L) d *= M;
lint ans = power(N, d, ok);
return ans;
};
// const int U = 30;
// FOR(N, 1, U) FOR(M, 1, U) FOR(L, 1, U) FOR(B, 1, U) {
// B_tmp = B;
// auto ans = Solve(N, M, L);
// auto phi = euler_phi(B);
// const lint e = power(M, L, phi);
// auto ans2 = B == 1 ? 0 : power(N, e, B);
// if (N == 0) {
// if (M == 0 and L == 0) ans2 = 0;
// else if (M) {
// ans2 = 0;
// } else if (L) {
// ans2 = 1 % B;
// }
// }
// if (ans != ans2) {
// dbg(make_tuple(N, M, L, B, ans, ans2, phi, e));
// exit(1);
// }
// // REP(_, L) e = e * M % phi;
// }
int N, M, L;
cin >> N >> M >> L;
auto ans = Solve(N, N, L);
cout << "! " << ans << '\n';
}