結果

問題 No.435 占い(Extra)
ユーザー char134217728
提出日時 2017-12-30 19:57:44
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 325 ms / 2,000 ms
コード長 2,370 bytes
コンパイル時間 2,141 ms
コンパイル使用メモリ 194,592 KB
最終ジャッジ日時 2025-01-05 06:42:19
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 32
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp:59:1: warning: ISO C++ forbids declaration of ‘main’ with no type [-Wreturn-type]
   59 | main(){
      | ^~~~

ソースコード

diff # 

#include <bits/stdc++.h>
#define FOR(i,a,b) for (int i=(a);i<(b);i++)
#define FORR(i,a,b) for (int i=(a);i>=(b);i--)
#define pb push_back
#define pcnt __builtin_popcount
#define show(x) cout<<#x<<" = "<<x<<endl;
#define maxs(x,y) x = max(x,y)
#define mins(x,y) x = min(x,y)
#define fi first
#define se second
#define rng(a) a.begin(),a.end()
#define each(it,c) for(__typeof((c).begin()) it=(c).begin();it!=(c).end();it++)
#define sz(x) (int)(x).size()
#define mp make_pair

using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef vector<pii> vpii;
typedef set<int> si;
typedef pair<ll,ll> pll;
typedef vector<ll> vl;
typedef vector<vl> vvl;
typedef vector<pll> vpll;
typedef set<ll> sl;
typedef __int128_t lll;
typedef pair<lll,lll> plll;
typedef vector<lll> vlll;
template<typename T>string join(vector<T>&v)
{stringstream s;FOR(i,0,sz(v))s<<' '<<v[i];return s.str().substr(1);}
ll gcd(ll a,ll b){if(a>b)swap(a,b);for(;a>0;b%=a,swap(a,b));return b;}
int modpow(ll a,ll n,int m){if(a==0)return a;ll p=1;for(;n>0;n/=2,a=a*a%m)if(n&1)p=p*a%m;return(int)p;}
void dout(double d){printf("%.12f\n",d);}

const int iinf = 1e9;
const ll linf = 1e18;
const int mod = 1e9+7;
const double pi = acos(-1);
const double eps = 1e-10;
int fac[3400000][2], inv[10] = {0, 1, 5, 0, 7, 2, 0, 4, 8, 0};
ll comb(int n, int m){
  int nm, _n, _m, _nm;
  nm = n-m;
  _n = n/3;
  _m = m/3;
  _nm = nm/3;
  ll k = fac[_n][1] - fac[_m][1] - fac[_nm][1];
  if(k > 1) return 0;
  ll ret = fac[_n][0] * inv[fac[_m][0]] * inv[fac[_nm][0]];
  FOR(i, _n*3+1, n+1)ret *= i%9;
  FOR(i, _m*3+1, m+1)ret *= inv[i%9];
  FOR(i, _nm*3+1, nm+1)ret *= inv[i%9];
  if(k > 0) ret *= 3;
  return ret%9;
}

main(){
  cin.tie(0);
  ios::sync_with_stdio(false);
  fac[0][0] = 1;
  FOR(i, 1, 3400000){
    int j = i*3, c = 0;
    while(j%3==0){
      j /= 3;
      c++;
    }
    fac[i][0] = fac[i-1][0] * (i*3-2) % 9 * (i*3-1) % 9 * j % 9;
    fac[i][1] = fac[i-1][1] + c;
  }
  int t;
  ll n, x, a, b, m;
  ll d, ans = 0;
  bool az = true;
  cin >> t;
  FOR(i, 0, t){
    cin >> n >> x >> a >> b >> m;
    ans = 0;
    az = true;
    FOR(j, 0, n){
      d = x % 10;
      if(d > 0) az = false;
      ans += comb(n-1, j) * d;
      x = ((x^a) + b ) % m;
    }
    ans %= 9;
    if(ans == 0 && !az) ans = 9;
    cout << ans << "\n";
  }
  return 0;
}
0