The OEIS mourns the passing of Jim Simons and is grateful to the Simons Foundation for its support of research in many branches of science, including the OEIS.
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!)
A353000 Quotients obtained when sigma(k) divides antisigma(k) with k = A076617(n), sigma (A000203) is the sum of divisors function and antisigma (A024816) is the sum of the non-divisors of n less than n function. 0
0, 0, 4, 4, 4, 37, 25, 68, 49, 122, 115, 340, 544, 487, 959, 2167, 1926, 4837, 3847, 6757, 6452, 3620, 11353, 13934, 9371, 16353, 9211, 30949, 49702, 17330, 32575, 72544, 62348, 109769, 145892, 51270, 173914, 130687, 61665, 102887, 351770, 446927, 504949, 258079 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,3
COMMENTS
Note that the quotient obtained when sigma(k) divides k*(k+1)/2 with k = A076617(n) is a(n) + 1.
LINKS
FORMULA
a(n) = A024816(A076617(n)) / A000203(A076617(n)).
EXAMPLE
A076617(6) = 95; sigma(95) = 120 and antisigma(95) = 4440, hence a(6) = 4440 / 120 = 37.
MATHEMATICA
Select[Table[(k*(k + 1)/2)/DivisorSigma[1, k] - 1, {k, 1, 10^6}], IntegerQ] (* Amiram Eldar, Apr 14 2022 *)
PROG
(PARI) is(n) = n*(n+1)/2%sigma(n) == 0; \\ A076617
f(n) = n*(n+1)/(2*sigma(n)) - 1;
lista(nn) = apply(f, select(is, [1..nn])); \\ Michel Marcus, Apr 15 2022
CROSSREFS
Sequence in context: A059812 A077725 A217940 * A372020 A182065 A325487
KEYWORD
nonn
AUTHOR
Bernard Schott, Apr 14 2022
EXTENSIONS
More terms from Amiram Eldar, Apr 14 2022
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 June 3 13:21 EDT 2024. Contains 373060 sequences. (Running on oeis4.)