結果

問題 No.435 占い(Extra)
ユーザー r_dream0r_dream0
提出日時 2017-02-21 21:08:38
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 484 ms / 2,000 ms
コード長 1,456 bytes
コンパイル時間 981 ms
コンパイル使用メモリ 74,268 KB
実行使用メモリ 26,068 KB
最終ジャッジ日時 2024-10-08 12:21:19
合計ジャッジ時間 7,438 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 79 ms
25,984 KB
testcase_01 AC 79 ms
26,004 KB
testcase_02 AC 79 ms
25,984 KB
testcase_03 AC 79 ms
25,984 KB
testcase_04 AC 78 ms
25,968 KB
testcase_05 AC 79 ms
26,068 KB
testcase_06 AC 83 ms
25,984 KB
testcase_07 AC 79 ms
25,984 KB
testcase_08 AC 80 ms
25,984 KB
testcase_09 AC 79 ms
25,984 KB
testcase_10 AC 79 ms
25,984 KB
testcase_11 AC 81 ms
25,984 KB
testcase_12 AC 91 ms
26,016 KB
testcase_13 AC 91 ms
25,984 KB
testcase_14 AC 91 ms
25,984 KB
testcase_15 AC 93 ms
25,984 KB
testcase_16 AC 122 ms
25,984 KB
testcase_17 AC 375 ms
25,984 KB
testcase_18 AC 90 ms
25,984 KB
testcase_19 AC 92 ms
25,984 KB
testcase_20 AC 190 ms
25,984 KB
testcase_21 AC 190 ms
25,984 KB
testcase_22 AC 192 ms
26,000 KB
testcase_23 AC 200 ms
25,956 KB
testcase_24 AC 222 ms
25,984 KB
testcase_25 AC 484 ms
25,948 KB
testcase_26 AC 192 ms
25,984 KB
testcase_27 AC 192 ms
25,984 KB
testcase_28 AC 78 ms
25,984 KB
testcase_29 AC 189 ms
25,984 KB
testcase_30 AC 191 ms
25,984 KB
testcase_31 AC 181 ms
25,984 KB
testcase_32 AC 184 ms
25,984 KB
testcase_33 AC 206 ms
25,984 KB
testcase_34 AC 308 ms
25,984 KB
testcase_35 AC 198 ms
25,984 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
using namespace std;

const int MAX_N = 10000001;
// (3n)! = 3^thress[n]k (kは3の倍数ではない)
vector<int> threes;
vector<char> fact;
vector<int> inv;
void init() {
  threes.resize((MAX_N + 2) / 3);
  threes[0] = 0;
  for(size_t i = 1; i < threes.size(); i++) {
    size_t j = i * 3;
    threes[i] = threes[i - 1];
    while(j % 3 == 0) {
      threes[i]++;
      j /= 3;
    }
  }
  fact.resize(MAX_N);
  fact[0] = 1;
  for(size_t i = 1; i < fact.size(); i++) {
    int j = i;
    while(j % 3 == 0) j /= 3;
    fact[i] = (fact[i - 1] * j) % 9;
  }
  inv.resize(9);
  for(int i = 0; i < 9; i++) {
    for(int j = 0; j < 9; j++) {
      if(i * j % 9 == 1) inv[i] = j;
    }
  }
}

int comb(int n, int k) {
  int t = threes[n/3] - threes[(n - k) / 3] - threes[k / 3];
  if(t >= 2) return 0;
  int ret = fact[n] * inv[fact[n - k]] * inv[fact[k]] % 9;
  if(t == 1) ret = ret * 3 % 9;
  return ret;
}

int main() {
  init();
  int T;
  cin >> T;
  for(int i = 0; i < T; i++) {
    int n, x, a, b, m;
    cin >> n >> x >> a >> b >> m;
    int lastS = x;
    bool allzero = (lastS % 10) == 0;
    int ans = (lastS % 10) * comb(n - 1, 0) % 9;
    for(int i = 1; i < n; i++) {
      int nextS = ((lastS ^ a) + b) % m;
      ans = (ans + (nextS % 10) * comb(n - 1, i)) % 9;
      if(nextS % 10)allzero = false;
      lastS = nextS;
    }
    ans %= 9;
    if(ans == 0 && !allzero) ans = 9;
    cout << ans << endl;
  }
}
0