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!)
A048166 Numbers k that are divisible by the number of unitary divisors of k (A034444). 6
1, 2, 4, 8, 12, 16, 20, 24, 28, 32, 36, 40, 44, 48, 52, 56, 64, 68, 72, 76, 80, 88, 92, 96, 100, 104, 108, 112, 116, 120, 124, 128, 136, 144, 148, 152, 160, 164, 168, 172, 176, 184, 188, 192, 196, 200, 208, 212, 216, 224, 232, 236, 240, 244, 248, 256, 264, 268 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

LINKS

Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms n = 1..5000 from G. C. Greubel)

FORMULA

Binomial transform of [1, 1, 1, 1, -3, 5, -7, 9, -11, 13, ...]. Binomial transform of this sequence is A048481. - Gary W. Adamson, Oct 23 2007

EXAMPLE

a(81) = 392 = 2^3*7^2 has 4 unitary divisors, {1, 392, 8, 49}, and 4 divides 392.

MATHEMATICA

okQ[n_] := Divisible[n, DivisorSum[n, Boole[CoprimeQ[#, n/#]]&]]; Select[ Range[300], okQ] (* Jean-Fran├žois Alcover, Dec 05 2015 *)

Select[Range[270], Divisible[#, 2^PrimeNu[#]] &] (* Amiram Eldar, Jul 16 2019 *)

PROG

(PARI) isok(n) = !(n % sumdiv(n, d, gcd(d, n/d)==1)); \\ Michel Marcus, Feb 25 2014

(PARI) isok(n) = !(n % 2^omega(n)); \\ Amiram Eldar, Jul 16 2019

CROSSREFS

Cf. A033950, A034444, A000005, A048167, A048481.

Sequence in context: A325326 A215459 A019442 * A335738 A010066 A180490

Adjacent sequences:  A048163 A048164 A048165 * A048167 A048168 A048169

KEYWORD

nonn

AUTHOR

Labos Elemer

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 8 19:29 EDT 2020. Contains 336298 sequences. (Running on oeis4.)