login
A258189
Pseudoprimes to base 3, written in base 3.
7
10101, 11111, 101121, 220212, 222001, 1022011, 1111221, 2010002, 2101001, 2121001, 10101022, 10122201, 10201001, 10212111, 11111112, 11121201, 12010221, 20210202, 21121121, 100001111, 101010101, 102112011, 110020001, 112112001, 120002211, 122000101, 201201201, 202212002
OFFSET
1,1
LINKS
FORMULA
a(n) = A007089(A005935(n)). - Michel Marcus, Sep 11 2015
MATHEMATICA
base = 3; t = {}; n = 1;
While[Length[t] < 40, n++;
If[! PrimeQ[n] && PowerMod[base, n - 1, n] == 1,
AppendTo[t, FromDigits@IntegerDigits[n, 3]]]]; t
PROG
(PARI) lista(nn, b=3) = {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
CROSSREFS
Cf. A005935 (pseudoprimes to base 3), A007089 (numbers in base 3).
Sequence in context: A043641 A210890 A210892 * A114385 A267680 A267876
KEYWORD
nonn,base
AUTHOR
Abdul Gaffar Khan, Sep 11 2015
EXTENSIONS
a(25) corrected by Georg Fischer, Dec 18 2020
STATUS
approved