OFFSET
0,2
COMMENTS
Let USigma denote the unitary sigma function, A034448.
As in A146891, let PF_p(n) denote the largest power of the prime p dividing n. PF_2 is A006519, and PF_3 is A038500. Furthermore define PF_1(n)=1.
Extension to multi-prime-indices is done by multiplying the corresponding functions: PF_{p,q,..}(n) = PF_p(n)*PF_q(n)*... An example of this is PF_{2,3} = A065331.
[How to compute c(m)]
Case of Base Primes = {2}{3}
c(0)=2^m, b(0)=2^m
c(n)=c(n-1)/PF_2[USigma[b(n-1)]]*PF_3[USigma[b(n-1)]]
b(n)=USigma[b(n-1)]/ PF_2,3[USigma[b(n-1)]]
IF b(k)=1 THEN END
a(m)=c(k)
Sequence gives a(m)
Factorization of term becomes 2^r*3^s.
MAPLE
CROSSREFS
KEYWORD
nonn,uned
AUTHOR
Yasutoshi Kohmoto, Apr 17 2009
EXTENSIONS
More terms from R. J. Mathar, Jun 24 2009
STATUS
approved