結果

問題 No.3187 Mingle
ユーザー The Forsaking
提出日時 2025-06-21 19:48:09
言語 C++17(clang)
(17.0.6 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 1,482 bytes
コンパイル時間 4,563 ms
コンパイル使用メモリ 170,204 KB
実行使用メモリ 175,204 KB
最終ジャッジ日時 2025-06-21 19:48:19
合計ジャッジ時間 9,001 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 2 TLE * 1 -- * 27
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp:21:25: warning: operator '>>' has lower precedence than '+'; '+' will be evaluated first [-Wshift-op-parentheses]
   21 |             int mid = l + r >> 1, v = i / mid < 175;
      |                       ~~^~~ ~~
main.cpp:21:25: note: place parentheses around the '+' expression to silence this warning
   21 |             int mid = l + r >> 1, v = i / mid < 175;
      |                         ^
      |                       (    )
1 warning generated.

ソースコード

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

int f[300086], p, s[300086][175], inv[300086];
int main() {
    map<int, int> ma;
    cin >> n >> p;
    inv[1] = 1;
    for (int i = 2; i <= 300000; i++) inv[i] = (ll)(p - p / i) * inv[p % i] % p;
    for (int i = 3; i <= n; i++) {
        ll c = 0;
        int l = 1, r = i;
        while (l < r) {
            int mid = l + r >> 1, v = i / mid < 175;
            r = v ? mid : r;
            l = v ? l : mid + 1;
        }
        for (int j = 1; j < l; j++) {
            int t = i - (i % j);
            c += t == i;
            f[i] = f[i] + (t != i ? f[t] : 0);
            f[i] = f[i] >= p ? f[i] - p : f[i];
        }
        while (l <= i) {
            int g = i / l, r = i / g, pr = i % l, pl = pr - (r - l) * g;
            c += pl == 0, pl = pl == 0 ? g : pl;
            int t = (pl <= pr ? s[i - pl][g] - s[i - pr][g] + f[i - pr] : 0);
            t += t >= p ? -p : t < 0 ? p : 0;
            f[i] = f[i] + t;
            f[i] = f[i] >= p ? f[i] - p : f[i];
            l = r + 1;
        }
        f[i] = ((ll)f[i] * inv[i] + 1) % p;
        f[i] = (ll)f[i] * inv[i - c] % p * i % p;
        for (int j = 1; j < 175; j++)
            s[i][j] = i >= j ? s[i - j][j] + f[i] >= p ? s[i - j][j] + f[i] - p : s[i - j][j] + f[i] : f[i];
    }
    printf("%d\n", f[n]);
    return 0;
}
0