結果

問題 No.3030 ミラー・ラビン素数判定法のテスト
ユーザー furafura
提出日時 2020-12-11 17:48:50
言語 C++17
(gcc 13.2.0 + boost 1.83.0)
結果
RE  
実行時間 -
コード長 828 bytes
コンパイル時間 2,265 ms
コンパイル使用メモリ 199,996 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-08-11 23:30:18
合計ジャッジ時間 4,264 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge12
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 AC 1 ms
4,380 KB
testcase_02 AC 1 ms
4,380 KB
testcase_03 AC 1 ms
4,380 KB
testcase_04 RE -
testcase_05 AC 328 ms
4,376 KB
testcase_06 AC 51 ms
4,376 KB
testcase_07 AC 48 ms
4,376 KB
testcase_08 AC 50 ms
4,380 KB
testcase_09 AC 705 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

#define rep(i,n) for(int i=0;i<(n);i++)

using namespace std;
using lint=long long;

random_device seed_gen;
mt19937 rng(seed_gen());

lint modpow(lint a,lint k,lint m){
	lint r=1%m;
	for(lint x=a%m;k>0;k>>=1,x=__int128(x)*x%m) if(k&1) r=__int128(r)*x%m;
	return r;
}

bool Miller_Rabin(lint n,int k=20){
	if(n<=1)   return false;
	if(n%2==0) return n==2;

	int s=0;
	lint d=n-1;
	while(~d&1) s++, d>>=1;

	while(k--){
		lint a=rng()%(n-3)+2;
		lint x=modpow(a,d,n);
		if(x==1 || x==n-1) continue;
		bool b=false;
		for(int r=1;r<s;r++){
			x=__int128(x)*x%n;
			if(x== 1 ) return false;
			if(x==n-1){ b=true; break; }
		}
		if(!b) return false;
	}

	return true;
}

int main(){
	int q; scanf("%d",&q);
	rep(_,q){
		lint x; scanf("%lld",&x);
		printf("%lld %d\n",x,!!Miller_Rabin(x));
	}
	return 0;
}
0