A328310 -id:A328310 - cp's OEIS frontend

A003415 a(n) = n' = arithmetic derivative of n: a(0) = a(1) = 0, a(prime) = 1, a(mn) = ma(n) + n*a(m).

0, 0, 1, 1, 4, 1, 5, 1, 12, 6, 7, 1, 16, 1, 9, 8, 32, 1, 21, 1, 24, 10, 13, 1, 44, 10, 15, 27, 32, 1, 31, 1, 80, 14, 19, 12, 60, 1, 21, 16, 68, 1, 41, 1, 48, 39, 25, 1, 112, 14, 45, 20, 56, 1, 81, 16, 92, 22, 31, 1, 92, 1, 33, 51, 192, 18, 61, 1, 72, 26, 59, 1, 156, 1, 39, 55, 80, 18, 71

Offset: 0

N. J. A. Sloane, R. K. Guy

Can be extended to negative numbers by defining a(-n) = -a(n).
Based on the product rule for differentiation of functions: for functions f(x) and g(x), (fg)' = f'g + fg'. So with numbers, (ab)' = a'b + ab'. This implies 1' = 0. - Kerry Mitchell, Mar 18 2004
The derivative of a number x with respect to a prime number p as being the number "dx/dp" = (x-x^p)/p, which is an integer due to Fermat's little theorem. - Alexandru Buium, Mar 18 2004
The relation (ab)' = a'b + ab' implies 1' = 0, but it does not imply p' = 1 for p a prime. In fact, any function f defined on the primes can be extended uniquely to a function on the integers satisfying this relation: f(Product_i p_i^e_i) = (Product_i p_i^e_i) * (Sum_i e_i*f(p_i)/p_i). - Franklin T. Adams-Watters, Nov 07 2006
See A131116 and A131117 for record values and where they occur. - Reinhard Zumkeller, Jun 17 2007
Let n be the product of a multiset P of k primes. Consider the k-dimensional box whose edges are the elements of P. Then the (k-1)-dimensional surface of this box is 2*a(n). For example, 2*a(25) = 20, the perimeter of a 5 X 5 square. Similarly, 2*a(18) = 42, the surface area of a 2 X 3 X 3 box. - David W. Wilson, Mar 11 2011
The arithmetic derivative n' was introduced, probably for the first time, by the Spanish mathematician José Mingot Shelly in June 1911 with "Una cuestión de la teoría de los números", work presented at the "Tercer Congreso Nacional para el Progreso de las Ciencias, Granada", cf. link to the abstract on Zentralblatt MATH, and L. E. Dickson, History of the Theory of Numbers. - Giorgio Balzarotti, Oct 19 2013
a(A235991(n)) odd; a(A235992(n)) even. - Reinhard Zumkeller, Mar 11 2014
Sequence A157037 lists numbers with prime arithmetic derivative, i.e., indices of primes in this sequence. - M. F. Hasler, Apr 07 2015
Maybe the simplest "natural extension" of the arithmetic derivative, in the spirit of the above remark by Franklin T. Adams-Watters (2006), is the "pi based" version where f(p) = primepi(p), see sequence A258851. When f is chosen to be the identity map (on primes), one gets A066959. - M. F. Hasler, Jul 13 2015
When n is composite, it appears that a(n) has lower bound 2*sqrt(n), with equality when n is the square of a prime, and a(n) has upper bound (n/2)*log_2(n), with equality when n is a power of 2. - Daniel Forgues, Jun 22 2016
If n = p1*p2*p3*... where p1, p2, p3, ... are all the prime factors of n (not necessarily distinct), and h is a real number (we assume h nonnegative and < 1), the arithmetic derivative of n is equivalent to n' = lim_{h->0} ((p1+h)*(p2+h)*(p3+h)*... - (p1*p2*p3*...))/h. It also follows that the arithmetic derivative of a prime is 1. We could assume h = 1/N, where N is an integer; then the limit becomes {N -> oo}. Note that n = 1 is not a prime and plays the role of constant. - Giorgio Balzarotti, May 01 2023

			6' = (2*3)' = 2'*3 + 2*3' = 1*3 + 2*1 = 5.
Note that, for example, 2' + 3' = 1 + 1 = 2, (2+3)' = 5' = 1. So ' is not linear.
G.f. = x^2 + x^3 + 4*x^4 + x^5 + 5*x^6 + x^7 + 12*x^8 + 6*x^9 + 7*x^10 + ...

