A342463 -id:A342463 - cp's OEIS frontend

A342001 Arithmetic derivative without its inherited divisor; the arithmetic derivative of n divided by A003557(n), which is a common divisor of both n and A003415(n).

0, 1, 1, 2, 1, 5, 1, 3, 2, 7, 1, 8, 1, 9, 8, 4, 1, 7, 1, 12, 10, 13, 1, 11, 2, 15, 3, 16, 1, 31, 1, 5, 14, 19, 12, 10, 1, 21, 16, 17, 1, 41, 1, 24, 13, 25, 1, 14, 2, 9, 20, 28, 1, 9, 16, 23, 22, 31, 1, 46, 1, 33, 17, 6, 18, 61, 1, 36, 26, 59, 1, 13, 1, 39, 11, 40, 18, 71, 1, 22, 4, 43, 1, 62, 22, 45, 32, 35, 1, 41, 20

Offset: 1

Antti Karttunen, Feb 28 2021

nonn

See also the scatter plot of A342002 that seems to reveal some interesting internal structure in this sequence, not fully explained by the regularity of primorial base expansion used in the latter sequence. - Antti Karttunen, May 09 2022

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, A003557, A341998, A342416, A342459.
Cf. A342002 [= a(A276086(n))], A342463 [= a(A342456(n))], A351945 [= a(A181819(n))], A353571 [= a(A003961(n))].
Cf. A346485 (Möbius transform), A347395 (convolution with Liouville's lambda), A347961 (with itself), and A347234, A347235, A347954, A347959, A347963, A349396, A349612 (for convolutions with other sequences).
Cf. A007947.

Array[#1/#2 & @@ {If[# < 2, 0, # Total[#2/#1 & @@@ FactorInteger[#]]] &@ Abs[#], #/Times @@ FactorInteger[#][[All, 1]]} &, 91] (* Michael De Vlieger, Mar 11 2021 *)

A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A003557(n) = (n/factorback(factorint(n)[, 1]));
A342001(n) = (A003415(n) / A003557(n));

from math import prod
from sympy import factorint
def A342001(n):
    q = prod(f:=factorint(n))
    return sum(q*e//p for p, e in f.items()) # Chai Wah Wu, Nov 04 2022

a(n) = A003415(n) / A003557(n).
For all n >= 0, a(A276086(n)) = A342002(n).
a(n) = A342414(n) * A342416(n) = A342459(n) * A342919(n). - Antti Karttunen, Apr 30 2022
Dirichlet g.f.: Dirichlet g.f. of A007947 * Sum_{p prime} p^s/((p^s-1)*(p^s+p-1)) = zeta(s) * Product_{p prime} (1+p^(1-s)-p^(-s)) * Sum_{p prime} p^s/((p^s-1)*(p^s+p-1)). - Sebastian Karlsson, May 05 2022
Sum_{k=1..n} a(k) ~ c * A065464 * Pi^2 * n^2 / 12, where c = Sum_{j>=2} (1/2 + (-1)^j * (Fibonacci(j) - 1/2))*PrimeZetaP(j) = 0.4526952873143153104685540856936425315834753528741817723313791528384... - Vaclav Kotesovec, May 09 2022

A342002 Čiurlionis sequence: Arithmetic derivative without its inherited divisor applied to the primorial base exp-function: a(n) = A342001(A276086(n)).

Original entry on oeis.org

0, 1, 1, 5, 2, 7, 1, 7, 8, 31, 13, 41, 2, 9, 11, 37, 16, 47, 3, 11, 14, 43, 19, 53, 4, 13, 17, 49, 22, 59, 1, 9, 10, 41, 17, 55, 12, 59, 71, 247, 106, 317, 19, 73, 92, 289, 127, 359, 26, 87, 113, 331, 148, 401, 33, 101, 134, 373, 169, 443, 2, 11, 13, 47, 20, 61, 17, 69, 86, 277, 121, 347, 24, 83, 107, 319, 142, 389, 31

Offset: 0

Antti Karttunen, Feb 28 2021

The scatter plot shows an interesting structure.
The terms are essentially the "wild" or "unherited" part of the arithmetic derivative (A003415) of those natural numbers (A048103) that are not immediately beyond all hope of reaching zero by iteration (as the terms of A100716 are), ordered by the primorial base expansion of n as in A276086. Sequence A342018 shows the positions of the terms here that have just moved to the "no hope" region, while A342019 shows how many prime powers in any term have breached the p^p limit. Note that the results are same as for A327860(n), as the division by "regular part", A328572(n) does not affect the "wild part" of the arithmetic derivative of A276086(n). - Antti Karttunen, Mar 12 2021
I decided to name this sequence in honor of Lithuanian artist Mikalojus Čiurlionis, 1875 - 1911, as the scatter plot of this sequence reminds me thematically of his work "Pyramid sonata", with similar elements: fractal repetition in different scales and high tension present, discharging as lightning. Like Čiurlionis's paintings, this sequence has many variations, see the Formula and Crossrefs sections. - Antti Karttunen, Apr 30 2022

Antti Karttunen, Table of n, a(n) for n = 0..11550
Victor Ufnarovski and Bo Åhlander, How to Differentiate a Number, J. Integer Seqs., Vol. 6, 2003, #03.3.4.
Wikimedia, Čiurlionis: Piramidziu sonata, Allegro (a medium resolution scan of the painting "Pyramid Sonata, Allegro")
Wikipedia, Mikalojus Konstantinas Čiurlionis
Index entries for sequences related to primorial base

Cf. A002110 (positions of 1's), A003415, A003557, A083345, A085731, A276086, A289234, A327860, A328571, A328572, A342001, A342005, A342006, A342016, A342022 (rgs-transform), A342417, A342419.
Cf. A342463 [= a(A329886(n))], A342920 [= a(A108951(n))], A342921 [= a(A276156(n))], A342017 [= A342007(a(n))], A342019 [= A129251(a(n))].
Cf. A166486 (a(n) mod 2, parity of terms, see comment in A327860), A353640 (a(n) mod 4).
Cf. A344760, A344761, A344762, A346252, A346253 and A345930, A353572, A353574 for permuted and other variants.
Cf. A351952 (similar definition, but using factorial base, with quite a different look).

A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A342002(n) = A342001(A276086(n)); \\ Uses also code from A342001.

A342002(n) = { my(s=0, m=1, p=2, e); while(n, e = (n%p); m *= p^(e>0); s += (e/p); n = n\p; p = nextprime(1+p)); (s*m); }; \\ Antti Karttunen, Mar 12 2021

A342002(n) = { my(s=0, p=2, e); while(n, e = (n%p); s += (e/p); n = n\p; p = nextprime(1+p)); numerator(s); }; \\ (Taking denominator instead would give A328571) - Antti Karttunen, Mar 12 2021

a(n) = A342001(A276086(n)) = A083345(A276086(n)).
a(n) = A327860(n) / A328572(n) = A003415(A276086(n)) / A003557(A276086(n)).
From Antti Karttunen, Jul 18 2021: (Start)
There are several permutations of this sequence. The following formulas show the relations:
a(n) = A344760(A289234(n)).
a(n) = A346252(A328623(n)) = A346253(A328622(n)).
a(n) = A344761(A328626(n)) = A344762(A328625(n)).
(End)

Sequence renamed as "Čiurlionis sequence" to honor Lithuanian artist Mikalojus Čiurlionis - Antti Karttunen, Apr 30 2022

A342920 a(n) = A342002(A108951(n)).

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 8, 12, 2, 1, 4, 1, 2, 6, 16, 1, 24, 1, 4, 6, 2, 1, 26, 50, 2, 16, 4, 1, 62, 1, 10, 6, 2, 126, 48, 1, 2, 6, 18, 1, 24, 1, 4, 46, 2, 1, 22, 1486, 100, 6, 4, 1, 32, 94, 8, 6, 2, 1, 54, 1, 2, 72, 20, 264, 12, 1, 4, 6, 120, 1, 376, 1, 2, 1142, 4, 242, 12, 1, 36, 342, 2, 1, 48, 272, 2, 6, 8, 1, 92, 318

Offset: 1

Antti Karttunen, Apr 06 2021

nonn

Cf. A108951, A324886, A327860, A328572, A342002, A342463.

A034386(n) = prod(i=1, primepi(n), prime(i));
A108951(n) = { my(f=factor(n)); prod(i=1, #f~, A034386(f[i, 1])^f[i, 2]) };
A342920(n) = A342002(A108951(n)); \\ Needs also code from A342002.

a(n) = A342002(A108951(n)) = A327860(A108951(n)) / A328572(A108951(n)).

A342456 A276086 applied to the primorial inflation of Doudna-tree, where A276086(n) is the prime product form of primorial base expansion of n.

Original entry on oeis.org

2, 3, 5, 9, 7, 25, 35, 15, 11, 49, 117649, 625, 717409, 1225, 55, 225, 13, 121, 1771561, 2401, 36226650889, 184877, 1127357, 875, 902613283, 514675673281, 3780549773, 1500625, 83852850675321384784127, 3025, 62004635, 21, 17, 169, 4826809, 14641, 8254129, 143, 2924207, 77, 8223741426987700773289, 59797108943, 546826709

Offset: 0

Antti Karttunen and Michael De Vlieger, Mar 14 2021

nonn

This sequence (which could be viewed as a binary tree, like the underlying A005940 and A329886) is similar to A324289, but unlike its underlying tree A283477 that generates only numbers that are products of distinct primorial numbers (i.e., terms of A129912), here the underlying tree A329886 generates all possible products of primorial numbers, i.e., terms of A025487, but in different order.

Cf. A005940, A025487, A108951, A129912, A276086, A283980, A324886, A342457 [= 2*A246277(a(n))], A342461 [= A001221(a(n))], A342462 [= A001222(a(n))], A342463 [= A342001(a(n))], A342464 [= A051903(a(n))].

Cf. A324289 (a subset of these terms, in different order).

Block[{a, f, r = MixedRadix[Reverse@ Prime@ Range@ 24]}, f[n_] :=
Times @@ MapIndexed[Prime[First[#2]]^#1 &, Reverse@ IntegerDigits[n, r]]; a[0] = 1; a[1] = 2; a[n_] := a[n] = If[EvenQ@ n, (Times @@ Map[Prime[PrimePi@ #1 + 1]^#2 & @@ # &, FactorInteger[#]] - Boole[# == 1])*2^IntegerExponent[#, 2] &[a[n/2]], 2 a[(n - 1)/2]]; Array[f@ a[#] &, 43, 0]] (* Michael De Vlieger, Mar 17 2021 *)

A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A283980(n) = {my(f=factor(n)); prod(i=1, #f~, my(p=f[i, 1], e=f[i, 2]); if(p==2, 6, nextprime(p+1))^e)};
A329886(n) = if(n<2,1+n,if(!(n%2),A283980(A329886(n/2)),2*A329886(n\2)));
A342456(n) = A276086(A329886(n));

a(n) = A276086(A329886(n)) = A324886(A005940(1+n)).
For all n >= 0, gcd(a(n), A329886(n)) = 1.
For all n >= 1, A055396(a(n))-1 = A061395(A329886(n)) = A290251(n) = 1+A080791(n).
For all n >= 0, a(2^n) = A000040(2+n).

A342462 Sum of digits when A329886(n) is written in primorial base, where A329886 is the primorial inflation of Doudna-tree.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 2, 2, 1, 2, 6, 4, 6, 4, 2, 4, 1, 2, 6, 4, 10, 6, 6, 4, 8, 12, 10, 8, 22, 4, 8, 2, 1, 2, 6, 4, 6, 2, 6, 2, 18, 10, 8, 6, 18, 12, 16, 4, 26, 16, 24, 8, 20, 14, 4, 6, 26, 16, 14, 8, 30, 6, 8, 4, 1, 2, 6, 4, 14, 12, 12, 8, 18, 12, 24, 4, 8, 12, 14, 4, 24, 20, 28, 20, 26, 16, 16, 12, 32, 26, 24, 14, 28, 16

Offset: 0

Antti Karttunen, Mar 15 2021

nonn

From David A. Corneth's Feb 27 2019 comment in A276150 follows that the only odd terms in this sequence are 1's occurring at 0 and at two's powers.

Subsequences starting at each n = 2^k are slowly converging towards A329886: 1, 2, 6, 4, 30, 12, 36, 8, 210, 60, 180, 24, etc.. Compare also to the behaviors of A324342 and A342463.

Antti Karttunen, Table of n, a(n) for n = 0..8192
Antti Karttunen, Data supplement: n, a(n) computed for n = 0..65537
Index entries for sequences related to binary expansion of n
Index entries for sequences related to primorial base

Cf. A001222, A005940, A108951, A276150, A329886, A342456, A342457, A342461, A342463, A342464.

Cf. also A324342, A324383, A324387, A324888.

A342462(n) = bigomega(A342456(n)); \\ Other code as in A342456.

A283980(n) = {my(f=factor(n)); prod(i=1, #f~, my(p=f[i, 1], e=f[i, 2]); if(p==2, 6, nextprime(p+1))^e)};
A329886(n) = if(n<2,1+n,if(!(n%2),A283980(A329886(n/2)),2*A329886(n\2)));
A276150(n) = { my(s=0, p=2, d); while(n, d = (n%p); s += d; n = (n-d)/p; p = nextprime(1+p)); (s); };
A342462(n) = A276150(A329886(n));

a(n) = A001222(A342456(n)) = A001222(A342457(n)).
a(n) = A276150(A329886(n)) = A324888(A005940(1+n)).
a(n) >= A342461(n).
For n >= 0, a(2^n) = 1.

cp's OEIS Frontend

A342001 Arithmetic derivative without its inherited divisor; the arithmetic derivative of n divided by A003557(n), which is a common divisor of both n and A003415(n).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A342002 Čiurlionis sequence: Arithmetic derivative without its inherited divisor applied to the primorial base exp-function: a(n) = A342001(A276086(n)).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

PARI

PARI

Formula

Extensions

A342920 a(n) = A342002(A108951(n)).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula

A342456 A276086 applied to the primorial inflation of Doudna-tree, where A276086(n) is the prime product form of primorial base expansion of n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A342462 Sum of digits when A329886(n) is written in primorial base, where A329886 is the primorial inflation of Doudna-tree.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

PARI

Formula