A327969 -id:A327969 - cp's OEIS frontend

A351071 Number of integers x in range A002110(n) .. A002110(1+n)-1 such that the k-th arithmetic derivative of A276086(x) is zero for some k, where A002110(n) is the n-th primorial.

1, 4, 8, 44, 216, 1474, 11130, 92489

Offset: 0

Antti Karttunen, Feb 02 2022

a(n) is the number of terms of A328116 in range A002110(n) .. A002110(1+n)-1.
a(n) is the number of terms in A351255 (and in A099308) whose largest prime factor (A006530) is A000040(1+n).
Ratio a(n) / A061720(n) develops as:
  0:     1 / 1       = 1.0
  1:     4 / 4       = 1.0
  2:     8 / 24      = 0.333...
  3:    44 / 180     = 0.244...
  4:   216 / 2100    = 0.1029...
  5:  1474 / 27720   = 0.05317...
  6: 11130 / 480480  = 0.02316...
  7: 92489 / 9189180 = 0.01006...
Computing term a(8) would need processing over 213393180 integers whose greatest prime factor is 23, from single A351255(105368) = 23 at start to product (2^1)*(3^2)*(5*4)*(7^6)*(11^10)*(13^12)*(17^16)*(19^18)*(23^22) at the end of the batch [number whose size in binary is 346 bits], and would required factoring integers of comparable size and more (see A351261), that might not all be easily factorable.

			There are eight terms [6, 7, 9, 12, 15, 20, 21, 28] that are >= A002110(2) and < A002110(3) in A328116 for which the corresponding terms [5, 10, 30, 25, 150, 375, 750, 5625] in A276086 (and A351255) are all in A099308, therefore a(2) = 8.

Cf. A002110, A003415, A061720, A099308, A328306, A328307, A328308, A328116, A351067, A351069, A351255, A351261, A351072 (partial sums).

Cf. also A327969.

