%I #21 Jun 06 2021 02:53:58
%S 6,6,6,8,10,12,12,12,15,15,18,18,20,20,20,24,24,24,24,24,28,30,30,30,
%T 30,35,35,35,36,36,40,40,40,40,40,45,45,45,45,48,48,48,54,54,54,56,56,
%U 56,56,60,60,60,60,60,60,63,70,70,70,70,70,70,70,72,72,72
%N a(n) = n + A332558(n) + 1.
%H Robert Israel, <a href="/A332559/b332559.txt">Table of n, a(n) for n = 1..10000</a>
%H J. S. Myers, R. Schroeppel, S. R. Shannon, N. J. A. Sloane, and P. Zimmermann, <a href="http://arxiv.org/abs/2004.14000">Three Cousins of Recaman's Sequence</a>, arXiv:2004:14000 [math.NT], April 2020.
%p f:= proc(n) local k,p;
%p p:= n;
%p for k from 1 do
%p p:= p*(n+k);
%p if (p/(n+k+1))::integer then return n+k+1 fi
%p od
%p end proc:
%p map(f, [$1..100]); # _Robert Israel_, Feb 25 2020
%t a[n_] := Module[{k, p = n}, For[k = 1, True, k++, p *= (n+k); If[Divisible[p, n+k+1], Return[n+k+1]]]];
%t Array[a, 100] (* _Jean-François Alcover_, Jul 18 2020, after Maple *)
%o (Python)
%o def a(n):
%o k, p = 1, n*(n+1)
%o while p%(n+k+1): k += 1; p *= (n+k)
%o return n + k + 1
%o print([a(n) for n in range(1, 67)]) # _Michael S. Branicky_, Jun 06 2021
%Y Cf. A061836, A332558, A332560, A332561.
%K nonn
%O 1,1
%A _Scott R. Shannon_ and _N. J. A. Sloane_, Feb 24 2020