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

A152454 Irregular triangle in which row n lists the numbers whose proper divisors sum to n.

Original entry on oeis.org

4, 9, 6, 25, 8, 10, 49, 15, 14, 21, 121, 27, 35, 22, 169, 16, 33, 12, 26, 39, 55, 289, 65, 77, 34, 361, 18, 51, 91, 20, 38, 57, 85, 529, 95, 119, 143, 46, 69, 133, 28, 115, 187, 841, 32, 125, 161, 209, 221, 58, 961, 45, 87, 247, 62, 93, 145, 253, 24, 155, 203, 299, 323, 1369
Offset: 2

Views

Author

T. D. Noe, Dec 05 2008

Keywords

Comments

In an aliquot sequence, all numbers in row n can be predecessors of n. This sequence is a permutation of the composite numbers; number k appears in row A001065(k). We start with n=2 because every prime would be in row 1. Note that row 2 is empty -- as are all the rows listed in A005114. Row n contains A048138(n) numbers. When n is prime, the largest number in row n+1 is n^2. When n>7 is odd, the largest number in row n is less than ((n-1)/2)^2 and (if a strong form of the Goldbach conjecture is true) has the form pq, with primes p
In row n, the first term is A070015(n), and the last term is A135244(n). - Michel Marcus, Nov 11 2014
The first row with several terms is row(6), where the difference between extreme terms is 25-6=19. The next row with a smaller difference is row(13) with a difference 35-27=8. And the next one is row(454) with a difference 602-596=6. Is there a next row with a smaller difference? - Michel Marcus, Nov 11 2014

Examples

			Irregular triangle starts:
  ; (empty row at n=2)
  4;
  9;
  ; (empty row at n=5)
  6, 25;
  8;
  10, 49;
  15;
  14;
  21;
  121;
  27, 35;
  22, 169;
  16, 33;
  12, 26;
  39, 55;
  289;
  ...
		

Crossrefs

Programs

  • Maple
    N:= 100: # for rows 2 to N, flattened
    for s from 2 to N do B[s]:= NULL od:
    for k from 1 to N^2 do
      if not isprime(k) then
        s:= numtheory:-sigma(k)-k;
        if s <= N then
           B[s]:= B[s],k;
        fi
      fi
    od:
    seq(B[s],s=2..N); # Robert Israel, Nov 11 2014
  • Mathematica
    nn=100; s=Table[{},{nn}]; Do[k=DivisorSigma[1,n]-n; If[1
    				
  • PARI
    row(n) = select(x->(sigma(x)-x)==n, [1..n^2]); \\ Michel Marcus, Feb 25 2025

A135244 Largest m such that the sum of the aliquot parts of m (A001065) equals n, or 0 if no such number exists.

Original entry on oeis.org

0, 4, 9, 0, 25, 8, 49, 15, 14, 21, 121, 35, 169, 33, 26, 55, 289, 77, 361, 91, 38, 85, 529, 143, 46, 133, 28, 187, 841, 221, 961, 247, 62, 253, 24, 323, 1369, 217, 81, 391, 1681, 437, 1849, 403, 86, 493, 2209, 551, 94, 589, 0, 667, 2809, 713, 106, 703, 68, 697, 3481
Offset: 2

Author

Ophir Spector (ospectoro(AT)yahoo.com), Nov 25 2007

Keywords

Comments

Previous name: Aliquot predecessors with the largest values.
Find each node's predecessors in aliquot sequences and choose the largest predecessor.
Climb the aliquot trees on shortest paths (see A135245 = Climb the aliquot trees on thickest branches).
The sequence starts at offset 2, since all primes satisfy sigma(n)-n = 1. - Michel Marcus, Nov 11 2014

Examples

			a(25) = 143 since 25 has 3 predecessors (95,119,143), 143 being the largest.
a(5) = 0 since it has no predecessors (see Untouchables - A005114).
		

Programs

  • Mathematica
    seq[max_] := Module[{s = Table[0, {n, 1, max}], i}, Do[If[(i = DivisorSigma[1, n] - n) <= max, s[[i]] = Max[s[[i]], n]], {n, 2, (max - 1)^2}]; Rest @ s]; seq[50]
  • PARI
    lista(nn) = {for (n=2, nn, k = (n-1)^2; while(k && (sigma(k)-k != n), k--); print1(k, ", "););} \\ Michel Marcus, Nov 11 2014

Extensions

a(1)=0 removed and offset set to 2 by Michel Marcus, Nov 11 2014
New name from Michel Marcus, Oct 31 2023
Showing 1-2 of 2 results.