Penta-power prime seeds: Difference between revisions
Content added Content deleted
(Added Go) |
(Realize in F#) |
||
Line 230: | Line 230: | ||
</pre> |
</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}}== |
=={{header|Factor}}== |
||
{{works with|Factor|0.99 2022-04-03}} |
{{works with|Factor|0.99 2022-04-03}} |