結果

問題 No.916 Encounter On A Tree
ユーザー ganmodokixganmodokix
提出日時 2019-10-25 22:55:59
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 10 ms / 2,000 ms
コード長 4,485 bytes
コンパイル時間 2,497 ms
コンパイル使用メモリ 185,932 KB
実行使用メモリ 11,516 KB
最終ジャッジ日時 2024-04-27 02:58:31
合計ジャッジ時間 3,450 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
6,812 KB
testcase_01 AC 1 ms
6,816 KB
testcase_02 AC 2 ms
6,940 KB
testcase_03 AC 10 ms
11,408 KB
testcase_04 AC 2 ms
6,940 KB
testcase_05 AC 9 ms
11,464 KB
testcase_06 AC 1 ms
6,940 KB
testcase_07 AC 10 ms
11,280 KB
testcase_08 AC 3 ms
6,944 KB
testcase_09 AC 4 ms
6,940 KB
testcase_10 AC 2 ms
6,940 KB
testcase_11 AC 2 ms
6,948 KB
testcase_12 AC 1 ms
6,940 KB
testcase_13 AC 9 ms
11,416 KB
testcase_14 AC 1 ms
6,940 KB
testcase_15 AC 9 ms
11,444 KB
testcase_16 AC 1 ms
6,944 KB
testcase_17 AC 9 ms
11,312 KB
testcase_18 AC 2 ms
6,940 KB
testcase_19 AC 9 ms
11,400 KB
testcase_20 AC 2 ms
6,944 KB
testcase_21 AC 9 ms
11,396 KB
testcase_22 AC 1 ms
6,944 KB
testcase_23 AC 9 ms
11,312 KB
testcase_24 AC 1 ms
6,940 KB
testcase_25 AC 8 ms
11,388 KB
testcase_26 AC 2 ms
6,944 KB
testcase_27 AC 6 ms
7,296 KB
testcase_28 AC 1 ms
6,940 KB
testcase_29 AC 9 ms
11,516 KB
testcase_30 AC 1 ms
6,940 KB
testcase_31 AC 2 ms
6,940 KB
testcase_32 AC 2 ms
6,944 KB
testcase_33 AC 1 ms
6,940 KB
testcase_34 AC 4 ms
6,940 KB
testcase_35 AC 3 ms
6,944 KB
testcase_36 AC 9 ms
11,344 KB
testcase_37 AC 9 ms
11,472 KB
testcase_38 AC 4 ms
6,940 KB
testcase_39 AC 2 ms
6,940 KB
testcase_40 AC 9 ms
11,384 KB
testcase_41 AC 2 ms
6,940 KB
testcase_42 AC 1 ms
6,944 KB
testcase_43 AC 2 ms
6,944 KB
testcase_44 AC 2 ms
6,944 KB
testcase_45 AC 3 ms
6,944 KB
testcase_46 AC 1 ms
6,940 KB
testcase_47 AC 1 ms
6,940 KB
testcase_48 AC 1 ms
6,944 KB
testcase_49 AC 2 ms
6,944 KB
testcase_50 AC 1 ms
6,940 KB
testcase_51 AC 1 ms
6,944 KB
testcase_52 AC 1 ms
6,944 KB
testcase_53 AC 1 ms
6,940 KB
testcase_54 AC 2 ms
6,940 KB
testcase_55 AC 2 ms
6,940 KB
testcase_56 AC 2 ms
6,940 KB
testcase_57 AC 2 ms
6,940 KB
testcase_58 AC 2 ms
6,940 KB
testcase_59 AC 2 ms
6,940 KB
testcase_60 AC 1 ms
6,944 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

// May this submission get accepted!

#pragma GCC optimize ("O3")
#pragma GCC target ("tune=native")
#pragma GCC target ("avx")

#include <bits/stdc++.h>

