A327934 Numbers k such that there is no prime p such that p^p divides k, but for its arithmetic derivative a positive finite number of such primes exist.
15, 35, 39, 51, 55, 87, 91, 95, 111, 115, 119, 123, 143, 155, 158, 159, 183, 187, 203, 215, 219, 225, 235, 247, 259, 267, 275, 287, 291, 295, 299, 303, 319, 323, 327, 329, 335, 339, 355, 371, 374, 391, 395, 403, 407, 410, 411, 415, 427, 441, 447, 451, 471, 473, 482, 511, 515, 519, 525, 527, 533, 535, 543, 551, 559, 579
Offset: 1
Keywords
Examples
From _Antti Karttunen_, Feb 01 2023: (Start) 1 has arithmetic derivative 1' = A003415(1) = 0. There are an infinite number of primes p such that p^p | 0, but because infinity is not a finite number, 1 is excluded from this sequence. 275 = 5^2 * 11 has no divisor of the form p^p, but its arithmetic derivative 275' = 135 = 3^3 * 5 has one divisor of the form p^p, therefore 275 is included in this sequence. 4125 = 3 * 5^3 * 11 has arithmetic derivative 4125' = A003415(4125) = 4225 = 5^2 * 13^2, that has no divisor of the form p^p, therefore 4125 is not included in this sequence. (End)
Links
- Antti Karttunen, Table of n, a(n) for n = 1..20000
Crossrefs
Programs
-
PARI
A003415(n) = {my(fac); if(n<1, 0, fac=factor(n); sum(i=1, matsize(fac)[1], n*fac[i, 2]/fac[i, 1]))}; \\ From A003415 A129251(n) = { my(f = factor(n)); sum(k=1, #f~, (f[k, 2]>=f[k, 1])); }; A327928(n) = if(n<=1,0,A129251(A003415(n))); isA327934(n) = ((0==A129251(n))&&(A327928(n)>0));
-
PARI
isA327934(n) = A360111(n); \\ Antti Karttunen, Feb 01 2023
Extensions
Definition clarified to exclude the term 1 by Antti Karttunen, Feb 01 2023
Comments