結果

問題 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
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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));
    }
}
0