結果

問題 No.97 最大の値を求めるくえり
コンテスト
ユーザー koba-e964
提出日時 2015-07-28 17:15:32
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
TLE  
実行時間 -
コード長 1,649 bytes
コンパイル時間 573 ms
コンパイル使用メモリ 57,824 KB
実行使用メモリ 11,808 KB
最終ジャッジ日時 2024-07-17 13:20:49
合計ジャッジ時間 14,394 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 17 TLE * 1 -- * 1
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>

#define REP(i,s,n) for(int i=(int)(s);i<(int)(n);i++)

using namespace std;
typedef long long int ll;

unsigned xor128_x = 123456789, xor128_y = 362436069, xor128_z = 521288629, xor128_w = 88675123;
unsigned xor128() {
	unsigned t = xor128_x ^ (xor128_x << 11);
	xor128_x = xor128_y; xor128_y = xor128_z; xor128_z = xor128_w;
	return xor128_w = xor128_w ^ (xor128_w >> 19) ^ (t ^ (t >> 8));
}
void generateA(int N, int A[]) {
    for(int i = 0; i < N; ++ i)
        A[i] = xor128() % 100003;
}

const int M =100005;
int a[M];
int b[M];
int aset[M] = {};
int memo[M];
int imod[M];
const ll mod = 100003;

ll invmod(ll x) {
  ll e = mod - 2;
  ll sum = 1;
  ll cur = x;
  while (e > 0) {
    if (e % 2) {
      sum = sum * cur % mod;
    }
    cur = cur * cur % mod;
    e /= 2;
  }
  return sum;
}

ll solve(const int *as, int q) {
  ll s = mod - 1;
  while (1) {
    if (aset[(s * imod[q]) % mod]) {
      return s;
    }
    s--;
  }
}

int main(void){
  int n, q;
  cin >> n >> q;
  REP(i, 0, q) {
    cin >> b[i];
  }
  generateA(n, a);
  REP(i, 0, mod) {
    memo[i] = -1;
  }
  REP(i, 1, mod) {
    imod[i] = invmod(i);
  }
  if (n >= 20000) {
    REP(i, 0, n) {
      aset[a[i]] = 1;
    }
    REP(i, 0, q) {
      if (memo[b[i]] >= 0) {
	cout << memo[b[i]] << endl;
	continue;
      }
      memo[b[i]] = solve(aset, b[i]);
      cout << memo[b[i]] << endl;
    }
    return 0;
  }
  REP(i, 0, q) {
    if (memo[b[i]] >= 0) {
      cout << memo[b[i]] << endl;
      continue;
    }
    ll ma = 0;
    REP(j, 0, n) {
      ma = max(ma, ((ll)a[j] * (ll)b[i]) % mod);
    }
    cout << ma << endl;
    memo[b[i]] = ma;
  }
}
0