G. Balzarotti, P. P. Lava, La derivata aritmetica, Editore U. Hoepli, Milano, 2013.
E. J. Barbeau, Problem, Canad. Math. Congress Notes, 5 (No. 8, April 1973), 6-7.
L. E. Dickson, History of the Theory of Numbers, Vol. 1, Chapter XIX, p. 451, Dover Edition, 2005. (Work originally published in 1919.)
A. M. Gleason et al., The William Lowell Putnam Mathematical Competition: Problems and Solutions 1938-1964, Math. Assoc. America, 1980, p. 295.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe, Table of n, a(n) for n = 0..10000
Krassimir T. Atanassov, A formula for the n-th prime number, Comptes rendus de l'Académie bulgare des Sciences, Tome 66, No 4, 2013.
E. J. Barbeau, Remark on an arithmetic derivative, Canad. Math. Bull. vol. 4, no. 2, May 1961.
A. Buium, Home Page
A. Buium, Differential characters of Abelian varieties over p-adic fields, Invent. Math. 122 (1995), no. 2, 309-340.
A. Buium, Geometry of p-jets, Duke Math. J. 82 (1996), no. 2, 349-367.
A. Buium, Arithmetic analogues of derivations, J. Algebra 198 (1997), no. 1, 290-299.
A. Buium, Differential modular forms, J. Reine Angew. Math. 520 (2000), 95-167.
Brad Emmons and Xiao Xiao, The Arithmetic Partial Derivative, arXiv:2201.12453 [math.NT], 2022.
José María Grau and Antonio M. Oller-Marcén, Giuga Numbers and the Arithmetic Derivative, Journal of Integer Sequences, Vol. 15 (2012), #12.4.1.
R. K. Guy, Letter to N. J. A. Sloane, Apr 1975
P. Haukkanen, M. Mattila, J. K. Merikoski and T. Tossavainen, Can the Arithmetic Derivative be Defined on a Non-Unique Factorization Domain?, Journal of Integer Sequences, 16 (2013), #13.1.2. - From _N. J. A. Sloane_, Feb 03 2013
P. Haukkanen, J. K. Merikoski and T. Tossavainen, Asymptotics of partial sums of the Dirichlet series of the arithmetic derivative, Mathematical Communications 25 (2020), 107-115.
Antti Karttunen, Program in LODA-assembly
J. Kovič, The Arithmetic Derivative and Antiderivative, Journal of Integer Sequences 15 (2012), Article 12.3.8.
Michael Penn, When the derivative of a number is not zero -- The arithmetic derivative., YouTube video, 2022.
Ivars Peterson, Deriving the Structure of Numbers, Science News, March 20, 2004.
D. J. M. Shelly, Una cuestión de la teoria de los numeros, Asociation Esp. Granada 1911, 1-12 S (1911). (Abstract of ref. JFM42.0209.02 on zbMATH.org)
T. Tossavainen, P. Haukkanen, J. K. Merikoski, and M. Mattila, We can differentiate numbers, too, The College Mathematics Journal 55 (2024), no. 2, 100-108.
Victor Ufnarovski and Bo Åhlander, How to Differentiate a Number, J. Integer Seqs., Vol. 6, 2003, #03.3.4.
Linda Westrick, Investigations of the Number Derivative, Siemens Foundation competition 2003 and Intel Science Talent Search 2004.
Wikipedia, Arithmetic derivative

