OFFSET
1,2
COMMENTS
This sequence is infinite as it contains the positive terms of A000225.
All terms are odd.
Run lengths in first differences appear to be regular and suggest a simple procedure to generate the sequence.
LINKS
EXAMPLE
The binary representation of 3^3 ("11011") ends with that of 3 ("11"), so 3 is a term.
MATHEMATICA
Select[Range[1200], Take[IntegerDigits[#^#, 2], -IntegerLength[ #, 2]] == IntegerDigits[ #, 2]&] (* Harvey P. Dale, Jan 12 2022 *)
PROG
(PARI) is(n, base=2) = Mod(n, base^#digits(n, base))^n==n
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Rémy Sigrist, Oct 11 2020
STATUS
approved