cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 26 results. Next

A333427 Numbers k such that k and k+1 are both primorial base Niven numbers (A333426).

Original entry on oeis.org

1, 8, 24, 32, 44, 64, 65, 132, 212, 224, 244, 245, 296, 368, 424, 425, 468, 560, 656, 720, 728, 737, 869, 1056, 1088, 1416, 1572, 1728, 2100, 2312, 2324, 2344, 2345, 2524, 2525, 2568, 2600, 2672, 2820, 2960, 3032, 3132, 3156, 3200, 3288, 3392, 3444, 4096, 4424
Offset: 1

Views

Author

Amiram Eldar, Mar 20 2020

Keywords

Examples

			1 is a term since 1 and 2 are both primorial base Niven numbers.
		

Crossrefs

Programs

  • Mathematica
    max = 6; bases = Prime @ Range[max, 1, -1]; nmax = Times @@ bases - 1; primNivenQ[n_] := Divisible[n, Plus @@ IntegerDigits[n, MixedRadix[bases]]]; q1 = primNivenQ[1]; seq = {}; Do[q2 = primNivenQ[n]; If[q1 && q2, AppendTo[seq, n - 1]]; q1 = q2, {n, 2, nmax}]; seq

A333428 Starts of runs of 3 consecutive primorial base Niven numbers (A333426).

Original entry on oeis.org

64, 244, 424, 2344, 2524, 4624, 16180, 30064, 30244, 32344, 43900, 60064, 71620, 91408, 99340, 127060, 154780, 182500, 210220, 250936, 338632, 365860, 477280, 510544, 510724, 512824, 513160, 540544, 540880, 790900, 842884, 876988, 1021024, 1021648, 1024000, 1051720
Offset: 1

Views

Author

Amiram Eldar, Mar 20 2020

Keywords

Examples

			64 is a term since 64, 65 and 66 are all primorial base Niven numbers.
		

Crossrefs

Programs

  • Mathematica
    max = 7; bases = Prime @ Range[max, 1, -1]; nmax = Times @@ bases - 1; primNivenQ[n_] := Divisible[n, Plus @@ IntegerDigits[n, MixedRadix[bases]]]; q1 = primNivenQ[1]; q2 = primNivenQ[2]; seq = {}; Do[q3 = primNivenQ[n]; If[q1 && q2 && q3, AppendTo[seq, n - 2]]; q1 = q2; q2 = q3, {n, 3, nmax}]; seq

A347496 Primorial base Niven numbers (A333426) with a record gap to the next primorial base Niven number.

Original entry on oeis.org

1, 2, 9, 12, 25, 50, 120, 344, 400, 770, 1120, 3920, 13566, 13734, 19845, 22748, 148148, 167854, 176220, 889896, 2946216, 3685416, 5072256, 7139280, 8521056, 9058900, 9625336, 17825857, 19392072, 27504848, 76952788, 106691001, 162789696, 198582784, 212847225
Offset: 1

Views

Author

Amiram Eldar, Sep 03 2021

Keywords

Comments

The corresponding gaps are 1, 2, 3, 4, 5, 10, 12, 16, 20, 34, 37, 48, 54, 66, 75, 121, 132, 146, 180, 238, 241, 248, 288, 302, 314, 332, 336, 343, 348, 400, 476, 479, 484, 496, 500, ...

Examples

			The first 8 primorial base Niven numbers are 1, 2, 4, 6, 8, 9, 12 and 16. The gaps between them are 1, 2, 2, 2, 1, 3 and 4. The record gaps, 1, 2, 3 and 4, occur after the terms 1, 2, 9 and 12.
		

Crossrefs

Programs

  • Mathematica
    max = 7; bases = Prime @ Range[max, 1, -1]; nmax = Times @@ bases - 1; sumdig[n_] := Plus @@ IntegerDigits[n, MixedRadix[bases]]; primoNivenQ[n_] := Divisible[n, sumdig[n]]; gapmax = 0; n1 = 1; s = {}; Do[If[primoNivenQ[n], gap = n - n1; If[gap > gapmax, gapmax = gap; AppendTo[s, n1]]; n1 = n], {n, 2, nmax}]; s

A276150 Sum of digits when n is written in primorial base (A049345); minimal number of primorials (A002110) that add to n.

Original entry on oeis.org

