Quad-power prime seeds: Difference between revisions
Content added Content deleted
Thundergnat (talk | contribs) (sigh) |
(Realize in F#) |
||
Line 117: | Line 117: | ||
</pre> |
</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}}== |
=={{header|Factor}}== |
||
{{works with|Factor|0.99 2022-04-03}} |
{{works with|Factor|0.99 2022-04-03}} |