結果
問題 | No.1573 Divisor Function |
ユーザー | どらら |
提出日時 | 2021-06-27 15:15:06 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 21 ms / 2,000 ms |
コード長 | 833 bytes |
コンパイル時間 | 4,031 ms |
コンパイル使用メモリ | 229,700 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-06-25 11:46:22 |
合計ジャッジ時間 | 5,683 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 43 |
ソースコード
#include <bits/stdc++.h> #include <atcoder/all> using namespace std; using namespace atcoder; #define REP(i,a,n) for(int i=(a); i<(int)(n); i++) #define rep(i,n) REP(i,0,n) #define FOR(it,c) for(__typeof((c).begin()) it=(c).begin(); it!=(c).end(); ++it) #define ALLOF(c) (c).begin(), (c).end() typedef long long ll; typedef unsigned long long ull; //using mint = modint1000000007; using mint = modint998244353; int main(){ int n, m; cin >> n >> m; if(n < m) m = n; mint ret(0); int lhs = 1; while(lhs <= m){ int rhs = n / (n/lhs); if(rhs >= m) rhs = m; mint a = (rhs+1-lhs); a *= (lhs+rhs); a /= 2; int x = n / lhs; mint mx(x); mint b = (mx*mx); b += mx * 3; b /= 2; ret += a * b; lhs = rhs + 1; } cout << ret.val() << endl; return 0; }