OFFSET
1,1
COMMENTS
Includes 2^(3*a) * 3^(4*b) if 3*a >= 4*b. - Robert Israel, Mar 30 2023
LINKS
Robert Israel, Table of n, a(n) for n = 1..155
EXAMPLE
Prime factors of 3960 are 2^3, 3^2, 5 and 11.
Sum_{i=1..7} (p(i)/(p(i)+1)) = 3*(2/(2+1)) + 2*(3/(3+1)) + 5/(5+1) + 11/(11+1) = 21/4.
Product_{i=1..7} (p(i)/(p(i)-1)) = (2/(2+1))^3 * (3/(3-1))^2 * 5/(5-1) * 11/(11-1) = 99/4.
Their sum is an integer: 21/4 + 99/4 = 30.
MAPLE
with(numtheory); P:=proc(i) local b, d, n, p;
for n from 2 to i do p:=ifactors(n)[2];
b:=add(op(2, d)*op(1, d)/(op(1, d)+1), d=p)+mul((op(1, d)/(op(1, d)-1))^op(2, d), d=p);
if trunc(b)=b then print(n); fi; od; end: P(10^7);
CROSSREFS
KEYWORD
nonn
AUTHOR
Paolo P. Lava, Oct 09 2013
STATUS
approved