結果

問題 No.3187 Mingle
ユーザー 蜜蜂
提出日時 2025-06-02 22:39:36
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 582 ms / 2,500 ms
コード長 638 bytes
コンパイル時間 1,000 ms
コンパイル使用メモリ 95,528 KB
実行使用メモリ 36,968 KB
最終ジャッジ日時 2025-06-13 14:38:04
合計ジャッジ時間 13,215 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

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);

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

  vector<mint> f(n + 1, 0);
  mint tot = 0;
  for(int i = 3; i <= n; i++){
    mint sum = tot;
    for(int j: d[i]) sum -= f[j];
    f[i] = (sum / i + 1) * i / (i - d[i].size());
    tot = sum + f[i] * d[i].size();
    for(int j: d[i]) f[j] = f[i];
  }
  cout << f[n].val() << endl;
}
0