A129152 -id:A129152 - cp's OEIS frontend

A129286 a(n) = A129152(n) / 5^5, where A129152 is the trajectory of 5^6 under A003415, the arithmetic derivative.

5, 6, 11, 12, 28, 60, 152, 388, 780, 2036, 4076, 8156, 16316, 32636, 66232, 169612, 339228, 1244808, 4856004, 14568120, 56648484, 190791072, 866874960, 4084880112, 16342519392, 73816544592, 305811402048, 1732931293824, 10145637246528, 52715820454848

Offset: 0

Reinhard Zumkeller, Apr 07 2007

nonn

Cf. A129284, A129285, A051674.

a129286 n = a129286_list !! n
a129286_list = iterate a129283 5  -- Reinhard Zumkeller, Nov 01 2013

a(n+1) = A129283(a(n)), a(0) = 5.

a(n) = A090635(n+3) / 4, i.e., the trajectory of 20 under A003415, divided by the common factor 4. - M. F. Hasler, Nov 27 2019

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

Original entry on oeis.org

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

A051674 a(n) = prime(n)^prime(n).

Original entry on oeis.org

4, 27, 3125, 823543, 285311670611, 302875106592253, 827240261886336764177, 1978419655660313589123979, 20880467999847912034355032910567, 2567686153161211134561828214731016126483469, 17069174130723235958610643029059314756044734431

Offset: 1

Asher Auel

Numbers k such that bigomega(k)^(bigomega(k)) = k, where bigomega = A001222. - Lekraj Beedassy, Aug 21 2004
Positive k such that k' = k, where k' is the arithmetic derivative of k. - T. D. Noe, Oct 12 2004
David Beckwith proposes (in the AMM reference): "Let n be a positive integer and let p be a prime number. Prove that (p^p) | n! implies that (p^(p + 1)) | n!". - Jonathan Vos Post, Feb 20 2006
Subsequence of A100716; A003415(m*a(n)) = A129283(m)*a(n), especially A003415(a(n)) = a(n). - Reinhard Zumkeller, Apr 07 2007
A168036(a(n)) = 0. - Reinhard Zumkeller, May 22 2015

			a(1) = 2^2 = 4.
a(2) = 3^3 = 27.
a(3) = 5^5 = 3125.

J.-M. De Koninck & A. Mercier, 1001 Problemes en Theorie Classique Des Nombres, Problem 740 pp. 95; 312, Ellipses Paris 2004.

T. D. Noe, Table of n, a(n) for n = 1..40
David Beckwith, Problem 11158, American Mathematical Monthly, Vol. 112, No. 5 (May 2005), p. 468.
Jurij Kovic, The Arithmetic Derivative and Antiderivative, Journal of Integer Sequences, Vol. 15 (2012), #12.3.8.

Cf. A000040, A000312, A003415 (arithmetic derivative of n), A129150, A129151, A129152, A048102, A072873 (multiplicative closure), A104126.
Subsequence of A100717; A203908(a(n)) = 0.
Subsequence of A097764.
Cf. A168036, A094289 (decimal expansion of Sum(1/p^p)).

a051674_list = map (\p -> p ^ p) a000040_list
-- Reinhard Zumkeller, Jan 21 2012

[p^p: p in PrimesUpTo(30)]; // Vincenzo Librandi, Mar 27 2014

A051674:=n->ithprime(n)^ithprime(n): seq(A051674(n), n=1..10); # Wesley Ivan Hurt, Jun 25 2016

