This documentation is automatically generated by online-judge-tools/verification-helper
View the Project on GitHub plasmatic1/templates
#define PROBLEM "https://judge.u-aizu.ac.jp/onlinejudge/description.jsp?id=DPL_5_A" #include "../../template.hpp" #include "../test_utils.hpp" #include "../../math/mod.hpp" using MI = ModInt<int, 1000000007>; int main() { fast_io(); int N = readi(), K = readi(); print(MI::pow(K, N)); }
#line 1 "tests/math/modint.test.cpp" #define PROBLEM "https://judge.u-aizu.ac.jp/onlinejudge/description.jsp?id=DPL_5_A" #line 2 "template.hpp" #include <bits/stdc++.h> #define DEBUG 1 using namespace std; // Defines #define fs first #define sn second #define pb push_back #define eb emplace_back #define mpr make_pair #define mtp make_tuple #define all(x) (x).begin(), (x).end() // Basic type definitions #if __cplusplus == 201703L // CPP17 only things template <typename T> using opt_ref = optional<reference_wrapper<T>>; // for some templates #endif using ll = long long; using ull = unsigned long long; using ld = long double; using pii = pair<int, int>; using pll = pair<long long, long long>; #ifdef __GNUG__ // PBDS order statistic tree #include <ext/pb_ds/assoc_container.hpp> // Common file #include <ext/pb_ds/tree_policy.hpp> using namespace __gnu_pbds; template <typename T, class comp = less<T>> using os_tree = tree<T, null_type, comp, rb_tree_tag, tree_order_statistics_node_update>; template <typename K, typename V, class comp = less<K>> using treemap = tree<K, V, comp, rb_tree_tag, tree_order_statistics_node_update>; // HashSet #include <ext/pb_ds/assoc_container.hpp> template <typename T, class Hash> using hashset = gp_hash_table<T, null_type, Hash>; template <typename K, typename V, class Hash> using hashmap = gp_hash_table<K, V, Hash>; const ll RANDOM = chrono::high_resolution_clock::now().time_since_epoch().count(); struct chash { ll operator()(ll x) const { return x ^ RANDOM; } }; #endif // More utilities int SZ(string &v) { return v.length(); } template <typename C> int SZ(C &v) { return v.size(); } template <typename C> void UNIQUE(vector<C> &v) { sort(v.begin(), v.end()); v.resize(unique(v.begin(), v.end()) - v.begin()); } template <typename T, typename U> void maxa(T &a, U b) { a = max(a, b); } template <typename T, typename U> void mina(T &a, U b) { a = min(a, b); } const ll INF = 0x3f3f3f3f, LLINF = 0x3f3f3f3f3f3f3f3f; #line 3 "tests/test_utils.hpp" // I/O template <typename T> void print(T v) { cout << v << '\n'; } template <typename T, typename... Rest> void print(T v, Rest... vs) { cout << v << ' '; print(vs...); } void fast_io() { ios_base::sync_with_stdio(false); cin.tie(NULL); } // Reading operators template <typename T, typename U> istream& operator>>(istream& in, pair<T, U> &o) { return in >> o.first >> o.second; } // Read helpers int readi() { int x; cin >> x; return x; } ll readl() { ll x; cin >> x; return x; } template <typename T> vector<T> readv(int n) { vector<T> res(n); for (auto &x : res) cin >> x; return res; } // Functional stuff template <typename T> vector<pair<int, T>> enumerate(vector<T> v, int start = 0) { vector<pair<int, T>> res; for (auto &x : v) res.emplace_back(start++, x); return res; } #line 3 "math/eea.hpp" /* * ax + by = gcd(a, b) * * we know * bx' + (a%b)y' = gcd(a, b) * * bx' + (a-b*(a//b))y' = gcd(a, b) * bx' + ay' - b*(a//b)y' = gcd(a, b) * ay' + b(x' - (a//b)y') = gcd(a, b) */ template <typename T> T extgcd(T a, T b, T &x, T &y) { if (b == 0) { x = 1; y = 0; return a; } T x0, y0, res = extgcd(b, a%b, x0, y0); x = y0; y = x0 - (a / b) * y0; return res; } #line 4 "math/mod.hpp" // based on Tourist modInt orz template <typename MD> struct _ModInt { using T = typename decay<decltype(MD::value)>::type; static_assert(sizeof(T) >= 4, "size of T must be at least 32 bits"); static_assert(sizeof(T) <= 8, "size of T must be at most 64 bits"); static_assert(is_integral<T>::value, "T must be an integral type"); #ifdef __SIZEOF_INT128__ using mul_t = typename conditional<sizeof(T) <= 4, int64_t, __int128>::type; #else using mul_t = int64_t; static_assert(sizeof(T) <= 4, "int128 not available, cannot use 64-bit size of T"); #endif constexpr static T mod() { return MD::value; } template <typename U> static T normalize(const U& x) { T res = x; res %= mod(); if (res < 0) res += mod(); return res; } T value; constexpr _ModInt() : value() {} template <typename U> _ModInt(const U& x) { value = normalize(x); } const T& operator()() const { return value; } template <typename U> operator U() const { return static_cast<U>(value); } // FastPow template <typename U> static _ModInt pow(_ModInt x, U y) { _ModInt res(1); for (; y; y /= 2) { if (y & 1) res *= x; x *= x; } return res; } static _ModInt inv(const _ModInt &x) { T inv, _; extgcd(x.value, mod(), inv, _); return _ModInt(inv); } // Arithmetic Operators w/ _ModInt // Assignment operators here _ModInt& operator+=(const _ModInt &o) { if ((value += o.value) >= mod()) value -= mod(); return *this; } template <typename U> _ModInt& operator+=(const U &o) { return *this += _ModInt(o); } _ModInt& operator-=(const _ModInt &o) { if ((value -= o.value) < 0) value += mod(); return *this; } template <typename U> _ModInt& operator-=(const U &o) { return *this -= _ModInt(o); } _ModInt& operator++() { return *this += 1; } _ModInt operator++(int) { _ModInt res(*this); *this += 1; return res; } _ModInt& operator--() { return *this -= 1; } _ModInt operator--(int) { _ModInt res(*this); *this -= 1; return res; } _ModInt& operator*=(const _ModInt &o) { value = (mul_t)value * o.value % mod(); if (value < 0) value += mod(); return *this; } // make sure cast to mul_t!!! template <typename U> _ModInt& operator*=(const U &o) { return *this *= _ModInt(o); } _ModInt& operator/=(const _ModInt &o) { return *this *= inv(o.value); } template <typename U> _ModInt& operator/=(const U &o) { return *this /= _ModInt(o); } _ModInt operator-() const { return _ModInt(value); } // Other Operators T& operator()() { return value; } // Definitions of some operators }; // Binary operators #define OP_CMP(op) template <typename T> bool operator op(const _ModInt<T> &lhs, const _ModInt<T> &rhs) { return lhs.value op rhs.value; } \ template <typename T, typename U> bool operator op(const _ModInt<T> &lhs, U rhs) { return lhs op _ModInt<T>(rhs); } \ template <typename T, typename U> bool operator op(U lhs, const _ModInt<T> &rhs) { return _ModInt<T>(lhs) op rhs; } #define OP_ARI(op) template <typename T> _ModInt<T> operator op(const _ModInt<T> &lhs, const _ModInt<T> &rhs) { return _ModInt<T>(lhs) op##= rhs; } \ template <typename T, typename U> _ModInt<T> operator op(U lhs, const _ModInt<T> &rhs) { return _ModInt<T>(lhs) op##= rhs; } \ template <typename T, typename U> _ModInt<T> operator op(const _ModInt<T> &lhs, U rhs) { return _ModInt<T>(lhs) op##= rhs; } OP_CMP(==) OP_CMP(!=) OP_CMP(<) OP_CMP(>) OP_CMP(<=) OP_CMP(>=) OP_ARI(+) OP_ARI(-) OP_ARI(*) OP_ARI(/) #undef OP_CMP #undef OP_ARI template <typename T> istream& operator>>(istream& in, _ModInt<T> &o) { return in >> o(); } template <typename T> ostream& operator<<(ostream& out, _ModInt<T> &o) { return out << o(); } // Definitions template <typename T, T mod> using ModInt = _ModInt<integral_constant<T, mod>>; template <typename T> struct VarMod { static T value; static void read(istream& in) { in >> value; } static void set(T v0) { value = v0; } }; template <typename T> using VarModInt = _ModInt<VarMod<T>>; #line 5 "tests/math/modint.test.cpp" using MI = ModInt<int, 1000000007>; int main() { fast_io(); int N = readi(), K = readi(); print(MI::pow(K, N)); }