OFFSET
1,2
COMMENTS
By convention 1 is included as the first term, in order to simplify the statement of certain results concerning this sequence.
LINKS
Douglas Latimer, Table of n, a(n) for n = 1..1000
EXAMPLE
3 is included, as its canonical prime factorization 3^1 contains only an odd number of exponents, all of which are congruent to 1 modulo 3.
81 is in the sequence, because its canonical prime factorization is 3^4, and that one exponent, 4, is congruent to 1 modulo 3.
MATHEMATICA
oneQ[n_]:=Module[{f=FactorInteger[n][[All, 2]]}, OddQ[Length[f]]&&Union[ Mod[ f, 3]]=={1}]; Select[Range[200], oneQ] (* Harvey P. Dale, Jul 03 2019 *)
PROG
(PARI) {plnt=1; k=1; print1(k, ", "); plnt++;
mxind=76 ; mxind++ ; for(k=2, 10^6,
M=factor(k); passes=1;
sz = matsize(M)[1];
for(k=1, sz, if(sz%2 != 1, passes=0; break()); if( M[k, 2] % 3 != 1, passes=0));
if( passes == 1 , print1(k, ", "); plnt++) ; if(mxind == plnt, break() ))}
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Douglas Latimer, Apr 12 2012
STATUS
approved