Array[Prime[ # ]^Prime[ # ] &, 12] (* Vladimir Joseph Stephan Orlovsky, May 01 2008 *)
#^#&/@Prime[Range[10]] (* Harvey P. Dale, May 17 2024 *)

a(n)=n=prime(n);n^n \\ Charles R Greathouse IV, Mar 20 2013

from gmpy2 import mpz
[mpz(prime(n))**mpz(prime(n)) for n in range(1,100)] # Chai Wah Wu, Jul 28 2014

a(n) = A000312(A000040(n)). - Altug Alkan, Sep 01 2016

Sum_{n>=1} 1/a(n) = A094289. - Amiram Eldar, Oct 13 2020

A100716 Numbers k such that p^p divides k for some prime p.

Original entry on oeis.org

4, 8, 12, 16, 20, 24, 27, 28, 32, 36, 40, 44, 48, 52, 54, 56, 60, 64, 68, 72, 76, 80, 81, 84, 88, 92, 96, 100, 104, 108, 112, 116, 120, 124, 128, 132, 135, 136, 140, 144, 148, 152, 156, 160, 162, 164, 168, 172, 176, 180, 184, 188, 189, 192, 196, 200, 204, 208, 212

Offset: 1

Leroy Quet, Dec 10 2004

nonn

Complement of A048103; A129251(a(n)) > 0; A051674 is a subsequence; A129254 = (terms a(k) such that a(k+1)=a(k)+1). - Reinhard Zumkeller, Apr 07 2007

A027748(a(n),k) <= A124010(a(n),k) for some k<=A001221(a(n)). - Reinhard Zumkeller, Apr 28 2012

			54 is included because 3^3 divides 54.

Reinhard Zumkeller, Table of n, a(n) for n = 1..1000

Complement: A048103.
Positions of nonzeros in A129251.
Cf. A100717, A129150, A129152.
Cf. A054744.
Cf. A051674 (a subsequence).
Subsequence of A276079 from which it differs for the first time at n=175, where a(175) = 628, while A276079(175) = 625, a value missing from here.

a100716 n = a100716_list !! (n-1)
a100716_list = filter (\x -> or $
   zipWith (<=) (a027748_row x) (map toInteger $ a124010_row x)) [1..]
-- Reinhard Zumkeller, Apr 28 2012
(Scheme, with Antti Karttunen's IntSeq-library)
(define A100716 (NONZERO-POS 1 1 A129251))
;; Antti Karttunen, Aug 18 2016

fQ[n_] := Union[ Table[ #[[1]] <= #[[2]]] & /@ FactorInteger[n]][[ -1]] == True; Select[ Range[2, 215], fQ[ # ] &] (* Robert G. Wilson v, Dec 14 2004 *)
f[n_] := Module[{aux=FactorInteger[n]}, Last@Union@Table[aux[[i,1]] <=  aux[[i,2]], {i,Length[aux]}] == True]; Select[Range[2,215], f] (* José María Grau Ribas, Jan 25 2012 *)
Rest@ Select[Range@ 216, Times @@ Boole@ Map[First@ # > Last@ # &, FactorInteger@ #] == 0 &] (* Michael De Vlieger, Aug 19 2016 *)

is(n)=forprime(p=2,default(primelimit),if(n%p^p==0,return(1));if(p^p>n,return(0))) \\ Charles R Greathouse IV, Jan 24 2012

from itertools import count, islice
from sympy import factorint
def A100716_gen(startvalue=1): # generator of terms >= startvalue
    return filter(lambda n:any(map(lambda d:d[1]>=d[0],factorint(n).items())),count(max(startvalue,1)))
A100716_list = list(islice(A100716_gen(),30)) # Chai Wah Wu, Jan 05 2023

a(n) ~ k*n with k = 1/(1 - Product(1 - p^-p)) = 3.5969959469... where the product is over all primes p. - Charles R Greathouse IV, Jan 24 2012

More terms from Robert G. Wilson v, Dec 14 2004

A258651 A(n,k) = n^(k) = k-th arithmetic derivative of n; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

0, 0, 1, 0, 0, 2, 0, 0, 1, 3, 0, 0, 0, 1, 4, 0, 0, 0, 0, 4, 5, 0, 0, 0, 0, 4, 1, 6, 0, 0, 0, 0, 4, 0, 5, 7, 0, 0, 0, 0, 4, 0, 1, 1, 8, 0, 0, 0, 0, 4, 0, 0, 0, 12, 9, 0, 0, 0, 0, 4, 0, 0, 0, 16, 6, 10, 0, 0, 0, 0, 4, 0, 0, 0, 32, 5, 7, 11, 0, 0, 0, 0, 4, 0, 0, 0, 80, 1, 1, 1, 12

Offset: 0

Alois P. Heinz, Jun 06 2015

			Square array A(n,k) begins:
  0,  0,  0,  0,  0,   0,   0,   0,    0,    0, ...
  1,  0,  0,  0,  0,   0,   0,   0,    0,    0, ...
  2,  1,  0,  0,  0,   0,   0,   0,    0,    0, ...
  3,  1,  0,  0,  0,   0,   0,   0,    0,    0, ...
  4,  4,  4,  4,  4,   4,   4,   4,    4,    4, ...
  5,  1,  0,  0,  0,   0,   0,   0,    0,    0, ...
  6,  5,  1,  0,  0,   0,   0,   0,    0,    0, ...
  7,  1,  0,  0,  0,   0,   0,   0,    0,    0, ...
  8, 12, 16, 32, 80, 176, 368, 752, 1520, 3424, ...
  9,  6,  5,  1,  0,   0,   0,   0,    0,    0, ...

Alois P. Heinz, Antidiagonals n = 0..115, flattened
J. Kovič, The Arithmetic Derivative and Antiderivative, Journal of Integer Sequences 15 (2012), Article 12.3.8
Wikipedia, Arithmetic derivative

Columns k=0-10 give: A001477, A003415, A068346, A099306, A258644, A258645, A258646, A258647, A258648, A258649, A258650.
Rows n=0,1,4,8 give: A000004, A000007, A010709, A129150.
Row 15: A090636, Row 28: A090637, Row 63: A090635, Row 81: A129151, Row 128: A369638, Row 1024: A214800, Row 15625: A129152.
Main diagonal gives A185232.
Antidiagonal sums give A258652.
Cf. also A328383.

d:= n-> n*add(i[2]/i[1], i=ifactors(n)[2]):
A:= proc(n, k) option remember; `if`(k=0, n, d(A(n, k-1))) end:
seq(seq(A(n, h-n), n=0..h), h=0..14);

d[n_] := n*Sum[i[[2]]/i[[1]], {i, FactorInteger[n]}]; d[0] = d[1] = 0;
A[n_, k_] := A[n, k] = If[k == 0, n, d[A[n, k-1]]];
Table[A[n, h-n], {h, 0, 14}, {n, 0, h}] // Flatten (* Jean-François Alcover, Apr 27 2017, translated from Maple *)

A(n,k) = A003415^k(n).

A129150 The n-th arithmetic derivative of 2^3.

Original entry on oeis.org

8, 12, 16, 32, 80, 176, 368, 752, 1520, 3424, 8592, 20096, 70464, 235072, 705280, 3023616, 13223680, 55540736, 278539264, 1392697344, 9541095424, 58609614848, 410267320320, 3397142953984, 24143851798528, 176071227916288, 1232666139967488, 9523075842834432

Offset: 0

Reinhard Zumkeller, Apr 01 2007

nonn

Conjecture: a strictly increasing sequence. - J. Lowell, Sep 10 2008
The sequence is strictly increasing because (4*n)' = 4*n + 4*n'. - David Radcliffe, Aug 19 2014
8 is the smallest integer that has a nontrivial trajectory (not going to 0 nor reduced to a fixed point as 4) under A003415, but 15 = A090636(1) has 8 as second term in its trajectory. 20 is the next larger such integer with a distinct trajectory, but has two larger predecessors, cf. A090635. - M. F. Hasler, Nov 27 2019
In general, the trajectory of p^(p+1) under A003415 has a common factor p^p, and divided by p^p it gives the trajectory of p under A129283: n -> n + n'. Here we have the case p = 2 (see A129284 for a(n)/2^2), cf. A129151 and A129152 for p = 3 and 5. - M. F. Hasler, Nov 28 2019

Charles R Greathouse IV, Table of n, a(n) for n = 0..100

Cf. A129151, A129152, A068327, A099309, A051674, A100716, A129284.

Row n = 8 of A258651.

a129150 n = a129150_list !! n
a129150_list = iterate a003415 8  -- Reinhard Zumkeller, Apr 29 2012

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]]]]]; s = 2^3; Join[{s}, Table[s = dn[s], {28}]] (* T. D. Noe, Mar 07 2013 *)

A129150(n,a=8)={if(n<0, vector(-n,n, if(n>1, a=A003415(a), a)), for(n=1,n, a=A003415(a)); a)}  \\ For n<0 return the vector a[0..-n-1]. - M. F. Hasler, Nov 27 2019

a(n+1) = A003415(a(n)), a(0) = 2^3 = 8.
a(n) = A090636(n+2).
A129251(a(n)) > 0. - Reinhard Zumkeller, Apr 07 2007
a(n) = 4*A129284(n). - M. F. Hasler, Nov 27 2019

a(21)-a(27) from Paolo P. Lava, Apr 16 2012

A129151 The n-th arithmetic derivative of 3^4.

Original entry on oeis.org

81, 108, 216, 540, 1188, 2484, 5076, 10260, 23112, 57996, 135648, 475632, 1586736, 4760640, 20409408, 89259840, 374899968, 1880140032, 9400707072, 64402394112, 395614900224, 2769304412160, 22930714939392, 162970999640064, 1188480788434944, 8320496444780544

Offset: 0

Reinhard Zumkeller, Apr 01 2007

nonn

In general, the trajectory of p^(p+1) under A003415 is equal to p^p times the trajectory of p under A129283: n -> n + n'. Here we have the case p = 3 (see A129285 for a(n)/3^3), see A129150 and A129152 for p = 2 and 5. - M. F. Hasler, Nov 28 2019

Cf. A129150, A129152, A068327.

Cf. A099309, A051674, A100716.

a129151 n = a129151_list !! n
a129151_list = iterate a003415 81  -- Reinhard Zumkeller, Apr 29 2012

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]]]]]; s = 3^4; Join[{s}, Table[s = dn[s], {25}]] (* T. D. Noe, Mar 07 2013 *)

