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

%I #20 Dec 15 2014 05:19:45

%S 4,9,6,25,8,10,49,15,14,21,121,27,35,22,169,16,33,12,26,39,55,289,65,

%T 77,34,361,18,51,91,20,38,57,85,529,95,119,143,46,69,133,28,115,187,

%U 841,32,125,161,209,221,58,961,45,87,247,62,93,145,253,24,155,203,299,323,1369

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

%C 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<q and p+q=n-1.

%C In row n, the first term is A070015(n), and the last term is A135244(n). - _Michel Marcus_, Nov 11 2014

%C 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

%H T. D. Noe, <a href="/A152454/b152454.txt">Rows n=2..1000 of triangle, flattened</a>

%H Michel Marcus, <a href="/A152454/a152454.txt">Rows n=2..1000 of triangle, not flattened</a>

%e Irregular triangle starts:

%e ; (empty row at n=2)

%e 4;

%e 9;

%e ; (empty row at n=5)

%e 6, 25;

%e 8;

%e 10, 49;

%e 15;

%e 14;

%e 21;

%e 121;

%e 27, 35;

%e 22, 169;

%e 16, 33;

%e 12, 26;

%e 39, 55;

%e 289;

%e ...

%p N:= 100: # for rows 2 to N, flattened

%p for s from 2 to N do B[s]:= NULL od:

%p for k from 1 to N^2 do

%p if not isprime(k) then

%p s:= numtheory:-sigma(k)-k;

%p if s <= N then

%p B[s]:= B[s],k;

%p fi

%p fi

%p od:

%p seq(B[s],s=2..N); # _Robert Israel_, Nov 11 2014

%t nn=100; s=Table[{},{nn}]; Do[k=DivisorSigma[1,n]-n; If[1<k<=nn, AppendTo[s[[k]],n]], {n,nn^2}]; Flatten[s]

%Y Cf. A001065, A005114, A048138.

%Y Cf. A070015, A135244, A135245.

%K nonn,tabf

%O 2,1

%A _T. D. Noe_, Dec 05 2008