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