0, 1, 1, 2, 2, 3, 1, 2, 2, 3, 3, 4, 2, 3, 3, 4, 4, 5, 3, 4, 4, 5, 5, 6, 4, 5, 5, 6, 6, 7, 1, 2, 2, 3, 3, 4, 2, 3, 3, 4, 4, 5, 3, 4, 4, 5, 5, 6, 4, 5, 5, 6, 6, 7, 5, 6, 6, 7, 7, 8, 2, 3, 3, 4, 4, 5, 3, 4, 4, 5, 5, 6, 4, 5, 5, 6, 6, 7, 5, 6, 6, 7, 7, 8, 6, 7, 7, 8, 8, 9, 3, 4, 4, 5, 5, 6, 4, 5, 5, 6, 6, 7, 5, 6, 6, 7, 7, 8, 6, 7, 7, 8, 8, 9, 7, 8, 8, 9, 9, 10, 4
Offset: 0

Views

Author

Antti Karttunen, Aug 22 2016

Keywords

Comments

The sum of digits of n in primorial base is odd if n is 1 or 2 (mod 4) and even if n is 0 or 3 (mod 4). Proof: primorials are 1 or 2 (mod 4) and a(n) can be constructed via the greedy algorithm. So if n = 4k + r where 0 <= r < 4, 4k needs an even number of primorials and r needs hammingweight(r) = A000120(r) primorials. Q.E.D. - David A. Corneth, Feb 27 2019

Examples

			For n=24, which is "400" in primorial base (as 24 = 4*(3*2*1) + 0*(2*1) + 0*1, see A049345), the sum of digits is 4, thus a(24) = 4.
		

Crossrefs

Cf. A333426 [k such that a(k)|k], A339215 [numbers not of the form x+a(x) for any x], A358977 [k such that gcd(k, a(k)) = 1].
Cf. A014601, A042963 (positions of even and odd terms), A343048 (positions of records).
Differs from analogous A034968 for the first time at n=24.

