OFFSET
1,1
COMMENTS
Subsequence of A144711.
LINKS
Daniel Lignon, Table of n, a(n) for n = 1..1000
MAPLE
filter:= proc(n)
local P, p;
P:= numtheory:-factorset(n);
nops(P) > 1 and issqr(add(p^2, p=P)/nops(P))
end proc:
select(filter, [$1..10000]); # Robert Israel, Feb 26 2015
MATHEMATICA
Complement[Select[Range[2, 5000], IntegerQ[RootMeanSquare[Select[Divisors[#], PrimeQ]]]&], Select[Range[2, 5000], Length[FactorInteger[#]]==1&]] (* Daniel Lignon, Feb 26 2015 *)
PROG
(PARI) isok(n) = ((nbp=omega(n)) > 1) && (f=factor(n)) && (x = sum(k=1, nbp, f[k, 1]^2)/nbp) && issquare(x) && (type(x) == "t_INT"); \\ Michel Marcus, Mar 03 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Daniel Lignon, Feb 26 2015
STATUS
approved