結果
問題 | No.2891 Mint |
ユーザー | srjywrdnprkt |
提出日時 | 2024-10-02 15:48:44 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,181 bytes |
コンパイル時間 | 2,400 ms |
コンパイル使用メモリ | 203,364 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-10-02 15:48:53 |
合計ジャッジ時間 | 4,699 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 1 ms
5,248 KB |
testcase_01 | AC | 2 ms
5,248 KB |
testcase_02 | AC | 2 ms
5,248 KB |
testcase_03 | AC | 2 ms
5,248 KB |
testcase_04 | WA | - |
testcase_05 | AC | 2 ms
5,248 KB |
testcase_06 | AC | 46 ms
5,248 KB |
testcase_07 | AC | 2 ms
5,248 KB |
testcase_08 | AC | 2 ms
5,248 KB |
testcase_09 | AC | 2 ms
5,248 KB |
testcase_10 | AC | 2 ms
5,248 KB |
testcase_11 | WA | - |
testcase_12 | AC | 2 ms
5,248 KB |
testcase_13 | AC | 2 ms
5,248 KB |
testcase_14 | AC | 2 ms
5,248 KB |
testcase_15 | AC | 2 ms
5,248 KB |
testcase_16 | AC | 2 ms
5,248 KB |
testcase_17 | AC | 35 ms
5,248 KB |
testcase_18 | AC | 20 ms
5,248 KB |
testcase_19 | AC | 19 ms
5,248 KB |
testcase_20 | AC | 30 ms
5,248 KB |
testcase_21 | AC | 20 ms
5,248 KB |
testcase_22 | AC | 39 ms
5,248 KB |
testcase_23 | AC | 15 ms
5,248 KB |
testcase_24 | AC | 31 ms
5,248 KB |
testcase_25 | AC | 41 ms
5,248 KB |
testcase_26 | AC | 41 ms
5,248 KB |
testcase_27 | AC | 39 ms
5,248 KB |
testcase_28 | AC | 44 ms
5,248 KB |
testcase_29 | AC | 19 ms
5,248 KB |
testcase_30 | AC | 33 ms
5,248 KB |
testcase_31 | AC | 21 ms
5,248 KB |
testcase_32 | AC | 36 ms
5,248 KB |
testcase_33 | AC | 25 ms
5,248 KB |
testcase_34 | AC | 29 ms
5,248 KB |
testcase_35 | AC | 32 ms
5,248 KB |
testcase_36 | AC | 33 ms
5,248 KB |
testcase_37 | WA | - |
testcase_38 | WA | - |
testcase_39 | WA | - |
testcase_40 | WA | - |
testcase_41 | WA | - |
testcase_42 | AC | 2 ms
5,248 KB |
testcase_43 | AC | 2 ms
5,248 KB |
testcase_44 | AC | 2 ms
5,248 KB |
testcase_45 | AC | 2 ms
5,248 KB |
testcase_46 | AC | 2 ms
5,248 KB |
testcase_47 | WA | - |
testcase_48 | WA | - |
testcase_49 | WA | - |
testcase_50 | WA | - |
testcase_51 | WA | - |
testcase_52 | WA | - |
testcase_53 | WA | - |
testcase_54 | WA | - |
testcase_55 | WA | - |
testcase_56 | WA | - |
ソースコード
#include <bits/stdc++.h> #include <atcoder/modint> using namespace std; using ll = long long; using namespace atcoder; using mint = modint998244353; //calculate floor(sqrt(s)) ll isqrt(ll s){ ll l=0, r=3e9, c; while(r-l>1){ c = (l+r)/2; if (c*c <= s) l = c; else r = c; } return l; } int main(){ cin.tie(nullptr); ios_base::sync_with_stdio(false); /* sum(i=1, N)(M-(M/k)k)=MN-sum(i=1,N)(M/k)k M/iがsqrt(M-1)以下のiはO(sqrt(M))個しかないので愚直にやる。 1<=j<=sqrt(M)に対して、M/k=jとなるkをまとめて j*(kの総和) j<=M/k<j+1 M/(j+1)<k<=M/j L = floor(M/(j+1))+1 R = max(floor(M/j), N) */ ll N, M, K, L, R; cin >> N >> M; mint ans=0, S, iv=mint(2).inv(); K = isqrt(M); for (int i=1; i<=1e9; i++){ if (M/i <= K) break; ans += (M/i)*i; } for (int j=K; j>=1; j--){ L = M/(j+1)+1; R = min(M/j, N); if (L>R) continue; //L~Rの和 S = mint(R)*(R+1)*iv-mint(L)*(L-1)*iv; ans += S*j; } cout << (mint(N)*M-ans).val() << endl; return 0; }