Programs

  • Mathematica
    nn = 120; b = MixedRadix[Reverse@ Prime@ NestWhileList[# + 1 &, 1, Times @@ Prime@ Range[# + 1] <= nn &]]; Table[Total@ IntegerDigits[n, b], {n, 0, nn}] (* Version 10.2, or *)
    nn = 120; f[n_] := Block[{a = {{0, n}}}, Do[AppendTo[a, {First@ #, Last@ #} &@ QuotientRemainder[a[[-1, -1]], Times @@ Prime@ Range[# - i]]], {i, 0, #}] &@ NestWhile[# + 1 &, 0, Times @@ Prime@ Range[# + 1] <= n &]; Rest[a][[All, 1]]]; Table[Total@ f@ n, {n, 0, 120}] (* Michael De Vlieger, Aug 26 2016 *)
  • PARI
    A276150(n) = { my(s=0, p=2, d); while(n, d = (n%p); s += d; n = (n-d)/p; p = nextprime(1+p)); (s); }; \\ Antti Karttunen, Feb 27 2019
  • Python
    from sympy import prime, primefactors
    def Omega(n): return 0 if n==1 else Omega(n//primefactors(n)[0]) + 1
    def a276086(n):
        i=0
        m=pr=1
        while n>0:
            i+=1
            N=prime(i)*pr
            if n%N!=0:
                m*=(prime(i)**((n%N)/pr))
                n-=n%N
            pr=N
        return m
    def a(n): return Omega(a276086(n))
    print([a(n) for n in range(201)]) # Indranil Ghosh, Jun 23 2017
    

Formula

a(n) = 1 + a(A276151(n)) = 1 + a(n-A002110(A276084(n))), a(0) = 0.
or for n >= 1: a(n) = 1 + a(n-A260188(n)).
Other identities and observations. For all n >= 0:
a(n) = A001222(A276086(n)) = A001222(A278226(n)).
a(n) >= A371091(n) >= A267263(n).
From Antti Karttunen, Feb 27 2019: (Start)
a(n) = A000120(A277022(n)).
a(A283477(n)) = A324342(n).
(End)
a(n) = A373606(n) + A373607(n). - Antti Karttunen, Jun 19 2024

A334308 Base phi Niven numbers: numbers divisible by the number of 1's in their base phi representation (A055778).

Original entry on oeis.org

1, 2, 6, 12, 15, 16, 18, 20, 30, 35, 36, 45, 48, 55, 60, 70, 72, 78, 84, 90, 91, 95, 96, 98, 104, 108, 132, 144, 147, 154, 168, 175, 184, 189, 208, 224, 231, 232, 245, 252, 256, 261, 264, 270, 275, 280, 282, 287, 294, 315, 322, 324, 330, 336, 340, 342, 351, 357
Offset: 1

Views

Author

Amiram Eldar, Apr 22 2020

Keywords

Examples

			6 is a term since its base phi representation is 1010.0001, and the number of 1's is 3, which is a divisor of 6.
		

Crossrefs

Programs

  • Mathematica
    phiDigSum[1] = 1; phiDigSum[n_] := Plus @@ RealDigits[n, GoldenRatio, 2*Ceiling[ Log[GoldenRatio, n]] ][[1]]; Select[Range[360], Divisible[#, phiDigSum[#]] &]

A342426 Niven numbers in base 3/2: numbers divisible by their sum of digits in fractional base 3/2 (A244040).

Original entry on oeis.org

1, 2, 6, 9, 14, 21, 40, 42, 56, 72, 84, 108, 110, 120, 126, 130, 143, 154, 156, 162, 165, 168, 169, 176, 180, 182, 189, 198, 220, 225, 231, 243, 252, 280, 288, 297, 306, 308, 320, 322, 330, 336, 348, 350, 364, 390, 423, 430, 432, 459, 460, 462, 480, 490, 504
Offset: 1

Views

Author

Amiram Eldar, Mar 11 2021

Keywords

Examples

			6 is a term since its representation in base 3/2 is 210 and 2 + 1 + 0 = 3 is a divisor of 6.
9 is a term since its representation in base 3/2 is 2100 and 2 + 1 + 0 + 0 = 3 is a divisor of 9.
		

Crossrefs

Subsequences: A342427, A342428, A342429.
Similar sequences: A005349 (decimal), A049445 (binary), A064150 (ternary), A064438 (quaternary), A064481 (base 5), A118363 (factorial), A328208 (Zeckendorf), A328212 (lazy Fibonacci), A331085 (negaFibonacci), A333426 (primorial), A334308 (base phi), A331728 (negabinary).

Programs

  • Mathematica
    s[0] = 0; s[n_] := s[n] = s[2*Floor[n/3]] + Mod[n, 3]; q[n_] := Divisible[n, s[n]]; Select[Range[500], q]

A342726 Niven numbers in base i-1: numbers that are divisible by the sum of their digits in base i-1.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 10, 12, 15, 16, 18, 20, 24, 25, 30, 32, 33, 35, 36, 40, 42, 44, 45, 48, 50, 54, 60, 64, 65, 66, 70, 77, 80, 88, 90, 96, 99, 100, 110, 112, 120, 124, 125, 126, 130, 140, 144, 145, 147, 150, 156, 160, 168, 170, 180, 182, 184, 185, 186, 190, 192
Offset: 1

Views

Author

Amiram Eldar, Mar 19 2021

Keywords

Comments

Numbers k that are divisible by A066323(k).
Equivalently, Niven numbers in base -4, since A066323(k) is also the sum of the digits of k in base -4.

Examples

			2 is a term since its representation in base i-1 is 1100 and 1+1+0+0 = 2 is a divisor of 2.
10 is a term since its representation in base i-1 is 111001100 and 1+1+1+0+0+1+1+0+0 = 5 is a divisor of 10.
		

Crossrefs

Similar sequences: A005349 (decimal), A049445 (binary), A064150 (ternary), A064438 (quaternary), A064481 (base 5), A118363 (factorial), A328208 (Zeckendorf), A328212 (lazy Fibonacci), A331085 (negaFibonacci), A333426 (primorial), A334308 (base phi), A331728 (negabinary), A342426 (base 3/2).

Programs

  • Mathematica
    v = {{0, 0, 0, 0}, {0, 0, 0, 1}, {1, 1, 0, 0}, {1, 1, 0, 1}}; q[n_] := Divisible[n, Total[Flatten @ v[[1 + Reverse @ Most[Mod[NestWhileList[(# - Mod[#, 4])/-4 &, n, # != 0 &], 4]]]]]]; Select[Range[200], q]

A344341 Gray-code Niven numbers: numbers divisible by the number of 1's in their binary reflected Gray code (A005811).

Original entry on oeis.org

1, 2, 3, 4, 6, 7, 8, 9, 12, 14, 15, 16, 20, 24, 27, 28, 30, 31, 32, 33, 36, 39, 40, 42, 44, 45, 48, 51, 52, 56, 57, 60, 62, 63, 64, 68, 72, 75, 76, 80, 84, 88, 90, 92, 96, 99, 100, 104, 105, 108, 111, 112, 116, 120, 123, 124, 126, 127, 128, 129, 132, 135, 136
Offset: 1

Views

Author

Amiram Eldar, May 15 2021

Keywords

Examples

			2 is a term since its Gray code is 11 and 1+1 = 2 is a divisor of 2.
6 is a term since its Gray code is 101 and 1+0+1 = 2 is a divisor of 6.
		

Crossrefs

Subsequences: A344342, A344343, A344344.
Similar sequences: A005349 (decimal), A049445 (binary), A064150 (ternary), A064438 (quaternary), A064481 (base 5), A118363 (factorial), A328208 (Zeckendorf), A328212 (lazy Fibonacci), A331085 (negaFibonacci), A333426 (primorial), A334308 (base phi), A331728 (negabinary), A342426 (base 3/2), A342726 (base i-1).

Programs

  • Mathematica
    gcNivenQ[n_] := Divisible[n, DigitCount[BitXor[n, Floor[n/2]], 2, 1]]; Select[Range[150], gcNivenQ]

A351714 Lucas-Niven numbers: numbers that are divisible by the number of terms in their minimal (or greedy) representation in terms of the Lucas numbers (A130310).

Original entry on oeis.org

1, 2, 3, 4, 6, 7, 8, 10, 11, 12, 14, 18, 20, 22, 24, 27, 29, 30, 32, 36, 39, 40, 42, 47, 48, 50, 54, 57, 58, 60, 64, 66, 69, 72, 76, 78, 80, 81, 84, 90, 92, 94, 96, 100, 104, 108, 120, 123, 124, 126, 129, 130, 132, 134, 135, 138, 140, 144, 152, 153, 156, 159, 160
Offset: 1

Views

Author

Amiram Eldar, Feb 17 2022

Keywords

Comments

Numbers k such that A116543(k) | k.

Examples

			6 is a term since its minimal Lucas representation, A130310(6) = 1001, has A116543(6) = 2 1's and 6 is divisible by 2.
		

Crossrefs

Programs

  • Mathematica
    lucasNivenQ[n_] := Module[{s = {}, m = n, k = 1}, While[m > 0, If[m == 1, k = 1; AppendTo[s, k]; m = 0, If[m == 2, k = 0; AppendTo[s, k]; m = 0, While[LucasL[k] <= m, k++]; k--; AppendTo[s, k]; m -= LucasL[k]; k = 1]]]; Divisible[n, Plus @@ IntegerDigits[Total[2^s], 2]]]; Select[Range[160], lucasNivenQ]

A351719 Lazy-Lucas-Niven numbers: numbers divisible by the number of terms in their maximal (or lazy) representation in terms of the Lucas numbers (A130311).

Original entry on oeis.org

1, 2, 4, 6, 9, 12, 16, 20, 25, 40, 42, 54, 60, 66, 78, 84, 91, 96, 104, 112, 120, 126, 144, 154, 161, 168, 175, 176, 180, 182, 184, 192, 203, 210, 216, 217, 224, 232, 234, 240, 243, 264, 270, 280, 288, 304, 306, 310, 315, 320, 322, 328, 336, 344, 350, 360, 378
Offset: 1

Views

Author

Amiram Eldar, Feb 17 2022

Keywords

Comments

Numbers k such that A131343(k) | k.

Examples

			6 is a term since its maximal Lucas representation, A130311(6) = 111, has A131343(6) = 3 1's and 6 is divisible by 3.
		

Crossrefs

Programs

  • Mathematica
    lazy = Select[IntegerDigits[Range[3000], 2], SequenceCount[#, {0, 0}] == 0 &]; t = Total[# * Reverse @ LucasL[Range[0, Length[#] - 1]]] & /@ lazy; s = FromDigits /@ lazy[[TakeWhile[Flatten[FirstPosition[t, #] & /@ Range[Max[t]]], NumberQ]]]; Position[Divisible[Range[Length[s]], Plus @@@ IntegerDigits[s]], True] // Flatten
Showing 1-10 of 26 results. Next