A351952 -id:A351952 - cp's OEIS frontend

A351953 a(n) = A351952(A225901(n)).

0, 1, 2, 7, 1, 5, 3, 11, 19, 53, 14, 43, 2, 9, 16, 47, 11, 37, 1, 7, 13, 41, 8, 31, 4, 15, 26, 73, 19, 59, 41, 117, 193, 491, 158, 421, 34, 103, 172, 449, 137, 379, 27, 89, 151, 407, 116, 337, 3, 13, 23, 67, 16, 53, 36, 107, 178, 461, 143, 391, 29, 93, 157, 419, 122, 349, 22, 79, 136, 377, 101, 307, 2, 11, 20, 61

Offset: 0

Antti Karttunen, Apr 02 2022

Antti Karttunen, Table of n, a(n) for n = 0..20160
Index entries for sequences related to factorial base representation

Cf. A003415, A003557, A225901, A276076, A342001, A351952.

Cf. also A344760.

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]));
A276076(n) = { my(i=0,m=1,f=1,nextf); while((n>0),i=i+1; nextf = (i+1)*f; if((n%nextf),m*=(prime(i)^((n%nextf)/f));n-=(n%nextf));f=nextf); m; };
A351952(n) = { my(u=A276076(n)); (A003415(u) / A003557(u)); };
A225901(n) = { my(s=0, d, k=2); while(n, d=n%k; n=n\k; if(d, s=s+(k-d)*(k-1)!); k=k+1); return(s); }; \\ From A225901.
A351953(n) = A351952(A225901(n));

a(n) = A351952(A225901(n)) = A342001(A276076(A225901(n))).

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

A276076 Factorial base exp-function: digits in factorial base representation of n become the exponents of successive prime factors whose product a(n) is.

Original entry on oeis.org

1, 2, 3, 6, 9, 18, 5, 10, 15, 30, 45, 90, 25, 50, 75, 150, 225, 450, 125, 250, 375, 750, 1125, 2250, 7, 14, 21, 42, 63, 126, 35, 70, 105, 210, 315, 630, 175, 350, 525, 1050, 1575, 3150, 875, 1750, 2625, 5250, 7875, 15750, 49, 98, 147, 294, 441, 882, 245, 490, 735, 1470, 2205, 4410, 1225, 2450, 3675, 7350, 11025, 22050, 6125, 12250, 18375, 36750, 55125, 110250, 343

Offset: 0

Antti Karttunen, Aug 18 2016

These are prime-factorization representations of single-variable polynomials where the coefficient of term x^(k-1) (encoded as the exponent of prime(k) in the factorization of n) is equal to the digit in one-based position k of the factorial base representation of n. See the examples.

			   n  A007623   polynomial     encoded as             a(n)
   -------------------------------------------------------
   0       0    0-polynomial   (empty product)        = 1
   1       1    1*x^0          prime(1)^1             = 2
   2      10    1*x^1          prime(2)^1             = 3
   3      11    1*x^1 + 1*x^0  prime(2) * prime(1)    = 6
   4      20    2*x^1          prime(2)^2             = 9
   5      21    2*x^1 + 1*x^0  prime(2)^2 * prime(1)  = 18
   6     100    1*x^2          prime(3)^1             = 5
   7     101    1*x^2 + 1*x^0  prime(3) * prime(1)    = 10
and:
  23     321  3*x^2 + 2*x + 1  prime(3)^3 * prime(2)^2 * prime(1)
                                      = 5^3 * 3^2 * 2 = 2250.

Antti Karttunen, Table of n, a(n) for n = 0..5040
Indranil Ghosh, Python program for computing this sequence.
Index entries for sequences related to factorial base representation.

Cf. A000040, A007623, A084558, A099563, A257687, A276009.
Cf. A276075 (a left inverse).
Cf. A276078 (same terms in ascending order).
Cf. also A000142, A001221, A001222, A002110, A007489, A008836, A019565, A033312, A034968, A048675, A051903, A059590, A060130, A076954, A246359, A248663, A262725, A276073, A276074, A351576, A351577, A351950, A351951, A351952, A351954.
Cf. also A275733, A275734, A275735, A275725 for other such encodings of factorial base related polynomials, and A276086 for a primorial base analog.

a[n_] := Module[{k = n, m = 2, r, p = 2, q = 1}, While[{k, r} = QuotientRemainder[k, m]; k != 0|| r != 0, q *= p^r; p = NextPrime[p]; m++]; q]; Array[a, 100, 0] (* Amiram Eldar, Feb 07 2024 *)

a(0) = 1, for n >= 1, a(n) = A275733(n) * a(A276009(n)).
Or: for n >= 1, a(n) = a(A257687(n)) * A000040(A084558(n))^A099563(n).
Other identities.
For all n >= 0:
A276075(a(n)) = n.
A001221(a(n)) = A060130(n).
A001222(a(n)) = A034968(n).
A051903(a(n)) = A246359(n).
A048675(a(n)) = A276073(n).
A248663(a(n)) = A276074(n).
a(A007489(n)) = A002110(n).
a(A059590(n)) = A019565(n).
For all n >= 1:
a(A000142(n)) = A000040(n).
a(A033312(n)) = A076954(n-1).
From Antti Karttunen, Apr 18 2022: (Start)
a(n) = A276086(A351576(n)).
A276085(a(n)) = A351576(n)
A003557(a(n)) = A351577(n).
A003415(a(n)) = A351950(n).
A069359(a(n)) = A351951(n).
A083345(a(n)) = A342001(a(n)) = A351952(n).
A351945(a(n)) = A351954(n).
A181819(a(n)) = A275735(n).
(End)
lambda(a(n)) = A262725(n+1), where lambda is Liouville's function, A008836. - Antti Karttunen and Peter Munn, Aug 09 2024

