結果
問題 |
No.1573 Divisor Function
|
ユーザー |
|
提出日時 | 2021-06-18 21:23:25 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 14 ms / 2,000 ms |
コード長 | 537 bytes |
コンパイル時間 | 3,808 ms |
コンパイル使用メモリ | 229,388 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-06-25 11:16:29 |
合計ジャッジ時間 | 5,400 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 43 |
ソースコード
#include <bits/stdc++.h> #include <atcoder/all> using namespace std; using namespace atcoder; using ll=long long; ll n,m; modint998244353 ans=0,j_sum; //解法 Σ[j=1→m]j×(n/j)×(n/j+3)/2 n/jの値で場合分け void clac(ll x) { // j_sum=n/jがxとなるjの和 j_sum=(min(m+1,n/(x+1)+1)+min(m,n/x))*(min(m,n/x)-min(m,n/(x+1)))/2; ans+=j_sum*x*(x+3)/2; } int main() { cin>>n>>m; for(ll i=1;i*i<=n;i++) { // n/j=iの場合 clac(i); // n/j=n/iの場合 if(n/i!=i)clac(n/i); } cout<<ans.val()<<endl; }