結果

問題 No.613 Solitude by the window
ユーザー rickythetarickytheta
提出日時 2017-12-13 11:33:09
言語 C++11
(gcc 11.4.0)
結果
TLE  
実行時間 -
コード長 1,632 bytes
コンパイル時間 1,324 ms
コンパイル使用メモリ 158,096 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-05-08 07:13:40
合計ジャッジ時間 13,694 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
5,248 KB
testcase_01 AC 1 ms
5,376 KB
testcase_02 TLE -
testcase_03 AC 5 ms
5,376 KB
testcase_04 AC 245 ms
5,376 KB
testcase_05 AC 4 ms
5,376 KB
testcase_06 AC 5 ms
5,376 KB
testcase_07 AC 67 ms
5,376 KB
testcase_08 AC 46 ms
5,376 KB
testcase_09 AC 516 ms
5,376 KB
testcase_10 AC 29 ms
5,376 KB
testcase_11 AC 1,710 ms
5,376 KB
testcase_12 TLE -
testcase_13 TLE -
testcase_14 AC 3 ms
5,376 KB
testcase_15 AC 4 ms
5,376 KB
testcase_16 AC 3 ms
5,376 KB
testcase_17 AC 1 ms
5,376 KB
testcase_18 AC 1 ms
5,376 KB
testcase_19 AC 2 ms
5,376 KB
testcase_20 AC 2 ms
5,376 KB
testcase_21 AC 476 ms
5,376 KB
testcase_22 AC 483 ms
5,376 KB
testcase_23 AC 357 ms
5,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef vector<int> vi;
typedef vector<ll> vl;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;

typedef int _loop_int;
#define REP(i,n) for(_loop_int i=0;i<(_loop_int)(n);++i)
#define FOR(i,a,b) for(_loop_int i=(_loop_int)(a);i<(_loop_int)(b);++i)
#define FORR(i,a,b) for(_loop_int i=(_loop_int)(b)-1;i>=(_loop_int)(a);--i)

#define DEBUG(x) cout<<#x<<": "<<x<<endl
#define DEBUG_VEC(v) cout<<#v<<":";REP(i,v.size())cout<<" "<<v[i];cout<<endl
#define ALL(a) (a).begin(),(a).end()

#define CHMIN(a,b) a=min((a),(b))
#define CHMAX(a,b) a=max((a),(b))

#define rsh 31
#define r (1ll<<31)
#define rmsk (r-1ll)

ll n,m;
ll r2, mm;
ll m4;

ll calc_mm(){
  ll res = 0;
  ll t = 0;
  ll rr = r;
  ll i = 1;
  while(rr > 1){
    if(!(t%2)){
      t += m;
      res += i;
    }
    t /= 2;
    rr /= 2;
    i *= 2;
  }
  return res;
}

inline ll MR(ll x){
  ll t = (x + ((x * mm) & rmsk)*m) >> rsh;
  return t<m ? t : t-m;
}

#define f(x) ((x)*((x)+4ll)%m)

int floyd_cycle_algorithm(ll init){
  int mu = 0;
  init = MR(r2 * init);
  ll x = init;
  do {
    mu++;
    x = MR(x*(x+m4));
  }while(x!=init);
  return mu;
}

ll naive(ll x, ll n){
  x = MR(r2 * x);
  while(n--){
    x = MR(x*(x+m4));
  }
  return MR(x);
}

int main(){
  cin>>n>>m;
  if(m == 2){
    cout<<0<<endl;
    return 0;
  }
  // montgomery
  r2 = r*r%m;
  mm = calc_mm();
  m4 = MR(r2 * 4);

  // solve
  if(n < 100){
    cout<<naive(2,n)<<endl;
    return 0;
  }
  ll x = naive(2,100);
  n -= 100;
  int mu = floyd_cycle_algorithm(x);
  n %= mu;
  x = naive(x,n);
  cout<<x<<endl;
  return 0;
}
0