cp's OEIS Frontend

A362128 Indices of records in A360179.

%I A362128 #8 Apr 14 2023 11:36:12
%S A362128 1,3,5,7,9,12,15,21,25,29,52,69,87,100,138,176,233,252,265,338,371,
%T A362128 384,400,447,459,492,515,594,606,677,779,794,850,896,1100,1144,1204,
%U A362128 1225,1243,1584,1674,1700,2214,2239,2265,2338,2568,3014,3621,3652,4027,4058,4437,4464,4584,4625,5111,5201,5281
%N A362128 Indices of records in A360179.
%H A362128 Michael De Vlieger, <a href="/A362128/b362128.txt">Table of n, a(n) for n = 1..8238</a>
%t A362128 nn = 6000; c[_] := False; h[_] := 0; f[n_] := DivisorSigma[0, n];
%t A362128 a[1] = j = u = 1; r = 0;
%t A362128 {1}~Join~Rest@ Reap[Do[
%t A362128       If[c[j],
%t A362128         k = j + f[u]; h[j]++; h[u]--,
%t A362128         k = f[j]; c[j] = True; h[j]++ ];
%t A362128       u = Min[u, j]; Set[{a[n], j}, {k, k}]; While[h[u] == 0, u++];
%t A362128       If[k > r, r = k; Sow[n]], {n, 2, nn}] ][[-1, -1]]
%Y A362128 Cf. A000005, A360179, A362127.
%K A362128 nonn
%O A362128 1,2
%A A362128 _Michael De Vlieger_, Apr 09 2023