login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A051448 Sum of prime divisors of n (with multiplicity) is a square. 2
1, 4, 14, 20, 24, 27, 39, 46, 55, 66, 94, 98, 140, 152, 155, 158, 168, 171, 183, 186, 189, 200, 203, 225, 240, 255, 256, 270, 272, 288, 290, 291, 295, 299, 306, 323, 324, 334, 344, 348, 354, 363, 387, 446, 455, 506, 539, 546, 578, 579, 583, 615, 650, 656, 695 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

Numbers k such that A001414(k) is a perfect square. - Michel Lagneau, May 28 2012

LINKS

Harvey P. Dale, Table of n, a(n) for n = 1..1000

EXAMPLE

a(2) = 14 because 14 = 2*7 and 2 + 7 = 3^2.

MAPLE

A:= proc(n) local e, j; e := ifactors(n)[2]: add (e[j][1]*e[j][2], j=1..nops(e)) end:

for m from 1 to 1000 do: m2:=sqrt(A(m)):if m2=floor(m2) then printf(`%d, `, m):else fi:od: # Michel Lagneau, May 28 2012

MATHEMATICA

lst = {}; Do[ww = Transpose[FactorInteger[k]]; w = ww[[1]].ww[[2]]; If[IntegerQ[Sqrt[w]], AppendTo[lst, k]], {k, 1, 1000}]; lst (* Michel Lagneau, May 28 2012 *)

Select[Range[700], IntegerQ[Sqrt[Total[Flatten[Table[#[[1]], #[[2]]]&/@ FactorInteger[ #]]]]]&] (* Harvey P. Dale, Dec 12 2018 *)

CROSSREFS

Sequence in context: A059007 A035401 A139330 * A075319 A030470 A326004

Adjacent sequences:  A051445 A051446 A051447 * A051449 A051450 A051451

KEYWORD

easy,nonn

AUTHOR

Joe K. Crump (joecr(AT)carolina.rr.com)

EXTENSIONS

More terms from James A. Sellers, Sep 08 2000

STATUS

approved

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 November 12 09:29 EST 2019. Contains 329054 sequences. (Running on oeis4.)