// 汎用マクロ
#define ALL_OF(x) (x).begin(), (x).end()
#define REP(i,n) for (long long i=0, i##_len=(n); i<i##_len; i++)
#define RANGE(i,is,ie) for (long long i=(is), i##_end=(ie); i<=i##_end; i++)
#define DSRNG(i,is,ie) for (long long i=(is), i##_end=(ie); i>=i##_end; i--)
#define STEP(i, is, ie, step) for (long long i=(is), i##_end=(ie), i##_step = (step); i<=i##_end; i+=i##_step)
#define UNIQUE(v) do { sort((v).begin(), (v).end()); (v).erase(unique((v).begin(), (v).end()), (v).end()); } while (false)
template<class T> bool chmax(T &a, const T &b) {if (a < b) {a = b; return true;} return false; }
template<class T> bool chmin(T &a, const T &b) {if (a > b) {a = b; return true;} return false; }
#define INF 0x7FFFFFFF
#define LINF 0x7FFFFFFFFFFFFFFFLL
#define Yes(q) ((q) ? "Yes" : "No")
#define YES(q) ((q) ? "YES" : "NO")
#define Possible(q) ((q) ? "Possible" : "Impossible")
#define POSSIBLE(q) ((q) ? "POSSIBLE" : "IMPOSSIBLE")
#define DUMP(q) cerr << "[DEBUG] " #q ": " << (q) << " at " __FILE__ ":" << __LINE__ << endl
#define DUMPALL(q) do { cerr << "[DEBUG] " #q ": ["; REP(i, (q).size()) { cerr << (q)[i] << (i == i_len-1 ? "" : ", "); } cerr << "] at " __FILE__ ":" << __LINE__ << endl; } while (false)
template<class T> T gcd(T a, T b) { if (a < b) std::swap(a, b); while (b) std::swap(a %= b, b); return a; }
template<class T> T lcm(const T a, const T b) { return a / gcd(a, b) * b; }

// gcc拡張マクロ
#define popcount __builtin_popcount
#define popcountll __builtin_popcountll

// エイリアス
#define DANCE_ long
#define ROBOT_ unsigned
#define HUMAN_ signed
#define CHOKUDAI_ const
using  ll = DANCE_ HUMAN_ DANCE_;
using ull = DANCE_ ROBOT_ DANCE_;
using cll = DANCE_ DANCE_ CHOKUDAI_;
using  ld = long double;
using namespace std;

// モジュール
// 剰余演算
constexpr ll pdiv = 1000000007LL; // 10^9+7
ll powll(ll a, ll n) {
    ll result = 1;
    ll base = a;
    ll mask = 1;
    while (mask <= n) {
        if (n & mask) (result *= base) %= pdiv;
        mask <<= 1;
        (base *= base) %= pdiv;
    }
    return result;
}

inline ll invll(ll a) { return powll(a, pdiv-2); }

vector<ll> fac_cache = {1}, invfac_cache = {1};
void make_fac_cache(ll a) {
    ll old_max = fac_cache.size() - 1;
    if (a > old_max) {
        fac_cache   .resize(a+1);
        invfac_cache.resize(a+1);
        for (ll i = old_max + 1; i <= a; i++) {
            fac_cache[i] = fac_cache[i-1] * i % pdiv;
        }
        invfac_cache[a] = invll(fac_cache[a]);
        for (ll i = a-1; i > old_max; i--) {
            invfac_cache[i] = invfac_cache[i+1] * (i + 1) % pdiv;
        }
    }
}

inline ll    facll(ll a) { make_fac_cache(a); return fac_cache[a]; }
inline ll invfacll(ll a) { make_fac_cache(a); return invfac_cache[a]; }
inline bool isoutll(ll n, ll r) { return n < 0 || r < 0 || n < r; }
inline ll nPr(ll n, ll r) { return isoutll(n, r) ? 0 : facll(n) * invfacll(n-r) % pdiv; }
inline ll nCr(ll n, ll r) { return isoutll(n, r) ? 0 : facll(n) * invfacll(n-r) % pdiv * invfacll(r) % pdiv; }
inline ll nHr(ll n, ll r) { return nCr(n+r-1, r); }

// 処理内容
int main() {
    
    // インタラクティブ問題では除去した方がいいかも
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    
    ll d, l, r, k; cin >> d >> l >> r >> k;

    ll ld = 0, rd = 0;
    while (l > 1) { l /= 2; ld++; }
    while (r > 1) { r /= 2; rd++; }

    if ((k + rd + ld) % 2) {
        cout << 0 << endl;
        return 0;
    }

    ll lcad = (ld + rd - k) / 2;
    if (lcad < 0 || rd <= lcad || ld < lcad) {
        cout << 0 << endl;
        return 0;
    }

    make_fac_cache(1LL << (d-1));

    ll ans = 1;
    REP(dc, d) {
        if (dc == rd) {
            ll plen = rd - lcad;
            ll n = 1LL << dc;
            if (ld != rd) {
                ans *= facll(n - 1);
                ans %= pdiv;
                ans *= 1LL << (lcad == ld ? plen : plen-1);
                ans %= pdiv;
            } else {
                ans *= n;
                ans %= pdiv;
                ans *= 1LL << (plen-1);
                ans %= pdiv;
                ans *= facll(n - 2);
                ans %= pdiv;
            }
        } else {
            (ans *= facll(1LL << dc)) %= pdiv;
        }
    }

    DUMP(lcad);
    DUMP(ld);
    DUMP(rd);

    cout << ans << endl;
    
}
0