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-3 of 3 results.

A337381 Numbers k for which A003973(k) >= 2*sigma(k).

Original entry on oeis.org

6, 8, 9, 12, 14, 15, 16, 18, 20, 21, 24, 27, 28, 30, 32, 35, 36, 40, 42, 44, 45, 48, 49, 50, 52, 54, 56, 60, 63, 64, 66, 68, 70, 72, 75, 76, 78, 80, 81, 84, 88, 90, 92, 96, 98, 99, 100, 102, 104, 105, 108, 110, 112, 114, 117, 120, 124, 125, 126, 128, 130, 132, 135, 136, 138, 140, 144, 147, 148, 150, 152, 153, 154
Offset: 1

Views

Author

Antti Karttunen, Aug 27 2020

Keywords

Comments

Note that A003973(n) >= sigma(n) for all n. See A336852.
Like the abundancy index (ratio A000203(n)/n), and ratio A003961(n)/n, the ratio A003973(n)/sigma(n) is also multiplicative and > 1 for all n > 1. Thus if the number has a proper divisor that is in this sequence, then the number itself is also. See A337543 for those terms included here, but which have no proper divisor in this sequence. - Antti Karttunen, Aug 31 2020
All terms are in A246282 because A341528(n) < A341529(n) for all n > 1. - Antti Karttunen, Feb 22 2021

Crossrefs

Cf. A337382 (complement), A337383 (characteristic function).
Subsequences: A337378, A337384, A337386, A337543 (primitive terms).
Subsequence of A246282.

Programs

  • PARI
    A003973(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); sigma(factorback(f)); };
    isA337381(n) = (A003973(n)>=2*sigma(n));

A337543 Numbers that are in A337381, but that have no proper divisor in A337381.

Original entry on oeis.org

6, 8, 9, 14, 15, 20, 21, 35, 44, 49, 50, 52, 68, 76, 92, 110, 124, 125, 130, 148, 172, 188, 190, 212, 230, 244, 275, 286, 292, 310, 325, 338, 356, 363, 425, 429, 452, 470, 475, 494, 507, 575, 598, 627, 663, 715, 741, 759, 775, 806, 845, 847, 874, 897, 925, 969, 1001, 1023, 1058, 1075, 1083, 1131, 1173, 1175, 1183, 1209
Offset: 1

Views

Author

Antti Karttunen, Aug 31 2020

Keywords

Comments

Numbers k that are in A337381 [i.e., for which A003973(k) >= 2*A000203(k)], but none of whose proper divisors are in A337381.

Crossrefs

Programs

  • PARI
    A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
    A337541(n) = sumdiv(n,d,sigma(A003961(d))>=2*sigma(d));
    isA337543(n) = (1==A337541(n));

Formula

{k: 1==A337541(k)}.

A337541 Number of divisors d of n for which sigma(A003961(d)) >= 2*sigma(d), where sigma is the sum of divisors, and A003961(x) shifts the prime factorization of x one step towards larger primes.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 2, 0, 1, 1, 2, 0, 3, 0, 1, 1, 0, 0, 4, 0, 0, 2, 2, 0, 3, 0, 3, 0, 0, 1, 5, 0, 0, 0, 3, 0, 4, 0, 1, 3, 0, 0, 6, 1, 1, 0, 1, 0, 5, 0, 4, 0, 0, 0, 6, 0, 0, 3, 4, 0, 2, 0, 1, 0, 3, 0, 8, 0, 0, 2, 1, 0, 2, 0, 5, 3, 0, 0, 7, 0, 0, 0, 3, 0, 7, 0, 1, 0, 0, 0, 8, 0, 3, 2, 3, 0, 2, 0, 3, 4
Offset: 1

Views

Author

Antti Karttunen, Aug 31 2020

Keywords

Comments

Number of terms of A337381 that divide n.

Crossrefs

Inverse Möbius transform of A337383.
Cf. A000203, A003961, A003973, A337381, A337542, A337543 (positions of ones).
Cf. also A337345.

Programs

  • PARI
    A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
    A337541(n) = sumdiv(n,d,sigma(A003961(d))>=2*sigma(d));

Formula

a(n) = Sum_{d|n} A337383(d).
a(n) = A337542(n) + A337383(n).
Showing 1-3 of 3 results.