結果

問題 No.1573 Divisor Function
コンテスト
ユーザー sushitoruna
提出日時 2021-06-18 21:52:21
言語 C++14
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=c++14 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 10 ms / 2,000 ms
コード長 537 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,721 ms
コンパイル使用メモリ 248,244 KB
実行使用メモリ 7,720 KB
最終ジャッジ日時 2026-03-11 19:55:23
合計ジャッジ時間 4,186 ms
ジャッジサーバーID
（参考情報） 		judge3_1 / judge1_1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 43
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#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 x=1;x*x<=n;x++)
  {
    // n/j=xの場合
    clac(x);
    // n/j=n/xの場合
    if(n/x!=x)clac(n/x);
  }
  cout<<ans.val()<<endl;
}
0