a(n+1) = A003415(a(n)), a(0) = 3^4 = 81.

a(n) = A129285(n)*3^3; A129251(a(n)) > 0. - Reinhard Zumkeller, Apr 07 2007

A129284 a(n) = A129150(n) / 4, where A129150(n) = n-th arithmetic derivative of 2^3.

Original entry on oeis.org

2, 3, 4, 8, 20, 44, 92, 188, 380, 856, 2148, 5024, 17616, 58768, 176320, 755904, 3305920, 13885184, 69634816, 348174336, 2385273856, 14652403712, 102566830080, 849285738496, 6035962949632, 44017806979072, 308166534991872, 2380768960708608, 23410894780694528

Offset: 0

Reinhard Zumkeller, Apr 07 2007

nonn

In general, the trajectory of p^(p+1) under A003415 has a common factor p^p, and divided by p^p it gives the trajectory of p under A129283: n -> n + n'. Here we have the case p = 2, see A129151 and A129152 for p = 3 and 5. - M. F. Hasler, Nov 28 2019

Paolo P. Lava, Table of n, a(n) for n = 0..75

Cf. A129285, A129286, A051674.

a129284 n = a129150 n `div` 4  -- Reinhard Zumkeller, Nov 01 2013, corrected by M. F. Hasler, Nov 29 2019

