A ソフィージェルマン素数が素数であるPよう2P + 1も素数です。たとえば、23も素数であるため、11はソフィージャーマン素数です。ソフィージェルマンの素数を昇順で計算する最短のプログラムを書く
回答:
17文字の場合、最大2 ^ 32の列挙が得られます。
G8#,{_mp*2*)mp},`
4文字以上の場合、範囲は2 ^ 32より大きいSG素数を含めるのに十分な大きさになります。
G8#K_*+,{_mp*2*)mp},`
4294967681以降= 2 ^ 32 + 385 <2 ^ 32 + 400。
もちろん、無料で範囲を等しく拡張することもできます。
C9#,{_mp*2*)mp},`
I,はI符号付き32ビット整数として扱われるため、の最大値はIです2 ** 31 - 1。
オンラインコンパイラ/エグゼキュータは無限ループであるため、出力を表示しないことに注意してください。
K1#~K1I&!tPK!tPhyKK
説明：
K1 K=1
# While true:
~K1 K+=1
I If
& logical AND
!tPK K is prime
!tPhyK 2*K+1 is prime (y is double, h is +1)
K Print K
PZ真実または偽の値を返しません。の素因数分解を返しますZ。素数のテストは!tPZであり、素因数分解に1つの因子のみが含まれているかどうかをチェックします。
!tP間違い0や1素数になることになります。簡単な修正は、すべてZを最初に置き換えてK割り当てることK2です。
K1代わりに割り当てK2、ifと増分を入れ替えます。このようにして、を削除できます)。そして+1*K2と同じですhyK。
Pythでの通常の素数チェックの方法を使用し!tPて、数値とその安全素に両方を適用します。両方を一度にチェックするためのちょっとしたトリックがあります。まで上がる10^10ので、ボーナスを獲得します。
f!+tPTtPhyTr2^TT
説明はすぐに来ます。
f r2^TT Filter from 2 till 10^10
! Logical not to detect empty lists
+ List concatenation
tP All but the firs element of the prime factorization
T The filter element
tP All but the firs element of the prime factorization
hyT 2n+1
#include<stdio.h>
#include<math.h>
int isprime(int);
int main(){
int check,n,secondcheck;
printf("enter how long you want to print\n");
scanf("%d",&n);
for(int i=2;i<n;i++){
check = isprime(i);
if(check==0){
secondcheck = isprime(2*i+1);
if(secondcheck==0){
printf("%d\t",i);
}
else
continue;
}
}
}
int isprime(int num){
int check = num,flag=0;
num = sqrt(num);
for(int i=2;i<=num;i++){
if(check%i==0){
flag=1;
return 1;
}
}
if(flag==0){
return 0;
}
}
i=3
q=[2]
while 1:
p=1
for x in range(2,round(i**.5)+1):p=min(p,i%x)
if p:
q+=[i];s=(i-1)/2
if s in q:print(s)
i+=2それはどのように機能しますか？
i=3 # Start at 3
q=[2] # Create list with first prime (2), to be list of primes.
while 1: # Loop forever
p=1 # Set p to 1 (true)
for x in range(2,round(i**0.5)+1): # Loop from 2 to the number's square root. x is the loop value
p=min(p,i%x) # Set p to the min of itself and the modulo of
# the number being tested and loop value (x).
# If p is 0 at the end, a modulo was 0, so it isn't prime.
if p: # Check if p is 0
q+=[i] # Add the current number (we know is prime) to list of primes (q)
s=(i-1)/2 # Generate s, the number that you would double and add 1 to make a prime.
if s in q:print(s) # If (i-1)/2 is a prime (in the list), then that prime satifies
# the condition 2p+1 is prime because i is 2p+1, and i is prime
i+=2 # Increment by 2 (no even numbers are prime, except 2)私のコンピューターによると、ソフィジェルマンの素数が2 ^ 32未満の0.023283％が生成されています。
それが終わったら、十分な行があればペーストビンに投稿します。あなたはそれらをすべて持っていることを確認するためにそれを使用できます。
.5より短い0.5