結果
問題 |
No.435 占い(Extra)
|
ユーザー |
![]() |
提出日時 | 2017-04-21 09:43:45 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 210 ms / 2,000 ms |
コード長 | 1,187 bytes |
コンパイル時間 | 1,421 ms |
コンパイル使用メモリ | 167,272 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-10-08 12:21:41 |
合計ジャッジ時間 | 5,485 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 32 |
ソースコード
#include<bits/stdc++.h> using namespace std; #define rep(i,a,b) for(int i=a;i<b;i++) const int inv[9] = { 0, 1, 5, 0, 7, 2, 0, 4, 8 }; int modpow(int a, int n, int mod = 9) { a %= mod; if (a == 0) return 0; int r = 1; while (n) r = r*((n % 2) ? a : 1) % mod, a = a*a%mod, n >>= 1; return r; } int sol(int N, int X, int A, int B, int M) { int a = 1, b = 1, t = 0; int v = X; int isAll0 = true; int res = 0; rep(i, 0, N) { int vv = v % 10; if (vv) isAll0 = false; int c = a * inv[b] * modpow(3, t) % 9; res = (res + vv * c) % 9; a *= (N - 1 - i); b *= (i + 1); while (a % 3 == 0 && 0 < a) a /= 3, t++; while (b % 3 == 0 && 0 < b) b /= 3, t--; a %= 9; b %= 9; v = ((v ^ A) + B) % M; } if (isAll0) return 0; if (res == 0) res = 9; return res; } //----------------------------------------------------------------------------------- int main() { cin.tie(0); ios::sync_with_stdio(false); int T; cin >> T; rep(t, 0, T) { int N, X, A, B, M; cin >> N >> X >> A >> B >> M; printf("%d\n", sol(N, X, A, B, M)); } }