結果

問題 No.681 Fractal Gravity Glue
コンテスト
ユーザー ei1333333
提出日時 2018-04-28 00:09:08
言語 C++17(gcc12)
(gcc 12.3.0 + boost 1.89.0)
結果
MLE  
実行時間 -
コード長 1,139 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 6,666 ms
コンパイル使用メモリ 254,412 KB
最終ジャッジ日時 2025-01-05 10:26:53
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 8 TLE * 5 MLE * 7
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<bits/stdc++.h>

using namespace std;

using int64 = long long;
const int mod = 1e9 + 7;

int64 N, B, D;

int64 f(int idx) {
  int64 base = 1;
  for(int i = 0; i < idx; i++) {
    base *= D + 1;
    if(base > mod) break;
  }
  return min(base - 1, 1LL << 30);
}


int64 mod_pow(int64 x, int64 n, int64 mod) {
  int64 ret = 1;
  while(n > 0) {
    if(n & 1) (ret *= x) %= mod;
    (x *= x) %= mod;
    n >>= 1;
  }
  return ret;
}

// あー人生ね　あー
int64 g(int idx) {
  int64 base = mod_pow(D + 1, idx + 1, mod);
  (base = mod - base) %= mod;
  base += D * idx % mod;
  base %= mod;
  base += D;
  base %= mod;
  base += 1;
  base %= mod;
  base *= mod_pow(D, mod - 2, mod);
  base %= mod;
  (base = mod - base) %= mod;
  return base;
}


int64 dfs(int idx, int pos) {
  if(pos == 0) return 0;
  if(idx == 1) return pos;
  int64 basesz = f(idx - 1) + 1;
  int64 sz = pos / basesz;
  int64 ret = 1LL * sz * idx % mod;
  ret += g(idx - 1) * sz % mod;
  ret %= mod;
  ret += dfs(idx - 1, pos - sz * basesz) % mod;
  return ret % mod;
}


int main() {
  cin >> N >> B >> D;
  cout << (g(B) + mod - dfs(B, N)) % mod << endl;
}
0