Penta-power prime seeds: Difference between revisions
Content added Content deleted
SqrtNegInf (talk | contribs) (Added Perl) |
|||
Line 147: | Line 147: | ||
1,121,189 1,158,975 1,483,691 1,490,475 1,512,321 |
1,121,189 1,158,975 1,483,691 1,490,475 1,512,321 |
||
</pre> |
</pre> |
||
=={{header|Perl}}== |
|||
{{libheader|ntheory}} |
|||
<syntaxhighlight lang="perl" line>use v5.36; |
|||
use bigint; |
|||
use ntheory 'is_prime'; |
|||
use List::Util 'max'; |
|||
sub comma { reverse ((reverse shift) =~ s/(.{3})/$1,/gr) =~ s/^,//r } |
|||
sub table ($c, @V) { my $t = $c * (my $w = 2 + max map {length} @V); ( sprintf( ('%'.$w.'s')x@V, @V) ) =~ s/.{1,$t}\K/\n/gr } |
|||
my($i,@ppps); |
|||
while (@ppps < 30) { |
|||
my $k = 1 + (my $n = 2 * $i++ + 1); |
|||
push @ppps, comma $n if |
|||
is_prime( 1 + $k) and |
|||
is_prime($n + $k) and |
|||
is_prime($n**2 + $k) and |
|||
is_prime($n**3 + $k) and |
|||
is_prime($n**4 + $k); |
|||
} |
|||
say 'First thirty penta-power prime seeds:'; |
|||
say table(10,@ppps);</syntaxhighlight> |
|||
{{out}} |
|||
<pre>First thirty penta-power prime seeds: |
|||
1 5 69 1,665 2,129 25,739 29,631 62,321 77,685 80,535 |
|||
82,655 126,489 207,285 211,091 234,359 256,719 366,675 407,945 414,099 628,859 |
|||
644,399 770,531 781,109 782,781 923,405 1,121,189 1,158,975 1,483,691 1,490,475 1,512,321</pre> |
|||
=={{header|Phix}}== |
=={{header|Phix}}== |