結果

問題 No.435 占い(Extra)
コンテスト
ユーザー はまやんはまやん
提出日時 2017-04-21 09:40:04
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
MLE  
実行時間 -
コード長 1,255 bytes
コンパイル時間 1,966 ms
コンパイル使用メモリ 169,900 KB
実行使用メモリ 68,832 KB
最終ジャッジ日時 2024-10-08 12:21:06
合計ジャッジ時間 7,746 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 29 MLE * 3
権限があれば一括ダウンロードができます

ソースコード

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) {
    vector<int> v;
    v.push_back(X);
    rep(i, 0, N - 1) v.push_back(((v[i] ^ A) + B) % M);
    rep(i, 0, N) v[i] %= 10;
    int sm = 0;
    for (int i : v) sm += i;
    if (sm == 0) return 0;

    int a = 1, b = 1, t = 0;
    int res = 0;
    int n = v.size();
    rep(i, 0, n) {
        int c = a * inv[b] * modpow(3, t) % 9;
        res = (res + v[i] * 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;
    }
    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