結果

問題 No.3187 Mingle
コンテスト
ユーザー The Forsaking
提出日時 2025-06-21 17:51:13
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 1,297 bytes
コンパイル時間 2,014 ms
コンパイル使用メモリ 194,364 KB
実行使用メモリ 136,864 KB
最終ジャッジ日時 2025-06-21 17:51:23
合計ジャッジ時間 8,282 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 2 TLE * 1 -- * 27
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

typedef pair<int, int> pii;
typedef long long ll;
const int N = 2000008, MOD = 998244353, INF = 0x3f3f3f3f;
ll res;
int n, m, cnt, w[N];

inline ll qmi(ll a, ll b, ll c) { ll res = 1; while (b) { if (b & 1) res = res * a % c; a = a * a % c; b >>= 1; } return res; }

ll f[300086], p;
int s[300086][333];
int main() {
    cin >> n >> p;
    for (int i = 3; i <= n; i++) {
        int l = 1;
        ll c = 0;
        while (l < i) {
            int r = i / (i / l);
            int g = i / l, pr = i % l, pl = pr - (r - l) * g;
            if (pl == 0) pl += g, c++;
            if (pl <= pr) {
                if (g < 333) {
                    f[i] = (f[i] + (s[i - pl][g] - s[i - pr][g] + f[i - pr] + p)) % p;
                } else {
                    for (int j = i - pr; j <= i - pl; j += g) {
                        f[i] = f[i] + f[j];
                        f[i] = f[i] > p ? f[i] - p : f[i];
                    }
                }
            }
            l = r + 1;
        }
        f[i] = (f[i] * qmi(i, p - 2, p) + 1) % p;
        f[i] = f[i] * qmi(i - c, p - 2, p) % p * i % p;
        for (int j = 1; j <= 332; j++)
            s[i][j] = i >= j ? (s[i - j][j] + f[i]) % p : f[i];
    }
    printf("%lld\n", f[n]);
    return 0;
}
0