# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/ Search: id:a141231 Showing 1-1 of 1 %I A141231 #14 Dec 11 2018 07:07:45 %S A141231 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 A141231 87,95,99,101,107,111,117,121,125,131,135,139,143,145,149,159,163,165, %U A141231 169,171,173,177,179,181,183,185,195,197,201,209,211,225,227,243,247,249 %N A141231 Positive integers whose prime factorizations contain only primes from A001122. %t A141231 aQ[n_] := Length[Select[FactorInteger[n][[;; , 1]], # > 1 && PrimitiveRoot@# != 2 &]] == 0; Select[Range[2, 250], aQ] (* _Amiram Eldar_, Dec 09 2018 *) %o A141231 (PARI) isokp(p) = znorder(Mod(2, p))==(p-1); %o A141231 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 A141231 Cf. A001122, A000040. %K A141231 nonn %O A141231 1,1 %A A141231 _Vladimir Shevelev_, Jun 16 2008 %E A141231 More terms from _Michel Marcus_, Dec 09 2018 # Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE