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!)
A349330 a(n) = Sum_{d|n} d^c(d), where c is the characteristic function of squares (A010052). 0
1, 2, 2, 6, 2, 4, 2, 7, 11, 4, 2, 9, 2, 4, 4, 23, 2, 14, 2, 9, 4, 4, 2, 11, 27, 4, 12, 9, 2, 8, 2, 24, 4, 4, 4, 55, 2, 4, 4, 11, 2, 8, 2, 9, 14, 4, 2, 28, 51, 30, 4, 9, 2, 16, 4, 11, 4, 4, 2, 15, 2, 4, 14, 88, 4, 8, 2, 9, 4, 8, 2, 58, 2, 4, 30, 9, 4, 8, 2, 28, 93, 4, 2, 15, 4, 4, 4 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
For each divisor d of n, add d if d is a square, otherwise add 1 [see example].
LINKS
FORMULA
a(p) = 2 iff p is prime. - Wesley Ivan Hurt, Nov 28 2021
EXAMPLE
The divisors of 12 are 1, 2, 3, 4, 6, and 12 with squares 1 and 4, so a(12) = 1 + 1 + 1 + 4 + 1 + 1 = 9 (respectively).
MATHEMATICA
a[n_] := DivisorSum[n, If[IntegerQ @ Sqrt[#], #, 1] &]; Array[a, 100] (* Amiram Eldar, Nov 15 2021 *)
PROG
(PARI) a(n) = sumdiv(n, d, if (issquare(d), d, 1)); \\ Michel Marcus, Nov 15 2021
(PARI) a(n) = {my(f = factor(n), cf = f, res); cf[, 2]\=2; res = numdiv(f)-prod(i = 1, #f~, cf[i, 2]+1); res+=prod(i = 1, #f~, ((f[i, 1]^(2*(cf[i, 2]+1))-1)/(f[i, 1]^2-1))); res } \\ David A. Corneth, Nov 16 2021
CROSSREFS
Cf. A010052.
Sequence in context: A247765 A129750 A278234 * A367203 A068976 A265392
KEYWORD
nonn
AUTHOR
Wesley Ivan Hurt, Nov 15 2021
STATUS
approved

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 March 29 06:57 EDT 2024. Contains 371265 sequences. (Running on oeis4.)