結果

問題 No.312 置換処理
ユーザー bal4ubal4u
提出日時 2019-04-19 13:51:38
言語 C
(gcc 12.3.0)
結果
WA  
実行時間 -
コード長 1,304 bytes
コンパイル時間 703 ms
コンパイル使用メモリ 32,000 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-09-22 11:26:43
合計ジャッジ時間 1,764 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
6,816 KB
testcase_01 AC 8 ms
6,940 KB
testcase_02 AC 11 ms
6,944 KB
testcase_03 AC 11 ms
6,944 KB
testcase_04 AC 1 ms
6,940 KB
testcase_05 AC 1 ms
6,944 KB
testcase_06 AC 1 ms
6,940 KB
testcase_07 AC 4 ms
6,940 KB
testcase_08 AC 1 ms
6,944 KB
testcase_09 AC 1 ms
6,944 KB
testcase_10 AC 1 ms
6,944 KB
testcase_11 AC 2 ms
6,944 KB
testcase_12 AC 1 ms
6,944 KB
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 AC 1 ms
6,940 KB
testcase_17 AC 1 ms
6,944 KB
testcase_18 AC 1 ms
6,944 KB
testcase_19 WA -
testcase_20 AC 1 ms
6,944 KB
testcase_21 AC 1 ms
6,944 KB
testcase_22 AC 1 ms
6,940 KB
testcase_23 AC 1 ms
6,944 KB
testcase_24 AC 1 ms
6,940 KB
testcase_25 AC 1 ms
6,944 KB
testcase_26 AC 1 ms
6,944 KB
testcase_27 AC 1 ms
6,944 KB
testcase_28 AC 1 ms
6,940 KB
testcase_29 AC 1 ms
6,940 KB
testcase_30 AC 3 ms
6,948 KB
testcase_31 WA -
testcase_32 AC 2 ms
6,940 KB
testcase_33 AC 4 ms
6,940 KB
testcase_34 AC 4 ms
6,940 KB
testcase_35 AC 1 ms
6,944 KB
testcase_36 AC 1 ms
6,944 KB
testcase_37 AC 1 ms
6,940 KB
testcase_38 WA -
testcase_39 AC 1 ms
6,944 KB
testcase_40 AC 1 ms
6,944 KB
testcase_41 AC 1 ms
6,940 KB
testcase_42 AC 1 ms
6,940 KB
testcase_43 AC 1 ms
6,940 KB
testcase_44 AC 1 ms
6,940 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

// yukicoder: No.312 置換処理
// 2019.4.19 bal4u

#include <stdio.h>
#include <math.h>

//// 素数テスト
long long calc_pow(int a, long long k, long long n)
{
	long long p;
	unsigned long long bit;

	p = 1;
	for (bit = 0x8000000000000000U; bit > 0; bit >>= 1) {
		if (p > 1) p = (p*p) % n;
		if (k & bit) p = (p*a) % n;
	}
	return p % n;
}

int suspect(int b, long long n)
{
	int i, t;
	long long u, n1;
	long long x0, x1;

	u = n1 = n - 1;
	t = 0;
	while ((u & 1) == 0) {
		t++;
		u >>= 1;
	}
	x0 = calc_pow(b, u, n);
	for (i = 1; i <= t; i++) {
		x1 = (x0*x0) % n;
		if (x1 == 1 && x0 != 1 && x0 != n1) return 0;
		x0 = x1;
	}
	return x1 == 1;
}

int prime_test(long long n)
{
	if (n == 1) return 0;          /* 1 */
	if (n == 2 || n == 3) return 1;          /* 2, 3 */
	if ((n & 1) == 0) return 0;  /* other even integers */
	if (!suspect(2, n)) return 0;
	if (n <= 1000000) {
		if (!suspect(3, n)) return 0;
	}
	else {
		if (!suspect(7, n)) return 0;
		if (!suspect(61, n)) return 0;
	}
	return 1;
}

int main()
{
	int b, a;
	long long N, ans;

	scanf("%lld", &N);
	ans = N, b = (int)sqrt((double)N);
	if ((N & 1) == 0 && N != 4) N >>= 1, ans = N;
	if (!prime_test(N)) {
		for (a = 3; a <= b; a++) {
			if (N % a == 0) break;
		}
		if (a <= b) ans = a;
	}
	printf("%lld\n", ans);
	return 0;
}
0