A349093 a(n) is the smallest nonprime number m (m = A018252(t)) such that n divides the product P(t) of all nonprime numbers up to and including m (P(t) = A036691(t-1)).

1, 4, 6, 4, 10, 6, 14, 6, 9, 10, 22, 6, 26, 14, 10, 8, 34, 9, 38, 10, 14, 22, 46, 6, 15, 26, 9, 14, 58, 10, 62, 8, 22, 34, 14, 9, 74, 38, 26, 10, 82, 14, 86, 22, 10, 46, 94, 8, 21, 15, 34, 26, 106, 9, 22, 14, 38, 58

Offset

1,2

Comments

a(n) >= 2*gpf(n) for n > 1, where gpf(n) denotes the greatest prime factor of n (A006530(n)).

Conjecture: the equation a(n) = a(n+1) has no solutions. This holds up to at least n = 10^7.

Consecutive solutions of the equation a(n) = 2*K(n) (where K(n) is the Kempner number A002034(n)) are consecutive terms of A048839.

Links

Table of n, a(n) for n=1..58.

J. Sondow and E. W. Weisstein, MathWorld: Smarandache Function

Formula

a(p) = 2*p for prime p.

a(p_1*p_2*...*p_u) = 2*p_u, where p_i's are distinct primes and p_1 < p_2 < ... < p_u.

a(n) where n is factored as n = p_1^k_1*p_2^k_2*...*p_u^k_u is given by a(n) = max( a(p_1^k_1), a(p_2^k_2), ..., a(p_u^k_u) ), where a(p_i^k_i) = w*p_i and w is the smallest m >= 2 satisfying the inequality:

-1 + Sum_{t=1..m} floor((m*p_i)/(p_i)^t) >= k_i.

Example

a(15) = 10 because:
15 does not divide 1=A036691(0)=1, 1*4=A036691(1)=4, 1*4*6=A036691(2)=24, 1*4*6*8=A036691(3)=192, 1*4*6*8*9=A036691(4)=1728 and does divide 1*4*6*8*9*10=A036691(5)=17280.

Prog

(Maxima)
f(p, k):=(z:2, for m:2 while -1+sum(floor((p*m)/(p^t)), t, 1, m)<k do (z:z+1), z*p);
a(n):=(for d:1 thru length(ifactors(n)) do (A[d]:f(ifactors(n)[d][1], ifactors(n)[d][2])), lmax(makelist(A[d], d, 1, length(ifactors(n)))));
makelist(a(n), n, 2, 200);

Crossrefs

Cf. A002034, A006530, A018252, A036691.

Sequence in context: A088738 A114595 A256415 * A143545 A338155 A328045

Adjacent sequences: A349090 A349091 A349092 * A349094 A349095 A349096

Keyword

nonn

Author

Lechoslaw Ratajczak, Mar 25 2022

Status

approved