Cf. A086134 (least prime factor of n').
Cf. A086131 (greatest prime factor of n').
Cf. A068719 (derivative of 2n).
Cf. A068720 (derivative of n^2).
Cf. A068721 (derivative of n^3).
Cf. A001787 (derivative of 2^n).
Cf. A027471 (derivative of 3^(n-1)).
Cf. A085708 (derivative of 10^n).
Cf. A068327 (derivative of n^n).
Cf. A024451 (derivative of p#).
Cf. A068237 (numerator of derivative of 1/n).
Cf. A068238 (denominator of derivative of 1/n).
Cf. A068328 (derivative of squarefree numbers).
Cf. A068311 (derivative of n!).
Cf. A168386 (derivative of n!!).
Cf. A260619 (derivative of hyperfactorial(n)).
Cf. A260620 (derivative of superfactorial(n)).
Cf. A068312 (derivative of triangular numbers).
Cf. A068329 (derivative of Fibonacci(n)).
Cf. A096371 (derivative of partition number).
Cf. A099301 (derivative of d(n)).
Cf. A099310 (derivative of phi(n)).
Cf. A342925 (derivative of sigma(n)).
Cf. A349905 (derivative of prime shift).
Cf. A327860 (derivative of primorial base exp-function).
Cf. A369252 (derivative of products of three odd primes), A369251 (same sorted).
Cf. A068346 (second  derivative of n).
Cf. A099306 (third   derivative of n).
Cf. A258644 (fourth  derivative of n).
Cf. A258645 (fifth   derivative of n).
Cf. A258646 (sixth   derivative of n).
Cf. A258647 (seventh derivative of n).
Cf. A258648 (eighth  derivative of n).
Cf. A258649 (ninth   derivative of n).
Cf. A258650 (tenth   derivative of n).
Cf. A185232 (n-th    derivative of n).
Cf. A258651 (A(n,k) = k-th arithmetic derivative of n).
Cf. A085731 (gcd(n,n')), A083345 (n'/gcd(n,n')), A057521 (gcd(n, (n')^k) for k>1).
Cf. A342014 (n' mod n), A369049 (n mod n').
Cf. A341998 (A003557(n')), A342001 (n'/A003557(n)).
Cf. A098699 (least x such that x' = n, antiderivative of n).
Cf. A098700 (n such that x' = n has no integer solution).
Cf. A099302 (number of solutions to x' = n).
Cf. A099303 (greatest x such that x' = n).
Cf. A051674 (n such that n' = n).
Cf. A083347 (n such that n' < n).
Cf. A083348 (n such that n' > n).
Cf. A099304 (least k such that (n+k)' = n' + k').
Cf. A099305 (number of solutions to (n+k)' = n' + k').
Cf. A328235 (least k > 0 such that (n+k)' = u * n' for some natural number u).
Cf. A328236 (least m > 1 such that (m*n)' = u * n' for some natural number u).
Cf. A099307 (least k such that the k-th arithmetic derivative of n is zero).
Cf. A099308 (k-th arithmetic derivative of n is zero for some k).
Cf. A099309 (k-th arithmetic derivative of n is nonzero for all k).
Cf. A129150 (n-th derivative of 2^3).
Cf. A129151 (n-th derivative of 3^4).
Cf. A129152 (n-th derivative of 5^6).
Cf. A189481 (x' = n has a unique solution).
Cf. A190121 (partial sums).
Cf. A258057 (first differences).
Cf. A229501 (n divides the n-th partial sum).
Cf. A165560 (parity).
Cf. A235991 (n' is odd), A235992 (n' is even).
Cf. A327863, A327864, A327865 (n' is a multiple of 3, 4, 5).
Cf. A157037 (n' is prime), A192192 (n'' is prime), A328239 (n''' is prime).
Cf. A328393 (n' is squarefree), A328234 (squarefree and > 1).
Cf. A328244 (n'' is squarefree), A328246 (n''' is squarefree).
Cf. A328303 (n' is not squarefree), A328252 (n' is squarefree, but n is not).
Cf. A328248 (least k such that the (k-1)-th derivative of n is squarefree).
Cf. A328251 (k-th arithmetic derivative is never squarefree for any k >= 0).
Cf. A256750 (least k such that the k-th derivative is either 0 or has a factor p^p).
Cf. A327928 (number of distinct primes p such that p^p divides n').
Cf. A342003 (max. exponent k for any prime power p^k that divides n').
Cf. A327929 (n' has at least one divisor of the form p^p).
Cf. A327978 (n' is primorial number > 1).
Cf. A328243 (n' is a partial sum of primorial numbers and larger than one).
Cf. A328310 (maximal prime exponent of n' minus maximal prime exponent of n).
Cf. A328320 (max. prime exponent of n' is less than that of n).
Cf. A328321 (max. prime exponent of n' is >= that of n).
Cf. A328383 (least k such that the k-th derivative of n is either a multiple or a divisor of n, but not both).
Cf. A263111 (the ordinal transform of a).
Cf. A300251, A319684 (Möbius and inverse Möbius transform).
Cf. A305809 (Dirichlet convolution square).
Cf. A349133, A349173, A349394, A349380, A349618, A349619, A349620, A349621 (for miscellaneous Dirichlet convolutions).
Cf. A069359 (similar formula which agrees on squarefree numbers).
Cf. A258851 (the pi-based arithmetic derivative of n).
Cf. A328768, A328769 (primorial-based arithmetic derivatives of n).
Cf. A328845, A328846 (Fibonacci-based arithmetic derivatives of n).
Cf. A302055, A327963, A327965, A328099 (for other variants and modifications).
Cf. A038554 (another sequence using "derivative" in its name, but involving binary expansion of n).
Cf. A322582, A348507 (lower and upper bounds), also A002620.

A003415:= Concatenation([0,0],List(List([2..10^3],Factors),
i->Product(i)*Sum(i,j->1/j))); # Muniru A Asiru, Aug 31 2017
(APL, Dyalog dialect) A003415 ← { ⍺←(0 1 2) ⋄ ⍵≤1:⊃⍺ ⋄ 0=(3⊃⍺)|⍵:((⊃⍺+(2⊃⍺)×(⍵÷3⊃⍺)) ((2⊃⍺)×(3⊃⍺)) (3⊃⍺)) ∇ ⍵÷3⊃⍺ ⋄ ((⊃⍺) (2⊃⍺) (1+(3⊃⍺))) ∇ ⍵} ⍝ Antti Karttunen, Feb 18 2024

a003415 0 = 0
a003415 n = ad n a000040_list where
  ad 1 _             = 0
  ad n ps'@(p:ps)
     | n < p * p     = 1
     | r > 0         = ad n ps
     | otherwise     = n' + p * ad n' ps' where
       (n',r) = divMod n p
-- Reinhard Zumkeller, May 09 2011

Ad:=func; [n le 1 select 0 else Ad(n): n in [0..80]]; // Bruno Berselli, Oct 22 2013

A003415 := proc(n) local B,m,i,t1,t2,t3; B := 1000000000039; if n<=1 then RETURN(0); fi; if isprime(n) then RETURN(1); fi; t1 := ifactor(B*n); m := nops(t1); t2 := 0; for i from 1 to m do t3 := op(i,t1); if nops(t3) = 1 then t2 := t2+1/op(t3); else t2 := t2+op(2,t3)/op(op(1,t3)); fi od: t2 := t2-1/B; n*t2; end;
A003415 := proc(n)
        local a,f;
        a := 0 ;
        for f in ifactors(n)[2] do
                a := a+ op(2,f)/op(1,f);
        end do;
        n*a ;
end proc: # R. J. Mathar, Apr 05 2012

a[ n_] := If[ Abs @ n < 2, 0, n Total[ #2 / #1 & @@@ FactorInteger[ Abs @ n]]]; (* Michael Somos, Apr 12 2011 *)
dn[0] = 0; dn[1] = 0; dn[n_?Negative] := -dn[-n]; dn[n_] := Module[{f = Transpose[FactorInteger[n]]}, If[PrimeQ[n], 1, Total[n*f[[2]]/f[[1]]]]]; Table[dn[n], {n, 0, 100}] (* T. D. Noe, Sep 28 2012 *)

A003415(n) = {local(fac);if(n<1,0,fac=factor(n);sum(i=1,matsize(fac)[1],n*fac[i,2]/fac[i,1]))} /* Michael B. Porter, Nov 25 2009 */

apply( A003415(n)=vecsum([n/f[1]*f[2]|f<-factor(n+!n)~]), [0..99]) \\ M. F. Hasler, Sep 25 2013, updated Nov 27 2019

A003415(n) = { my(s=0, m=1, spf); while(n>1, spf = A020639(n); n /= spf; s += m*n; m *= spf); (s); }; \\ Antti Karttunen, Mar 10 2021

a(n) = my(f=factor(n), r=[1/(e+!e)|e<-f[,1]], c=f[,2]); n*r*c; \\ Ruud H.G. van Tol, Sep 03 2023

from sympy import factorint
def A003415(n):
    return sum([int(n*e/p) for p,e in factorint(n).items()]) if n > 1 else 0
# Chai Wah Wu, Aug 21 2014

def A003415(n):
    F = [] if n == 0 else factor(n)
    return n * sum(g / f for f, g in F)
[A003415(n) for n in range(79)] # Peter Luschny, Aug 23 2014

If n = Product p_i^e_i, a(n) = n * Sum (e_i/p_i).
a(m*p^p) = (m + a(m))*p^p, p prime: a(m*A051674(k))=A129283(m)*A051674(k). - Reinhard Zumkeller, Apr 07 2007
For n > 1: a(n) = a(A032742(n)) * A020639(n) + A032742(n). - Reinhard Zumkeller, May 09 2011
a(n) = n * Sum_{p|n} v_p(n)/p, where v_p(n) is the largest power of the prime p dividing n. - Wesley Ivan Hurt, Jul 12 2015
For n >= 2, Sum_{k=2..n} floor(1/a(k)) = pi(n) = A000720(n) (see K. T. Atanassov article). - Ivan N. Ianakiev, Mar 22 2019
From A.H.M. Smeets, Jan 17 2020: (Start)
Limit_{n -> oo} (1/n^2)*Sum_{i=1..n} a(i) = A136141/2.
Limit_{n -> oo} (1/n)*Sum_{i=1..n} a(i)/i = A136141.
a(n) = n if and only if n = p^p, where p is a prime number. (End)
Dirichlet g.f.: zeta(s-1)*Sum_{p prime} 1/(p^s-p), see A136141 (s=2), A369632 (s=3) [Haukkanen, Merikoski and Tossavainen]. - Sebastian Karlsson, Nov 25 2021
From Antti Karttunen, Nov 25 2021: (Start)
a(n) = Sum_{d|n} d * A349394(n/d).
For all n >= 1, A322582(n) <= a(n) <= A348507(n).
If n is not a prime, then a(n) >= 2*sqrt(n), or in other words, for all k >= 1 for which A002620(n)+k is not a prime, we have a(A002620(n)+k) > n. [See Ufnarovski and Åhlander, Theorem 9, point (3).]
(End)

More terms from Michel ten Voorde, Apr 11 2001

A051903 Maximum exponent in the prime factorization of n.

Original entry on oeis.org

0, 1, 1, 2, 1, 1, 1, 3, 2, 1, 1, 2, 1, 1, 1, 4, 1, 2, 1, 2, 1, 1, 1, 3, 2, 1, 3, 2, 1, 1, 1, 5, 1, 1, 1, 2, 1, 1, 1, 3, 1, 1, 1, 2, 2, 1, 1, 4, 2, 2, 1, 2, 1, 3, 1, 3, 1, 1, 1, 2, 1, 1, 2, 6, 1, 1, 1, 2, 1, 1, 1, 3, 1, 1, 2, 2, 1, 1, 1, 4, 4, 1, 1, 2, 1, 1, 1, 3, 1, 2, 1, 2, 1, 1, 1, 5, 1, 2, 2, 2, 1, 1, 1, 3, 1

Offset: 1

Labos Elemer, Dec 16 1999

Smallest number of factors of all factorizations of n into squarefree numbers, see also A128651, A001055. - Reinhard Zumkeller, Mar 30 2007
Maximum number of invariant factors among abelian groups of order n. - Álvar Ibeas, Nov 01 2014
a(n) is the highest of the frequencies of the parts of the partition having Heinz number n. We define the Heinz number of a partition p = [p_1, p_2, ..., p_r] as Product(p_j-th prime, j=1..r) (concept used by Alois P. Heinz in A215366 as an "encoding" of a partition). For example, for the partition [1, 1, 2, 4, 10] we get 2*2*3*7*29 = 2436. Example: a(24) = 3; indeed, the partition having Heinz number 24 = 2*2*2*3 is [1,1,1,2], where the distinct parts 1 and 2 have frequencies 3 and 1, respectively. - Emeric Deutsch, Jun 04 2015
From Thomas Ordowski, Dec 02 2019: (Start)
a(n) is the smallest k such that b^(phi(n)+k) == b^k (mod n) for all b.
The Euler phi function can be replaced by the Carmichael lambda function.
Problems:
(*) Are there composite numbers n > 4 such that n == a(n) (mod phi(n))? By Lehmer's totient conjecture, there are no such squarefree numbers.
(**) Are there odd numbers n such that a(n) > 1 and n == a(n) (mod lambda(n))? These are odd numbers n such that a(n) > 1 and b^n == b^a(n) (mod n) for all b.
(***) Are there odd numbers n such that a(n) > 1 and n == a(n) (mod ord_{n}(2))? These are odd numbers n such that a(n) > 1 and 2^n == 2^a(n) (mod n).
Note: if (***) do not exist, then (**) do not exist. (End)
Niven (1969) proved that the asymptotic mean of this sequence is 1 + Sum_{j>=2} 1 - (1/zeta(j)) (A033150). - Amiram Eldar, Jul 10 2020

			For n = 72 = 2^3*3^2, a(72) = max(exponents) = max(3,2) = 3.

Daniel Forgues, Table of n, a(n) for n = 1..100000 (first 10000 terms from T. D. Noe)
Benjamin Merlin Bumpus and Zoltan A. Kocsis, Spined categories: generalizing tree-width beyond graphs, arXiv:2104.01841 [math.CO], 2021.
Cao Hui-Zhong, The Asymptotic Formulas Related to Exponents in Factoring Integers, Math. Balkanica, Vol. 5 (1991), Fasc. 2.
Ivan Niven, Averages of Exponents in Factoring Integers, Proc. Amer. Math. Soc., Vol. 22, No. 2 (1969), pp. 356-360.
Eric Weisstein's World of Mathematics, Niven's Constant.
Index entries for sequences computed from exponents in factorization of n

Average value is A033150 = 1.7052....

Cf. A002322, A005361, A008479, A028234, A051904, A052409, A067029, A091050, A129132, A327295, A328310, A329885.

a051903 1 = 0
a051903 n = maximum $ a124010_row n -- Reinhard Zumkeller, May 27 2012

A051903 := proc(n)
        a := 0 ;
        for f in ifactors(n)[2] do
                a := max(a,op(2,f)) ;
        end do:
        a ;
end proc: # R. J. Mathar, Apr 03 2012
# second Maple program:
a:= n-> max(0, seq(i[2], i=ifactors(n)[2])):
seq(a(n), n=1..120);  # Alois P. Heinz, May 09 2020

Table[If[n == 1, 0, Max @@ Last /@ FactorInteger[n]], {n, 100}] (* Ray Chandler, Jan 24 2006 *)

a(n)=if(n>1,vecmax(factor(n)[,2]),0) \\ Charles R Greathouse IV, Oct 30 2012

from sympy import factorint
def A051903(n):
    return max(factorint(n).values()) if n > 1 else 0
# Chai Wah Wu, Jan 03 2015

;; With memoization-macro definec.
(definec (A051903 n) (if (= 1 n) 0 (max (A067029 n) (A051903 (A028234 n))))) ;; Antti Karttunen, Aug 08 2016

a(n) = max_{k=1..A001221(n)} A124010(n,k). - Reinhard Zumkeller, Aug 27 2011
a(1) = 0; for n > 1, a(n) = max(A067029(n), a(A028234(n))). - Antti Karttunen, Aug 08 2016
Conjecture: a(n) = a(A003557(n)) + 1. This relation together with a(1) = 0 defines the sequence. - Velin Yanev, Sep 02 2017
Comment from David J. Seal, Sep 18 2017: (Start)
This conjecture seems very easily provable to me: if the factorization of n is p1^k1 * p2^k2 * ... * pm^km, then the factorization of the largest squarefree divisor of n is p1 * p2 * ... * pm. So the factorization of A003557(n) is p1^(k1-1) * p2^(k2-1) * ... * pm^(km-1) if exponents of zero are allowed, or with the product terms that have an exponent of zero removed if they're not (if that results in an empty product, consider it to be 1 as usual).
The formula then follows from the fact that provided all ki >= 1, Max(k1, k2, ..., km) = Max(k1-1, k2-1, ..., km-1) + 1, and Max(k1-1, k2-1, ..., km-1) is not altered by removing the ki-1 values that are 0, provided we treat the empty Max() as being 0. That proves the formula and the provisos about empty products and Max() correspond to a(1) = 0.
Also, for any n, applying the formula Max(k1, k2, ..., km) times to n = p1^k1 * p2^k2 * ... * pm^km reduces all the exponents to zero, i.e., to the case a(1) = 0, so that case and the formula generate the sequence. (End)
Sum_{k=1..n} (-1)^k * a(k) ~ c * n, where c = Sum_{k>=2} 1/((2^k-1)*zeta(k)) = 0.44541445377638761933... . - Amiram Eldar, Jul 28 2024
a(n) <= log(n)/log(2). - Hal M. Switkay, Jul 03 2025

A328321 Numbers n for which A328311(n) = 1 + A051903(A003415(n)) - A051903(n) is strictly positive.

Original entry on oeis.org

4, 6, 10, 12, 14, 15, 16, 20, 21, 22, 26, 27, 28, 30, 33, 34, 35, 36, 38, 39, 42, 44, 46, 48, 50, 51, 52, 54, 55, 57, 58, 60, 62, 64, 65, 66, 68, 69, 70, 74, 76, 77, 78, 80, 82, 84, 85, 86, 87, 91, 92, 93, 94, 95, 99, 100, 102, 105, 106, 108, 110, 111, 112, 114, 115, 116, 118, 119, 122, 123, 124, 129, 130, 132, 133

Offset: 1

Antti Karttunen, Oct 13 2019

nonn

Numbers n for which A051903(A003415(n)) >= A051903(n), i.e., numbers such that taking their arithmetic derivative does not decrease their "degree", A051903, the maximal exponent in prime factorization.

			10 = 2*5 has maximal exponent (A051903) 1, and its arithmetic derivative A003415(10) = 2+5 = 7 also has maximal exponent 1, thus 10 is included in this sequence.
15 = 3*5 has maximal exponent 1, and its arithmetic derivative A003415(15) = 3+5 = 8 = 2^3 has maximal exponent 3, thus 15 is included in this sequence.
For 8 = 2^3, its arithmetic derivative A003415(8) = 12 = 2^2 * 3, and as 2 < 3 (highest exponent of 12 is less than that of 8), 8 is NOT included here, and from this we also see that A100716 is not a subsequence of this sequence.

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

Cf. A003415, A051903, A100716, A328302, A328310, A328311.

Cf. A328320 (complement), A051674, A157037, A328304, A328305 (subsequences).

A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A051903(n) = if((1==n),0,vecmax(factor(n)[, 2]));
A328311(n) = if(n<=1,0,1+(A051903(A003415(n)) - A051903(n)));
isA328321(n) = (A328311(n)>0);

A328320 Numbers for which A328311(n) = 1 + A051903(A003415(n)) - A051903(n) is zero (including 1 as the initial term).

Original entry on oeis.org

1, 2, 3, 5, 7, 8, 9, 11, 13, 17, 18, 19, 23, 24, 25, 29, 31, 32, 37, 40, 41, 43, 45, 47, 49, 53, 56, 59, 61, 63, 67, 71, 72, 73, 75, 79, 81, 83, 88, 89, 90, 96, 97, 98, 101, 103, 104, 107, 109, 113, 117, 120, 121, 125, 126, 127, 128, 131, 136, 137, 139, 147, 149, 150, 151, 152, 153, 157, 160, 162, 163, 167, 168, 169

Offset: 1

Antti Karttunen, Oct 13 2019

nonn

After 1, the numbers whose "degree" (maximal exponent, A051903) is decremented by one when arithmetic derivative (A003415) is applied to them.

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

Cf. A003415, A051903, A328302, A328310.
Indices of zeros in A328311.
Cf. A328321 (complement), A328252 (a subsequence).

A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A051903(n) = if((1==n),0,vecmax(factor(n)[, 2]));
A328311(n) = if(n<=1,0,1+(A051903(A003415(n)) - A051903(n)));
isA328320(n) = (0==A328311(n));

A328311 a(n) = 1 + A051903(A003415(n)) - A051903(n), a(1) = 0 by convention.

Original entry on oeis.org

0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 3, 0, 2, 3, 2, 0, 0, 0, 2, 1, 1, 0, 0, 0, 1, 1, 4, 0, 1, 0, 0, 1, 1, 2, 1, 0, 1, 4, 0, 0, 1, 0, 3, 0, 2, 0, 1, 0, 1, 2, 2, 0, 2, 4, 0, 1, 1, 0, 1, 0, 1, 0, 1, 2, 1, 0, 2, 1, 1, 0, 0, 0, 1, 0, 3, 2, 1, 0, 1, 0, 1, 0, 1, 1, 2, 5, 0, 0, 0, 2, 4, 1, 2, 3, 0, 0, 0, 1, 1, 0, 1, 0, 0, 1

Offset: 1

Antti Karttunen, Oct 13 2019

nonn

All terms are nonnegative because taking the arithmetic derivative (A003415) of n may decrease its "degree" (i.e., its maximal exponent, A051903) by at most one, and in many cases may also increase it, or keep it same.

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

One more than A328310.
Cf. A003415, A051903, A328312.
Cf. A328320 (indices of zeros), A328321 (of nonzero terms).

A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A051903(n) = if((1==n),0,vecmax(factor(n)[, 2]));
A328311(n) = if(n<=1,0,1+(A051903(A003415(n)) - A051903(n)));

a(1) = 0, for n > 1, a(n) = 1 + A051903(A003415(n)) - A051903(n).

For n > 1, a(n) = 1 + A328310(n).

A328384 If n is of the form p^p, a(n) = 0, otherwise a(n) gives the number of iterations of x -> A003415(x) needed to reach the first number different from n which is either a prime, or whose degree (A051903) differs from the degree of n. If the degree of the final number is <= that of n, then a(n) = -1 * iteration count.

Original entry on oeis.org

-1, -1, -1, 0, -1, -1, -1, -1, -1, -1, -1, 1, -1, 1, 1, 1, -1, -1, -1, 1, -2, -1, -1, -1, -1, 2, 0, 1, -1, -1, -1, -1, 2, -1, 1, 3, -1, -3, 1, -1, -1, -1, -1, 1, -1, 1, -1, 3, -1, -2, 1, 1, -1, 1, 1, -1, -2, -1, -1, 2, -1, 3, -1, 2, 1, -1, -1, 1, 3, -1, -1, -1, -1, 2, -1, 1, 1, -1, -1, 5, -1, -1, -1, 2, -2, 1, 1, -1, -1, -1, 1, 1, -2, 1

Offset: 1

Antti Karttunen, Oct 15 2019

sign

The records -1, 0, 1, 2, 3, 5, 8, 10, 11, 13, ... occur at n = 1, 4, 12, 26, 36, 80, 108, 4887, 18688, 22384, ...

			For n = 21 = 3*7, A051903(21) = 1, A003415(21) = 10 = 2*5, is of the same degree as A051903(10) = 1, but then A003415(10) = 7, which is a prime, having degree <= of the starting value (as we have A051903(7) <= A051903(21)), thus a(21) = -1 * 2 = -2.
For n = 33 = 3*11, A051903(33) = 1, A003415(33) = 14 = 2*7, is of the same degree, but on the second iteration, A003415(14) = 9 = 3^2, with A051903(9) = 2, which is larger than the initial degree, thus a(33) = +2.

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

Cf. A003415, A051903, A051674, A157037, A258651, A328252, A328310, A328320.
Cf. A328385 (the number found in the iteration).
Cf. A256750, A328248, A328383 for similar counts.

A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A051903(n) = if((n<=1),n-1,vecmax(factor(n)[, 2]));
A328384(n) = { my(d=A051903(n), u=A003415(n), k=1); if(u==n,return(0)); while(u && (u!=n) && !isprime(u) && A051903(u)==d, k++; n = u; u = A003415(u)); if(A051903(u)<=d,-k,k); };

a(1) = -1 as 0 is here considered having a smaller degree than 1.
a(p) = -1 for all primes.
a(A051674(n)) = 0.
a(A157037(n)) = -1.
a(A328252(n)) = -1.
a(A328320(n)) = -1.

A328385 If n is of the form p^p, a(n) = n, otherwise a(n) is the first number found by iterating the map x -> A003415(x) that is different from n and either a prime, or whose degree (A051903) differs from the degree of n.

Original entry on oeis.org

0, 1, 1, 4, 1, 5, 1, 12, 6, 7, 1, 16, 1, 9, 8, 32, 1, 21, 1, 24, 7, 13, 1, 44, 10, 8, 27, 32, 1, 31, 1, 80, 9, 19, 12, 96, 1, 7, 16, 68, 1, 41, 1, 48, 39, 25, 1, 608, 14, 39, 20, 56, 1, 81, 16, 92, 13, 31, 1, 96, 1, 9, 51, 640, 18, 61, 1, 72, 8, 59, 1, 156, 1, 16, 55, 80, 18, 71, 1, 3424, 108, 43, 1, 128, 13, 45, 32, 140, 1, 123, 20, 96, 19

Offset: 1

Antti Karttunen, Oct 14 2019

nonn

			For n = 3, 3 is a prime, thus a(3) = 1.
For n = 4, A003415(4) = 4, thus as it is among the fixed points of A003415 and a(4) = 4.
For n = 8 = 2^3, its "degree" is A051903(33) = 3, but A003415(8) = 12 = 2^2 * 3, with degree 2, thus a(8) = 12.
For n = 21 = 3*7, A051903(21) = 1, the first derivative A003415(21) = 10 = 2*5 is of the same degree as A051903(10) = 1, but then continuing, we have A003415(10) = 7, which is a prime, thus a(21) = 7.
For n = 33 = 3*11, A051903(33) = 1, A003415(33) = 14 = 2*7, is of the same degree, but on the second iteration, A003415(14) = 9 = 3^2, with A051903(9) = 2, different from the initial degree, thus a(33) = 9.

Antti Karttunen, Table of n, a(n) for n = 1..16384
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..65537

Cf. A003415, A051674, A051903, A157037.
Cf. A328384 (the number of iterations needed to reach such a number).
Cf. also A327968, A327975, A327977, A328252, A328305, A328310, A328311, A328320, A328321.

A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A051903(n) = if((1==n),0,vecmax(factor(n)[, 2]));
A328385(n) = { my(d=A051903(n), u=A003415(n)); while(u && (u!=n) && !isprime(u) && A051903(u)==d, n = u; u = A003415(u)); (u); };

a(1) = 0 [as here the degrees of 0 and 1 are considered different].
a(p) = 1 for all primes.
a(A051674(n)) = A051674(n).
a(A157037(n)) = A003415(A157037(n)), a prime.
a(A328252(n)) = A003415(A328252(n)), a squarefree number.
a(n) = A003415^(k)(n), when k = abs(A328384(n)). [Taking the abs(A328384(n))-th arithmetic derivative of n gives a(n)]

A342003 Maximal exponent in the prime factorization of the arithmetic derivative of n: a(n) = A051903(A003415(n)).

Original entry on oeis.org

0, 0, 2, 0, 1, 0, 2, 1, 1, 0, 4, 0, 2, 3, 5, 0, 1, 0, 3, 1, 1, 0, 2, 1, 1, 3, 5, 0, 1, 0, 4, 1, 1, 2, 2, 0, 1, 4, 2, 0, 1, 0, 4, 1, 2, 0, 4, 1, 2, 2, 3, 0, 4, 4, 2, 1, 1, 0, 2, 0, 1, 1, 6, 2, 1, 0, 3, 1, 1, 0, 2, 0, 1, 1, 4, 2, 1, 0, 4, 3, 1, 0, 2, 1, 2, 5, 2, 0, 1, 2, 5, 1, 2, 3, 4, 0, 1, 2, 2, 0, 1, 0, 2, 1

Offset: 2

Antti Karttunen, Mar 01 2021

nonn

Antti Karttunen, Table of n, a(n) for n = 2..65537

Cf. A003415, A051903, A328310, A328311, A328320, A328321, A341994, A341995, A342004.

Cf. A000040 (indices of zeros), A328234 (of ones), A328393 (of the terms < 2).

A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A051903(n) = if((1==n),0,vecmax(factor(n)[, 2]));
A342003(n) = A051903(A003415(n));

a(n) = A051903(A003415(n)).
a(n) = A051903(n) + A328310(n).
a(n) = 1 iff A341994(n) = 1.

A342026 Difference in the maximal prime exponents between the arithmetic derivative of A276086(n) and A276086(n) itself, which is the prime product form of primorial base expansion of n.

Original entry on oeis.org

-1, -1, 0, -1, -1, -1, 0, 2, 0, -1, -1, -1, 0, -1, -1, 2, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 1, 0, 0, -1, -1, 1, 0, 0, 0, -1, -1, -1, -1, 0, 0, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 0, -1, -1, -1, -1, -1, 0, -1, 1, -1, -1, -1, -1, -1, -1, -1, 4, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 1

Offset: 1

Antti Karttunen, Mar 13 2021

sign

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

Cf. A003415, A276086, A342005, A327860, A328114, A328310, A328391, A342006 (positions of nonnegative terms).

Cf. also A342016, A342019.

A051903(n) = if((1==n), 0, vecmax(factor(n)[, 2]));
A327860(n) = { my(s=0, m=1, p=2, e); while(n, e = (n%p); m *= (p^e); s += (e/p); n = n\p; p = nextprime(1+p)); (s*m); };
A328391(n) = if(!n,n,A051903(A327860(n)));
A328114(n) = { my(s=0, p=2); while(n, s = max(s, (n%p)); n = n\p; p = nextprime(1+p)); (s); };
A342026(n) = (A328391(n) - A328114(n));

a(n) = A328310(A276086(n)) = A328391(n) - A328114(n).

a(n) = -1 iff A342005(n) = 1.

cp's OEIS Frontend

A003415 a(n) = n' = arithmetic derivative of n: a(0) = a(1) = 0, a(prime) = 1, a(m*n) = m*a(n) + n*a(m).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

PARI

Python

Sage

Formula

Extensions

A051903 Maximum exponent in the prime factorization of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Scheme

Formula

A328321 Numbers n for which A328311(n) = 1 + A051903(A003415(n)) - A051903(n) is strictly positive.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

A328320 Numbers for which A328311(n) = 1 + A051903(A003415(n)) - A051903(n) is zero (including 1 as the initial term).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

A328311 a(n) = 1 + A051903(A003415(n)) - A051903(n), a(1) = 0 by convention.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Formula

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

Formula

A328385 If n is of the form p^p, a(n) = n, otherwise a(n) is the first number found by iterating the map x -> A003415(x) that is different from n and either a prime, or whose degree (A051903) differs from the degree of n.

A003415 a(n) = n' = arithmetic derivative of n: a(0) = a(1) = 0, a(prime) = 1, a(mn) = ma(n) + n*a(m).