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 11 results. Next

A325499 Difference sequence of A325424.

Original entry on oeis.org

1, 5, 2, 2, 2, 1, 3, 3, 1, 4, 1, 5, 1, 1, 4, 1, 1, 6, 2, 2, 1, 5, 1, 1, 2, 2, 7, 1, 2, 2, 1, 7, 2, 2, 1, 1, 2, 3, 1, 4, 6, 1, 1, 2, 2, 1, 7, 4, 1, 3, 2, 1, 1, 2, 2, 1, 1, 5, 1, 4, 1, 5, 2, 2, 2, 1, 1, 2, 3, 1, 4, 7, 1, 4, 1, 1, 5, 1, 2, 2, 1, 3, 2, 1, 1, 2
Offset: 1

Views

Author

Clark Kimberling, May 05 2019

Keywords

Comments

See A325417 for a guide to related sequences.
Conjecture: every term is in {1,2,3,4,5,6,7}.

Links

Crossrefs

Cf. A325424, A325417, A036668, A325498.

Programs

  • Mathematica
    a = {1}; Do[AppendTo[a, NestWhile[# + 1 &, Last[a] + 1,
Apply[Or, Map[MemberQ[a, #] &, Select[Flatten[{#/3, #/2}],
IntegerQ]]] &]], {2000}]; a ;       (* A036668 *)
c = Complement[Range[Last[a]], a] ; (* A325424 *)
Differences[a]  (* A325498 *)
Differences[c]  (* A325499 *)
(* Peter J. C. Moses, Apr 23 2019 *)

Formula

Asymptotic mean: lim_{n->oo} (1/n) * Sum_{k=1..n} a(k) = 12/5. - Amiram Eldar, Nov 26 2020

A325417 a(n) is the least number not 2*a(m) or 3*a(m)+1 for any m < n.

Original entry on oeis.org

1, 3, 5, 7, 8, 9, 11, 12, 13, 15, 17, 19, 20, 21, 23, 27, 29, 31, 32, 33, 35, 36, 39, 41, 43, 44, 45, 47, 48, 49, 50, 51, 53, 55, 56, 57, 59, 60, 63, 65, 67, 68, 69, 71, 73, 74, 75, 76, 77, 79, 80, 81, 83, 84, 85, 87, 89, 91, 92, 93, 95, 99, 101, 103, 104
Offset: 1

Views

Author

Clark Kimberling, Apr 24 2019

Keywords

Comments

In column 1 of the following guide to related sequences, disallowed terms are indicated by the variable x representing a(m) for m < n.
Disallowed Sequence(a) Complement(c) Differences(a) Differences(c)
2x, 3x+1 A325417 A325418 A325444 A325445
3x, 2x+1 A325419 A325420 A325494 A325495
2x+1, 3x+1 A077477 A325422 A325496 A325497
2x, 3x A036668 A325424 A325498 A325499
[3x/2], 2x A325425 A325426 A325518 A325519
[3x/2], 2x+1 A325427 A325428 A325520 A325521
[3x/2], 3x A325429 A325430 A325522 A325523
3x, 4x A325431 A325432 A325525 A325526
2x, 3x-1 A325462 A325463 A325526 A325527
2x, 3x-2 A325464 A325465 A325528 A325529
2x-1, 3x-1 A325440 A325441 A325530 A325531
2x-1, 3x A325442 A325443 A325532 A325533
2x+1, 3x+2 A325539 A325540 A325541 A325542

Examples

			The sequence necessarily starts with 1.  The next 2 terms are determined as follows:  because a(1) = 1, the numbers 2 and 4 are disallowed, so that a(2) = 3, whence the numbers 6 and 10 are disallowed, so that a(3) = 5.

Links

Crossrefs

Cf. A325418, A325444.

Programs

  • Mathematica
    a = {1}; Do[AppendTo[a, NestWhile[# + 1 &, Last[a] + 1,
Apply[Or, Map[MemberQ[a, #] &, Select[Flatten[{#/2, (# - 1)/3}],
IntegerQ]]] &]], {150}]; a     (* A325417 *)
Complement[Range[Last[a]], a]  (* A325418 *)
(* Peter J. C. Moses, Apr 23 2019 *)

