A144100 Numbers k such that k is strictly greater than f(k), where f(k) = 1 if k is prime, 2 * rad(k) if 4 divides k and rad(k) otherwise.

2, 3, 5, 7, 8, 9, 11, 13, 16, 17, 18, 19, 23, 24, 25, 27, 29, 31, 32, 36, 37, 40, 41, 43, 45, 47, 48, 49, 50, 53, 54, 56, 59, 61, 63, 64, 67, 71, 72, 73, 75, 79, 80, 81, 83, 88, 89, 90, 96, 97, 98, 99, 100, 101, 103, 104, 107, 108, 109, 112, 113, 117, 120, 121, 125, 126

Offset

1,1

Comments

This is the set of all integers k such that there exists a full period linear congruential pseudorandom number generator x -> bx + c (mod k), where b is not a multiple of k, b - 1 is a multiple of f(k) and c is a positive integer relatively prime to k.

4 is the only prime power not a member of the set: f(4) = 2 * rad(4) = 4.

This sequence consists of the primes and 2*A013929. - Charlie Neder, Jan 28 2019

Links

T. D. Noe, Table of n, a(n) for n = 1..1000

Formula

A144907(a(n)) < a(n). - Reinhard Zumkeller, Mar 12 2014

Example

2 is a member: f(2) = 1 and the sequence (0, 1, 0, ...) given by x -> x + 1 (mod 2) has period 2.
8 is a member: f(8) = 4 and the sequence (0, 1, 6, 7, 4, 5, 2, 3, 0, ...) given by x -> 5x + 1 (mod 8) has period 8.
18 is a member: f(18) = 6 and the sequence (0, 1, 14, 3, 4, 17, 6, 7, 2, 9, 10, 5, 12, 13, 8, 15, 16, 11, 0, ...) given by x -> 13x + 1 (mod 18) has period 18.

Prog

(PARI) rad(n) = local(p); p=factor(n)[, 1]; prod(i=1, length(p), p[i]) ;
f(n) = if (isprime(n), 1, if ((n % 4)==0 , 2*rad(n), rad(n))); isok(n) = n > f(n); \\ Michel Marcus, Aug 09 2013
(Haskell)
a144100 n = a144100_list !! (n-1)
a144100_list = filter (\x -> a144907 x < x) [1..]
-- Reinhard Zumkeller, Mar 12 2014

Crossrefs

Cf. A000040, A007947, A071139, A082377, A086486, A089352, A133810, A133811.

Sequence in context: A358978 A063743 A353968 * A359059 A328320 A086539

Adjacent sequences: A144097 A144098 A144099 * A144101 A144102 A144103

Keyword

easy,nice,nonn

Author

Reikku Kulon, Sep 10 2008

Status

approved