結果
問題 | No.673 カブトムシ |
ユーザー |
![]() |
提出日時 | 2020-09-04 01:53:52 |
言語 | C++17(gcc12) (gcc 12.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 2,749 bytes |
コンパイル時間 | 2,862 ms |
コンパイル使用メモリ | 91,776 KB |
最終ジャッジ日時 | 2025-01-14 04:20:45 |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 12 WA * 2 |
ソースコード
#pragma GCC target("avx2")#pragma GCC optimize("O3")#pragma GCC optimize("unroll-loops")#include<iostream>using namespace std;using i32 = int_fast32_t;using i64 = int_fast64_t;#define rep(i, n) for (i32 i = 0; i < (i32)(n); i++)#define all(a) (a).begin(),(a).end()#define rall(a) (a).rbegin(),(a).rend()using P = pair<i64,i64>;using i64 = int_fast64_t;constexpr i64 MAX = 10000000;constexpr i64 MOD = 1000000007;i64 fac[MAX], finv[MAX], inv[MAX];template <i64 modulus>class modcal{public:i64 a;constexpr modcal(const i64 x = 0) noexcept : a(x % modulus) {}constexpr i64 &value() noexcept { return a; }constexpr const i64 &value() const noexcept { return a; }constexpr modcal operator+(const modcal rhs) const noexcept{return modcal(*this) += rhs;}constexpr modcal operator-(const modcal rhs) const noexcept{return modcal(*this) -= rhs;}constexpr modcal operator*(const modcal rhs) const noexcept{return modcal(*this) *= rhs;}constexpr modcal operator/(const modcal rhs) noexcept{return modcal(*this) /= rhs;}constexpr modcal &operator+=(const modcal rhs) noexcept{a += rhs.a;if (a >= modulus){a -= modulus;}return *this;}constexpr modcal &operator-=(const modcal rhs) noexcept{if (a < rhs.a){a += modulus;}a -= rhs.a;return *this;}constexpr modcal &operator*=(const modcal rhs) noexcept{a = a * rhs.a % modulus;return *this;}constexpr modcal &operator/=(modcal rhs) noexcept{i64 exp = modulus - 2;while (exp){if (exp % 2){*this *= rhs;}rhs *= rhs;exp /= 2;}return *this;}void COMninit(){fac[0] = fac[1] = 1;finv[0] = finv[1] = 1;inv[1] = 1;for (int i = 2; i < MAX; i++){fac[i] = fac[i - 1] * i % modulus;inv[i] = MOD - inv[modulus % i] * (modulus / i) % modulus;finv[i] = finv[i - 1] * inv[i] % modulus;}}i64 COMn(i64 n, i64 k){if (n < k)return 0;if (n < 0 || k < 0)return 0;return fac[n] * (finv[k] * finv[n - k] % modulus) % modulus;}constexpr modcal<1000000007> modpow(const modcal<1000000007> &a, i64 n){if (n == 0)return 1;auto t = modpow(a, n / 2);t = t * t;if (n & 1)t = t * a;return t;}};using modc = modcal<1000000007>;int main(){ios::sync_with_stdio(false);std::cin.tie(nullptr);i64 b,c,d;cin >> b >> c >> d;if(c == 1){modc ans = b;ans *= 1;ans *= d;cout << ans.a << endl;return 0;}modc ans = c;ans = ans.modpow(c,d);ans -= 1;ans *= b;ans /= (c - 1);ans *= c;cout << ans.a << endl;}