A036668 Hati numbers: of form 2^i*3^j*k, i+j even, (k,6)=1.

Original entry on oeis.org

1, 4, 5, 6, 7, 9, 11, 13, 16, 17, 19, 20, 23, 24, 25, 28, 29, 30, 31, 35, 36, 37, 41, 42, 43, 44, 45, 47, 49, 52, 53, 54, 55, 59, 61, 63, 64, 65, 66, 67, 68, 71, 73, 76, 77, 78, 79, 80, 81, 83, 85, 89, 91, 92, 95, 96, 97, 99, 100, 101, 102, 103, 107
Offset: 1

Views

Author

N. J. A. Sloane, Antreas P. Hatzipolakis (xpolakis(AT)hol.gr)

Keywords

Comments

If n appears then 2n and 3n do not. - Benoit Cloitre, Jun 13 2002
Closed under multiplication. Each term is a product of a unique subset of {6} U A050376 \ {2,3}. - Peter Munn, Sep 14 2019

Links

Crossrefs

Cf. A003159, A007310, A014601, A036667, A050376, A052330, A325424 (complement), A325498 (first differences), A373136 (characteristic function).
Positions of 0's in A182582.
Subsequences: A084087, A339690, A352272, A352273.
Cf. also A121537, A145204, A225837, A307150, A329575, A334747, A334748, A339746.

Programs

  • Maple
    N:= 1000: # to get all terms up to N
A:= {seq(2^i,i=0..ilog2(N))}:
Ae,Ao:= selectremove(issqr,A):
Be:= map(t -> seq(t*9^j, j=0 .. floor(log[9](N/t))),Ae):
Bo:= map(t -> seq(t*3*9^j,j=0..floor(log[9](N/(3*t)))),Ao):
B:= Be union Bo:
C1:= map(t -> seq(t*(6*i+1),i=0..floor((N/t -1)/6)),B):
C2:= map(t -> seq(t*(6*i+5),i=0..floor((N/t - 5)/6)),B):
A036668:= C1 union C2; # Robert Israel, May 09 2014
  • Mathematica
    a = {1}; Do[AppendTo[a, NestWhile[# + 1 &, Last[a] + 1,
Apply[Or, Map[MemberQ[a, #] &, Select[Flatten[{#/3, #/2}],
IntegerQ]]] &]], {150}]; a  (* A036668 *)
(* Peter J. C. Moses, Apr 23 2019 *)
  • PARI
    twos(n) = {local(r,m);r=0;m=n;while(m%2==0,m=m/2;r++);r}
threes(n) = {local(r,m);r=0;m=n;while(m%3==0,m=m/3;r++);r}
isA036668(n) = (twos(n)+threes(n))%2==0 \\ Michael B. Porter, Mar 16 2010
  • PARI
    is(n)=(valuation(n,2)+valuation(n,3))%2==0 \\ Charles R Greathouse IV, Sep 10 2015
  • PARI
    list(lim)=my(v=List(),N);for(n=0,logint(lim\=1,3),N=if(n%2,2*3^n,3^n); while(N<=lim, forstep(k=N,lim,[4*N,2*N], listput(v,k)); N<<=2)); Set(v) \\ Charles R Greathouse IV, Sep 10 2015
  • Python
    from itertools import count
def A036668(n):
    def f(x):
        c = n+x
        for i in range(x.bit_length()+1):
            i2 = 1<x:
                    break
                m = x//k
                c -= (m-1)//6+(m-5)//6+2
        return c
    m, k = n, f(n)
    while m != k: m, k = k, f(k)
    return m # Chai Wah Wu, Jan 28 2025

Formula

a(n) = 12/7 * n + O(log^2 n). - Charles R Greathouse IV, Sep 10 2015
{a(n)} = A052330({A014601(n)}), where {a(n)} denotes the set of integers in the sequence. - Peter Munn, Sep 14 2019

