結果
| 問題 | 
                            No.1573 Divisor Function
                             | 
                    
| コンテスト | |
| ユーザー | 
                             | 
                    
| 提出日時 | 2021-06-17 20:38:46 | 
| 言語 | C++14  (gcc 13.3.0 + boost 1.87.0)  | 
                    
| 結果 | 
                             
                                WA
                                 
                             
                            
                         | 
                    
| 実行時間 | - | 
| コード長 | 527 bytes | 
| コンパイル時間 | 3,930 ms | 
| コンパイル使用メモリ | 229,856 KB | 
| 実行使用メモリ | 6,948 KB | 
| 最終ジャッジ日時 | 2024-06-25 11:15:52 | 
| 合計ジャッジ時間 | 5,610 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge4 / judge2 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 3 | 
| other | AC * 32 WA * 11 | 
ソースコード
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;
using ll=long long;
ll n,m;
modint998244353 ans=0,sum;
void f(ll x)
{
  sum=(min(m,n/x)+min(m,n/(x+1)+1))*(min(m,n/x)-min(m,n/(x+1)))/2;// n/j=iとなるjの和
  ans+=sum*x*(x+3)/2;
}
int main()
{
  /*解法
  	Σ[j=1→m]j×(n/j)×(n/j+3)/2
    
    n=20,m=10
    20 10 6 5 4 3 2 2 2 2
  */
  cin>>n>>m;
  for(ll i=1;i*i<=n;i++)
  {
    // n/j=iの場合
    f(i);
    // n/j=n/iの場合
    f(n/i);
  }
  cout<<ans.val()<<endl;
}