OFFSET
1,2
COMMENTS
If rad(s(n-1)) > rad(s(n-2)), then a(n) = rad(s(n-1)).
Let r = rad(s(n-1)). a(n) = r implies a(n+1) = 2*r for n >= 3.
Let R(r) = {m : rad(m) | r} sorted such that R(r,i) denotes the i-th smallest term.
The sequence can be interpreted as an irregular table T(i,j) = R(r,j)*r, where row i shares the same squarefree kernel r. Since T(i,j) == 0 (mod r), we may express the sequence of partial sums also as an irregular table S(i,j) = m*r that continues so long as rad(m) | r.
When rad(m) does not divide r, we then have S(i+1,1) = m'*r' where m' = m*r/r' and r' > r for r' = rad(S(i+1,1)).
Since the squarefree kernels r are nondecreasing as n increases, this sequence is not a permutation of natural numbers and only a finite number of smallest terms that have the same squarefree kernel appear in the sequence. For example, the intersection of A033845 and this sequence is {6, 12, 18}. Furthermore, not all squarefree numbers appear in this sequence.
Intersection of this sequence and A001694 is {1}.
Intersection of this sequence and the primes is {2, 3}.
The only primorials (in A002110) in this sequence appear to be {1, 2, 6}.
The behavior of this sequence regarding its partial sums is similar to that of A392975 regarding its partial sums.
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
EXAMPLE
Table of n, a(n), s(n) for n = 1..16:
n a(n) = r * m s(n) rad(s(n)) = r
----------------------------------------------
1 1 = 1 * 1 1 1
2 2 = 2 * 1 3 3
3 3 = 3 * 1 6 2 * 3
4 6 = 6 * 1 12 2 * 3
5 12 = 6 * 2 24 2 * 3
6 18 = 6 * 3 42 2 * 3 * 7
7 42 = 42 * 1 84 2 * 3 * 7
8 84 = 42 * 2 168 2 * 3 * 7
9 126 = 42 * 3 294 2 * 3 * 7
10 168 = 42 * 4 462 2 * 3 * 7 * 11
11 462 = 462 * 1 924 2 * 3 * 7 * 11
12 924 = 462 * 2 1848 2 * 3 * 7 * 11
13 1386 = 462 * 3 3234 2 * 3 * 7 * 11
14 1848 = 462 * 4 5082 2 * 3 * 7 * 11
15 2772 = 462 * 6 7854 2 * 3 * 7 * 11 * 17
16 7854 = 7854 * 1 15708 2 * 3 * 7 * 11 * 17
This sequence as an irregular table T(i,j) where row i has the same squarefree kernel:
i\j | 1 2 3 4 5 6 7 8
----+--------------------------------------------------------
0 | 1;
1 | 2;
2 | 3;
3 | 6, 12, 18;
4 | 42, 84, 126, 168;
5 | 462, 924, 1386, 1848, 2772;
6 | 7854, 15708, 23562, 31416, 47124, 54978, 62832, 70686;
...
In this table, T(i,1) is squarefree and all the rest of the terms are nonsquarefree.
MATHEMATICA
Block[{a, c, f, k, r, s, u, nn}, nn = 40; c[_] := False; u[_] := 1; f[x_] := Times @@ FactorInteger[x][[;; , 1]]; a = Range[1, 2]; Map[Set[c[#], True] &, a]; s = Total[a]; a~Join~Reap[Do[r = f[s]; While[Set[k, u[r]*r]; Nand[! c[k], Divisible[r, f@ u[r] ] ], u[r]++]; Sow[k]; s += k; Set[c[k], True], {n, Length[a], nn}] ][[-1, 1]] ]
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Michael De Vlieger and James C. McMahon, Mar 31 2026
STATUS
approved