OFFSET
1,3
COMMENTS
A046311 is a subset of this sequence.
Some numbers admit more than one couple of divisors a, b: 3456 = 8 * 432 = 54 * 64 and 3456' = 15552 = 8' * 432' = 12 * 1296 = 54' * 64' = 81 * 192.
LINKS
Paolo P. Lava, Table of n, a(n) for n = 1..100
EXAMPLE
a(0) = 0 because 0 = 0 * b and 0' = 0' * b' = 0;
a(1) = 1 because 1 = 1 * 1 and 1' = 1' * 1' = 0;
a(2) = 256 because 256 = 16 * 16 and 256' = 16' * 16' = 32 * 32 = 1024;
a(3) = 512 because 512 = 8 * 64 and 512' = 8' * 64' = 12 * 192 = 2304.
MAPLE
with(numtheory): P:=proc(q) local a, b, c, j, k, n, p;
for n from 1 to q do j:=sort([op(divisors(n))]);
for k from 2 to trunc((nops(j)+1)/2) do
a:=j[k]*add(op(2, p)/op(1, p), p=ifactors(j[k])[2]);
b:=(n/j[k])*add(op(2, p)/op(1, p), p=ifactors(n/j[k])[2]);
c:=n*add(op(2, p)/op(1, p), p=ifactors(n)[2]);
if c=a*b then print(n); break; fi; od; od; end: P(10^6);
CROSSREFS
KEYWORD
nonn
AUTHOR
Paolo P. Lava, Oct 24 2017
STATUS
approved