結果
| 問題 |
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;
}