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.

A339938 Odd non-coreful abundant numbers: the odd terms of A308127.

Original entry on oeis.org

15015, 19635, 21945, 23205, 25935, 26565, 31395, 33495, 33915, 35805, 39585, 41055, 42315, 42735, 45885, 47355, 49665, 50505, 51765, 54285, 55965, 58695, 61215, 64155, 68145, 70455, 72345, 75075, 77385, 80535, 82005, 83265, 84315, 91245, 95865, 102795, 105105
Offset: 1

Views

Author

Amiram Eldar, Dec 23 2020

Keywords

Comments

First differs from A112643, A129485 and A249263 at n = 28.

Examples

			15015 is a term since it is odd and the sum of its non-coreful divisors is A308135(15015) = 17241 > 15015.
		

Crossrefs

Intersection of A005408 and A308127.
Cf. A308135.

Programs

  • Mathematica
    f[p_, e_] := (p^(e + 1) - 1)/(p - 1); fc[p_, e_] := f[p, e] - 1; s[1] = 0; s[n_] := Times @@ (f @@@ FactorInteger[n]) - Times @@ (fc @@@ FactorInteger[n]); Select[Range[1, 10^5, 2], s[#] > # &]

A308135 Sum of non-coreful divisors of n.

Original entry on oeis.org

0, 1, 1, 1, 1, 6, 1, 1, 1, 8, 1, 10, 1, 10, 9, 1, 1, 15, 1, 12, 11, 14, 1, 18, 1, 16, 1, 14, 1, 42, 1, 1, 15, 20, 13, 19, 1, 22, 17, 20, 1, 54, 1, 18, 18, 26, 1, 34, 1, 33, 21, 20, 1, 42, 17, 22, 23, 32, 1, 78, 1, 34, 20, 1, 19, 78, 1, 24, 27, 74, 1, 27, 1, 40
Offset: 1

Views

Author

Amiram Eldar and Paolo P. Lava, May 14 2019

Keywords

Comments

Non-coreful divisor d of a number k is a divisor such that rad(d) != rad(k), where rad(k) is the largest squarefree divisor of k (A007947).

Examples

			a(15) = 9. Prime factors of 15 are 3, 5 and its divisors are 1, 3, 5, 15. The non-coreful divisors are 1, 3, 5 and their sum is 9.
		

Crossrefs

Programs

  • Maple
    with(numtheory): P:=proc(k) local a,n; a:=mul(n,n=factorset(k));
    sigma(k)-a*sigma(k/a); end: seq(P(i),i=1..74);
  • Mathematica
    f[p_, e_] := (p^(e + 1) - 1)/(p - 1); fc[p_, e_] := f[p, e] - 1; a[1] = 0; a[n_] := Times @@ (f @@@ FactorInteger[n]) - Times @@ (fc @@@ FactorInteger[n]); Array[a, 100]

Formula

a(n) = A000203(n) - A057723(n).
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = A013661 - A065487 = 0.413642... . - Amiram Eldar, Dec 08 2023

A360525 Numbers k such that A360522(k) > 2*k.

Original entry on oeis.org

30, 42, 60, 66, 70, 78, 84, 90, 102, 114, 120, 126, 132, 138, 140, 150, 156, 168, 174, 180, 186, 204, 210, 222, 228, 246, 252, 258, 276, 282, 294, 300, 318, 330, 348, 354, 360, 366, 372, 390, 402, 420, 426, 438, 444, 462, 474, 492, 498, 510, 516, 534, 546, 564
Offset: 1

Views

Author

Amiram Eldar, Feb 10 2023

Keywords

Comments

First differs from A308127 at n = 15.
Analogous to abundant numbers (A005101) with A360522 instead of A000203.
Subsequence of A005101 because A360522(n) <= A000203(n) for all n.
The least odd term is a(1698) = A360526(1) = 15015.
The numbers of terms not exceeding 10^k, for k = 1, 2, ..., are 0, 8, 95, 1135, 10890, 110867, 1104596, 11048123, 110534517, 1105167384, 11051009278, ... . Apparently, the asymptotic density of this sequence exists and equals 0.1105...

Examples

			30 is a term since A360522(30) = 72 > 2*30.
		

Crossrefs

Subsequence of A005101.

Programs

  • Mathematica
    f[p_, e_] := p^e + e; s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; q[n_] := s[n] > 2*n; Select[Range[1000], q]
  • PARI
    is(n) = {my(f = factor(n)); prod(i = 1, #f~, f[i,1]^f[i,2] + f[i,2]) > 2*n;}
Showing 1-3 of 3 results.