結果
問題 | No.2891 Mint |
ユーザー |
👑 ![]() |
提出日時 | 2024-09-13 22:14:04 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 55 ms / 2,000 ms |
コード長 | 1,518 bytes |
コンパイル時間 | 833 ms |
コンパイル使用メモリ | 80,152 KB |
最終ジャッジ日時 | 2025-02-24 07:46:14 |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 54 |
ソースコード
#ifdef NACHIA#define _GLIBCXX_DEBUG#else#define NDEBUG#endif#include <iostream>#include <string>#include <vector>#include <algorithm>using i64 = long long;using u64 = unsigned long long;#define rep(i,n) for(int i=0; i<int(n); i++)const i64 INF = 1001001001001001001;template<typename A> void chmin(A& l, const A& r){ if(r < l) l = r; }template<typename A> void chmax(A& l, const A& r){ if(l < r) l = r; }#include <atcoder/modint>namespace nachia{template<class Int = long long>std::vector<Int> EnumerateQuotients(Int N, bool doInsertZero){std::vector<Int> res;if(doInsertZero) res.push_back(0);Int f = 1;for( ; f*f<N; f++) res.push_back(f);int qp1 = res.size();for(Int k=1; k*f<=N; k++){res.push_back((k&1) ? (N/k) : (res[qp1+k/2-1]/2));}std::reverse(res.begin() + qp1, res.end());return res;}}using Modint = atcoder::static_modint<998244353>;using namespace std;void testcase(){i64 N, M; cin >> N >> M;Modint ans = Modint(M) * Modint(N);auto q = nachia::EnumerateQuotients(M, true);Modint inv2 = Modint(2).inv();rep(i,q.size()-1){i64 l = min(N, q[i]);i64 r = min(N, q[i+1]);i64 f = q[q.size() - 1 - i];Modint d = 0;d += Modint(r) * Modint(r+1) * inv2;d -= Modint(l) * Modint(l+1) * inv2;d *= f;ans -= d;}cout << ans.val() << endl;}int main(){ios::sync_with_stdio(false); cin.tie(nullptr);testcase();return 0;}