結果

問題 No.1573 Divisor Function
ユーザー sushitorunasushitoruna
提出日時 2021-06-17 20:55:22
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 15 ms / 2,000 ms
コード長 516 bytes
コンパイル時間 4,224 ms
コンパイル使用メモリ 228,940 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-06-25 11:15:58
合計ジャッジ時間 5,628 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 43
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;
using ll=long long;
ll n,m;
modint998244353 ans=0,sum;
//解法	Σ[j=1→m]j×(n/j)×(n/j+3)/2
void f(ll x)
{
  ll j_min=min(m+1,n/(x+1)+1),j_max=min(m,n/x);
  sum=(j_min+j_max)*(j_max-j_min+1)/2;// n/j=iとなるjの和
  ans+=sum*x*(x+3)/2;
}
int main()
{
  cin>>n>>m;
  for(ll i=1;i*i<=n;i++)
  {
    // n/j=iの場合
    f(i);
    // n/j=n/iの場合
    if(n/i==i)continue;
    f(n/i);
  }
  cout<<ans.val()<<endl;
}
0