結果

問題 No.312 置換処理
ユーザー bal4ubal4u
提出日時 2019-04-19 13:59:54
言語 C
(gcc 12.3.0)
結果
AC  
実行時間 13 ms / 2,000 ms
コード長 1,375 bytes
コンパイル時間 150 ms
コンパイル使用メモリ 32,000 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-04-27 11:01:30
合計ジャッジ時間 1,351 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
6,812 KB
testcase_01 AC 6 ms
6,940 KB
testcase_02 AC 13 ms
6,940 KB
testcase_03 AC 8 ms
6,944 KB
testcase_04 AC 1 ms
6,940 KB
testcase_05 AC 1 ms
6,940 KB
testcase_06 AC 1 ms
6,944 KB
testcase_07 AC 4 ms
6,944 KB
testcase_08 AC 1 ms
6,940 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 AC 1 ms
6,944 KB
testcase_14 AC 1 ms
6,940 KB
testcase_15 AC 1 ms
6,944 KB
testcase_16 AC 1 ms
6,944 KB
testcase_17 AC 1 ms
6,940 KB
testcase_18 AC 1 ms
6,940 KB
testcase_19 AC 1 ms
6,944 KB
testcase_20 AC 1 ms
6,940 KB
testcase_21 AC 1 ms
6,940 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,940 KB
testcase_26 AC 1 ms
6,940 KB
testcase_27 AC 1 ms
6,944 KB
testcase_28 AC 1 ms
6,944 KB
testcase_29 AC 1 ms
6,940 KB
testcase_30 AC 2 ms
6,940 KB
testcase_31 AC 1 ms
6,940 KB
testcase_32 AC 2 ms
6,940 KB
testcase_33 AC 5 ms
6,940 KB
testcase_34 AC 4 ms
6,944 KB
testcase_35 AC 1 ms
6,944 KB
testcase_36 AC 1 ms
6,940 KB
testcase_37 AC 1 ms
6,944 KB
testcase_38 AC 1 ms
6,940 KB
testcase_39 AC 1 ms
6,940 KB
testcase_40 AC 1 ms
6,944 KB
testcase_41 AC 1 ms
6,944 KB
testcase_42 AC 1 ms
6,944 KB
testcase_43 AC 1 ms
6,944 KB
testcase_44 AC 1 ms
6,944 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);
	if (N % 3 == 0) ans = 3;
	else if ((N & 3) == 0) ans = 4;
	else {
		ans = N;
		if ((N & 1) == 0) N >>= 1, ans = N;
		if (!prime_test(N)) {
			b = (int)sqrt((double)N);
			for (a = 5; a <= b; a++) {
				if (N % a == 0) break;
			}
			if (a <= b) ans = a;
		}
	}
	printf("%lld\n", ans);
	return 0;
}
0