Quad-power prime seeds: Difference between revisions

NB: The source of the ALGOL 68-primes library is on a Rosetta Code code page linked from the above.<br>

Note that to run this with ALGOL 68G under Windows (and probably Linux) a large heap size must be specified on the command line, e.g. <code>-heap 1024M</code>.

# n^p + n + 1 is prime for p = 1, 2, 3, 4 #

PR read "primes.incl.a68" PR # include prime utilities #

FI

OD

{{out}}

<pre>

 

=={{header|F_Sharp|F#}}==

// Quad-power prime seeds. Nigel Galloway: August 22nd., 2022

let fG n g=let n=bigint(n:int) in let n=n**g+n+1I in Open.Numeric.Primes.MillerRabin.IsProbablePrime &n

Seq.initInfinite((+)1)|>Seq.filter(fun n->let g=fG n in g 1&&g 2&&g 3&&g 4)|>Seq.take 50|>Seq.iter(printf "%d "); printfn "\n"

[1000000..1000000..10000000]|>Seq.scan(fun(n,g,x) l->let i,e=fN(g,l) in (n+i,e,l))(0,0,0)|>Seq.skip 1|>Seq.iter(fun(n,g,l)->printfn $"First element over %8d{l} is %9d{g} at index %3d{n}")

{{out}}

<pre>

=={{header|Factor}}==

{{works with|Factor|0.99 2022-04-03}}

math.primes prettyprint sequences tools.memory.private ;

 

 

"First fifty quad-power prime seeds:" print

{{out}}

<pre>

 

=={{header|Phix}}==

<span style="color: #008080;">with</span> <span style="color: #008080;">javascript_semantics</span>

<span style="color: #008080;">include</span> <span style="color: #004080;">mpfr</span><span style="color: #0000FF;">.</span><span style="color: #000000;">e</span>

<span style="color: #000000;">n</span> <span style="color: #0000FF;">+=</span> <span style="color: #000000;">1</span>

<span style="color: #008080;">end</span> <span style="color: #008080;">while</span>

{{out}}

<pre>

 

=={{header|Raku}}==

 

my @qpps = lazy (1..*).hyper(:5000batch).grep: -> \n { my \k = n + 1; (n+k).is-prime && (n²+k).is-prime && (n³+k).is-prime && (n⁴+k).is-prime }

my $key = @qpps.first: * > $threshold, :k;

say "{$threshold.&cardinal.fmt: '%13s'} is the {ordinal-digit $key + 1}: {@qpps[$key].&comma}";

{{out}}

<pre>First fifty quad-power prime seeds:

{{libheader|Wren-fmt}}

GMP allows us to stretch a little more.

import "./fmt" for Fmt

 

}

n = n + 1

 

{{out}}