結果

問題 No.2994 べき内積
ユーザー 👑 ygussany
提出日時 2024-12-19 00:50:40
言語 C
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 2,494 bytes
コンパイル時間 523 ms
コンパイル使用メモリ 37,760 KB
実行使用メモリ 19,232 KB
最終ジャッジ日時 2024-12-19 00:50:52
合計ジャッジ時間 11,295 ms
ジャッジサーバーID
(参考情報)
judge3 / judge4
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 16 WA * 7
権限があれば一括ダウンロードができます

ソースコード

diff #

#pragma GCC target("avx2")
#pragma GCC optimize("Ofast,unroll-loops")
#include <stdio.h>

const int p = 1009;

long long pow_mod(int n, long long k)
{
	long long N, ans = 1;
	for (N = n; k > 0; k >>= 1, N = N * N % p) if (k & 1) ans = ans * N % p;
	return ans;
}

#define MT_N 624
#define MT_M 397
#define MT_MATRIX_A 0x9908b0dfUL
#define MT_UPPER_MASK 0x80000000UL
#define MT_LOWER_MASK 0x7fffffffUL

static unsigned int mt[MT_N];
static int mti = MT_N + 1;

void init_genrand(unsigned int s)
{
    mt[0] = s & 0xffffffffUL;
    for (mti = 1; mti < MT_N; mti++) {
        mt[mti] = (1812433253UL * (mt[mti-1] ^ (mt[mti-1] >> 30)) + mti); 
        mt[mti] &= 0xffffffffUL;
    }
}

unsigned int genrand()
{
    unsigned int y;
    static unsigned int mag01[2] = {0x0UL, MT_MATRIX_A};

    if (mti >= MT_N) {
        int kk;
        if (mti == MT_N + 1) init_genrand(5489UL);
		
        for (kk = 0; kk < MT_N - MT_M; kk++) {
            y = (mt[kk] & MT_UPPER_MASK) | (mt[kk+1] & MT_LOWER_MASK);
            mt[kk] = mt[kk+MT_M] ^ (y >> 1) ^ mag01[y&0x1UL];
        }
        for (; kk < MT_N - 1; kk++) {
            y = (mt[kk] & MT_UPPER_MASK) | (mt[kk+1] & MT_LOWER_MASK);
            mt[kk] = mt[kk+(MT_M-MT_N)] ^ (y >> 1) ^ mag01[y&0x1UL];
        }
        y = (mt[MT_N-1] & MT_UPPER_MASK) | (mt[0] & MT_LOWER_MASK);
        mt[MT_N-1] = mt[MT_M-1] ^ (y >> 1) ^ mag01[y&0x1UL];

        mti = 0;
    }
  
    y = mt[mti++];

    y ^= (y >> 11);
    y ^= (y << 7) & 0x9d2c5680UL;
    y ^= (y << 15) & 0xefc60000UL;
    y ^= (y >> 18);

    return y;
}

int main()
{
	int i, M, N, K[400001], A[2001];
	scanf("%d %d", &M, &N);
	for (i = 0; i <= M; i++) scanf("%d", &(K[i]));
	for (i = 0; i <= N; i++) scanf("%d", &(A[i]));
	// for (i = 0, K[0] = p - 1; i <= N; i++) A[i] = genrand() % p;
	
	int j, k;
	long long B[1009][2001];
	for (i = 1, B[0][0] = 1; i <= N; i++) B[0][i] = 0;
	for (i = 0; i <= N; i++) B[1][i] = A[i];
	for (k = 2; k < p; k++) {
		for (i = 0; i <= N; i++) for (j = 0, B[k][i] = 0; j <= i; j++) B[k][i] += B[k-1][j] * A[i-j];
		for (i = 0; i <= N; i++) B[k][i] %= p;
	}
	
	if (A[0] == 0) {
		for (i = M, k = 0; i >= 0 && k < p; i--) k = k * p + K[i];
		if (k < p) for (i = 0; i <= N; i++) printf("%lld ", B[k][i]);
		else for (i = 0; i <= N; i++) printf("0 ");
	} else {
		for (i = M, k = 0; i >= 1; i--) k = (k * p + K[i]) % (p - 1);
		k = pow_mod(A[0], k);
		for (i = 0; i <= N; i++) printf("%lld ", B[K[0]%p][i] * k % p);
	}
	printf("\n");
	fflush(stdout);
	return 0;
}
0