結果

問題 No.8030 ミラー・ラビン素数判定法のテスト
ユーザー 👑 testestest
提出日時 2017-10-22 00:03:01
言語 C
(gcc 13.3.0)
結果
AC  
実行時間 6,841 ms / 9,973 ms
コード長 1,127 bytes
コンパイル時間 1,114 ms
コンパイル使用メモリ 30,208 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-11-16 22:34:45
合計ジャッジ時間 19,161 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 10
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.c:39:1: warning: return type defaults to 'int' [-Wimplicit-int]
   39 | main(){
      | ^~~~
main.c: In function 'main':
main.c:40:9: warning: 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: warning: 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 # 

#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;
}

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