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.

Showing 1-5 of 5 results.

A099308 Numbers m whose k-th arithmetic derivative is zero for some k. Complement of A099309.

Original entry on oeis.org

0, 1, 2, 3, 5, 6, 7, 9, 10, 11, 13, 14, 17, 18, 19, 21, 22, 23, 25, 29, 30, 31, 33, 34, 37, 38, 41, 42, 43, 46, 47, 49, 53, 57, 58, 59, 61, 62, 65, 66, 67, 70, 71, 73, 77, 78, 79, 82, 83, 85, 89, 93, 94, 97, 98, 101, 103, 105, 107, 109, 113, 114, 118, 121, 126, 127, 129, 130
Offset: 1

Views

Author

T. D. Noe, Oct 12 2004

Keywords

Comments

The first derivative of 0 and 1 is 0. The second derivative of a prime number is 0.
For all n, A003415(a(n)) is also a term of the sequence. A351255 gives the nonzero terms as ordered by their position in A276086. - Antti Karttunen, Feb 14 2022

Examples

			18 is on this list because the first through fifth derivatives are 21, 10, 7, 1, 0.

References

Links

Crossrefs

Cf. A003415 (arithmetic derivative of n), A099307 (least k such that the k-th arithmetic derivative of n is zero), A099309 (complement, numbers whose k-th arithmetic derivative is nonzero for all k), A351078 (first noncomposite reached when iterating the derivative from these numbers), A351079 (the largest term on such paths).
Cf. A328308, A328309 (characteristic function and their partial sums), A341999 (1 - charfun).
Cf. A276086, A328116, A351255 (permutation of nonzero terms), A351257, A351259, A351261, A351072 (number of prime(k)-smooth terms > 1).
Cf. also A256750 (number of iterations needed to reach either 0 or a number with a factor of the form p^p), A327969, A351088.
Union of A359544 and A359545.

Programs

  • Mathematica
    dn[0]=0; dn[1]=0; dn[n_]:=Module[{f=Transpose[FactorInteger[n]]}, If[PrimeQ[n], 1, Plus@@(n*f[[2]]/f[[1]])]]; d1=Table[dn[n], {n, 40000}]; nLim=200; lst={1}; i=1; While[i<=Length[lst], currN=lst[[i]]; pre=Intersection[Flatten[Position[d1, currN]], Range[nLim]]; pre=Complement[pre, lst]; lst=Join[lst, pre]; i++ ]; Union[lst]
  • PARI
    \\ The following program would get stuck in nontrivial loops. However, we assume that the conjecture 3 in Ufnarovski & Åhlander paper holds ("The differential equation n^(k) = n has only trivial solutions p^p for primes p").
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));
isA099308(n) = if(!n, 1, while(n>1, n = A003415checked(n)); (n)); \\ Antti Karttunen, Feb 14 2022

Formula

For all n >= 0, A328309(a(n)) = n. - Antti Karttunen, Feb 14 2022

A099309 Numbers n whose k-th arithmetic derivative is nonzero for all k. Complement of A099308.

Original entry on oeis.org

4, 8, 12, 15, 16, 20, 24, 26, 27, 28, 32, 35, 36, 39, 40, 44, 45, 48, 50, 51, 52, 54, 55, 56, 60, 63, 64, 68, 69, 72, 74, 75, 76, 80, 81, 84, 86, 87, 88, 90, 91, 92, 95, 96, 99, 100, 102, 104, 106, 108, 110, 111, 112, 115, 116, 117, 119, 120, 122, 123, 124, 125, 128, 132
Offset: 1

Views

Author

T. D. Noe, Oct 12 2004

Keywords

Comments

Numbers of the form n = m*p^p (where p is prime), i.e., multiples of some term in A051674, have n' = (m + m')*p^p, which is again of the same form, but strictly larger iff m > 1. Therefore successive derivatives grow to infinity in this case, and they are constant when m = 1. There are other terms in this sequence, but I conjecture that they all eventually lead to a term of this form, e.g., 26 -> 15 -> 8 etc. - M. F. Hasler, Apr 09 2015

References

Links

Crossrefs

Cf. A003415 (arithmetic derivative of n), A099307 (least k such that the k-th arithmetic derivative of n is zero), A099308 (numbers whose k-th arithmetic derivative is zero for some k).
Cf. A341999 (characteristic function),
Positions of zeros in A256750, A351078, A351079 (after their initial zeros), also in A328308, A328312.
Subsequences include: A100716, A327929, A327934, A328251, A359547 (intersection with A048103).

Programs

  • PARI
    is(n)=until(4>n=factorback(n~)*sum(i=1,#n,n[2,i]/n[1,i]), for(i=1,#n=factor(n)~,n[1,i]>n[2,i]||return(1))) \\ M. F. Hasler, Apr 09 2015

A351255 Numbers whose k-th arithmetic derivative is zero for some k>0, ordered by their position in A276086.

Original entry on oeis.org

1, 2, 3, 6, 9, 18, 5, 10, 30, 25, 150, 375, 750, 5625, 7, 14, 21, 42, 126, 70, 105, 315, 350, 1575, 3150, 1750, 2625, 49, 98, 882, 490, 735, 4410, 2450, 3675, 11025, 12250, 30625, 61250, 183750, 686, 3430, 5145, 25725, 77175, 385875, 1929375, 3858750, 4802, 72030, 120050, 180075, 33614, 100842, 117649, 705894, 26471025
Offset: 1

Views

Author

Antti Karttunen, Feb 10 2022

Keywords

Comments

Equal to nonzero terms of A099308 when sorted into ascending order. In this order, which is dictated by the primorial base expansion of n (A049345) and mapped to products of prime powers by A276086, all terms of A099308 that are prime(k)-smooth appear before the terms that are not prime(k)-smooth.
Number of terms whose greatest prime factor (A006530) is prime(n) [in other words, that are prime(n)-smooth but not prime(n-1)-smooth] is given by A351071(n): 1, 4, 8, 44, 216, 1474, 11130, ...
For all n > 1, A003415(a(n)) is also a term of the sequence.
Note that only 451 of the first 105367 terms (all 19-smooth terms) are such that there occurs a 19-smooth number (A080682) larger than 1 on the path before 1 is encountered, when starting from x = a(n) and iterating with map x -> A003415(x).

Links

Crossrefs

Cf. A003415, A049345, A099307, A099308, A276086, A328116, A351071, A351072 (number of prime(n)-smooth terms).
Cf. A351256 [= A051903(a(n))], A351257 [= A099307(a(n))], A351258, A351259 [= A351078(a(n))], A351261 [= A351079(a(n))].

Programs

  • PARI
    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));