Name changed by Antti Karttunen, Apr 18 2022

A351950 Arithmetic derivative of the factorial base exp-function: a(n) = A003415(A276076(n)).

Original entry on oeis.org

0, 1, 1, 5, 6, 21, 1, 7, 8, 31, 39, 123, 10, 45, 55, 185, 240, 705, 75, 275, 350, 1075, 1425, 3975, 1, 9, 10, 41, 51, 165, 12, 59, 71, 247, 318, 951, 95, 365, 460, 1445, 1905, 5385, 650, 2175, 2825, 8275, 11100, 30075, 14, 77, 91, 329, 420, 1281, 119, 483, 602, 1939, 2541, 7287, 840, 2905, 3745, 11165, 14910, 40845

Offset: 0

Antti Karttunen, Apr 01 2022

Antti Karttunen, Table of n, a(n) for n = 0..10080
Index entries for sequences related to factorial base representation

Cf. A003415, A276076, A351576, A351952.
Differs from a similarly defined A327860 for the first time at n=24.
Coincides with A351951 on n given by A059590.

A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A276076(n) = { my(i=0,m=1,f=1,nextf); while((n>0),i=i+1; nextf = (i+1)*f; if((n%nextf),m*=(prime(i)^((n%nextf)/f));n-=(n%nextf));f=nextf); m; };
A351950(n) = A003415(A276076(n));

a(n) = A003415(A276076(n)).

a(n) = A327860(A351576(n)).

A351954 Arithmetic derivative without its inherited divisor applied to the prime shadow of the factorial base exp-function: a(n) = A342001(A181819(A276076(n))).

Original entry on oeis.org

0, 1, 1, 2, 1, 5, 1, 2, 2, 3, 5, 8, 1, 5, 5, 8, 2, 7, 1, 7, 7, 12, 8, 31, 1, 2, 2, 3, 5, 8, 2, 3, 3, 4, 8, 11, 5, 8, 8, 11, 7, 10, 7, 12, 12, 17, 31, 46, 1, 5, 5, 8, 2, 7, 5, 8, 8, 11, 7, 10, 2, 7, 7, 10, 3, 9, 8, 31, 31, 46, 13, 41, 1, 7, 7, 12, 8, 31, 7, 12, 12, 17, 31, 46, 8, 31, 31, 46, 13, 41, 2, 9, 9, 14, 11

Offset: 0

Antti Karttunen, Apr 02 2022

Compare the scatter plot to those of A275735, A353575 and of A353577. - Antti Karttunen, Apr 30 2022

Antti Karttunen, Table of n, a(n) for n = 0..20160
Index entries for sequences related to factorial base representation

Cf. A003415, A003557, A181819, A275735, A276076, A342001, A351576, A351945, A351952, A353575, A353577.

Cf. A051683 (positions of 1's).

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]));
A181819(n) = factorback(apply(e->prime(e),(factor(n)[,2])));
A276076(n) = { my(i=0,m=1,f=1,nextf); while((n>0),i=i+1; nextf = (i+1)*f; if((n%nextf),m*=(prime(i)^((n%nextf)/f));n-=(n%nextf));f=nextf); m; };
A342001(n) = (A003415(n) / A003557(n));
A351945(n) = A342001(A181819(n));
A351954(n) = A351945(A276076(n));

a(n) = A342001(A275735(n)) = A351945(A276076(n)).

a(n) = A353577(A351576(n)). - Antti Karttunen, Apr 30 2022

Verbal description added to the definition by Antti Karttunen, Apr 30 2022

A351577 a(n) = A003557(A276076(n)).

Original entry on oeis.org

1, 1, 1, 1, 3, 3, 1, 1, 1, 1, 3, 3, 5, 5, 5, 5, 15, 15, 25, 25, 25, 25, 75, 75, 1, 1, 1, 1, 3, 3, 1, 1, 1, 1, 3, 3, 5, 5, 5, 5, 15, 15, 25, 25, 25, 25, 75, 75, 7, 7, 7, 7, 21, 21, 7, 7, 7, 7, 21, 21, 35, 35, 35, 35, 105, 105, 175, 175, 175, 175, 525, 525, 49, 49, 49, 49, 147, 147, 49, 49, 49, 49, 147, 147, 245

Offset: 0

Antti Karttunen, Apr 01 2022

Antti Karttunen, Table of n, a(n) for n = 0..10080
Index entries for sequences related to factorial base representation

Cf. A003557, A085731, A276009, A276076, A328572, A351576, A351950, A351952.

A003557(n) = (n/factorback(factorint(n)[, 1]));
A276076(n) = { my(i=0,m=1,f=1,nextf); while((n>0),i=i+1; nextf = (i+1)*f; if((n%nextf),m*=(prime(i)^((n%nextf)/f));n-=(n%nextf));f=nextf); m; };
A351577(n) = A003557(A276076(n));

a(n) = A003557(A276076(n)).
a(n) = A276076(A276009(n)).
a(n) = A328572(A351576(n)).
a(n) = A085731(A276076(n)) = gcd(A276076(n), A351950(n)).

cp's OEIS Frontend

A351953 a(n) = A351952(A225901(n)).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

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

A276076 Factorial base exp-function: digits in factorial base representation of n become the exponents of successive prime factors whose product a(n) is.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A351950 Arithmetic derivative of the factorial base exp-function: a(n) = A003415(A276076(n)).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula

A351954 Arithmetic derivative without its inherited divisor applied to the prime shadow of the factorial base exp-function: a(n) = A342001(A181819(A276076(n))).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Formula

Extensions

A351577 a(n) = A003557(A276076(n)).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula