A141231 - OEIS

%I #14 Dec 11 2018 07:07:45

%S 3,5,9,11,13,15,19,25,27,29,33,37,39,45,53,55,57,59,61,65,67,75,81,83,

%T 87,95,99,101,107,111,117,121,125,131,135,139,143,145,149,159,163,165,

%U 169,171,173,177,179,181,183,185,195,197,201,209,211,225,227,243,247,249

%N Positive integers whose prime factorizations contain only primes from A001122.

%t aQ[n_] := Length[Select[FactorInteger[n][[;; , 1]], # > 1 && PrimitiveRoot@# != 2 &]] == 0; Select[Range[2, 250], aQ] (* _Amiram Eldar_, Dec 09 2018 *)

%o (PARI) isokp(p) = znorder(Mod(2, p))==(p-1);

%o isok(n) = {if ((n>1) && (n % 2), my(f=factor(n)); #select(x->isokp(x), f[, 1]) == #f~; , 0); } \\ _Michel Marcus_, Dec 09 2018

%Y Cf. A001122, A000040.

%K nonn

%O 1,1

%A _Vladimir Shevelev_, Jun 16 2008

%E More terms from _Michel Marcus_, Dec 09 2018