結果

問題 No.3187 Mingle
ユーザー 蜜蜂
提出日時 2025-06-13 14:25:24
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 967 bytes
コンパイル時間 907 ms
コンパイル使用メモリ 96,068 KB
実行使用メモリ 418,176 KB
最終ジャッジ日時 2025-06-13 14:38:43
合計ジャッジ時間 46,052 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 9 WA * 21
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
using namespace std;

#include <atcoder/modint>
using namespace atcoder;

using mint = modint;

int main(){
  int n;
  cin >> n;
  long long p;
  cin >> p;

  mint::set_mod(p);

  int m = 350;
  vector<int> d(n + 1, 0);
  vector<mint> f(n + 1, 0);
  vector<vector<mint>> f_sum(m + 1, vector<mint>(n + 1, 0));

  for(int i = 1; i <= n; i++){
    for(int j = i; j <= n; j += i){
      d[j]++;
    }
  }

  for(int i = 3; i <= n; i++){
    mint sum = 0;

    for(int j = 1; j <= m; j++){
      if(i % j != 0) sum += f[i - i % j];
    }
    for(int j = 1; j <= m; j++){
      int l = i / (j + 1) + 1;
      int r = i / j;
      if(r * j == i) r--;
      l = max(l, m + 1);
      if(l > r) break;
      sum += f_sum[j][r * j] - f_sum[j][(l - 1) * j];
    }
    f[i] = (sum / i + 1) * i / (i - d[i]);
    for(int j = 1; j <= m; j++){
      if(i % j == 0) f_sum[j][i] = f_sum[j][i - j] + f[i];
    }
  }

  cout << f[n].val() << endl;
}
0