A262101 Pseudoprimes to base 4, written in base 4.
33, 1111, 1123, 11111, 12303, 13003, 20301, 22011, 22333, 101101, 103133, 103313, 111223, 120231, 122133, 123001, 131203, 131301, 133333, 200113, 212201, 222031, 230011, 300331, 303031, 310213, 321203, 333001, 1010101, 1010103, 1021021, 1022323, 1023323, 1111111, 1112233, 1213021, 1213303, 1330111, 2002001, 2010201, 2013313, 2023033, 2031211, 2032223
Offset: 1
Links
- Amiram Eldar, Table of n, a(n) for n = 1..10000
Crossrefs
Programs
-
Mathematica
BaseForm[Select[Range[4096], Not[PrimeQ[#]] && PowerMod[4, # - 1, #] == 1 &], 4]
-
PARI
lista(nn, b=4) = {for (n=1, nn, if (Mod(b, n)^(n-1)==1 && !ispseudoprime(n) && n>1, print1(subst(Pol(digits(n,b), x), x, 10), ", ");););} \\ Michel Marcus, Sep 30 2015