結果

問題 No.916 Encounter On A Tree
ユーザー otamay6otamay6
提出日時 2019-10-25 22:46:43
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 1,864 bytes
コンパイル時間 1,618 ms
コンパイル使用メモリ 172,636 KB
実行使用メモリ 21,732 KB
最終ジャッジ日時 2024-09-13 06:20:24
合計ジャッジ時間 3,687 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 AC 2 ms
6,940 KB
testcase_02 AC 2 ms
6,940 KB
testcase_03 AC 27 ms
20,324 KB
testcase_04 AC 2 ms
6,940 KB
testcase_05 AC 28 ms
21,068 KB
testcase_06 WA -
testcase_07 AC 28 ms
20,880 KB
testcase_08 AC 6 ms
6,944 KB
testcase_09 AC 10 ms
7,632 KB
testcase_10 AC 2 ms
6,944 KB
testcase_11 AC 2 ms
6,944 KB
testcase_12 AC 2 ms
6,940 KB
testcase_13 AC 27 ms
20,436 KB
testcase_14 AC 2 ms
6,940 KB
testcase_15 AC 28 ms
20,816 KB
testcase_16 AC 2 ms
6,944 KB
testcase_17 AC 27 ms
19,932 KB
testcase_18 AC 2 ms
6,940 KB
testcase_19 AC 28 ms
20,112 KB
testcase_20 WA -
testcase_21 AC 28 ms
20,396 KB
testcase_22 AC 2 ms
6,940 KB
testcase_23 AC 28 ms
20,056 KB
testcase_24 AC 2 ms
6,940 KB
testcase_25 AC 28 ms
19,824 KB
testcase_26 AC 2 ms
6,940 KB
testcase_27 AC 15 ms
11,992 KB
testcase_28 AC 2 ms
6,940 KB
testcase_29 AC 27 ms
20,300 KB
testcase_30 AC 2 ms
6,944 KB
testcase_31 AC 2 ms
6,940 KB
testcase_32 AC 2 ms
6,944 KB
testcase_33 AC 2 ms
6,940 KB
testcase_34 AC 10 ms
7,764 KB
testcase_35 AC 7 ms
6,944 KB
testcase_36 AC 28 ms
21,432 KB
testcase_37 AC 27 ms
20,244 KB
testcase_38 AC 10 ms
8,016 KB
testcase_39 AC 1 ms
6,944 KB
testcase_40 AC 30 ms
21,732 KB
testcase_41 AC 2 ms
6,940 KB
testcase_42 AC 2 ms
6,944 KB
testcase_43 AC 2 ms
6,940 KB
testcase_44 AC 2 ms
6,940 KB
testcase_45 AC 6 ms
6,940 KB
testcase_46 AC 2 ms
6,944 KB
testcase_47 AC 2 ms
6,940 KB
testcase_48 AC 2 ms
6,940 KB
testcase_49 AC 2 ms
6,944 KB
testcase_50 AC 2 ms
6,944 KB
testcase_51 AC 2 ms
6,940 KB
testcase_52 AC 2 ms
6,944 KB
testcase_53 AC 2 ms
6,944 KB
testcase_54 AC 2 ms
6,944 KB
testcase_55 AC 2 ms
6,944 KB
testcase_56 AC 2 ms
6,944 KB
testcase_57 AC 2 ms
6,944 KB
testcase_58 AC 2 ms
6,940 KB
testcase_59 AC 3 ms
6,940 KB
testcase_60 AC 3 ms
6,940 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function 'int main()':
main.cpp:56:13: warning: 'Rd' may be used uninitialized [-Wmaybe-uninitialized]
   56 |     ll p=(Ld+Rd-k);
      |           ~~^~~
main.cpp:49:11: note: 'Rd' was declared here
   49 |     ll Ld,Rd;
      |           ^~
main.cpp:67:9: warning: 'Ld' may be used uninitialized [-Wmaybe-uninitialized]
   67 |         if(i==Ld){
      |         ^~
main.cpp:49:8: note: 'Ld' was declared here
   49 |     ll Ld,Rd;
      |        ^~

ソースコード

diff # 

#include<bits/stdc++.h>
#define REP(i,n) for(int i=0,i##_len=(n);i<i##_len;++i)
#define rep(i,a,b) for(int i=int(a);i<int(b);++i)
#define All(x) (x).begin(),(x).end()
using namespace std;
typedef long long ll;
constexpr ll mod=1e9+7;
ll mpow(ll x, ll n){
    ll ans = 1; x%=mod;
    while(n != 0){
        if(n&1) ans = ans*x % mod;
        x = x*x % mod;
        n = n >> 1;
    }
    return ans;
}
ll inv_mod(ll a){return mpow(a,mod-2);}
class Factorial{//階乗とその逆元を求めて計算に利用するクラス
 private:
    vector<ll> fac;
    vector<ll> ifac;
 public:
    
    Factorial(ll N){
        fac.push_back(1);
        REP(i,N) fac.push_back(fac[i]*(i+1)%mod);
        ifac.resize(N+1);
        ifac[N]=inv_mod(fac[N]);
        REP(i,N) ifac[N-1-i]=(ifac[N-i]*(N-i))%mod;
    }

    ll fact(ll a){return fac[a];}
    ll ifact(ll a){return ifac[a];}

    ll cmb(ll a,ll b){
        if(a==0&&b==0) return 1;
        if(a<b||a<0||b<0) return 0;
        ll tmp =ifact(a-b)*ifact(b)%mod;
        return tmp*fac[a]%mod;
    }
    ll per(ll a,ll b){
        if(a==0&&b==0) return 1;
        if(a<b||a<0||b<0) return 0;
        return fac[a]*ifac[a-b]%mod;
    }
};
int main(){
    ll d,l,r,k;cin>>d>>l>>r>>k;
    ll Ld,Rd;
    ll tmp=1;
    REP(i,d){
        if(tmp<=l&&l<2*tmp) Ld=i;
        if(tmp<=r&&r<2*tmp) Rd=i;
        tmp*=2;
    }
    ll p=(Ld+Rd-k);
    if(p&1||p<0||d-1<=p/2){
        cout<<0<<endl;
        return 0;
    }
    p/=2;
    Factorial get(tmp);
    ll ans=1;
    ll v=1;
    REP(i,d){
        ll cnt=0;
        if(i==Ld){
            ans*=v;
            ans%=mod;
            cnt++;
        }
        if(i==Rd){
            ans=ans*mpow(2,Rd-p-1)%mod;
            if(p==Ld) ans*=2;
            ans%=mod;
            cnt++;
        }
        ans*=get.fact(v-cnt);
        ans%=mod;
        v*=2;
    }
    cout<<ans<<endl;
}
0