[go: up one dir, main page]

login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
Search: a066214 -id:a066214
Displaying 1-3 of 3 results found. page 1
     Sort: relevance | references | number | modified | created      Format: long | short | data
A066213 Numbers which are sums of squares of some subset of divisors. +10
4
1, 4, 9, 16, 20, 25, 30, 36, 49, 64, 80, 81, 90, 100, 120, 121, 126, 130, 144, 150, 169, 180, 195, 196, 210, 225, 252, 256, 264, 270, 272, 280, 289, 294, 300, 315, 320, 324, 330, 336, 350, 360, 361, 378, 390, 396, 400, 414, 420, 441, 450, 468, 480, 484, 500 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
If m is in the sequence then so is m*k^2 for k >= 1. - David A. Corneth, Jan 22 2024
LINKS
David A. Corneth, PARI program
EXAMPLE
20 is in the list since 20 = 2^2 + 4^2 and 2 and 4 are divisors of 20
MAPLE
isA066213 := proc(n)
local S, els;
S:=subsets(numtheory[divisors](n));
while not S[finished] do
els:=S[nextvalue]() ;
if add(d^2, d=els) = n then
return true ;
end if ;
end do;
false
end proc:
for n from 1 do
if isA066213(n) then
print(n) ;
end if;
end do: # R. J. Mathar, Oct 09 2023
MATHEMATICA
okQ[k_] := AnyTrue[Subsets[Select[Divisors[k]^2, # <= k&]], Total[#]==k&];
Reap[For[k = 1, k <= 5000, k++, If[okQ[k], Print[k]; Sow[k]]]][[2, 1]] (* Jean-François Alcover, May 27 2024 *)
PROG
(PARI) \\ See PARI link
CROSSREFS
KEYWORD
nonn
AUTHOR
Erich Friedman, Dec 17 2001
EXTENSIONS
Offset 1 from David A. Corneth, Jan 22 2024
STATUS
approved
A066215 Numbers which are sums of cubes of some subset of divisors. +10
4
1, 8, 27, 36, 64, 72, 125, 216, 252, 288, 343, 378, 512, 520, 576, 584, 729, 738, 756, 792, 828, 855, 954, 972, 1000, 1044, 1331, 1350, 1440, 1520, 1540, 1728, 1764, 1800, 1890, 1944, 1980, 2016, 2070, 2160, 2197, 2304, 2352, 2376, 2400, 2484, 2520, 2548 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
There are cubes that have not a single, trivial representation but more than one. These start with 27000 =+8^3+9^3+10^3+12^3+15^3+18^3+24^3 =+1^3+2^3+4^3+6^3+8^3+10^3+15^3+20^3+24^3 = +2^3+4^3+9^3+10^3+15^3+20^3+24^3 =+1^3+2^3+4^3+12^3+15^3+20^3+24^3 =+2^3+3^3+4^3+5^3+6^3+12^3+15^3+18^3+25^3 =+1^3+4^3+5^3+6^3+8^3+9^3+12^3+20^3+25^3 =+15^3+20^3+25^3 = +1^3+3^3+6^3+8^3+9^3+18^3+27^3 =+30^3 and 46656 =+1^3+2^3+3^3+6^3+8^3+9^3+12^3+16^3+18^3+24^3+27^3 =+4^3+24^3+32^3 =+36^3 and 74088 =+2^3+6^3+7^3+8^3+9^3+12^3+18^3+21^3+24^3+27^3+28^3 =+4^3+6^3+8^3+14^3+18^3+21^3+24^3+27^3+28^3 =+42^3. - R. J. Mathar, Jan 21 2024
If m is in the sequence then so is m*k^3 for k >= 1. - David A. Corneth, Jan 21 2024
LINKS
David A. Corneth, Table of n, a(n) for n = 1..10286 (first 649 terms from R. J. Mathar)
EXAMPLE
72 is in the list since 72 = 2^3 + 4^3 and 2 and 4 are divisors of 72
MATHEMATICA
okQ[k_] := AnyTrue[Subsets[Select[Divisors[k]^3, # <= k&]], Total[#]==k&];
Reap[For[k = 1, k <= 10000, k++, If[okQ[k], Print[k]; Sow[k]]]][[2, 1]] (* Jean-François Alcover, May 27 2024 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Erich Friedman, Dec 17 2001
STATUS
approved
A066216 Noncube numbers which are sums of cubes of some subset of divisors. +10
4
36, 72, 252, 288, 378, 520, 576, 584, 738, 756, 792, 828, 855, 954, 972, 1044, 1350, 1440, 1520, 1540, 1764, 1800, 1890, 1944, 1980, 2016, 2070, 2160, 2304, 2352, 2376, 2400, 2484, 2520, 2548, 2556, 2700, 2772, 2808, 2820, 2870, 2880, 3024, 3220, 3240 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
LINKS
Jean-François Alcover, Table of n, a(n) for n = 1..1600
EXAMPLE
72 is in the list since 72 = 2^3 + 4^3 and 2 and 4 are divisors of 72
MATHEMATICA
okQ[k_] := AnyTrue[Subsets[Select[Divisors[k]^3, # <= k&]], Total[#]==k&];
Reap[For[k = 1, k <= 100000, k++, If[!IntegerQ[k^(1/3)] && okQ[k], Print[k]; Sow[k]]]][[2, 1]] (* Jean-François Alcover, May 27 2024 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Erich Friedman, Dec 17 2001
STATUS
approved
page 1

Search completed in 0.007 seconds

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified August 29 21:13 EDT 2024. Contains 375518 sequences. (Running on oeis4.)