login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

a(1)=2, a(n) = 1 + (smallest integer > a(n-1) with the same number of divisors as a(n-1)).
1

%I #13 Aug 03 2016 18:18:43

%S 2,4,10,15,22,27,34,36,101,104,106,112,163,168,211,224,229,234,261,

%T 269,272,305,310,319,322,329,335,340,343,347,350,353,360,421,432,529,

%U 842,844,846,851,863,878,880,913,915,936,991,998

%N a(1)=2, a(n) = 1 + (smallest integer > a(n-1) with the same number of divisors as a(n-1)).

%H Harvey P. Dale, <a href="/A128513/b128513.txt">Table of n, a(n) for n = 1..1001</a>

%e a(6) = 27, which has 4 positive divisors. The next larger

%e integer with 4 divisors is 33. So a(7) = 1 + 33 = 34.

%t a = {2}; n = 3; While[n < 1000, If[Length[Divisors[n]] == Length[Divisors[a[[ -1]]]], AppendTo[a, n + 1]; n++;]; n++ ]; a (* _Stefan Steinerberger_, Jun 15 2007 *)

%t nxt[n_]:=Module[{k=n+1,ds=DivisorSigma[0,n]},While[DivisorSigma[ 0,k] != ds, k++]; k+1]; NestList[nxt,2,50] (* _Harvey P. Dale_, Aug 03 2016 *)

%K nonn

%O 1,1

%A _Leroy Quet_, Jun 12 2007

%E More terms from _Stefan Steinerberger_, Jun 15 2007