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.

A325651 a(n) = greatest k such that sigma(k) = sigma(n).

This page as a plain text file.
%I A325651 #17 Jul 26 2025 08:07:55
%S A325651 1,2,3,4,5,11,7,8,9,17,11,12,13,23,23,25,17,18,19,41,31,22,23,59,25,
%T A325651 41,27,39,29,71,31,32,47,53,47,36,37,59,39,89,41,77,43,83,45,71,47,75,
%U A325651 49,50,71,97,53,95,71,95,79,89,59,167,61,77,103,64,83,119
%N A325651 a(n) = greatest k such that sigma(k) = sigma(n).
%H A325651 Amiram Eldar, <a href="/A325651/b325651.txt">Table of n, a(n) for n = 1..10000</a>
%e A325651 a(6) = 11 because sigma(6) = sigma(11) = 12.
%t A325651 a[n_] := Block[{s = DivisorSigma[1, n], k}, k=s; While[ DivisorSigma[1, k] != s, k--]; k]; Array[a, 66] (* _Giovanni Resta_, May 20 2019 *)
%o A325651 (Magma) [Max([k: k in[1..10000] | SumOfDivisors(k) eq SumOfDivisors(n)]): n in [1..100]];
%o A325651 (PARI) a(n) = {my(s=sigma(n)); forstep(i=s, 1, -1, if (sigma(i) == s, return(i)););} \\ _Michel Marcus_, May 12 2019
%Y A325651 Cf. A000203, A275987 (least k such that sigma(k) = sigma(n)).
%Y A325651 See A070242, A325652 and A325653 for number, sum and product of such numbers k.
%K A325651 nonn
%O A325651 1,2
%A A325651 _Jaroslav Krizek_, May 12 2019