結果

問題 No.435 占い(Extra)
コンテスト
ユーザー tossy
提出日時 2016-10-24 13:44:19
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,688 bytes
コンパイル時間 905 ms
コンパイル使用メモリ 92,976 KB
実行使用メモリ 23,040 KB
最終ジャッジ日時 2024-10-08 12:11:54
合計ジャッジ時間 35,931 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 5 WA * 17 TLE * 10
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>
#include <cstring>
#include <string>
#include <cmath>
#include <cassert>
#include <iostream>
#include <algorithm>
#include <stack>
#include <queue>
#include <vector>
#include <set>
#include <map>
#include <bitset>
using namespace std;

#define repl(i,a,b) for(int i=(int)(a);i<(int)(b);i++)
#define rep(i,n) repl(i,0,n)
#define mp(a,b) make_pair(a,b)
#define pb(a) push_back(a)
#define all(x) (x).begin(),(x).end()
#define dbg(x) cout<<#x"="<<x<<endl
#define fi first
#define se second

#define INF 2147483600

vector<int> primes;

// x^n mod
long mod_pow(long x, long n, long mod){
  long res=1;
  while(n>0){
    if(n&1) res=res*x%mod;
    x=x*x%mod;
    n>>=1;
  }
  return res;
}

int main(){
  #define SZ 10000000
  vector<short> divs(SZ+1, 0); //divs[i] := iを割り切る素数のひとつ, 素数なら0
  for(int i=2; i<=SZ; i++) if(divs[i]==0){
    for(long j=(long)i*i; j<=SZ; j+=i) divs[j]=i;
  }

  int tc;
  scanf("%d\n", &tc);
  rep(loops, tc){
    int n,x,a,b,m;
    scanf("%d %d %d %d %d", &n ,&x, &a, &b, &m);

    bool isZero=true;
    vector<int> comb(9,0); //combinationの素因数分解
    int res = 0;

    rep(i,n){
      int r = x%10;
      if(r!=0) isZero=false;

      long cc=1;
      repl(j,2,8) cc *= mod_pow(j, comb[j]%9, 9); //2,3,5,7だけで十分
      res = (res + r * cc) % 9;

      int p = n-1-i;
      int q = i+1;
      while(divs[p]>0){
        comb[divs[p]%9]++;
        p /= divs[p];
      }
      comb[p%9]++;
      while(divs[q]>0){
        comb[divs[q]%9]--;
        q /= divs[q];
      }
      comb[q%9]--;

      x = ((x^a)%m + b)%m;
    }

    if(!isZero && res==0) res=9;
    printf("%d\n", res);
  }

  return 0;
}
0