OFFSET
1,1
COMMENTS
This is a second-order nonlinear ADE. It is known that many linear second-order ADEs have infinitely many solutions (A334261), but nonlinear cases haven't been studied.
LINKS
Nathan Mabey, C Script
EXAMPLE
k = 12 is in the sequence, since for m = 4, we have m' = m'' = 4, so m'm - m'' = 16 - 4 = 12 = k.
PROG
(C) See Link
(MATLAB)
function a = A353872( max_pow_2 )
a = [];
maxad2 = ad(ad(2^max_pow_2));
for m = 1:2^max_pow_2
if length(factor(m)) == 2
d = ad(m); b = ad(d); c = d*m;
k(m) = b - c;
end
end
for n = 1:length(k)
if k(n) > -maxad2;
if isempty(find(a == k(n), 1))
if 1 == length(find(k == k(n)))
a = [a k(n)];
end
end
end
end
a = sort(-a);
end
function y = ad( x )
y = 0;
if(x > 1)
p = factor(x); pu = unique(p);
for n = 1:length(pu);
y = y + (x*length(find(p == pu(n))))/pu(n);
end
end
end % Thomas Scheuerle, Jun 15 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Nathan Mabey, May 08 2022
EXTENSIONS
More terms from Jinyuan Wang, Jun 15 2022
STATUS
approved