\\ Memoization would work quite badly here. (See comments in A351255. In practice sequence A328306 was computed first, up to its term a(9699690). Same data is available in A328116.)
A002110(n) = prod(i=1,n,prime(i));
A003415checked(n) = if(n<=1, 0, my(f=factor(n), s=0); for(i=1, #f~, if(f[i, 2]>=f[i, 1], return(0), s += f[i, 2]/f[i, 1])); (n*s));
A328308(n) = if(!n, 1, while(n>1, n = A003415checked(n)); (n));
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A328306(n) = A328308(A276086(n));
A351071(n) = sum(k=A002110(n),A002110(1+n)-1,A328306(k));

a(n) = Sum_{k=A002110(n) .. A002110(1+n)-1} A328306(k).

a(n) = A328307(A002110(1+n)) - A328307(A002110(n)).

A328249 Numbers k such that the third arithmetic derivative of A276086(k) is prime.

Original entry on oeis.org

5, 15, 21, 31, 43, 79, 91, 98, 104, 106, 223, 229, 231, 254, 255, 289, 291, 301, 305, 423, 453, 456, 487, 670, 674, 677, 692, 702, 730, 736, 2343, 2378, 2380, 2400, 2409, 2534, 2537, 2543, 2552, 2562, 2585, 2602, 2618, 2629, 2767, 2804, 2821, 2831, 2839, 2942, 2943, 2957, 2962, 2963, 2974, 4621, 4669, 4672, 4687, 4717, 4841, 4844

Offset: 1

Antti Karttunen, Oct 11 2019

nonn

Numbers k such that A003415(A003415(A327860(k))) = A099306(A276086(k)) is a prime.
Numbers k such that A276086(k) is in A328239.
For all n, A327969(a(n)) <= 6. This is sharp for example with a(7) = 91.

Antti Karttunen, Table of n, a(n) for n = 1..1315

Cf. A003415, A099306, A276086, A327860, A327969, A328239, A328233, A328240.

A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A327860(n) = { my(m=1, i=0, s=0, pr=1, nextpr); while((n>0), i=i+1; nextpr = prime(i)*pr; if((n%nextpr), my(e=((n%nextpr)/pr)); m *= (prime(i)^e); s += (e / prime(i)); n-=(n%nextpr)); pr=nextpr); (s*m); };
isA328249(n) = isprime(A003415(A003415(A327860(n))));

A353515 The length of the shortest path from n to 1 when using the transitions x -> A003415(x) and x -> A003961(x), or -1 if no 1 can ever be reached from n.

Original entry on oeis.org

0, 1, 1, 4, 1, 2, 1, 7, 3, 2, 1, 6, 1, 4, 7, 8, 1, 4, 1, 6, 3, 2, 1, 7, 3, 6, 7, 8, 1, 2, 1, 10, 5, 2, 6, 5, 1, 4, 4, 6, 1, 2, 1, 6, 5, 4, 1, 9, 4, 5, 5, 8, 1, 6, 8, 8, 3, 2, 1, 4, 1, 6, 5, 10, 5, 2, 1, 5, 4, 2, 1, 8, 1, 5, 6, 6, 5, 2, 1, 9, 7, 2, 1, 4, 3, 5, 7, 8, 1, 5, 7, 7, 3, 5, 4, 9, 1, 6, 7, 6, 1, 3, 1, 7, 2

Offset: 1

Antti Karttunen, Apr 23 2022

nonn

This is a variant of A327969 that seems to be less in need of an escape clause. Note that enough prime shifts with A003961 will eventually transform every term of A100716 (which is a subsequence of A099309) to a term of A048103, and that A051903(A003961(n)) = A051903(n). See also the array A344027.

Records 0, 1, 4, 7, 8, 10, 12, 13, 14, 15, 16, 19, ... occur at 1, 2, 4, 8, 16, 32, 128, 256, 768, 1024, 2048, 4096, ..., etc.

			From n = 4, we can reach 1 with just four steps as A003961(4) = 9, A003415(9) = 6, A003415(6) = 5 and A003415(5) = 1, and because there are no shorter paths we have a(4) = 4.
From n = 8, we can reach 1 with seven steps, as A003415(8) = 12, A003961(12) = 45, A003415(45) = 39, A003961(39) = 85, A003415(85) = 22, A003415(22) = 13, A003415(13) = 1, and because there are no shorter paths we have a(8) = 7.
For n = 15, as A003415(15) = 8, we know that a(15) is at most a(8)+1, i.e., 8. But we can do better, as A003961(15) = 35, A003961(35) = 77, A003415(77) = 18, A003415(18) = 21, A003415(21) = 10, A003415(10) = 7, A003415(7) = 1, and because there are no shorter paths we have a(15) = 7.
From n = 49, we can reach 1 in four steps, as A003961(49) = 121, A003415(121) = 22, A003415(22) = 13, A003415(13) = 1. Note that this is less than A099307(49)-1, as it would take 5 steps to reach 1 if using the arithmetic derivative only, 49 -> 14 -> 9 -> 6 -> 5 -> 1.

Antti Karttunen, Table of n, a(n) for n = 1..16383
Index entries for sequences computed from indices in prime factorization

Cf. A003415, A003961, A048103, A051903, A099307, A099308, A099309, A100716, A349905.

Cf. also A327969 and A343221, A344027, A353484.

A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A003961(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); }; \\ From A003961
A353515(n) = if(1==n,0,my(xs=Set([n]),newxs,a,b,u); for(k=1,oo, newxs=Set([]); for(i=1,#xs,u = xs[i]; a = A003415(u); if(1==a, return(k)); if(isprime(a), return(k+1)); b = A003961(u); newxs = setunion([a],newxs); newxs = setunion([b],newxs)); xs = newxs));

a(1) = 0, a(p^p) = 1 + a(A003961(p^p)) for primes p, and for other numbers, a(n) = 1 + min(a(A003415(n)), a(A003961(n))).
a(p) = 1 for all primes p.
a(n) < A099307(n), unless A099307(n) = 0. [I.e., for all n in A099308]

A369652 Maximal digit in the primorial base representation of the n-th arithmetic derivative of 128.

Original entry on oeis.org

4, 4, 6, 2, 6, 11, 8, 7, 11, 11, 8, 10, 15, 15, 12, 18, 17, 30, 28, 22, 21, 37, 28, 38, 42, 33, 22, 52, 43, 56, 47, 49, 57, 60, 55, 54, 58, 70, 52, 66, 68, 57, 63, 86, 58, 88, 92, 66, 78, 95, 85, 52, 102, 70, 111, 57, 117, 99, 136, 104, 129, 110, 146, 127, 135, 132, 131, 129, 126, 145, 112, 150, 128, 129, 154, 161, 145

Offset: 0

Antti Karttunen, Feb 04 2024

nonn

This sequence relates to the question whether A327969(128) has a positive integer value, or whether it is -1 by the escape clause.

Note that when iterating the map k -> k' from A276086(A369638(4)) = A276086(15232) = 3299611946113357875 onward, the maximal exponent in the prime factorization (A051903) keeps on decreasing until it reaches 1 at the fifth iteration, and then stays as 1 for three more iterations (with k then 38863666759992439 = 643*60441161368573), but then, alas, on the next iteration, k' = 60441161369216 = 2^7 * 472196573197.

			The third arithmetic derivative (A099306) of 128 is 5056, which in primorial base (A049345) is written as 220220, therefore a(3) = 2.
The fourth arithmetic derivative (A258644) of 128 is 15232, which in primorial base is written as 663320, therefore a(4) = 6.

Sean A. Irvine, Table of n, a(n) for n = 0..142
Index entries for sequences related to primorial base

Cf. A003415, A049345, A099306, A258644, A328114, A369638.

Cf. also A327969.

a(n) = A328114(A369638(n)).

A378209 Antiderivatives of 334406399, numbers k for which A003415(k) = A024451(9) = A003415(A002110(9)).

Original entry on oeis.org

223092870, 975351895, 1527890095, 1885679383, 2189118743, 2329696457, 2338611863, 3485765789, 4586671213, 5453593183, 5472849253, 5674340053, 8071055747, 8931775397, 9332889127, 9453996491, 9601098443, 10293819917, 12717530039, 17343441881, 18636581773, 19498393573, 20167656703, 23244839627, 23515890737, 23556538969

Offset: 1

Antti Karttunen, Nov 20 2024

Apart from the initial term A002110(9), all other terms are products of three distinct odd primes, A046389. Compare to the comments in A369239.
Note that A024451(9) = 334406399 = 43 * 163 * 47711 == -1 (mod 12). Compare the sequences A369450, A369451 and A369452 to see why there is such a sudden peak in A377993 at n=9, when compared to the nearby terms before and after.
For all n=1..330, A327969(a(n)) <= 7 = A099307(a(n)), because, when we apply A003415 successively, we get: A003415(334406399) -> 9835475 [= A369651(9)] -> 4893565 -> 978718 -> 564671 (which is a prime) -> 1 -> 0.

Antti Karttunen, Table of n, a(n) for n = 1..330
Index entries for sequences related to primorial numbers

Row 9 of irregular triangle A377992.
Subsequence of A099308, and after the initial term, subsequence of A046389.
Cf. A002110, A003415, A024451, A099307, A327969, A369239, A369651.
Cf. A369450, A369451, A369452.

cp's OEIS Frontend

A351071 Number of integers x in range A002110(n) .. A002110(1+n)-1 such that the k-th arithmetic derivative of A276086(x) is zero for some k, where A002110(n) is the n-th primorial.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

PARI

Formula

A328249 Numbers k such that the third arithmetic derivative of A276086(k) is prime.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

A353515 The length of the shortest path from n to 1 when using the transitions x -> A003415(x) and x -> A003961(x), or -1 if no 1 can ever be reached from n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

Formula

A369652 Maximal digit in the primorial base representation of the n-th arithmetic derivative of 128.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Formula

A378209 Antiderivatives of 334406399, numbers k for which A003415(k) = A024451(9) = A003415(A002110(9)).

Views

Author

Keywords

Comments

Links

Crossrefs