結果

問題 No.2940 Sigma Sigma Div Floor Problem
ユーザー 👑 Nachia
提出日時 2024-10-18 22:27:49
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2 ms / 6,000 ms
コード長 1,001 bytes
コンパイル時間 793 ms
コンパイル使用メモリ 77,384 KB
最終ジャッジ日時 2025-02-24 20:57:17
ジャッジサーバーID
（参考情報） 		judge5 / judge
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <string>
#include <vector>
#include <algorithm>
using i64 = long long;
using u64 = unsigned long long;
#define rep(i,n) for(int i=0; i<int(n); i++)
const i64 INF = 1001001001001001001;
template<typename A> void chmin(A& l, const A& r){ if(r < l) l = r; }
template<typename A> void chmax(A& l, const A& r){ if(l < r) l = r; }
#include <atcoder/modint>
using Modint = atcoder::static_modint<998244353>;
using namespace std;

void testcase(){
    i64 n; cin >> n;
    i64 m = 1;
    while((m+1)*(m+1) <= n) m += 1;
    Modint f1 = 0;
    for(i64 a=1; a<=m; a++){
        i64 b = n / a;
        f1 += Modint(a) * Modint(b) * Modint(b+1);
    }
    f1 -= (Modint(m) * Modint(m+1) / 2).pow(2);
    Modint f2 = 0;
    for(i64 a=1; a<=m; a++){
        f2 += Modint(n / a) * 2;
    }
    f2 -= Modint(m).pow(2);
    Modint ans = f2 * Modint(n+1) - f1;
    cout << ans.val() << endl;
}

int main(){
    ios::sync_with_stdio(false); cin.tie(nullptr);
    testcase();
    return 0;
}
0