Sexy primes: Difference between revisions

{{Wikipedia|Sexy_prime}}

 

In mathematics, '''sexy primes''' are prime numbers that differ from each other by six.

 

 

The term "sexy prime" is a pun stemming from the Latin word for six: ''sex''.<br><br>

 

::*For each of pairs, triplets, quadruplets and quintuplets, Find and display the count of each group type of sexy primes less than one million thirty-five ('''1,000,035''').

::*Find and display the count of the unsexy primes less than one million thirty-five.

 

 

 

=={{header|Factor}}==

IN: rosetta-code.sexy-primes

 

CONSTANT: tuplet-names qw{ pair triplet quadruplet quintuplet }

 

: tuplet ( m n -- seq ) dupd 1 - 6 * + 6 <range> ;

 

: sexy-tuplets ( n -- seq ) [ primes ] dip '[

] { } make ;

 

: show-tuplets ( n -- )

 

: unsexy-primes ( -- seq ) primes [

 

: show-unsexy ( -- )

 

: main ( -- ) 2 5 [a,b] [ show-tuplets ] each show-unsexy ;

 

{{out}}

<pre>

 

 

 

 

</pre>

 

=={{header|Go}}==

 

import "fmt"

 

limit++

// True denotes composite, false denotes prime.

c[1] = true

// no need to bother with even numbers over 2 for this task

for {

p2 := p * p

}

 

}

}

 

func main() {

var pairs [][2]int

var trips [][3]int

var quins [][5]int

var unsexy = []int{2, 3}

unsexy = append(unsexy, i)

continue

}

pair := [2]int{i, i + 6}

pairs = append(pairs, pair)

continue

}

trip := [3]int{i, i + 6, i + 12}

trips = append(trips, trip)

continue

}

quad := [4]int{i, i + 6, i + 12, i + 18}

quads = append(quads, quad)

continue

}

quin := [5]int{i, i + 6, i + 12, i + 18, i + 24}

quins = append(quins, quin)

}

}

 

 

 

 

 

{{out}}

<pre>

 

 

 

 

</pre>

 

 

 

 

 

 

}

 

 

}

 

{{out}}

 

 

 

 

 

=={{header|Python}}==

def primes2(limit=LIMIT):

if limit < 2: return []

return [2] + [i + i + 3 for i, v in enumerate(buf) if v]

 

primeset = set(primes)

#%%

s = [[] for x in range(4)]

unsexy = []

for p in primes:

s[0].append((p, p+6))

else:

if p - 6 not in primeset:

unsexy.append(p)

continue

s[1].append((p, p+6, p+12))

else:

continue

s[2].append((p, p+6, p+12, p+18))

else:

continue

s[3].append((p, p+6, p+12, p+18, p+24))

#%%

print('"SEXY" PRIME GROUPINGS:')

for sexy, name in zip(s, 'pairs triplets quadruplets quintuplets'.split()):

for sx in sexy[-5:]:

print(' ',sx)

print(f'\nThere are {len(unsexy)} unsexy primes ending with ...')

for usx in unsexy[-10:]:

 

{{out}}

(5, 11, 17, 23, 29)

 

999853

999863

999961

999979

 

=={{header|REXX}}==

parse arg N endU end2 end3 end4 end5 . /*obtain optional argument from the CL.*/

if endU=='' | endU=="," then endU= 10 /* " " " " " " */

if end2=='' | end2=="," then end2= 5 /* " " " " " " */

if end4=='' | end4=="," then end4= 5 /* " " " " " " */

if end5=='' | end5=="," then end4= 5 /* " " " " " " */

say 'There are ' commas(w2%2) @sexy "pairs less than " Nc

say

say 'There are ' commas(w3%3) @sexy "triplets less than " Nc

say

say 'There are ' commas(w4%4) @sexy "quadruplets less than " Nc

say

say

exit /*stick a fork in it, we're all done. */

/*──────────────────────────────────────────────────────────────────────────────────────*/

commas: procedure; parse arg _; n= _'.9'; #= 123456789; b= verify(n, #, "M")

e= verify(n, #'0', , verify(n, #"0.", 'M') ) - 4

{{out|output|text=&nbsp; when using the default inputs:}}

 

(Shown at &nbsp; <big>'''⁵/₆'''</big> &nbsp; size.)

<pre style="font-size:83%">

 

The last 5 sexy prime pairs are:

999371~999377 999431~999437 999721~999727 999763~999769 999953~999959

 

The last 5 sexy prime triplets are:

997427~997433~997439 997541~997547~997553 998071~998077~998083 998617~998623~998629 998737~998743~998749

 

The last 5 sexy prime quadruplets are:

977351~977357~977363~977369 983771~983777~983783~983789 986131~986137~986143~986149 990371~990377~990383~990389 997091~997097~997103~997109

 

5~11~17~23~29

 

</pre>

 

 

[[Extensible prime generator#zkl]] could be used instead.

while(p.nextPrime()<=M){ ps[p]=1 } // bitmap of primes

 

 

 

 

 

 

{{out}}

The last 5 are:

(999371,999377) (999431,999437) (999721,999727) (999763,999769) (999953,999959)

 

The last 5 are:

(997427,997433,997439) (997541,997547,997553) (998071,998077,998083) (998617,998623,998629) (998737,998743,998749)

 

The last 5 are:

(977351,977357,977363,977369) (983771,983777,983783,983789) (986131,986137,986143,986149) (990371,990377,990383,990389) (997091,997097,997103,997109)

 

(5,11,17,23,29)

 

The last 10 are:

</pre>