cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 24 results. Next

A017666 Denominator of sum of reciprocals of divisors of n.

Original entry on oeis.org

1, 2, 3, 4, 5, 1, 7, 8, 9, 5, 11, 3, 13, 7, 5, 16, 17, 6, 19, 10, 21, 11, 23, 2, 25, 13, 27, 1, 29, 5, 31, 32, 11, 17, 35, 36, 37, 19, 39, 4, 41, 7, 43, 11, 15, 23, 47, 12, 49, 50, 17, 26, 53, 9, 55, 7, 57, 29, 59, 5, 61, 31, 63, 64, 65, 11, 67, 34, 23, 35, 71, 24, 73, 37, 75, 19
Offset: 1

Views

Author

Keywords

Comments

Sum_{ d divides n } 1/d^k is equal to sigma_k(n)/n^k. So sequences A017665-A017712 also give the numerators and denominators of sigma_k(n)/n^k for k = 1..24. The power sums sigma_k(n) are in sequences A000203 (k=1), A001157-A001160 (k=2,3,4,5), A013954-A013972 for k = 6,7,...,24. - Ahmed Fares (ahmedfares(AT)my-deja.com), Apr 05 2001
Denominators of coefficients in expansion of Sum_{n >= 1} x^n/(n*(1-x^n)) = Sum_{n >= 1} log(1/(1-x^n)).
Also n/gcd(n, sigma(n)) = n/A009194(n); also n/lcm(all common divisors of n and sigma(n)). Equals 1 if 6,28,120,496,672,8128,..., i.e., if n is from A007691. - Labos Elemer, Aug 14 2002
a(A007691(n)) = 1. - Reinhard Zumkeller, Apr 06 2012
Denominator of sigma(n)/n = A000203(n)/n. a(n) = 1 for numbers n in A007691 (multiply-perfect numbers), a(n) = 2 for numbers n in A159907 (numbers n with half-integral abundancy index), a(n) = 3 for numbers n in A245775, a(n) = n for numbers n in A014567 (numbers n such that n and sigma(n) are relatively prime). See A162657 (n) - the smallest number k such that a(k) = n. - Jaroslav Krizek, Sep 23 2014
For all n, a(n) <= n, and thus records are obtained for terms of A014567. - Michel Marcus, Sep 25 2014
Conjecture: If a(n) is in A005153, then n is in A005153. In particular, if n has dyadic rational abundancy index, i.e., a(n) is in A000079 (such as A007691 and A159907), then n is in A005153. Since every term of A005153 greater than 1 is even, any odd n such that a(n) in A005153 must be in A007691. It is natural to ask if there exists a generalization of the indicator function for A005153, call it m(n), such that m(n) = 1 for n in A005153, 0 < m(n) < 1 otherwise, and m(a(n)) <= m(n) for all n. See also A050972. - Jaycob Coleman, Sep 27 2014

Examples

			1, 3/2, 4/3, 7/4, 6/5, 2, 8/7, 15/8, 13/9, 9/5, 12/11, 7/3, 14/13, 12/7, 8/5, 31/16, ...
		

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 162, #16, (6), 4th formula.

Crossrefs

Programs

  • Haskell
    import Data.Ratio ((%), denominator)
    a017666 = denominator . sum . map (1 %) . a027750_row
    -- Reinhard Zumkeller, Apr 06 2012
    
  • Magma
    [Denominator(DivisorSigma(1,n)/n): n in [1..50]]; // G. C. Greubel, Nov 08 2018
    
  • Maple
    with(numtheory): seq(denom(sigma(n)/n), n=1..76) ; # Zerinvary Lajos, Jun 04 2008
  • Mathematica
    Table[Denominator[DivisorSigma[-1, n]], {n, 100}] (* Vladimir Joseph Stephan Orlovsky, Jul 21 2011 *)
    Table[Denominator[DivisorSigma[1, n]/n], {n, 1, 50}] (* G. C. Greubel, Nov 08 2018 *)
  • PARI
    a(n) = denominator(sigma(n)/n); \\ Michel Marcus, Sep 23 2014
    
  • Python
    from math import gcd
    from sympy import divisor_sigma
    def A017666(n): return n//gcd(divisor_sigma(n),n) # Chai Wah Wu, Mar 21 2023

