|
MAPLE
|
with(numtheory); P:= proc(q) local a, b, c, n, p;
for n from 1 to q do
a:=(n+2)*add(op(2, p)/op(1, p), p=ifactors(n+2)[2]);
b:=(n+1)*add(op(2, p)/op(1, p), p=ifactors(n+1)[2]);
c:=n*add(op(2, p)/op(1, p), p=ifactors(n)[2]);
if a=b+c then print(n); fi; od; end: P(10^9);
|