Extensions

Offset changed by Chai Wah Wu, Jan 28 2025

A334747 Let p be the smallest prime not dividing the squarefree part of n. Multiply n by p and divide by the product of all smaller primes.

Original entry on oeis.org

2, 3, 6, 8, 10, 5, 14, 12, 18, 15, 22, 24, 26, 21, 30, 32, 34, 27, 38, 40, 42, 33, 46, 20, 50, 39, 54, 56, 58, 7, 62, 48, 66, 51, 70, 72, 74, 57, 78, 60, 82, 35, 86, 88, 90, 69, 94, 96, 98, 75, 102, 104, 106, 45, 110, 84, 114, 87, 118, 120, 122, 93, 126, 128, 130, 55
Offset: 1

Views

Author

Peter Munn, May 09 2020

Keywords

Comments

A bijection from the positive integers to the nonsquares, A000037.
A003159 (which has asymptotic density 2/3) lists index n such that a(n) = 2n. The sequence maps the terms of A003159 1:1 onto A036554, defining a bijection between them.
Similarly, bijections are defined from A007417 to A325424, from A325424 to A145204\{0}, and from the first in each of the following pairs to the nonsquare integers in the second: (A145204\{0}, A036668), (A036668, A007417), (A036554, A003159), (A332820, A332821), (A332821, A332822), (A332822, A332820). Note that many of these are between sets where membership depends on whether a number's squarefree part divides by 2 and/or 3.
Starting from 1, and iterating the sequence as a(1) = 2, a(2) = 3, a(3) = 6, a(6) = 5, a(5) = 10, etc., runs through the squarefree numbers in the order they appear in A019565. - Antti Karttunen, Jun 08 2020

Examples

			168 = 42*4 has squarefree part 42 (and square part 4). The smallest prime absent from 42 = 2*3*7 is 5 and the product of all smaller primes is 2*3 = 6. So a(168) = 168*5/6 = 140.

Crossrefs

Permutation of A000037.
Row 2, and therefore column 2, of A331590. Cf. A334748 (row 3).
A007913, A034386, A053669, A225546 are used in formulas defining the sequence.
The formula section details how the sequence maps the terms of A002110, A003961, A019565; and how f(a(n)) relates to f(n) for f = A008833, A048675, A267116; making use of A003986.
Subsequences: A016825 (odd bisection), A036554, A329575.
Bijections are defined that relate to A003159, A007417, A036668, A145204, A325424, A332820, A332821, A332822.
Cf. also binary trees A334860, A334866 and A334870 (a left inverse).

Programs

  • PARI
    a(n) = {my(c=core(n), m=n); forprime(p=2, , if(c % p, m*=p; break, m/=p)); m;} \\ Michel Marcus, May 22 2020

Formula

a(n) = n * m / A034386(m-1), where m = A053669(A007913(n)).
a(n) = A331590(2, n) = A225546(2 * A225546(n)).
a(A019565(n)) = A019565(n+1).
a(k * m^2) = a(k) * m^2.
a(A003961(n)) = 2 * A003961(n).
a(2 * A003961(n)) = A003961(a(n)).
a(A002110(n)) = prime(n+1).
A048675(a(n)) = A048675(n) + 1.
A008833(a(n)) = A008833(n).
A267116(a(n)) = A267116(n) OR 1, where OR denotes the bitwise operation A003986.
a(A003159(n)) = A036554(n) = 2 * A003159(n).
A334870(a(n)) = n. - Antti Karttunen, Jun 08 2020

A329575 Numbers whose smallest Fermi-Dirac factor is 3.

Original entry on oeis.org

3, 12, 15, 21, 27, 33, 39, 48, 51, 57, 60, 69, 75, 84, 87, 93, 105, 108, 111, 123, 129, 132, 135, 141, 147, 156, 159, 165, 177, 183, 189, 192, 195, 201, 204, 213, 219, 228, 231, 237, 240, 243, 249, 255, 267, 273, 276, 285, 291, 297, 300, 303, 309, 321, 327, 336, 339, 345
Offset: 1

