結果
| 問題 | No.1573 Divisor Function |
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2021-06-17 20:59:28 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.89.0) |
| 結果 |
TLE
|
| 実行時間 | - |
| コード長 | 605 bytes |
| 記録 | |
| コンパイル時間 | 3,916 ms |
| コンパイル使用メモリ | 229,920 KB |
| 実行使用メモリ | 10,752 KB |
| 最終ジャッジ日時 | 2024-06-25 11:16:19 |
| 合計ジャッジ時間 | 7,396 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 TLE * 1 |
| other | -- * 43 |
ソースコード
#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);
}*/
for(ll i=1;i<=m;i++)
{
ans+=i*(n/i)*(n/i+3)/2;
}
cout<<ans.val()<<endl;
}
/*
5 2 1 1
20 10 6 8
*/