Quadrat special primes: Difference between revisions
Content added Content deleted
(Add Factor) |
Catskill549 (talk | contribs) (added AWK) |
||
Line 108: | Line 108: | ||
10667 11243 11279 11423 12323 12647 12791 13367 13691 14591 14627 14771 |
10667 11243 11279 11423 12323 12647 12791 13367 13691 14591 14627 14771 |
||
15671 |
15671 |
||
</pre> |
|||
=={{header|AWK}}== |
|||
<lang AWK> |
|||
# syntax: GAWK -f QUADRAT_SPECIAL_PRIMES.AWK |
|||
# converted from FreeBASIC |
|||
BEGIN { |
|||
stop = 15999 |
|||
p = 2 |
|||
j = 1 |
|||
printf("%5d ",p) |
|||
count++ |
|||
while (1) { |
|||
while (1) { |
|||
if (is_prime(p+j*j)) { break } |
|||
j++ |
|||
} |
|||
p += j*j |
|||
if (p > stop) { break } |
|||
printf("%5d%1s",p,++count%10?"":"\n") |
|||
j = 1 |
|||
} |
|||
printf("\nQuadrat special primes 1-%d: %d\n",stop,count) |
|||
exit(0) |
|||
} |
|||
function is_prime(x, i) { |
|||
if (x <= 1) { |
|||
return(0) |
|||
} |
|||
for (i=2; i<=int(sqrt(x)); i++) { |
|||
if (x % i == 0) { |
|||
return(0) |
|||
} |
|||
} |
|||
return(1) |
|||
} |
|||
</lang> |
|||
{{out}} |
|||
<pre> |
|||
2 3 7 11 47 83 227 263 587 911 |
|||
947 983 1019 1163 1307 1451 1487 1523 1559 2459 |
|||
3359 4259 4583 5483 5519 5843 5879 6203 6779 7103 |
|||
7247 7283 7607 7643 8219 8363 10667 11243 11279 11423 |
|||
12323 12647 12791 13367 13691 14591 14627 14771 15671 |
|||
Quadrat special primes 1-15999: 49 |
|||
</pre> |
</pre> |
||