結果

問題 No.1853 Many Operations
コンテスト
ユーザー 蜜蜂
提出日時 2021-12-28 18:10:16
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 3 ms / 2,000 ms
コード長 1,884 bytes
コンパイル時間 3,947 ms
コンパイル使用メモリ 237,572 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-11-17 10:47:22
合計ジャッジ時間 4,924 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 26
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function 'int main()':
main.cpp:85:10: warning: structured bindings only available with '-std=c++17' or '-std=gnu++17' [-Wc++17-extensions]
   85 |     auto [l,r]=ave.back();
      |          ^
main.cpp:99:12: warning: structured bindings only available with '-std=c++17' or '-std=gnu++17' [-Wc++17-extensions]
   99 |   for(auto [l,r]:ave){
      |            ^

ソースコード

diff # 

//g++ 1.cpp -std=c++14 -O2 -I .
#include <bits/stdc++.h>
using namespace std;

#include <atcoder/all>
using namespace atcoder;

using ll = long long;
using ld = long double;

using vi = vector<int>;
using vvi = vector<vi>;
using vll = vector<ll>;
using vvll = vector<vll>;
using vld = vector<ld>;
using vvld = vector<vld>;
using vst = vector<string>;
using vvst = vector<vst>;

#define fi first
#define se second
#define pb push_back
#define pq_big(T) priority_queue<T,vector<T>,less<T>>
#define pq_small(T) priority_queue<T,vector<T>,greater<T>>
#define all(a) a.begin(),a.end()
#define rep(i,start,end) for(ll i=start;i<(ll)(end);i++)
#define per(i,start,end) for(ll i=start;i>=(ll)(end);i--)
#define uniq(a) sort(all(a));a.erase(unique(all(a)),a.end())

constexpr ll mod = 998244353;

map<ll,ll> value;

// sum g[0] + ... + g[x]
ll h(ll x){
  if(x<=1){
    return x;
  }
  if(value[x]!=0){
    return value[x];
  }
  ll res=0;
  ll n;

  if(x>=2){
    //even g[2]+g[4]+...+g[2n] = g[1]+g[2]+...+g[n]+n
    n=x/2;
    res+=h(n)+n;
  }

  if(x>=1){
    //1 mod 4 g[1]+g[5]+...+g[4n+1] = g[0]+g[4]+...+g[4n]+n+1 = g[4]+...+g[4n]+n+1 = g[1]+...+g[n]+3n+1
    n=(x-1)/4;
    res+=h(n)+3*n+1;
  }

  if(x>=3){
    //3 mod 4 g[3]+g[7]+...+g[4n+3] = g[4]+g[8]+...+g[4n+4]+n+1 = g[1]+...+g[n+1]+3n+3
    n=(x-3)/4;
    res+=h(n+1)+3*n+3;
  }

  res%=mod;
  value[x]=res;

  return res;
}

// [a,b] [c,d]
ll range(ll a,ll b,ll c,ll d){
  ll l=max(a,c);
  ll r=min(b,d);

  return max(r-l+1,(ll)0);
}

int main(){
  ios::sync_with_stdio(false);
  cin.tie(nullptr);

  vector<pair<ll,ll>> ave;
  ave.emplace_back(3,3);

  while(true){
    auto [l,r]=ave.back();
    l=l*2-1+l%2;
    r=r*2+1-r%2;
    ave.emplace_back(l,r);
    if(l>=2e18){
      break;
    }
  }

  ll n;
  cin>>n;

  ll ans=h(n);

  for(auto [l,r]:ave){
    ans+=mod-range(1,n,l,r)%mod;
  }

  ans%=mod;

  cout<<ans<<endl;
}
0