結果

問題 No.916 Encounter On A Tree
ユーザー east1016
提出日時 2019-10-25 23:07:42
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 3,115 bytes
コンパイル時間 1,612 ms
コンパイル使用メモリ 170,988 KB
実行使用メモリ 11,916 KB
最終ジャッジ日時 2024-09-13 07:10:55
合計ジャッジ時間 3,817 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 5
other AC * 55 WA * 1
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp:45:12: warning: use of 'auto' in parameter declaration only available with '-std=c++20' or '-fconcepts'
   45 | void press(auto &v){v.erase(unique(ALL(v)),v.end());}  //  圧縮
      |            ^~~~

ソースコード

diff # 

#include <bits/stdc++.h>  // ver2.2.4
#define int long long
#define endl "\n"
#define ALL(v) (v).begin(),(v).end()
#define Vi vector<int>
#define VVi vector<Vi>
#define Vs vector<string>
#define Pii pair<int,int>
#define VPii vector<Pii>
#define Tiii tuple<int,int,int>
#define PQi priority_queue<int>
#define PQir priority_queue<int,vector<int>,greater<int>>
#define pb push_back
#define mp make_pair
#define mt make_tuple
#define itos to_string
#define stoi stoll
#define FI first
#define SE second
#define cYES cout<<"YES"<<endl
#define cNO cout<<"NO"<<endl
#define cYes cout<<"Yes"<<endl
#define cNo cout<<"No"<<endl
#define _ <<' '<<
#define sortr(v) sort(v,greater<>())
#define rep(i,a,b) for(int i=a;i<b;i++)
#define repeq(i,a,b) for(int i=a;i<=b;i++)
#define repreq(i,a,b) for(int i=a;i>=b;i--)
#define leng(n) (int)(log10(n)+1)
#define dem(a,b) ((a+b-1)/(b))
#define Vin(a) rep(iI,0,a.size())cin>>a[iI]
#define Vout(a) rep(lZ,0,a.size()-1)cout<<a[lZ]<<' ';cout<<a.back()<<endl
#define VVout(a) rep(lY,0,a.size()){Vout(a[lY]);}
#define VPout(a) rep(lX,0,a.size())cout<<a[lX].FI<<' '<<a[lX].SE<<endl
#define Verr(a) rep(iZ,0,a.size()-1)cerr<<a[iZ]<<' ';cerr<<a.back()<<endl
#define VVerr(a) rep(iY,0,a.size()){Verr(a[iY]);}
#define VPerr(a) rep(iX,0,a.size())cerr<<a[iX].FI<<' '<<a[iX].SE<<endl
#define INF 3000000000000000000  //  3.0*10^18(MAXの1/3くらい)
#define MAX LLONG_MAX
#define MOD 1000000007  //  10^9 + 7
using namespace std;

int gcd(int a,int b){return b?gcd(b,a%b):a;}  //  最大公約数gcd
int lcm(int a,int b){return a/gcd(a,b)*b;}  //  最小公倍数lcm
void press(auto &v){v.erase(unique(ALL(v)),v.end());}  //  圧縮
int mypow(int x, int n, int m){  //  累乗x^n(mod m)  O(log(n))
if(n==0)return 1;if(n%2==0)return mypow(x*x%m,n/2,m);else return x*mypow(x,n-1,m)%m;}

int stpow(int n){
  if(n==0)return 1;
  int a=1;
  rep(i,1,n+1)a=a*i%MOD;
  return a;
}
Vi stpowV(){
  Vi a(1100001);
  a[0]=1;
  repeq(i,1,1100000)a[i]=a[i-1]*i%MOD;
  return a;
}

signed main() {
  cin.tie(0);cout.tie(0);ios::sync_with_stdio(false);cout<<fixed<<setprecision(12);
  int n,l,r,k;
  cin >> n >> l >> r >> k;
  Vi stp = stpowV();
  
  int qw,er,ty;
  int depl=0,depr=0;
  qw = 1;
  while(1){
    if(qw*2 <= l) {
      qw *= 2;
      depl++;
    }
    else break;
  }
  depl++;
  cerr << depl <<endl;
  qw = 1;
  while(1){
    if(qw*2 <= r) {
      qw *= 2;
      depr++;
    }
    else break;
  }
  depr++;
  cerr << depr << endl;
  
  //ここまでok
  
  int minn = abs(depl-depr), maxx = depl + depr;
  
  if(k < minn || k > maxx) {cout<<0<<endl;return 0;}
  if(k%2 != maxx%2) {cout<<0<<endl;return 0;}
  
  int upper = (k-minn)/2, stand = mypow(2,minn,MOD);
  if(upper==0) upper++;
  int deci = stand * mypow(2,upper-1,MOD) % MOD;
  int ue = min(depr,depl);
  deci %= MOD;
  int ans = mypow(2,ue-1,MOD) * deci % MOD;
  ans %= MOD;
  cerr<<"nowans="<<ans<<endl;
  
  repeq(i,1,n){
    qw = mypow(2,i-1,MOD);
    
    if(i==depl) qw--;
    if(i==depr) qw--;
    
    ans = ans * stp[qw] % MOD;
  }
  ans %= MOD;
  cout << ans << endl;
  
  cerr<<'#'<<endl;
  cerr<<deci<<endl;
  
  return 0;
}
0