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!)
A286652 Unitary practical numbers: numbers n such that every 1 <= k <= usigma(n) is a sum of distinct unitary divisors of n. 6
1, 2, 6, 30, 42, 66, 78, 210, 330, 390, 462, 510, 546, 570, 690, 714, 798, 858, 870, 930, 966, 1050, 1110, 1122, 1218, 1230, 1254, 1290, 1302, 1326, 1410, 1470, 1482, 1518, 1554, 1590, 1650, 1722, 1770, 1794, 1806, 1830, 1914, 1950, 1974, 2010, 2046, 2130 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

The unitary version of A005153. The squarefree terms of both sequences are the same, A265501. The nonsquarefree terms of this sequence are in A287173.

LINKS

Amiram Eldar, Table of n, a(n) for n = 1..1000

MATHEMATICA

usigma[n_] :=  Block[{d = Divisors[n]}, Plus @@ Select[d, GCD[#, n/#] == 1 &]]; uPracticalQ[n_] :=  Module[{f, p, e, prod = 1, ok = True}, If[n < 1 || (n > 1 && OddQ[n]), False, If[n == 1, True, f = FactorInteger[n]; {p, e} = Transpose[f]; r = Sort[p^e]; Do[If[r[[i]] > 1 + usigma[prod], ok = False; Break[]]; prod = prod*r[[i]], {i, Length[p]}]; ok]]]; Select[ Range[100000], uPracticalQ]

CROSSREFS

Cf. A005153, A034448, A265501, A287173.

Sequence in context: A346106 A330326 A006954 * A265501 A090801 A166062

Adjacent sequences:  A286649 A286650 A286651 * A286653 A286654 A286655

KEYWORD

nonn

AUTHOR

Amiram Eldar, May 27 2017

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 August 3 12:51 EDT 2021. Contains 346438 sequences. (Running on oeis4.)