結果

問題 No.312 置換処理
コンテスト
ユーザー bal4u
提出日時 2019-04-19 13:59:54
言語 C
(gcc 13.3.0)
結果
AC  
実行時間 12 ms / 2,000 ms
コード長 1,375 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 252 ms
コンパイル使用メモリ 32,000 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-11-15 12:27:11
合計ジャッジ時間 1,508 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 45
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

// 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