%I #32 Jun 04 2016 06:59:27
%S 10,65,140,420,2100,2210,20737,32045,200725,207370,1204350,1347905,
%T 1762645,16502850,31427800,37741340,107671200,130643100,200728169,
%U 239719720,357491225,417225900,430085380,766750575,1088692500,1132409168,1328204850,1788379460
%N Numbers that divide the average of the squares of their aliquot parts.
%C Ratio: 1, 1, 5, 10, 78, 1, 109, 565,...
%C If the ratio is equal to 1 we have 10, 65, 20737 (A140362).
%H Giovanni Resta, <a href="/A255245/b255245.txt">Table of n, a(n) for n = 1..59</a> (terms < 10^11)
%e Aliquot parts of 10 are 1, 2, 5. The average of their squares is (1^2 + 2^2 + 5^2) / 3 = (1 + 4 + 25) / 3 = 30 / 3 = 10 and 10 / 10 = 1.
%p with(numtheory); P:=proc(q) local a,b,k,n;
%p for n from 2 to q do a:=sort([op(divisors(n))]);
%p b:=add(a[k]^2,k=1..nops(a)-1)/(nops(a)-1);
%p if type(b/n,integer) then lprint(n);
%p fi; od; end: P(10^6);
%t Select[Range[10^6],Mod[Mean[Most[Divisors[#]^2]],#]==0&] (* _Ivan N. Ianakiev_, Mar 03 2015 *)
%o (PARI) isok(n) = (q=(sumdiv(n, d, (d!=n)*d^2)/(numdiv(n)-1))) && (type(q)=="t_INT") && ((q % n) == 0); \\ _Michel Marcus_, Feb 20 2015
%o (Python)
%o from __future__ import division
%o from sympy import factorint
%o A255245_list = []
%o for n in range(2,10**9):
%o ....s0 = s2 = 1
%o ....for p,e in factorint(n).items():
%o ........s0 *= e+1
%o ........s2 *= (p**(2*(e+1))-1)//(p**2-1)
%o ....q, r = divmod(s2-n**2,s0-1)
%o ....if not (r or q % n):
%o ........A255245_list.append(n) # _Chai Wah Wu_, Mar 08 2015
%Y Cf. A001065, A255244.
%K nonn
%O 1,1
%A _Paolo P. Lava_, Feb 20 2015
%E More terms from _Michel Marcus_, Feb 20 2015
%E a(17)-a(24) from _Chai Wah Wu_, Mar 08 2015
%E a(25)-a(28) from _Giovanni Resta_, May 30 2016
|