結果
問題 |
No.916 Encounter On A Tree
|
ユーザー |
|
提出日時 | 2019-10-21 02:57:57 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 12 ms / 2,000 ms |
コード長 | 1,655 bytes |
コンパイル時間 | 1,573 ms |
コンパイル使用メモリ | 171,160 KB |
実行使用メモリ | 11,648 KB |
最終ジャッジ日時 | 2024-07-02 17:39:05 |
合計ジャッジ時間 | 3,325 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 5 |
other | AC * 56 |
ソースコード
#include <bits/stdc++.h> using namespace std; using ll = long long; const ll MOD = (ll)1e9 + 7; int main() { int d, l, r, k, MAX_d = 20; cin >> d >> l >> r >> k; assert(2<=d&&d<=20); assert(1<=l&&l<r&&r<pow(2,d)); assert(1<=k&&k<=100000); vector<ll> pow2(MAX_d + 1, 1), sum(MAX_d + 1, 0), fac((1<<MAX_d) + 1,1); //2冪を前計算 for (int i = 2; i < MAX_d + 1; ++i) pow2[i] = pow2[i - 1] * 2; //2冪の和を前計算 for (int i = 1; i < MAX_d + 1; ++i) sum[i] = sum[i - 1] + pow2[i]; //階乗を前計算 for (ll i = 1; i < (1 << MAX_d) + 1; ++i) fac[i] = (fac[i-1] * i) % MOD; //l,rを深さに変換 l = lower_bound(sum.begin(), sum.end(), l) - sum.begin(); r = lower_bound(sum.begin(), sum.end(), r) - sum.begin(); int lca = -1; //lcaが存在するならそのlcaの深さを計算 if ((l + r - k) > 1 && (l + r - k) % 2 == 0) lca = (l + r - k) / 2; //lcaが条件を満たしていない場合コーナー if(lca == -1 || lca > l || lca > r){ cout << 0 << endl; return 0; } ll ans = 1; //上の段からl,r以外の頂点の順列の数え上げ for (int i = 1; i <= d; ++i) { ll cnt = pow2[i]; if (i == l) cnt--; if (i == r) cnt--; (ans *= fac[cnt]) %= MOD; } //lcaとなるような頂点の位置についての数え上げ (ans *= pow2[lca]) %= MOD; //lcaを決め打ちした時のlの位置についての数え上げ (ans *= pow2[l - lca]) %= MOD; //lcaを決め打ちした時のrの位置についての数え上げ (ans *= pow2[r - lca]) %= MOD; //lcaについて、子は左右について2通りのパターンがあるため (ans *= 2LL) %= MOD; cout << ans << endl; return 0; }