Quad-power prime seeds: Difference between revisions
Realize in F#
Thundergnat (talk | contribs) (sigh) |
(Realize in F#) |
||
Line 117:
</pre>
=={{header|F_Sharp|F#}}==
<lang fsharp>
// 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
let fN(n,g)=Seq.initInfinite((+)n)|>Seq.filter(fun n->let g=fG n in g 1&&g 2&&g 3&&g 4)|>Seq.mapi(fun n g->(n,g))|>Seq.find(snd>>(<)g)
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}")
</lang>
{{out}}
<pre>
1 2 5 6 69 131 426 1665 2129 2696 5250 7929 9689 13545 14154 14286 16434 19760 25739 27809 29631 36821 41819 46619 48321 59030 60500 61955 62321 73610 77685 79646 80535 82655 85251 86996 91014 96566 97739 105939 108240 108681 119754 122436 123164 126489 140636 150480 153179 163070
First element over 1000000 is 1009286 at index 140
First element over 2000000 is 2015496 at index 233
First element over 3000000 is 3005316 at index 318
First element over 4000000 is 4004726 at index 382
First element over 5000000 is 5023880 at index 451
First element over 6000000 is 6000554 at index 513
First element over 7000000 is 7047129 at index 566
First element over 8000000 is 8005710 at index 600
First element over 9000000 is 9055151 at index 644
First element over 10000000 is 10023600 at index 700
</pre>
=={{header|Factor}}==
{{works with|Factor|0.99 2022-04-03}}
| |||