結果

問題 No.1968 Distance
ユーザー 👑 Nachia
提出日時 2022-06-03 22:41:31
言語 C++17(gcc12)
(gcc 12.3.0 + boost 1.87.0)
結果
AC  
実行時間 6 ms / 2,000 ms
コード長 1,404 bytes
コンパイル時間 2,434 ms
コンパイル使用メモリ 111,364 KB
最終ジャッジ日時 2025-01-29 17:52:30
ジャッジサーバーID
(参考情報)
judge1 / judge5
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 25
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <iostream>
#include <string>
#include <vector>
#include <algorithm>
#include <atcoder/modint>
using namespace std;
using i32 = int32_t;
using u32 = uint32_t;
using i64 = int64_t;
using u64 = uint64_t;
#define rep(i,n) for(int i=0; i<(int)(n); i++)


const i64 INF = 1001001001001001001;
using modint = atcoder::static_modint<998244353>;


modint triangleArea(i64 h){
    return (modint(h) * modint(h) * modint(2)).inv();
}


int main(){
    i64 N; cin >> N;
    if(N == 0){ cout << "748683265\n"; return 0; }
    modint ans = 0;
    for(i64 sh = 0; sh < 50; sh++){
        modint len = modint(((i64)1 << sh) + 1).inv() - modint(((i64)2 << sh) + 1).inv();
        modint area = triangleArea(((i64)1 << sh) + 1) - triangleArea(((i64)2 << sh) + 1);
        i64 givel = 1;
        i64 giver = ((i64)4 << sh) - 1;
        if(givel / 2 + giver / 2 >= N){ giver /= 2; givel /= 2; }
        while(givel + giver < N){ giver = giver * 2 + 1; givel = givel * 2 + 1; }
        if(givel + giver / 2 >= N){
            ans += area / modint(givel + 1);
        }
        else{
            ans += (len - area) / modint(giver + 1);
        }
    }
    ans += triangleArea(((i64)1 << 50) + 1);
    ans *= 2;
    cout << ans.val() << '\n';
    return 0;
}


struct ios_do_not_sync{
    ios_do_not_sync(){
        std::ios::sync_with_stdio(false);
        std::cin.tie(nullptr);
    }
} ios_do_not_sync_instance;


0