Views

Author

Peter Munn, Apr 27 2020

Keywords

Comments

Every positive integer is the product of a unique subset of the terms of A050376 (sometimes called Fermi-Dirac primes). This sequence lists the numbers where the relevant subset includes 3 but not 2.
Numbers whose squarefree part is divisible by 3 but not 2.
Positive multiples of 3 that survive sieving by the rule: if m appears then 2m, 3m and 6m do not. Asymptotic density is 1/6.

Examples

			6 is the product of the following terms of A050376: 2, 3. These terms include 2, so 6 is not in the sequence.
12 is the product of the following terms of A050376: 3, 4. These terms include 3, but not 2, so 12 is in the sequence.
20 is the product of the following terms of A050376: 4, 5. These terms do not include 3, so 20 is not in the sequence.

Links

Crossrefs

Cf. A007913, A036668, A050376, A059897, A223490, A307150, A339690.
Intersection of any 2 of A003159, A145204 and A325424; also subsequence of A028983.
Ordered 3rd quadrisection of A052330.

Programs

  • Mathematica
    f[p_, e_] := p^(2^IntegerExponent[e, 2]); fdmin[n_] := Min @@ f @@@ FactorInteger[n]; Select[Range[350], fdmin[#] == 3 &] (* Amiram Eldar, Nov 27 2020 *)
  • PARI
    isok(m) = core(m) % 6 == 3; \\ Michel Marcus, May 01 2020
  • Python
    from itertools import count
from sympy import integer_log
def A329575(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x):
        c = n+x
        for i in range(integer_log(x,9)[0]+1):
            i2 = 9**i
            for j in count(0,2):
                k = i2<x:
                    break
                m = x//k
                c -= (m-1)//6+(m-5)//6+2
        return c
    return 3*bisection(f,n,n) # Chai Wah Wu, Apr 10 2025

Formula

A223490(a(n)) = 3.
A007913(a(n)) == 3 (mod 6).
A059897(2, a(n)) = 2 * a(n).
A059897(3, a(n)) * 3 = a(n).
{a(n) : n >= 1} = {k : 3 * A307150(k) = 2 * k}.
A003159 = {a(n) / 3 : n >= 1} U {a(n) : n >= 1}.
A036668 = {a(n) / 3 : n >= 1} U {a(n) * 2 : n >= 1}.
A145204 \ {0} = {a(n) : n >= 1} U {a(n) * 2 : n >= 1}.
a(n) = 3*A339690(n). - Chai Wah Wu, Apr 10 2025

A225838 Numbers of form 2^i*3^j*(6k+5), i, j, k >= 0.

Original entry on oeis.org

5, 10, 11, 15, 17, 20, 22, 23, 29, 30, 33, 34, 35, 40, 41, 44, 45, 46, 47, 51, 53, 58, 59, 60, 65, 66, 68, 69, 70, 71, 77, 80, 82, 83, 87, 88, 89, 90, 92, 94, 95, 99, 101, 102, 105, 106, 107, 113, 116, 118, 119, 120, 123, 125, 130, 131, 132, 135, 136, 137, 138
Offset: 1

Views

Author

Ralf Stephan, May 16 2013

Keywords

Comments

Are a(n) > A225837(n) for all n? - Zak Seidov, May 17 2013
Yes. Imagine every 3-smooth number, m, visits you regularly, depositing a gold coin for safe keeping at each epoch (6k+1)*m and collecting it at epoch (6k+5)*m. If you run out of coins, you are doing something other than keeping them in a vault! - Peter Munn, Nov 13 2023
The asymptotic density of this sequence is 1/2. - Amiram Eldar, Apr 03 2022

Links

Crossrefs

Complement of A225837.
Symmetric difference of A003159 and A026225; of A189716 and A325424.

Programs

  • Magma
    [n: n in [1..200] | d mod 6 eq 5 where d is n div (2^Valuation(n,2)*3^Valuation(n,3))]; // Bruno Berselli, May 16 2013
  • Mathematica
    mx = 153; t = {}; Do[n = 2^i*3^j (6 k + 5); If[n <= mx, AppendTo[t, n]], {i, 0, Log[2, mx]}, {j, 0, Log[3, mx]}, {k, 0, mx/6}]; Union[t] (* T. D. Noe, May 16 2013 *)
  • PARI
    for(n=1,200,t=n/(2^valuation(n,2)*3^valuation(n,3));if((t%6==5),print1(n,",")))
  • Python
    from sympy import integer_log
def A225838(n):
    def f(x): return n+sum(((x//3**i>>j)+5)//6 for i in range(integer_log(x,3)[0]+1) for j in range((x//3**i).bit_length()))
    m, k = n, f(n)
    while m != k: m, k = k, f(k)
    return m # Chai Wah Wu, Feb 02 2025

A307150 Row 6 of array in A059897.

Original entry on oeis.org

6, 3, 2, 24, 30, 1, 42, 12, 54, 15, 66, 8, 78, 21, 10, 96, 102, 27, 114, 120, 14, 33, 138, 4, 150, 39, 18, 168, 174, 5, 186, 48, 22, 51, 210, 216, 222, 57, 26, 60, 246, 7, 258, 264, 270, 69, 282, 32, 294, 75, 34, 312, 318, 9, 330, 84, 38, 87, 354, 40, 366, 93
Offset: 1

Views

Author

N. J. A. Sloane, Mar 29 2019

Keywords

Comments

From Peter Munn, Apr 02 2019: (Start)
Also column 6 of A059897.
A self-inverse permutation of the positive integers with no fixed points; A073675 composed with A120229.
The permutation swaps pairs of integers whose ratio is 1:6 or 2:3, these ratios corresponding to the factorizations 1*6 = 2*3 = 6. Row 6 is the first row of A059897 to exhibit more than 1 such ratio.
(End)
The integers in the pairs with ratio 1:6 are listed in A036668, the integers in the pairs with ratio 2:3 are listed in A325424. - Peter Munn, Mar 05 2020

Links

Crossrefs

Cf. A036668, A059897, A073675, A120229, A325424.

Formula

From Peter Munn, Apr 02 2019: (Start)
a(n) = A059897(6,n) = A059897(n,6).
a(n) = A073675(A120229(n)) = A120229(A073675(n)) = A073675(n) * A120229(n) / n.
(End)

Extensions

More terms from Alois P. Heinz, Mar 31 2019

A334748 Let p be the smallest odd prime not dividing the squarefree part of n. Multiply n by p and divide by the product of all smaller odd primes.

Original entry on oeis.org

3, 6, 5, 12, 15, 10, 21, 24, 27, 30, 33, 20, 39, 42, 7, 48, 51, 54, 57, 60, 35, 66, 69, 40, 75, 78, 45, 84, 87, 14, 93, 96, 55, 102, 105, 108, 111, 114, 65, 120, 123, 70, 129, 132, 135, 138, 141, 80, 147, 150, 85, 156, 159, 90, 165, 168, 95, 174, 177, 28, 183, 186, 189
Offset: 1

Views

Author

Peter Munn, May 09 2020

Keywords

Comments

A permutation of A028983.
A007417 (which has asymptotic density 3/4) lists index n such that a(n) = 3n. The sequence maps the terms of A007417 1:1 onto A145204\{0}, defining a bijection between them.
Similarly, bijections are defined from the odd numbers (A005408) to the nonsquare odd numbers (A088828), from the positive even numbers (A299174) to A088829, from A003159 to the nonsquares in A003159, and from A325424 to the nonsquares in A036668. The latter two bijections are between sets where membership depends on whether a number's squarefree part divides by 2 and/or 3.

Examples

			84 = 21*4 has squarefree part 21 (and square part 4). The smallest odd prime absent from 21 = 3*7 is 5 and the product of all smaller odd primes is 3. So a(84) = 84*5/3 = 140.

Crossrefs

Permutation of A028983.
Row 3, and therefore column 3, of A331590. Cf. A334747 (row 2).
A007913, A034386, A225546, A284723 are used in formulas defining the sequence.
The formula section details how the sequence maps the terms of A003961, A019565, A070826; and how f(a(n)) relates to f(n) for f = A008833, A048675, A267116; making use of A003986.
Subsequences: A016051, A145204\{0}, A329575.
Bijections are defined that relate to A003159, A005408, A007417, A036668, A088828, A088829, A299174, A325424.

Programs

  • PARI
    a(n) = {my(c=core(n), m=n); forprime(p=3, , if(c % p, m*=p; break, m/=p)); m;} \\ Michel Marcus, May 22 2020

Formula

a(n) = n * p / (A034386(p-1)/2), where p = A284723(A007913(n)).
a(n) = A334747(A334747(n)).
a(n) = A331590(3, n) = A225546(4 * A225546(n)).
a(2*n) = 2 * a(n).
a(A019565(n)) = A019565(n+2).
a(k * m^2) = a(k) * m^2.
a(A003961(n)) = A003961(A334747(n)).
a(A070826(n)) = prime(n+1).
A048675(a(n)) = A048675(n) + 2.
A008833(a(n)) = A008833(n).
A267116(a(n)) = A267116(n) OR 1, where OR denotes the bitwise operation A003986.
a(A007417(n)) = A145204(n+1) = 3 * A007417(n).

A352273 Numbers whose squarefree part is congruent to 5 modulo 6.

Original entry on oeis.org

5, 11, 17, 20, 23, 29, 35, 41, 44, 45, 47, 53, 59, 65, 68, 71, 77, 80, 83, 89, 92, 95, 99, 101, 107, 113, 116, 119, 125, 131, 137, 140, 143, 149, 153, 155, 161, 164, 167, 173, 176, 179, 180, 185, 188, 191, 197, 203, 207, 209, 212, 215, 221, 227, 233, 236, 239, 245, 251
Offset: 1

Views

Author

Peter Munn, Mar 10 2022

Keywords

Comments

Numbers of the form 4^i * 9^j * (6k+5), i, j, k >= 0.
1/5 of each multiple of 5 in A352272.
The product of any two terms is in A352272.
The product of a term of this sequence and a term of A352272 is a term of this sequence.
The positive integers are usefully partitioned as {A352272, 2*A352272, 3*A352272, 6*A352272, {a(n)}, 2*{a(n)}, 3*{a(n)}, 6*{a(n)}}. There is a table in the example section giving sequences formed from unions of the parts.
The parts correspond to the cosets of A352272 considered as a subgroup of the positive integers under the operation A059897(.,.). Viewed another way, the parts correspond to the intersection of the integers with the cosets of the multiplicative subgroup of the positive rationals generated by the terms of A352272.
The asymptotic density of this sequence is 1/4. - Amiram Eldar, Apr 03 2022

Examples

			The squarefree part of 11 is 11, which is congruent to 5 (mod 6), so 11 is in the sequence.
The squarefree part of 15 is 15, which is congruent to 3 (mod 6), so 15 is not in the sequence.
The squarefree part of 20 = 2^2 * 5 is 5, which is congruent to 5 (mod 6), so 20 is in the sequence.
The table below lists OEIS sequences that are unions of the cosets described in the initial comments, and indicates the cosets included in each sequence. A352272 (as a subgroup) is denoted H, and this sequence (as a coset) is denoted H/5, in view of its terms being one fifth of the multiples of 5 in A352272.
             H    2H    3H    6H    H/5  2H/5  3H/5  6H/5
A003159      X           X           X           X
A036554            X           X           X           X
.
A007417      X     X                 X     X
A145204\{0}              X     X                 X     X
.
A026225      X           X                 X           X
A026179\{1}        X           X     X           X
.
A036668      X                 X     X                 X
A325424            X     X                 X     X
.
A055047      X                             X
A055048            X                 X
A055041                  X                             X
A055040                        X                 X
.
A189715      X                 X           X     X
A189716            X     X           X                 X
.
A225837      X     X     X     X
A225838                              X     X     X     X
.
A339690      X                       X
A329575                  X                       X
.
A352274      X           X
(The sequence groupings in the table start with the subgroup of the quotient group of H, followed by its cosets.)

Links

Crossrefs

Intersection of any three of A003159, A007417, A189716 and A225838.
Intersection of A036668 and A055048.
Complement within A339690 of A352272.
Cf. A007913, A059897, A307151, A334832.
Closure of A084088 under multiplication by 9.
Other subsequences: A033429\{0}, A016969.
Other sequences in the example table: A036554, A145204, A026179, A026225, A325424, A055040, A055041, A055047, A189715, A225837, A329575, A352274.

Programs

  • Mathematica
    q[n_] := Module[{e2, e3}, {e2, e3} = IntegerExponent[n, {2, 3}]; EvenQ[e2] && EvenQ[e3] && Mod[n/2^e2/3^e3, 6] == 5]; Select[Range[250], q] (* Amiram Eldar, Apr 03 2022 *)
  • PARI
    isok(m) = core(m) % 6 == 5;
  • Python
    from itertools import count
def A352273(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x):
        c = n+x
        for i in count(0):
            i2 = 9**i
            if i2>x: break
            for j in count(0,2):
                k = i2<x: break
                c -= (x//k-5)//6+1
        return c
    return bisection(f,n,n) # Chai Wah Wu, Feb 14 2025

Formula

{a(n) : n >= 1} = {m >= 1 : A007913(m) == 5 (mod 6)}.
{a(n) : n >= 1} = A334832/5 U A334832/11 U A334832/17 U A334832/23 where A334832/k denotes {A334832(m)/k : m >= 1, k divides A334832(m)}.
Using the same notation, {a(n) : n >= 1} = A352272/5 = {A307151(A352272(m)) : m >= 1}.
{A225838(n) : n >= 1} = {m : m = a(j)*k, j >= 1, k divides 6}.

A325498 Difference sequence of A036668.

Original entry on oeis.org

3, 1, 1, 1, 2, 2, 2, 3, 1, 2, 1, 3, 1, 1, 3, 1, 1, 1, 4, 1, 1, 4, 1, 1, 1, 1, 2, 2, 3, 1, 1, 1, 4, 2, 2, 1, 1, 1, 1, 1, 3, 2, 3, 1, 1, 1, 1, 1, 2, 2, 4, 2, 1, 3, 1, 1, 2, 1, 1, 1, 1, 4, 2, 3, 1, 1, 1, 1, 1, 2, 1, 1, 3, 1, 2, 4, 2, 4, 1, 1, 1, 3, 1, 1, 3, 1
Offset: 1

Views

Author

Clark Kimberling, May 05 2019

Keywords

Comments

See A325417 for a guide to related sequences.
Conjecture: every term is in {1,2,3,4}.

Examples

			A036668 is given by A(n) = least number not 2*A(m) or 3*A(m) for any m < n, so that A = (1,4,5,6,7,9,11,...), with differences (3,1,1,1,2,2,...).

Links

Crossrefs

Cf. A325417, A036668, A325499.

Programs

  • Mathematica
    a = {1}; Do[AppendTo[a, NestWhile[# + 1 &, Last[a] + 1,
Apply[Or, Map[MemberQ[a, #] &, Select[Flatten[{#/3, #/2}],
IntegerQ]]] &]], {2000}]; a ;       (* A036668 *)
c = Complement[Range[Last[a]], a] ; (* A325424 *)
Differences[a]  (* A325498 *)
Differences[c]  (* A325499 *)
(* Peter J. C. Moses, Apr 23 2019 *)

Formula

Asymptotic mean: lim_{n->oo} (1/n) * Sum_{k=1..n} a(k) = 12/7. - Amiram Eldar, Nov 26 2020
Showing 1-10 of 11 results. Next

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.