Extensions

More terms from Labos Elemer, Aug 14 2002

A054024 Sum of the divisors of n reduced modulo n.

Original entry on oeis.org

0, 1, 1, 3, 1, 0, 1, 7, 4, 8, 1, 4, 1, 10, 9, 15, 1, 3, 1, 2, 11, 14, 1, 12, 6, 16, 13, 0, 1, 12, 1, 31, 15, 20, 13, 19, 1, 22, 17, 10, 1, 12, 1, 40, 33, 26, 1, 28, 8, 43, 21, 46, 1, 12, 17, 8, 23, 32, 1, 48, 1, 34, 41, 63, 19, 12, 1, 58, 27, 4, 1, 51, 1, 40, 49, 64, 19, 12, 1, 26, 40
Offset: 1

Views

Author

Asher Auel, Jan 19 2000

Keywords

Comments

If a(n) = 0, then n is a multiply-perfect number (A007691). - Alonso del Arte, Mar 30 2014

Examples

			a(12) = 4 because sigma(12) = 28 and 28 == 4 (mod 12).
a(13) = 1 because 13 is prime.
a(14) = 10 because sigma(14) = 24 and 24 == 10 (mod 14).
		

Crossrefs

Cf. A000203 (sigma), A005114 (untouchable numbers), A007691 (positions of 0's), A045768, A045769, A088834, A045770, A076496, A159907.

Programs

Formula

a(n) = sigma(n) mod n.
a(p) = 1 for p prime.

A347391 Numbers k such that sigma(k) is either their sibling in Doudna tree (A005940), or one of the sibling's descendants.

Original entry on oeis.org

3, 4, 5, 15, 20, 189, 945, 2125, 6375, 9261, 46305, 401625, 19679625
Offset: 1

Views

Author

Antti Karttunen, Aug 30 2021

Keywords

Comments

Numbers k > 1 such that nearest common ancestor of k and sigma(k) in Doudna tree is the parent of k, and sigma(k) is not a descendant of k.
Any hypothetical odd term x in A005820 (triperfect numbers) would also be a member of this sequence. This is illustrated in the following diagram which shows how the neighborhood of such x would look like in the Doudna tree (A005940). If m (the parent of x, x = A003961(m), m = A064989(x)) is even, then x is a multiple of 3, while if m is odd, then 3 does not divide x. Because the abundancy index decreases when traversing leftwards in the Doudna tree, m must be a term of A068403. Both x and m would also need to be squares, by necessity.
.
<--A003961-- m ---(*2)--->
.............../ \...............
/ \
/ \
x 2m
/ \ / \
etc.../ \.....2x sigma(x) = 3x..../ \.....4m
/ \ / \ / \
etc. etc. etc. \ / etc.
\ /
6x 9x = sigma(2x)
/ \ / \
etc. \ etc. etc.
\
12x = sigma(3x) if m odd.
.
From the diagram we also see that 2x would then need to be a term of A347392 (as well as that of A159907 and also in A074388, thus sqrt(x) should be a term of A097023), and furthermore, if x is not a multiple of 3 (i.e., when m is odd), then sigma(3*x) = 4*sigma(x) = 4*(3*x), thus 3*x = sigma(x) would be a term of A336702 (particularly, in A027687) and x would be a term of A323653.
Moreover, any odd square x in this sequence (for which sigma(x) would also be odd), would have an abundancy index of at least three (sigma(x)/x >= 3). See comments in A347383.
Note how 401625 = 6375 * 63 = 945 * 425, 46305 = 945 * 49, 9261 = 189 * 49, 6375 = 2125 * 3, 945 = 189 * 5 = 15 * 63 and 9261*2125 = 19679625. It seems that when the multiplicands are coprime, then they are both terms of this sequence, e.g. 2125 and 3, 189 and 5, 2125 and 9261.
From Antti Karttunen, Jul 10 2024: (Start)
Regarding the observation above, for two coprime odd numbers x, y, if both are included here because sigma(x) = 2^a * A064989(x) and sigma(y) = 2^b * A064989(y), then also their product x*y is included because in that case sigma(x*y) = 2^(a+b) * A064989(x*y).
Also, for two coprime odd numbers x, y, if both are included here because sigma(x) = A065119(i) * x and sigma(y) = A065119(j) * y, then also their product x*y is included because sigma(x*y) = A065119(k) * x*y, where A065119(k) = A065119(i)*A065119(j). The existence of such numbers (that would include odd triperfect and odd 6-perfect numbers, see A046061) is so far hypothetical, none is known.
It is not possible that the odd x is in this sequence if sigma(x) = k*A003961^e(x) and e = A061395(k)-2 >= 1.
Note that all odd terms < 2^33 here are some of the exponentially odd divisors of 19679625 (see A374199, also A374463 and A374464).
(End)
Question: from a(6) = 189 onward, are the rest of terms all in A347390?
Conjecture: sequence is finite.
If it exists, a(14) > 2^33.

Examples

			Sigma(3) = 4 is located as the sibling of 3 in the Doudna-tree (see the illustration in A005940), thus 3 is included in this sequence.
Sigma(4) = 7 is located as a grandchild of 3 (which is the sibling of 4) in the Doudna-tree, thus 4 is included in this sequence.
Sigma(5) = 6 is located as the sibling of 5 in the Doudna-tree, thus 5 is included in this sequence.
189 (= 3^3 * 7) is a term, as sigma(189) = 320, and 320 occurs as a descendant of 80 (which is the right sibling of 189) in the Doudna tree, as illustrated below:
.
             40
            /  \
   A003961 /    \ *2
          /      \
        189       80
        / \      / \
     etc   etc etc  160
                   / \
                 etc  320
                     / \
                   etc. etc.
.
945 (= 3^3 * 5 * 7) is a term, as sigma(945) = 1920, and 1920 occurs as a descendant of 240, which is the right sibling of 945 in the Doudna tree, as illustrated below:
            120
            /  \
   A003961 /    \ *2
          /      \
        945       240
        / \      / \
     etc   etc  etc  480
                   / \
                 etc  960
                     / \
                   etc. 1920
                        / \
                     etc. etc.
		

Crossrefs

Programs

  • PARI
    isA347391(n) = (1==A347381(n));
    
  • PARI
    A064989(n) = {my(f); f = factor(n); if((n>1 && f[1,1]==2), f[1,2] = 0); for (i=1, #f~, f[i,1] = precprime(f[i,1]-1)); factorback(f)};
    A252463(n) = if(!(n%2),n/2,A064989(n));
    isA347391(n) = if(1==n,0,my(m=A252463(n), s=sigma(n)); while(s>m, if(s==n, return(0)); s = A252463(s)); (s==m));

A229110 Sum of non-divisors of n reduced modulo n.

Original entry on oeis.org

0, 0, 2, 3, 4, 3, 6, 5, 5, 7, 10, 2, 12, 11, 6, 9, 16, 6, 18, 8, 10, 19, 22, 0, 19, 23, 14, 14, 28, 3, 30, 17, 18, 31, 22, 35, 36, 35, 22, 10, 40, 9, 42, 26, 12, 43, 46, 44, 41, 32, 30, 32, 52, 15, 38, 20, 34, 55, 58, 42, 60, 59, 22, 33, 46, 21, 66, 44, 42
Offset: 1

Views

Author

Jaroslav Krizek, Sep 22 2013

Keywords

Comments

Numbers n such that a(n)=0 are: 1, 2, 24, 4320, 4680, ... (see A159907, conjecture by Jaroslav Krizek and further comments). - Michel Marcus, Sep 23 2013
Numbers n such that a(n)=n/2 are: 6, 28, 120, 496, 672, ... = A007691 \ {1}. - Michel Marcus, Sep 25 2013

Crossrefs

Programs

  • Haskell
    a229110 n = mod (a024816 n) n
    
  • PARI
    a(n) = lift(sum(i=1, n, if (n % i, Mod(i, n), 0))); \\ Michel Marcus, Sep 23 2013
    
  • PARI
    a(n)=(n*(n+1)/2-sigma(n))%n \\ Charles R Greathouse IV, Sep 23 2013

Formula

a(n) = A024816(n) mod n.

A141645 Numbers n such that sigma(n)/n = 9/2.

Original entry on oeis.org

8910720, 17428320, 8583644160, 57629644800, 206166804480, 1416963251404800, 15338300494970880, 6275163455171297280, 200286975596707184640, 215594611071909888000, 5997579964837140234240, 39887491844324122951680
Offset: 1

Views

Author

Yasutoshi Kohmoto, May 10 2008

Keywords

Comments

200286975596707184640 belongs to this sequence. - Gerard P. Michon, May 10 2009
81703797123392614369698250752 is in this sequence. - Gerard P. Michon, May 11 2009

Crossrefs

Cf. A159907 = (half-integer abundancy). - Gerard P. Michon, May 10 2009

Programs

Extensions

Definition rewritten by M. F. Hasler, May 10 2008
Is it certain that there are no other terms below the highest value shown? - N. J. A. Sloane, Sep 07 2008
a(1) corrected and a(3) confirmed by Ray Chandler, Sep 18 2008
Missing term a(2)=17428320 communicated by Walter Nissen, Apr 17 2009. There are no further terms through 2^34.
Missing term a(4) inserted and a(6) and a(8) added by Grzegorz Lach (137.036(AT)gmail.com), Apr 18 2009. a(4) was also sent by Avinoam Kalma (a.kalma(AT)gmail.com), Apr 20 2009.
a(7) = 15338300494970880 added by Gerard P. Michon, May 10 2009.
a(10)=215594611071909888000 [and above] from Michel Marcus, added by Gerard P. Michon, Jun 04 2009

A160320 Numbers n for which sigma(n)/n = k+1/3 with integer k.

Original entry on oeis.org

3, 12, 234, 1080, 6048, 6552, 435708, 4713984, 18506880, 36197280, 275890944, 299980800, 4138364160, 13286744064, 17827568640, 30600708096, 459010621440, 2198278051200, 2567400675840, 3321402084000, 4959751305600, 295185198672000, 1217745604232064
Offset: 1

Views

Author

J.C. Klein (hklein(AT)planet.nl), May 08 2009

Keywords

Examples

			a(1)=3 since sigma(3)/3=(1+3)/3=1+1/3.
a(2)=12 since sigma(12)/12=(1+2+3+4+6+12)/12=28/12=7/3=2+1/3.
		

Crossrefs

Cf. A160321 (sigma(n)/n=k+2/3), A159907 (sigma(n)/n=k+1/2).

Programs

Extensions

a(11)-a(16) from Donovan Johnson, May 26 2009
a(17)-a(23) from Michel Marcus, Sep 21 2012

A160321 Numbers n for which sigma(n)/n=k+2/3 with integer k.

Original entry on oeis.org

84, 270, 1488, 1638, 24384, 35640, 199584, 2142720, 12999168, 100651008, 208565280, 240589440, 470564640, 3899750400, 6039429120, 25769607168, 36639203328, 53798734080, 231758392320, 314039721600, 412316073984, 5566503720960, 5967138078720, 302512616524800
Offset: 1

Views

Author

J.C.Klein (hklein(AT)planet.nl), May 08 2009

Keywords

Examples

			84 is in the sequence because sigma(84)/84=224/84=2+2/3
270 is in the sequence because sigma(270)/270=720/270=2+2/3
		

Crossrefs

Cf. A160320 (sigma(n)/n=k+1/3), A159907 (sigma(n)/n=k+1/2).

Extensions

a(10)-a(18) from Donovan Johnson, May 26 2009
a(19)-a(24) from Michel Marcus, Sep 21 2012

A160678 Numbers n whose abundancy is equal to 13/2; sigma(n)/n = 13/2.

Original entry on oeis.org

170974031122008628879954060917200710847692800, 1893010442758976546037991125738431754692198400, 54361481238923605327597493185154939181072384000
Offset: 1

Views

Author

Gerard P. Michon, Jun 06 2009

Keywords

Comments

This sequence includes many terms but it is conjectured to be finite.

Examples

			a(1) = 2^23 3^9 5^2 7^5 11^5 13^2 17 19^3 31 37 43 61^2 97 181 241.
As the "sum of divisors" function (sigma) is a multiplicative function, sigma(a(1)) is the product of the values of sigma at the above prime powers, respectively given as follows, in factorized form:
sigma(a(1)) = (3^2 5 7 13 17 241) (2^2 11^2 61) (31) (2^3 3 19 43) (2^2 3^2 7 19 37) (3 61) (2 3^2) (2^3 5 181) (2^5) (2 19) (2^2 11) (3 13 97) (2 7 13) (2 7^2) (2 11^2).
a(1) belongs to the sequence because the latter product boils down to 13/2 times the former.
		

Crossrefs

Cf. A000203 (sigma function, sum of divisors), A141643 (abundancy = 5/2), A055153 (abundancy = 7/2), A141645 (abundancy = 9/2), A159271 (abundancy = 11/2), A159907 (half-integral abundancy, "hemiperfect numbers"), A088912 (least numbers of given half-integer abundancy). A007691 (multiperfect numbers, abundancy is an integer), A000396 (perfect numbers, abundancy = 2), A005101 (abundant numbers, abundancy is greater than 2), A005100 (deficient numbers, abundancy is less than 2).

Programs

A227302 Numbers m such that m divides sigma(2*m).

Original entry on oeis.org

1, 3, 12, 14, 60, 248, 336, 2160, 2340, 4064, 13104, 15120, 16380, 261888, 1089270, 4455360, 8714160, 10213632, 11784960, 16775168, 22766400, 45981824, 71495424, 98532480, 229909120, 689727360, 738152448, 4291822080, 4294934528, 5100118016, 7091219520
Offset: 1

Views

Author

Alex Ratushnyak, Jul 05 2013

Keywords

Comments

If m belongs to the sequence, then sigma(2*m)/m is an integer, so sigma(2*m)/(2*m) is either an integer or half of an integer, so 2*m is either perfect, multiperfect or hemiperfect. - Michel Marcus, Jul 09 2013

Crossrefs

Cf. A141643, A055153, A141645, A159271, A160678. (hemiperfect numbers)

Programs

A242484 Numbers n such that antisigma(n) mod n = 0, where antisigma(n) = A024816(n) = sum of numbers less than n which do not divide n.

Original entry on oeis.org

1, 2, 24, 4320, 4680, 26208, 8910720, 17428320, 20427264, 91963648, 197064960, 8583644160, 10200236032, 21857648640, 57575890944, 57629644800, 206166804480, 17116004505600, 1416963251404800, 15338300494970880
Offset: 1

Views

Author

Jaroslav Krizek, May 16 2014

Keywords

Comments

Numbers n such that antisigma(n) mod n = A229110(n) = 0.
If there are any odd multiply-perfect numbers, they are members of this sequence.
If there is no odd multiply-perfect number, then a(n) = A159907(n-1) for n >= 2.

Examples

			24 is in sequence because antisigma(24) mod 24 = 240 mod 24 = 0.
		

Crossrefs

Programs

  • Magma
    [n: n in [1..1000000] | 0 eq ((n*(n+1))div 2 - SumOfDivisors(n)) mod n];
Showing 1-10 of 24 results. Next