login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A252424 Numbers n such that sum of odd divisors of n equals sum of squares of primes dividing n. 1

%I

%S 18,36,72,144,234,288,468,576,936,1152,1872,2304,3744,4608,7488,9216,

%T 14976,18432,29952,36864,59904,73728,119808,147456,239616,294912,

%U 479232,589824,958464,1179648,1916928,2359296,3833856,4718592,7667712,9437184,15335424,18874368

%N Numbers n such that sum of odd divisors of n equals sum of squares of primes dividing n.

%C Numbers n such that A000593(n) = A005063(n).

%C a(n)==0 mod 18, and the numbers 18*2^m, m=0,1,... are in the sequence because the odd divisors are {1, 3, 9}, the prime factors are {2, 3} => 2^2 + 3^2 = 1 + 3 + 9 = 13.

%C The numbers of the form 18*13*2^m are in the sequence because the odd divisors are {1, 3, 9, 13, 39, 117}, the prime factors are {2, 3, 13} => 2^2 + 3^2 + 13^2 = 1 + 3 + 9 + 13 + 39 + 117 = 182.

%H Robert G. Wilson v, <a href="/A252424/b252424.txt">Table of n, a(n) for n = 1..56</a>

%e 18 is in the sequence because the prime factors of 18 are {2, 3}, the odd divisors of 18 are {1, 3, 9} => 2^2 + 3^2 = 1 + 3 + 9 = 13.

%e Or 18 => A000593(18) = A005063(18)= 13.

%p with(numtheory):nn:=10^5:

%p for n from 2 to nn do:

%p x:=factorset(n):n0:=nops(x):

%p s0:=sum('x[i]^2','i'=1..n0):

%p y:=divisors(n):n1:=nops(y):

%p s :=0 :

%p for j from 1 to n1 do :

%p if irem (y[j],2)=1 then s:=s+y[j]:

%p else

%p fi:

%p od:

%p if s=s0

%p then

%p printf(`%d, `,n):

%p else

%p fi:

%p od:

%t a252424[n_Integer] := Module[{f, g},

%t f[x_] := Plus @@ Select[Divisors[x], OddQ[#] &];

%t g[x_] := Plus @@ (First@Transpose@FactorInteger[x]^2);

%t Rest@Select[Range[n], f[#] == g[#] &]]; a252424[10^6] (* _Michael De Vlieger_, Dec 17 2014 *)

%t Select[Range[19*10^6],Total[Select[Divisors[#],OddQ]]==Total[ FactorInteger[ #][[All,1]]^2]&] (* _Harvey P. Dale_, May 11 2020 *)

%o (PARI) isok(n) = my(f = factor(n)); sum(i=1, #f~, f[i,1]^2) == sumdiv(n, d, d*(d%2)); \\ _Michel Marcus_, Dec 17 2014

%Y Cf. A000593, A005063.

%K nonn

%O 1,1

%A _Michel Lagneau_, Dec 17 2014

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified July 11 22:31 EDT 2020. Contains 335652 sequences. (Running on oeis4.)