cp's OEIS Frontend

A353125 a(1)=2. If a(n) is a novel term, a(n+1) = sopfr(a(n)), else if there are k occurrences of a(j)=a(n), (1<=j<=n), a(n+1)=k*a(n).

%I A353125 #32 Jul 01 2025 23:33:40
%S A353125 2,2,4,4,8,6,5,5,10,7,7,14,9,6,12,7,21,10,20,9,18,8,16,8,24,9,27,9,36,
%T A353125 10,30,10,40,11,11,22,13,13,26,15,8,32,10,50,12,24,48,11,33,14,28,11,
%U A353125 44,15,30,60,12,36,72,12,48,96,13,39,16,32,64,12,60,120
%N A353125 a(1)=2. If a(n) is a novel term, a(n+1) = sopfr(a(n)), else if there are k occurrences of a(j)=a(n), (1<=j<=n), a(n+1)=k*a(n).
%C A353125 1 and 3 cannot be terms. Let conditions 1, 2 refer to the effects of a novel and repeat term respectively, and let M(m) be the multiplicity of m in the sequence. M(m) >= 2 for all m in N \ {1,3}. All first occurrences of a prime p > 2 follow a novel term in A046363 (by condition 1).
%C A353125 If the first occurrence of composite m arises from condition 1, then so does the second. Proof: Suppose not, then the second m must be a multiple w*t of a prior term t (condition 2); w >= 2. The term following the second m must be (w+1)*t, and this must equal 2*m (condition 2). Thus 2*w*t = (w+1)*t, so then w=1. But w >= 2; contradiction. Corollary: Once composite m has occurred by condition 1, then all subsequent occurrences of m occur the same way.
%C A353125 Most composite terms appear first by condition 2 (e.g., 12 as 2 copies of 6), and then subsequently by condition 1. Thus for all k in the sequence M(k) >= A000607(k), with equality when k is prime > 2, or certain composite numbers.
%C A353125 Conjecture: The 10 composite numbers 6,9,15,25,35,49,77,121,143,169 behave as primes in this sequence, namely for any such m, M(m) = A000607(m). For all other composite m, M(m) > A000607(m), i.e., at least one (up front) copy by condition 2.
%C A353125 The first occurrences of primes appear in natural order initially, but this is not sustained (e.g., 61 appears before 59).
%H A353125 Michael De Vlieger, <a href="/A353125/b353125.txt">Table of n, a(n) for n = 1..10000</a>
%H A353125 Michael De Vlieger, <a href="/A353125/a353125.png">Annotated log-log scatterplot of a(n)</a>, n = 1..2^12, records in red and likely local minima in blue.
%H A353125 Michael De Vlieger, <a href="/A353125/a353125_1.png">Log-log scatterplot of a(n)</a>, n = 1..2^16.
%e A353125 a(1)=2, a novel term, so a(2)=sopfr(2)=2. 2 occurs twice and is the only prime p whose multiplicity is not A000607(p), simply because it is the seed term.
%e A353125 Since 2 has now appeared twice, a(3)=2*2=4, a novel term, so a(4)=sopfr(4)=4.
%e A353125 a(25)=24 (3 occurrences of 8), a(46)=24 (2 occurrences of 12). Subsequently all occurrences of 24 are from condition 1. Therefore M(24) = 2 + A000607(24) = 48.
%t A353125 Block[{a, c, j, k, nn}, nn = 120; c[_] = 0; j = a[1] = 2; c[2]++; Do[If[c[j] == 1, Set[k, Total@ Flatten[ConstantArray[#1, #2] & @@@ FactorInteger[j]]], Set[k, c[j] j]]; j = a[i] = k; c[k]++, {i, 2, nn}]; Array[a, nn] ] (* _Michael De Vlieger_, Apr 24 2022 *)
%o A353125 (PARI) sopfr(n) = (n=factor(n))[, 1]~*n[, 2]; \\ A001414
%o A353125 lista(nn) = {my(v=vector(nn), k); v[1] = 2; for (n=2, nn, if ((k=#select(x->(x==v[n-1]), Vec(v, n-1))) == 1, v[n] = sopfr(v[n-1]), v[n] = k*v[n-1]);); v;} \\ _Michel Marcus_, May 16 2022
%Y A353125 Cf. A001414 (sopfr), A046363, A000607, A000040.
%K A353125 nonn
%O A353125 1,1
%A A353125 _David James Sycamore_ and _Michael De Vlieger_, Apr 24 2022
%E A353125 More terms from _Michael De Vlieger_, Apr 24 2022