A129284_upto(n)=A129150_upto(n)\4 \\ M. F. Hasler, Nov 29 2019

a(n+1) = A129283(a(n)), a(0) = 2.

a(18)-a(28) from Paolo P. Lava, Apr 16 2012

Edited by M. F. Hasler, Nov 27 2019

A214800 The n-th arithmetic derivative of 4^5.

Original entry on oeis.org

1024, 5120, 26624, 148480, 777216, 4249600, 25123840, 134268928, 741121024, 3811486720, 19819736064, 112733967360, 644174359552, 3220871798784, 22009290649600, 142032743280640, 881678229782528, 6174273912848384, 44101956520460288, 312415523110846464

Offset: 0

Franz Tolosa, Mar 07 2013

nonn

See A003415 for the definition of the arithmetic derivative.

Cf. A003415, A129150, A129151, A129152.

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]]]]]; s = 4^5; Join[{s}, Table[s = dn[s], {19}]] (* T. D. Noe, Mar 07 2013 *)

cp's OEIS Frontend

A129286 a(n) = A129152(n) / 5^5, where A129152 is the trajectory of 5^6 under A003415, the arithmetic derivative.

Views

Author

Keywords

Crossrefs

Programs

Haskell

Formula

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

A051674 a(n) = prime(n)^prime(n).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Python

Formula

A100716 Numbers k such that p^p divides k for some prime p.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Python

Formula

Extensions

A258651 A(n,k) = n^(k) = k-th arithmetic derivative of n; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A129150 The n-th arithmetic derivative of 2^3.

Views

Author

Keywords

Comments

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