結果

問題 No.14 最小公倍数ソート
コンテスト
ユーザー bal4u
提出日時 2021-08-14 17:02:57
言語 C
(gcc 13.3.0)
結果
AC  
実行時間 11 ms / 5,000 ms
コード長 2,005 bytes
コンパイル時間 972 ms
コンパイル使用メモリ 30,976 KB
実行使用メモリ 9,472 KB
最終ジャッジ日時 2024-10-05 13:21:06
合計ジャッジ時間 2,137 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 20
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.c: In function 'in':
main.c:7:14: warning: implicit declaration of function 'getchar_unlocked' [-Wimplicit-function-declaration]
    7 | #define gc() getchar_unlocked()
      |              ^~~~~~~~~~~~~~~~
main.c:11:24: note: in expansion of macro 'gc'
   11 |         int n = 0, c = gc();
      |                        ^~
main.c: In function 'out':
main.c:8:15: warning: implicit declaration of function 'putchar_unlocked' [-Wimplicit-function-declaration]
    8 | #define pc(c) putchar_unlocked(c)
      |               ^~~~~~~~~~~~~~~~
main.c:29:35: note: in expansion of macro 'pc'
   29 |         prev = n, f[n]--; if (sp) pc(' ');
      |                                   ^~

ソースコード

diff # 

// yukicoder: 14 最小公倍数ソート
// 2021.8.14

#include <stdio.h>
#include <stdlib.h>

#define gc() getchar_unlocked()
#define pc(c) putchar_unlocked(c)

int in() {   // 非負整数の入力
	int n = 0, c = gc();
	do n = 10 * n + (c & 0xf), c = gc(); while (c >= '0');
	return n;
}

#define MAX 10000
#define MA 128
int a[MAX+5], f[MAX+5]; int sz;
int *qi[MAX+5], hi[MAX+5], ti[MAX+5], mi[MAX+5];
int *qj[MAX+5], hj[MAX+5], mj[MAX+5];
int prev;

int cmp(const void *a, const void *b) {
	return *(int *)a - *(int *)b;
}

void out(int n, int sp) { // 非負整数の表示（出力）
	char b[20];
	prev = n, f[n]--; if (sp) pc(' '); 
	int i = 0; while (n) b[i++] = n % 10 + '0', n /= 10;
	while (i--) pc(b[i]);
}

void mktbl(int ma) {
	int i, j;
	for (i = 2; i <= ma; ++i) for (j = i; j <= ma; j += i) if (f[j]) {
		if (hi[i] >= mi[i]) {
			mi[i] += MA, qi[i] = realloc(qi[i], sizeof(int)*mi[i]);
		}
		qi[i][hi[i]++] = j;
		if (hj[j] >= mj[j]) {
			mj[j] += MA, qj[j] = realloc(qj[j], sizeof(int)*mj[j]);
		}
		qj[j][hj[j]++] = i;
	}
}

int next(int x) {
	int i, u, v, g, t;
	g = prev * x;
	for (i = hj[prev]-1; i >= 0; --i) {
		u = qj[prev][i];
		while (ti[u] < hi[u]) {
			v = qi[u][ti[u]];
			if (f[v]) break;
			ti[u]++;
		}
		if (ti[u] >= hi[u]) continue; 
		if ((t = prev/u * v) < g ||
			(t == g && v < x)) g = t, x = v;
	}
	return x;
}

int main() {
	int i, k, N, one = 0;
	N = in();
	for (i = 0; i < N; i++) {
		k = in();
		if (i && k == 1) f[1]++, one++;
		else f[k]++, a[sz++] = k;
	}
	qsort(a+1, sz-1, sizeof(int), cmp);
	mktbl(a[sz-1]);
/*
for (i = 2; i < 20; ++i) {
	printf("(%d) sz=%d: ", i, hi[i]);
	for (k = 0; k < hi[i]; ++k) printf("%d ", qi[i][k]);
	printf("\n");
}
printf("\n");
for (i = 2; i < 20; ++i) {
	printf("(%d) sz=%d: ", i, hj[i]);
	for (k = 0; k < hj[i]; ++k) printf("%d ", qj[i][k]);
	printf("\n");
}
*/
	out(a[0], 0);
	while (one--) out(1, 1);
	i = 1; while (i < sz) {
		if (!f[a[i]]) ++i;
		else {
			k = next(a[i]);
			out(k, 1);
		}
	}
	pc('\n');
	return 0;
}
0