結果
| 問題 |
No.1573 Divisor Function
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2021-06-18 21:52:21 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 14 ms / 2,000 ms |
| コード長 | 537 bytes |
| コンパイル時間 | 4,254 ms |
| コンパイル使用メモリ | 228,856 KB |
| 実行使用メモリ | 5,376 KB |
| 最終ジャッジ日時 | 2024-06-25 11:16:35 |
| 合計ジャッジ時間 | 5,420 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| 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 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;
}