結果

問題 No.97 最大の値を求めるくえり
コンテスト
ユーザー paruki
提出日時 2016-10-05 07:46:00
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,804 bytes
コンパイル時間 1,745 ms
コンパイル使用メモリ 167,204 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-11-21 17:12:20
合計ジャッジ時間 4,757 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 18 WA * 1
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include "bits/stdc++.h"
using namespace std;
#define FOR(i,j,k) for(int (i)=(j);(i)<(int)(k);++(i))
#define rep(i,j) FOR(i,0,j)
#define each(x,y) for(auto &(x):(y))
#define mp make_pair
#define all(x) (x).begin(),(x).end()
#define debug(x) cout<<#x<<": "<<(x)<<endl
#define smax(x,y) (x)=max((x),(y))
#define smin(x,y) (x)=min((x),(y))
#define MEM(x,y) memset((x),(y),sizeof (x))
#define sz(x) (int)(x).size()
typedef long long ll;
typedef pair<int, int> pii;
typedef vector<int> vi;
typedef vector<ll> vll;

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;
}

int inverse(int x, int mod){
    long long a = x, b = mod, u = 1, v = 0;
    while(b) {
        long long t = a / b;
        a -= t * b;
        swap(a, b);
        u -= t * v;
        swap(u, v);
    }
    if(u < 0)u += mod;
    return (int)u;
}

const int MOD = 100003;
int N, Q, q[MOD], ans[MOD], a[MOD];
bool b[MOD];

void f(){
    rep(i, Q){
        int ma = -1;
        rep(j, N)smax(ma, (int)((ll)q[i] * a[j] % MOD));
        ans[i] = ma;
    }
}

void g(){
    rep(i, N)b[a[i]] = 1;
    rep(i, Q){
        for(int j = MOD - 1; j >= 0; --j){
            int x = (ll)j*inverse(q[i], MOD)%MOD;
            if(b[x]){
                ans[i] = j;
                break;
            }
        }
    }
}

int main(){
    cin >> N >> Q;
    rep(i, Q)scanf("%d", &q[i]);
    generateA(N, a);
    if((ll)N*Q <= 500000000){
        f();
    } else{
        g();
    }
    rep(i, Q)printf("%d\n", ans[i]);
}
0