A099307(n) = { my(s=1); while(n>1, n = A003415checked(n); s++); if(n,s,0); };
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
for(n=0, 2^9, u=A276086(n); c = A099307(u); if(c>0,print1(u, ", ")));

Formula

a(n) = A276086(A328116(n)).

A351079 a(n) is the largest term encountered on the path from n to 0 when iterating the map x -> x', or 0 if 0 cannot be reached from n (or if n is 0). Here x' is the arithmetic derivative of x, A003415.

Original entry on oeis.org

0, 1, 2, 3, 0, 5, 6, 7, 0, 9, 10, 11, 0, 13, 14, 0, 0, 17, 21, 19, 0, 21, 22, 23, 0, 25, 0, 0, 0, 29, 31, 31, 0, 33, 34, 0, 0, 37, 38, 0, 0, 41, 42, 43, 0, 0, 46, 47, 0, 49, 0, 0, 0, 53, 0, 0, 0, 57, 58, 59, 0, 61, 62, 0, 0, 65, 66, 67, 0, 0, 70, 71, 0, 73, 0, 0, 0, 77, 78, 79, 0, 0, 82, 83, 0, 85, 0, 0, 0, 89
Offset: 0

Views

Author

Antti Karttunen, Feb 11 2022

Keywords

Comments

Question: Is there any good upper bound for ratio a(n)/n? See also comments in A351261.

Examples

			For n = 15, if we iterate with A003415, we get a path 15 -> 8 -> 12 -> 16 -> 32 -> 80 -> 176 -> 368 -> ..., where the terms just keep on growing without ever reaching zero, therefore a(15) = 0.
For n = 18, its path down to zero, when iterating A003415 is: 18 -> 21 -> 10 -> 7 -> 1 -> 0, and the largest term is 21, therefore a(18) = 21.

Links

Crossrefs

Cf. A003415, A099307, A099309, A351078, A351261 [= a(A351255(n))].

Programs

  • PARI
    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));
A351079(n) = { my(m=n); while(n>1, n = A003415checked(n); m = max(m,n)); if(n,m); };

Formula

For n > 0, a(n) = 0 if A099307(n) = 0, otherwise a(n) = max(n, a(A003415(n))).
a(0) = 0 and a(A099309(n)) = 0 for all n.

A351259 First noncomposite number reached when iterating the map x -> x', when starting from x = A351255(n). Here x' is the arithmetic derivative of x, A003415.

Original entry on oeis.org

1, 2, 3, 5, 5, 7, 5, 7, 31, 7, 41, 71, 191, 2711, 7, 5, 7, 41, 103, 59, 71, 271, 71, 1031, 2887, 439, 5, 5, 7, 631, 251, 401, 3491, 1031, 1319, 17747, 9733, 1931, 16319, 77351, 131, 5, 419, 7079, 22343, 971, 5981, 6861581, 419, 18731, 11903, 33937, 7079, 15287, 15287, 6143, 6944111, 1415651, 11, 13, 5, 61, 103, 401, 631
Offset: 1

Views

Author

Antti Karttunen, Feb 11 2022

Keywords

Comments

For the initial 105367 19-smooth terms of A351255, the last 7 occurs here at a(54796), with A351255(54796) = 289993286583 = 3^2 * 7 * 11 * 13^2 * 19^5, and the last 5 occurs here at a(65777), with A351255(65777) = 391899820830375516750 = 2 * 3^2 * 5^3 * 7^3 * 13^3 * 17^3 * 19^6, already a moderately high starting value, in whose vicinity most ending primes for successful iterations are much larger. This observation motivates a conjecture: Even from large numbers with high exponents in their prime factorization it is sometimes possible to reach a small prime. Compare to the conjecture 8 in Ufnarovski & Åhlander paper.

Examples

			From A351255(27) = 2625 it takes 12 iterations of the map x -> A003415(x) to reach zero: 2625 -> 2825 -> 1155 -> 886 -> 445 -> 94 -> 49 -> 14 -> 9 -> 6 -> 5 -> 1 -> 0. Two steps before the final zero is the first and only prime on the path, 5, therefore a(27) = 5.

Links

Crossrefs

Cf. A000040, A003415, A276086, A327975, A327977, A351078, A351255, A351257, A351261.

Programs

  • PARI
    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));
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A351078(n) = { while(n>1&&!isprime(n), n = A003415checked(n)); (n); };
for(n=0, 2^9, u=A276086(n); p = A351078(u); if(p>0,print1(p, ", ")));

Formula

a(n) = A351078(A351255(n)).
a(1) = 1, and for n > 1, a(n) = A003415^[A351257(n)-2](A351255(n)). [This means: take the (A351257(n)-2)-th arithmetic derivative of A351255(n)].
Showing 1-5 of 5 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.