[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!)
Revision History for A348507 (Underlined text is an addition; strikethrough text is a deletion.)

Showing entries 1-10 | older changes
A348507 a(n) = A003959(n) - n, where A003959 is multiplicative with a(p^e) = (p+1)^e.
(history; published version)
#44 by Michael De Vlieger at Thu Nov 11 19:30:45 EST 2021
STATUS

reviewed

approved

#43 by Michel Marcus at Thu Nov 11 02:52:48 EST 2021
STATUS

proposed

reviewed

#42 by Antti Karttunen at Wed Nov 10 17:16:28 EST 2021
STATUS

editing

proposed

#41 by Antti Karttunen at Wed Nov 10 17:15:39 EST 2021
PROG

A348507(n) = { my(s=0, m=1, spf); while(n>1, spf = A020639(n); n /= spf; s += m*n; m *= (1+spf)); (s); }; \\ (Compare this with my similar programs programgiven forin A003415 and forin A322582) - Antti Karttunen, Nov 06 2021

#40 by Antti Karttunen at Wed Nov 10 17:14:43 EST 2021
CROSSREFS

Cf. also A168066A168065 (the arithmetic mean of this and A322582).), A168066.

#39 by Antti Karttunen at Wed Nov 10 17:12:56 EST 2021
COMMENTS

a(p*(n/p)) - (n/p) = (p+1)*a(n/p) holds for all prime divisors p of n, which can be seen by expanding the left hand side as (A003959(p*(n/p)) - (p*(n/p))) - (n/p) = (p+1)*A003959(n/p)-((p+1)*(n/p)) = (p+1)*(A003959(n/p)-(n/p)) = (p+1)*a(n/p). This implies that a(n) >= A003415(n) for all n. (See also comments in A348970). - Antti Karttunen, Nov 06 2021

FORMULA

For all n>= >= 1, A322582(n) <= A003415(n) <= a(n).

For n > 1, a(n) = a(A032742(n))*(1+A020639(n)) + A032742(n). [See the comments above, and compare this with Reinhard Zumkeller's May 09 2011 recursive formula for A003415 !]]

#38 by Antti Karttunen at Wed Nov 10 16:56:40 EST 2021
CROSSREFS

Cf. A348971 (Möbius transform), ) and A349139, A349140, A349141, A349142, A349143 (other Dirichlet convolutions).

#37 by Antti Karttunen at Wed Nov 10 16:54:01 EST 2021
CROSSREFS

Cf. A001065, A003415, A003959, A020639, A032742, A348029, A348508, A348929 [= gcd(n,a(n))], A348950, A348970, A348971 (Möbius transform)..

Cf. A348971 (Möbius transform), A349139, A349140, A349141, A349142, A349143 (other Dirichlet convolutions).

#36 by Antti Karttunen at Wed Nov 10 16:50:44 EST 2021
CROSSREFS

Cf. A001065, A003415, A003959, A020639, A032742, A348029, A348508, A348929 [= gcd(n,a(n))], A348950, A348970, A348971 (Möbius transform).

#35 by Antti Karttunen at Sun Nov 07 06:42:31 EST 2021
CROSSREFS

Cf. A001065, A003415, A003959, A020639, A032742, A348029, A348508, A348929 [= gcd(n,a(n))], A348970, A348971 (Möbius transform).

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 06:09 EDT 2024. Contains 375510 sequences. (Running on oeis4.)