結果
問題 | No.681 Fractal Gravity Glue |
ユーザー | fine |
提出日時 | 2018-04-28 01:08:22 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 1,430 bytes |
コンパイル時間 | 1,535 ms |
コンパイル使用メモリ | 168,188 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-06-27 23:10:04 |
合計ジャッジ時間 | 2,219 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
5,248 KB |
testcase_01 | AC | 2 ms
5,376 KB |
testcase_02 | AC | 2 ms
5,376 KB |
testcase_03 | AC | 2 ms
5,376 KB |
testcase_04 | AC | 2 ms
5,376 KB |
testcase_05 | AC | 2 ms
5,376 KB |
testcase_06 | AC | 2 ms
5,376 KB |
testcase_07 | AC | 2 ms
5,376 KB |
testcase_08 | AC | 2 ms
5,376 KB |
testcase_09 | AC | 2 ms
5,376 KB |
testcase_10 | AC | 2 ms
5,376 KB |
testcase_11 | AC | 2 ms
5,376 KB |
testcase_12 | AC | 2 ms
5,376 KB |
testcase_13 | AC | 2 ms
5,376 KB |
testcase_14 | AC | 2 ms
5,376 KB |
testcase_15 | AC | 2 ms
5,376 KB |
testcase_16 | AC | 2 ms
5,376 KB |
testcase_17 | AC | 2 ms
5,376 KB |
testcase_18 | AC | 2 ms
5,376 KB |
testcase_19 | AC | 1 ms
5,376 KB |
ソースコード
#include <bits/stdc++.h>using namespace std;using ll = long long;const ll MOD = 1000000007;ll modpow(ll x, ll n, ll mod = MOD) {ll res = 1;while (n > 0) {if (n & 1) res = res * x % mod;x = x * x % mod;n >>= 1;}return res;}ll calc(int i, int d) {ll tmp = modpow(d + 1, i - 1);ll res = (tmp + MOD - 1) % MOD * modpow(d, MOD - 2) % MOD;(res += tmp * (d + 2) % MOD) %= MOD;(res += (MOD - i - 1)) %= MOD;return res;}ll G[50], N[50];ll solve(int idx, ll num) {if (idx == 0) return 0;ll mod = num % (N[idx - 1] + 1);if (mod == 0) {return (G[idx - 1] + idx) * (num / (N[idx - 1] + 1)) % MOD;} else if (mod == N[idx - 1]) {return ((G[idx - 1] + idx) * (num / (N[idx - 1] + 1)) % MOD + G[idx - 1]) % MOD;} else {return ((G[idx - 1] + idx) * (num / (N[idx - 1] + 1)) % MOD + solve(idx - 1, mod)) % MOD;}}int main() {cin.tie(0);ios::sync_with_stdio(false);int n, b, d;cin >> n >> b >> d;int lim = 0;N[0] = 1;for (int i = 0; ; i++) {if (N[i] - 1 > n) {N[i]--;break;}lim++;N[i + 1] = N[i] * (d + 1);N[i]--;}for (int i = 1; i < lim; i++) {G[i] = calc(i, d);}ll ans = calc(b, d);(ans += MOD - solve(lim, n)) %= MOD;cout << ans << endl;return 0;}