結果
問題 |
No.681 Fractal Gravity Glue
|
ユーザー |
![]() |
提出日時 | 2020-03-31 18:57:06 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 32 ms / 2,000 ms |
コード長 | 809 bytes |
コンパイル時間 | 286 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 10,752 KB |
最終ジャッジ日時 | 2024-06-24 23:26:37 |
合計ジャッジ時間 | 1,607 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 20 |
ソースコード
#!/usr/bin/ python3.8 import sys read = sys.stdin.buffer.read readline = sys.stdin.buffer.readline readlines = sys.stdin.buffer.readlines from functools import lru_cache MOD = 10**9 + 7 N, B, D = map(int, read().split()) counts = [0] * 64 sizes = [0] * 64 for i in range(1, 64): counts[i] = counts[i - 1] * (D + 1) + D sizes[i] = sizes[i - 1] * (D + 1) + i * D @lru_cache(None) def size(n): return ((pow(D + 1, n, MOD) - 1) * (1 + pow(D, MOD - 2, MOD)) - n) % MOD for i in range(64): assert sizes[i] % MOD == size(i) def f(N): """下 $N$ 個の大きさの合計""" ret = 0 B = 63 while N: q, N = divmod(N, counts[B] + 1) ret += q * (sizes[B] + B + 1) B -= 1 return ret full = size(B) done = f(N) answer = (full - done) % MOD print(answer)