結果

問題 No.8030 ミラー・ラビン素数判定法のテスト
ユーザー testestest
提出日時 2017-10-22 00:00:09
言語 C(gnu17)
(gcc 15.2.0)
コンパイル:
gcc-15 -O2 -std=gnu17 -Wno-error=implicit-function-declaration -Wno-error=implicit-int -Wno-error=incompatible-pointer-types -Wno-error=int-conversion -DONLINE_JUDGE -o a.out _filename_ -lm
実行:
./a.out
結果
CE  
(最新)
AC  
(最初)
実行時間 -
コード長 1,128 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 121 ms
コンパイル使用メモリ 28,052 KB
最終ジャッジ日時 2026-02-22 00:21:25
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
コンパイルエラー時のメッセージ・ソースコードは、提出者また管理者しか表示できないようにしております。（リジャッジ後のコンパイルエラーは公開されます）
ただし、clay言語の場合は開発者のデバッグのため、公開されます。

コンパイルメッセージ
main.c:38:1: error: return type defaults to 'int' [-Wimplicit-int]
   38 | main(){
      | ^~~~
main.c: In function 'main':
main.c:40:9: error: implicit declaration of function 'scanf' [-Wimplicit-function-declaration]
   40 |         scanf("%d",&n);
      |         ^~~~~
main.c:1:1: note: include '<stdio.h>' or provide a declaration of 'scanf'
  +++ |+#include <stdio.h>
    1 | #define ll long long
main.c:40:9: warning: incompatible implicit declaration of built-in function 'scanf' [-Wbuiltin-declaration-mismatch]
   40 |         scanf("%d",&n);
      |         ^~~~~
main.c:40:9: note: include '<stdio.h>' or provide a declaration of 'scanf'
main.c:43:17: error: implicit declaration of function 'printf' [-Wimplicit-function-declaration]
   43 |                 printf("%lld %d\n",x,isp(x));
      |                 ^~~~~~
main.c:43:17: note: include '<stdio.h>' or provide a declaration of 'printf'
main.c:43:17: warning: incompatible implicit declaration of built-in function 'printf' [-Wbuiltin-declaration-mismatch]
main.c:43:17: note: include '<stdio.h>' or provide a declaration of 'printf'

ソースコード

diff #
raw source code 

#define ll long long
#define ull unsigned long long
#define bool int
#define T 1
#define F 0
ll pom(ll a,ll n,int m){ll x=1;for(a%=m;n;n/=2)n%2?x=x*a%m:0,a=a*a%m;return x;}
ll mumll(ll p,ll q,ll m){ull x=0;for(;q;q/=2)q&1?x=(x+p)%m:0,p=((ull)p+p)%m;return(ll)x;}
ll pomll(ll a,ll n,ll m){ll x=1;for(;n;n/=2)n&1?x=mumll(x,a,m):0,a=mumll(a,a,m);return x;}

//ミラーラビン
ll sus[]={2,7,61,0};
//ll susll[]={2,3,5,7,11,13,17,19,23,29,31,37,0};
ll susll[]={2,325,9375,28178,450775,9780504,1795265022,0};
bool isp(ll n){
	if(n==2)return T;
	if(n<2||n%2==0)return F;
	ll d=(n-1)/(n-1&1-n);
	if(n<(1LL<<31)){
		for(int i=0;sus[i];i++){
			if(sus[i]>=n)return T;
			ll t=d,a=pom(sus[i],t,n);
			if(a==1||a==n-1)continue;
			while((t*=2)!=n-1&&(a=a*a%n)!=n-1&&a!=1);
			if(t==n-1||a==1)return F;
		}
		return T;
	}
	for(int i=0;susll[i];i++){
		if(susll[i]%n==0)continue;
		ll t=d,a=pomll(susll[i],t,n);
		if(a==1||a==n-1)continue;
		while((t*=2)!=n-1&&(a=mumll(a,a,n))!=n-1&&a!=1);
		if(t==n-1||a==1)return F;
	}
	return T;
}

main(){
	ll n,x;
	scanf("%d",&n);
	while(n--){
		scanf("%lld",&x);
		printf("%lld %d\n",x,isp(x));
	}
}
0