結果

問題 No.1286 Stone Skipping
コンテスト
ユーザー mkmkmkmk
提出日時 2022-05-14 16:18:13
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 2,564 bytes
コンパイル時間 3,980 ms
コンパイル使用メモリ 230,852 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-07-23 02:06:17
合計ジャッジ時間 4,610 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <atcoder/all>
using namespace atcoder;
using mint = modint998244353;
//using mint = modint1000000007;
using namespace std;
typedef long long lint;
#define rep(i,n) for(lint i=0;i<n;i++)
#define repp(i,m,n) for(lint i=m;i<n;i++)
#define repm(i,n) for(lint i=n-1;i>=0;i--)
#define INF (1ll<<60)
#define all(x) (x).begin(),(x).end()
//const lint MOD =1000000007;
const lint MOD=998244353;
const lint MAX = 4000000;
using Graph =vector<vector<lint>>;
typedef pair<lint,lint> P;
typedef map<lint,lint> M;
#pragma GCC optimize("Ofast")
#define chmax(x,y) x=max(x,y)
#define chmin(x,y) x=min(x,y)

 
 
 
lint fac[MAX], finv[MAX], inv[MAX];
 
void COMinit() 
{
    fac[0] = fac[1] = 1;
    finv[0] = finv[1] = 1;
    inv[1] = 1;
    for (lint i = 2; i < MAX; i++)
    {
        fac[i] = fac[i - 1] * i % MOD;
        inv[i] = MOD - inv[MOD % i] * (MOD / i) % MOD;
        finv[i] = finv[i - 1] * inv[i] % MOD;
    }
}
 
long long COM(lint n, lint k)
{
    if (n < k)
        return 0;
    if (n < 0 || k < 0)
        return 0;
    return fac[n] * (finv[k] * finv[n - k] % MOD) % MOD;
}
 
lint primary(lint num)
{
    if (num < 2) return 0;
    else if (num == 2) return 1;
    else if (num % 2 == 0) return 0;
 
    double sqrtNum = sqrt(num);
    for (int i = 3; i <= sqrtNum; i += 2)
    {
        if (num % i == 0)
        {
            return 0;
        }
    }
 
    return 1;
}
   long long modpow(long long a, long long n, long long mod) {
    long long res = 1;
    while (n > 0) {
        if (n & 1) res = res * a % mod;
        a = a * a % mod;
        n >>= 1;
    }
    return res;
}
    lint lcm(lint a,lint b){
        return a/__gcd(a,b)*b;
    }
     lint gcd(lint a,lint b){
        return __gcd(a,b);
    }

    lint popcount(lint n){
        lint ret=0;
        while(n!=0){
            if(n%2)ret++;
            n/=2;
        }
        return ret;
    }

    lint calc(lint t,lint num){
        lint ret=t;
        rep(i,num){
            t/=2;
            ret+=t;
        }
        return ret;
    }

    int main(){
        lint d;
        cin>>d;
        lint ans=d;
        rep(h,65){
        lint l=0,r=d+1;
        while((r-l)>1){
            lint mid=(l+r)/2;
            lint cmp=calc(mid,h);
            if(cmp<d)l=mid;
            else r=mid;
        }
        if(calc(r,h)==d)chmin(ans,r);

        }
        cout<<ans<<endl;

    }
0