cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A333518 a(n) = A000720(A006530(A334468(n))).

This page as a plain text file.
%I A333518 #48 May 09 2020 00:45:48
%S A333518 1,2,1,2,3,1,2,2,2,3,1,2,3,3,2,2,3,4,1,4,2,3,3,2,3,2,3,4,2,3,3,1,3,4,
%T A333518 2,3,3,2,4,4,3,4,2,3,4,2,4,3,6,3,2,3,1,3,4,2,4,3,5,6,4,3,2,4,4,4,3,7,
%U A333518 3,4,2,4,3,3,6,4,6,2,4,4,3,5,6,8,7,3,2,5,3,4,1,4,5,3,5,4,4,2,4,5,3,5,6,3,4
%N A333518 a(n) = A000720(A006530(A334468(n))).
%C A333518 Indices of the greatest prime factor of A334468(n).
%C A333518 Consider A334468, a list of numbers m = n+j such that j > 0 is also the smallest number such that n+j has no prime factor > j for some n and j = A217287(n).
%C A333518 Since prime q always contributes a novel prime divisor (i.e., q itself) to the set of distinct primes that divide at least 1 number i the range n + i (1 <= i <= j), the numbers m in A334468 are composite, and given the above, m is a product of relatively small prime factors.
%e A333518 Start with n = 1, the empty product. Incrementing n and storing the distinct prime factors each time, we encounter 2, which does not divide any previous number n. Therefore we proceed to n = 3, which is prime and its distinct prime divisor again does not divide any previous number. Finally, at 4, we have the distinct prime divisor 2, since 2 divides the product of the previous range {1, 2, 3}, we end the chain. Therefore 4 is the first term of this sequence.
%e A333518 We list row n of A217438 below, starting with n aligned in columns:
%e A333518 1  2  3
%e A333518    2  3
%e A333518       3  4  5
%e A333518          4  5  6  7
%e A333518             5  6  7
%e A333518                6  7
%e A333518                   7  8  9  10  11
%e A333518                      8  9  10  11
%e A333518                         9  10  11
%e A333518                            10  11  12  13  14
%e A333518                                11  12  13  14  15
%e A333518                                    12  13  14  15
%e A333518                                        13  14  15
%e A333518                                            14  15
%e A333518                                                ...
%e A333518 Adding 1 to the last numbers seen in all the rows, we generate the sequence A334468: {4, 6, 8, 12, 15, 16, ...}. Of these, we have greatest prime factors {2, 3, 2, 3, 5, 2, ...} with indices {1, 2, 1, 2, 3, 1, ...}.
%e A333518 Least indices of prime(k) in a(n):
%e A333518    i  p(i)    n    a(n)
%e A333518   ---------------------
%e A333518    1    2     1      4
%e A333518    2    3     2      6
%e A333518    3    5     5     15
%e A333518    4    7    18     63
%e A333518    5   11    59    308
%e A333518    6   13    49    234
%e A333518    7   17    68    374
%e A333518    8   19    84    475
%e A333518    9   23   292   2392
%e A333518   10   29   401   3625
%e A333518   11   31   518   4991
%e A333518   12   37   791   8547
%e A333518   ...
%t A333518 Block[{nn = 2^10, r}, r = Array[If[# == 1, 0, Total[2^(PrimePi /@ FactorInteger[#][[All, 1]] - 1)]] &, nn]; Map[PrimePi@ FactorInteger[#][[-1, 1]] &, #] &@ Union@ Array[Block[{k = # + 1, s = r[[#]]}, While[UnsameQ[s, Set[s, BitOr[s, r[[k]] ] ] ], k++]; k] &, nn - Ceiling@ Sqrt@ nn] ]
%Y A333518 Cf. A000720, A006530, A217287, A217438, A334468.
%K A333518 nonn
%O A333518 1,2
%A A333518 _Michael De Vlieger_, May 05 2020