A357735 - OEIS

%I #10 Oct 12 2022 09:49:31

%S 1,2,6,3,4,8,9,11,10,12,14,5,15,16,18,20,7,21,13,24,27,22,26,28,23,25,

%T 30,32,33,31,34,35,40,44,55,36,37,39,17,42,45,38,46,48,50,19,57,52,41,

%U 62,43,54,56,58,60,64,51,63,66,49,70,77,68,69,72,75,74,76

%N a(1)=1, a(2)=2. Thereafter a(n+1) is least k != partial sum s(n) which has not occurred earlier, such that gcd(k, s(n)) > 1.

%C It follows from the definition that if s(n) is prime then a(n+1) = 2*s(n). This happens only once in the sequence, when a(3)=6, following s(2)=3. For all n > 2 s(n) is composite. Conjectured to be a permutation of the positive integers (primes not in natural order).

%H Michael De Vlieger, <a href="/A357735/a357735.png">Log-log scatterplot of a(n)</a> n = 1..2^17, labeling records in red, local minima in blue, and highlighting prime terms in green.

%e Since a(1)=1 and a(2)=2, we have s(2)=3, then a(3) is 6, the smallest unused term sharing a divisor with 3.

%t nn = 2^16; c[_] = False; Array[Set[{a[#], c[#]}, {#, True}] &, 2]; u = s = 3; Do[k = u; While[Nand[! c[k], ! CoprimeQ[k, s], k != s], k++]; Set[{a[n], c[k]}, {k, True}]; s += k; If[k == u, While[c[u], u++]], {n, 3, nn}]; Array[a, nn] (* _Michael De Vlieger_, Oct 11 2022 *)

%Y Cf. A064413, A084385, A347113, A351743.

%K nonn

%O 1,2

%A _David James Sycamore_, Oct 11 2022