結果

問題 No.2940 Sigma Sigma Div Floor Problem
ユーザー 👑 NachiaNachia
提出日時 2024-10-18 22:27:49
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 2 ms / 6,000 ms
コード長 1,001 bytes
コンパイル時間 719 ms
コンパイル使用メモリ 79,288 KB
実行使用メモリ 6,824 KB
最終ジャッジ日時 2024-10-18 22:51:19
合計ジャッジ時間 2,747 ms
ジャッジサーバーID
（参考情報） 		judge8 / judge
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
6,816 KB
testcase_01 AC 2 ms
6,816 KB
testcase_02 AC 2 ms
6,820 KB
testcase_03 AC 2 ms
6,816 KB
testcase_04 AC 2 ms
6,820 KB
testcase_05 AC 2 ms
6,816 KB
testcase_06 AC 2 ms
6,816 KB
testcase_07 AC 2 ms
6,816 KB
testcase_08 AC 1 ms
6,820 KB
testcase_09 AC 2 ms
6,816 KB
testcase_10 AC 1 ms
6,820 KB
testcase_11 AC 1 ms
6,816 KB
testcase_12 AC 2 ms
6,816 KB
testcase_13 AC 2 ms
6,820 KB
testcase_14 AC 2 ms
6,820 KB
testcase_15 AC 2 ms
6,816 KB
testcase_16 AC 1 ms
6,820 KB
testcase_17 AC 2 ms
6,824 KB
testcase_18 AC 2 ms
6,820 KB
testcase_19 AC 2 ms
6,820 KB
testcase_20 AC 2 ms
6,820 KB
testcase_21 AC 2 ms
6,824 KB
testcase_22 AC 2 ms
6,820 KB
testcase_23 AC 1 ms
6,816 KB
testcase_24 AC 2 ms
6,820 KB
testcase_25 AC 1 ms
6,816 KB
testcase_26 AC 2 ms
6,820 KB
testcase_27 AC 2 ms
6,820 KB
testcase_28 AC 2 ms
6,816 KB
testcase_29 AC 2 ms
6,820 KB
testcase_30 AC 2 ms
6,820 KB
testcase_31 AC 2 ms
6,816 KB
testcase_32 AC 2 ms
6,820 KB
08_evil_01.txt AC 1,019 ms
6,820 KB
権限があれば一括ダウンロードができます

ソースコード

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