Jump to content

Penta-power prime seeds: Difference between revisions

Realize in F#
(Added Go)
(Realize in F#)
Line 230:
</pre>
 
=={{header|F_Sharp|F#}}==
<syntaxhighlight lang="fsharp">
// Penta-power prime seeds. Nigel Galloway: April 5th., 2023
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 0&&g 1&&g 2&&g 3&&g 4)|>Seq.mapi(fun n g->(n,g))|>Seq.find(snd>>(<)g)
Seq.initInfinite((*)2>>(+)1)|>Seq.filter(fun n->let g=fG n in g 0&&g 1&&g 2&&g 3&&g 4)|>Seq.take 30|>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}")
</syntaxhighlight>
{{out}}
<pre>
1 5 69 1665 2129 25739 29631 62321 77685 80535 82655 126489 207285 211091 234359 256719 366675 407945 414099 628859 644399 770531 781109 782781 923405 1121189 1158975 1483691 1490475 1512321
 
First element over 1000000 is 1121189 at index 25
First element over 2000000 is 2066079 at index 38
First element over 3000000 is 3127011 at index 46
First element over 4000000 is 4059525 at index 50
First element over 5000000 is 5279175 at index 58
First element over 6000000 is 6320601 at index 62
First element over 7000000 is 7291361 at index 67
First element over 8000000 is 8334915 at index 68
First element over 9000000 is 9100671 at index 70
First element over 10000000 is 10347035 at index 71
</pre>
=={{header|Factor}}==
{{works with|Factor|0.99 2022-04-03}}
2,172

edits

Cookies help us deliver our services. By using our services, you agree to our use of cookies.