A342921 -id:A342921 - cp's OEIS frontend

A276156 Numbers obtained by reinterpreting base-2 representation of n in primorial base: a(0) = 0, a(2n) = A276154(a(n)), a(2n+1) = 1 + A276154(a(n)).

0, 1, 2, 3, 6, 7, 8, 9, 30, 31, 32, 33, 36, 37, 38, 39, 210, 211, 212, 213, 216, 217, 218, 219, 240, 241, 242, 243, 246, 247, 248, 249, 2310, 2311, 2312, 2313, 2316, 2317, 2318, 2319, 2340, 2341, 2342, 2343, 2346, 2347, 2348, 2349, 2520, 2521, 2522, 2523, 2526, 2527, 2528, 2529, 2550, 2551, 2552, 2553, 2556, 2557, 2558, 2559, 30030, 30031

Offset: 0

Antti Karttunen, Aug 24 2016

Numbers that are sums of distinct primorial numbers, A002110.

Numbers with no digits larger than one in primorial base, A049345.

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

Cf. A000040, A001511, A002110, A007088, A007814, A019565, A049345, A257993, A276084, A276085, A276154, A351073, A328461, A328473, A328474, A328571, A328831, A328836, A341518, A344591.
Complement of A177711.
Subsequences: A328233, A328832, A328462 (odd bisection).
Conjectured subsequences: A328110, A380527.
Fixed points of A328841, positions of zeros in A328828, A328842, and A329032, positions of ones in A328581, A328582, and A381032.
Positions of terms < 2 in A328114.
Indices where A327860 and A329029 coincide.
Cf. also table A328464 (and its rows).
Cf. also A059590, A283985, A290249, A342921.

nn = 65; b = MixedRadix[Reverse@ Prime@ Range[IntegerLength[nn, 2] - 1]]; Table[FromDigits[IntegerDigits[n, 2], b], {n, 0, 65}] (* Version 10.2, or *)
Table[Total[Times @@@ Transpose@ {Map[Times @@ # &, Prime@ Range@ Range[0, Length@ # - 1]], Reverse@ #}] &@ IntegerDigits[n, 2], {n, 0, 65}] (* Michael De Vlieger, Aug 26 2016 *)

A276156(n) = { my(s=0, p=1, r=1); while(n, if(n%2, s += r); n>>=1; p = nextprime(1+p); r *= p); (s); }; \\ Antti Karttunen, Feb 03 2022

from sympy import prime, primorial, primepi, factorint
from operator import mul
def a002110(n): return 1 if n<1 else primorial(n)
def a276085(n):
    f=factorint(n)
    return sum([f[i]*a002110(primepi(i) - 1) for i in f])
def a019565(n): return reduce(mul, (prime(i+1) for i, v in enumerate(bin(n)[:1:-1]) if v == '1')) # after Chai Wah Wu
def a(n): return 0 if n==0 else a276085(a019565(n))
print([a(n) for n in range(101)]) # Indranil Ghosh, Jun 23 2017

a(0) = 0, a(2n) = A276154(a(n)), a(2n+1) = 1+A276154(a(n)).
Other identities. For all n >= 0:
a(n) = A276085(A019565(n)).
A049345(a(n)) = A007088(n).
A257993(a(n)) = A001511(n).
A276084(a(n)) = A007814(n).
A051903(a(n)) = A351073(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

cp's OEIS Frontend

A276156 Numbers obtained by reinterpreting base-2 representation of n in primorial base: a(0) = 0, a(2n) = A276154(a(n)), a(2n+1) = 1 + A276154(a(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