Wilson primes of order n: Difference between revisions

<lang racket>#lang racket

(require math/number-theory)

(define ((wilson-prime? n) p)

(and (>= p n)

(prime? p)

(divides? (sqr p)

(- (* (factorial (- n 1))

(factorial (- p n)))

(expt -1 n)))))

(define primes<11000 (filter prime? (range 1 11000)))

(for ((n (in-range 1 (add1 11))))

(printf "~a: ~a~%" n (filter (wilson-prime? n) primes<11000)))</lang>

<pre>1: (5 13 563)

2: (2 3 11 107 4931)

3: (7)

4: (10429)

5: (5 7 47)

6: (11)

7: (17)

8: ()

9: (541)

10: (11 1109)

11: (17 2713)</pre>