A003596 -id:A003596 - cp's OEIS frontend

A107788 Numbers of the form (8^i)*(11^j), with i, j >= 0.

1, 8, 11, 64, 88, 121, 512, 704, 968, 1331, 4096, 5632, 7744, 10648, 14641, 32768, 45056, 61952, 85184, 117128, 161051, 262144, 360448, 495616, 681472, 937024, 1288408, 1771561, 2097152, 2883584, 3964928, 5451776, 7496192, 10307264

Offset: 1

Douglas Winston (douglas.winston(AT)srupc.com), Jun 14 2005

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

Cf. A025633, A025634, A107764, A003596, A003597, A107988, A003598, A003599, A108090.

import Data.Set (singleton, deleteFindMin, insert)
a107788 n = a107788_list !! (n-1)
a107788_list = f $ singleton (1,0,0) where
   f s = y : f (insert (8 * y, i + 1, j) $ insert (11 * y, i, j + 1) s')
         where ((y, i, j), s') = deleteFindMin s
-- Reinhard Zumkeller, May 15 2015

Take[Union[8^First[#]*11^Last[#]&/@Tuples[Range[0,20],2]],40] (* Harvey P. Dale, Jan 17 2015 *)
n = 10^6; Flatten[Table[8^i*11^j, {i, 0, Log[8, n]}, {j, 0, Log[11, n/8^i]}]] // Sort (* Amiram Eldar, Oct 07 2020 *)

Sum_{n>=1} 1/a(n) = (8*11)/((8-1)*(11-1)) = 44/35. - Amiram Eldar, Oct 07 2020

a(n) ~ exp(sqrt(2*log(8)*log(11)*n)) / sqrt(88). - Vaclav Kotesovec, Oct 07 2020

A108687 Numbers of the form (9^i)*(11^j), with i, j >= 0.

Original entry on oeis.org

1, 9, 11, 81, 99, 121, 729, 891, 1089, 1331, 6561, 8019, 9801, 11979, 14641, 59049, 72171, 88209, 107811, 131769, 161051, 531441, 649539, 793881, 970299, 1185921, 1449459, 1771561, 4782969, 5845851, 7144929, 8732691, 10673289, 13045131

Offset: 1

Douglas Winston (douglas.winston(AT)srupc.com), Jun 17 2005

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

Subsequence of A003597.

Cf. A025633, A025634, A107788, A107764, A003596, A107988, A003598, A003599, A108090.

import Data.Set (singleton, deleteFindMin, insert)
a108687 n = a108687_list !! (n-1)
a108687_list = f $ singleton (1,0,0) where
   f s = y : f (insert (9 * y, i + 1, j) $ insert (11 * y, i, j + 1) s')
         where ((y, i, j), s') = deleteFindMin s
-- Reinhard Zumkeller, May 15 2015

f[upto_]:=With[{max9=Floor[Log[9,upto]],max11=Floor[Log[11,upto]]}, Select[Union[Times@@{9^First[#],11^Last[#]}&/@Tuples[{Range[0, max9], Range[0, max11]}]], #<=upto&]]; f[14000000]  (* Harvey P. Dale, Mar 11 2011 *)

from sympy import integer_log
def A108687(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): return n+x-sum(integer_log(x//11**i,9)[0]+1 for i in range(integer_log(x,11)[0]+1))
    return bisection(f,n,n) # Chai Wah Wu, Mar 25 2025

Sum_{n>=1} 1/a(n) = (9*11)/((9-1)*(11-1)) = 99/80. - Amiram Eldar, Sep 24 2020

a(n) ~ exp(sqrt(2*log(9)*log(11)*n)) / sqrt(99). - Vaclav Kotesovec, Sep 24 2020

A025635 Numbers of form 9^i*10^j, with i, j >= 0.

Original entry on oeis.org

1, 9, 10, 81, 90, 100, 729, 810, 900, 1000, 6561, 7290, 8100, 9000, 10000, 59049, 65610, 72900, 81000, 90000, 100000, 531441, 590490, 656100, 729000, 810000, 900000, 1000000, 4782969, 5314410, 5904900, 6561000, 7290000, 8100000, 9000000

Offset: 1

David W. Wilson

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

Cf. A025612, A025616, A025621, A025625, A025629, A025632, A025634, A108761, A003596, A003597, A107988, A003598, A108698, A003599, A107788, A108687, A108779, A108090.

import Data.Set (singleton, deleteFindMin, insert)
a025635 n = a025635_list !! (n-1)
a025635_list = f $ singleton (1,0,0) where
   f s = y : f (insert (9 * y, i + 1, j) $ insert (10 * y, i, j + 1) s')
         where ((y, i, j), s') = deleteFindMin s
-- Reinhard Zumkeller, May 15 2015

With[{max = 10^7}, Flatten[Table[9^i*10^j, {i, 0, Log[9, max]}, {j, 0, Log[10, max/9^i]}]] // Sort] (* Amiram Eldar, Mar 29 2025 *)

list(lim)=my(v=List(), N); for(n=0, logint(lim\=1, 10), N=10^n; while(N<=lim, listput(v, N); N*=9)); Set(v) \\ Charles R Greathouse IV, Jan 10 2018

from sympy import integer_log
def A025635(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): return n+x-sum(integer_log(x//10**i,9)[0]+1 for i in range(integer_log(x,10)[0]+1))
    return bisection(f,n,n) # Chai Wah Wu, Mar 25 2025

Sum_{n>=1} 1/a(n) = 5/4. - Amiram Eldar, Mar 29 2025

a(n) = 9^A025683(n) * 10^A025691(n). - R. J. Mathar, Jul 06 2025

A107466 Numbers of the form (5^i)*(13^j).

Original entry on oeis.org

1, 5, 13, 25, 65, 125, 169, 325, 625, 845, 1625, 2197, 3125, 4225, 8125, 10985, 15625, 21125, 28561, 40625, 54925, 78125, 105625, 142805, 203125, 274625, 371293, 390625, 528125, 714025, 1015625, 1373125, 1856465, 1953125, 2640625

Offset: 1

Douglas Winston (douglas.winston(AT)srupc.com), May 27 2005

nonn

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A003586, A003592, A003593, A003591, A003594, A003595, A003596, A003597, A003598, A003599, A107326, A107364.

mx = 2700000; Sort@ Flatten@ Table[5^i*13^j, {i, 0, Log[5, mx]}, {j, 0, Log[13, mx/5^i]}] (* Robert G. Wilson v, Aug 17 2012 *)

list(lim)=my(v=List(),N);for(n=0,log(lim)\log(13),N=13^n;while(N<=lim,listput(v,N);N*=5));vecsort(Vec(v)) \\ Charles R Greathouse IV, Jun 28 2011

from sympy import integer_log
def A107466(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): return n+x-sum(integer_log(x//13**i,5)[0]+1 for i in range(integer_log(x,13)[0]+1))
    return bisection(f,n,n) # Chai Wah Wu, Mar 25 2025

Sum_{n>=1} 1/a(n) = (5*13)/((5-1)*(13-1)) = 65/48. - Amiram Eldar, Sep 23 2020

a(n) ~ exp(sqrt(2*log(5)*log(13)*n)) / sqrt(65). - Vaclav Kotesovec, Sep 23 2020

A108698 Numbers of the form (6^i)*(11^j), with i, j >= 0.

Original entry on oeis.org

1, 6, 11, 36, 66, 121, 216, 396, 726, 1296, 1331, 2376, 4356, 7776, 7986, 14256, 14641, 26136, 46656, 47916, 85536, 87846, 156816, 161051, 279936, 287496, 513216, 527076, 940896, 966306, 1679616, 1724976, 1771561, 3079296, 3162456

Offset: 1

Douglas Winston (douglas.winston(AT)srupc.com), Jun 19 2005

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

Cf. A025626, A025627, A025628, A025629, A064476, A107710, A003596, A003597, A107988, A003598, A108687, A003599, A107788, A108090.

import Data.Set (singleton, deleteFindMin, insert)
a108698 n = a108698_list !! (n-1)
a108698_list = f $ singleton (1,0,0) where
   f s = y : f (insert (6 * y, i + 1, j) $ insert (11 * y, i, j + 1) s')
         where ((y, i, j), s') = deleteFindMin s
-- Reinhard Zumkeller, May 15 2015

n = 10^6; Flatten[Table[6^i*11^j, {i, 0, Log[6, n]}, {j, 0, Log[11, n/6^i]}]] // Sort (* Amiram Eldar, Oct 07 2020 *)

Sum_{n>=1} 1/a(n) = (6*11)/((6-1)*(11-1)) = 33/25. - Amiram Eldar, Oct 07 2020

a(n) ~ exp(sqrt(2*log(6)*log(11)*n)) / sqrt(66). - Vaclav Kotesovec, Oct 07 2020

A108056 Numbers of the form (7^i)*(13^j).

Original entry on oeis.org

1, 7, 13, 49, 91, 169, 343, 637, 1183, 2197, 2401, 4459, 8281, 15379, 16807, 28561, 31213, 57967, 107653, 117649, 199927, 218491, 371293, 405769, 753571, 823543, 1399489, 1529437, 2599051, 2840383, 4826809, 5274997, 5764801, 9796423, 10706059, 18193357, 19882681

Offset: 1

Douglas Winston (douglas.winston(AT)srupc.com), Jun 02 2005

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A003586, A003592, A003593, A003591, A003594, A003595, A003596, A003597, A003598, A003599, A107326, A107364, A107466.

n = 10^7; Flatten[Table[7^i*13^j, {i, 0, Log[7, n]}, {j, 0, Log[13, n/7^i]}]] // Sort (* Amiram Eldar, Sep 23 2020 *)

list(lim)=my(v=List(),N);for(n=0,log(lim)\log(13),N=13^n;while(N<=lim,listput(v,N);N*=7));vecsort(Vec(v)) \\ Charles R Greathouse IV, Jun 28 2011

from sympy import integer_log
def A108056(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return n+x-sum(integer_log(x//13**i,7)[0]+1 for i in range(integer_log(x,13)[0]+1))
    return bisection(f,n,n) # Chai Wah Wu, Oct 22 2024

Sum_{n>=1} 1/a(n) = (7*13)/((7-1)*(13-1)) = 91/72. - Amiram Eldar, Sep 23 2020

a(n) ~ exp(sqrt(2*log(7)*log(13)*n)) / sqrt(91). - Vaclav Kotesovec, Sep 23 2020

More terms from Amiram Eldar, Sep 23 2020

A108779 Numbers of the form (10^i)*(11^j), with i, j >= 0.

Original entry on oeis.org

1, 10, 11, 100, 110, 121, 1000, 1100, 1210, 1331, 10000, 11000, 12100, 13310, 14641, 100000, 110000, 121000, 133100, 146410, 161051, 1000000, 1100000, 1210000, 1331000, 1464100, 1610510, 1771561, 10000000, 11000000, 12100000, 13310000

Offset: 1

Douglas Winston (douglas.winston(AT)srupc.com), Jun 26 2005

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

Cf. A025612, A025616, A025621, A025625, A025629, A025632, A025634, A025635, A108761, A003596, A003597, A107988, A003598, A108698, A003599, A107788, A108687, A108090.

import Data.Set (singleton, deleteFindMin, insert)
a108779 n = a108779_list !! (n-1)
a108779_list = f $ singleton (1,0,0) where
   f s = y : f (insert (10 * y, i + 1, j) $ insert (11 * y, i, j + 1) s')
         where ((y, i, j), s') = deleteFindMin s
-- Reinhard Zumkeller, May 15 2015

n = 10^7; Flatten[Table[10^i*11^j, {i, 0, Log10[n]}, {j, 0, Log[11, n/10^i]}]] // Sort (* Amiram Eldar, Sep 25 2020 *)

Sum_{n>=1} 1/a(n) = (10*11)/((10-1)*(11-1)) = 11/9. - Amiram Eldar, Sep 25 2020

a(n) ~ exp(sqrt(2*log(10)*log(11)*n)) / sqrt(110). - Vaclav Kotesovec, Sep 25 2020

A317804 Numbers of form 2^i*12^j, with i, j >= 0.

Original entry on oeis.org

1, 2, 4, 8, 12, 16, 24, 32, 48, 64, 96, 128, 144, 192, 256, 288, 384, 512, 576, 768, 1024, 1152, 1536, 1728, 2048, 2304, 3072, 3456, 4096, 4608, 6144, 6912, 8192, 9216, 12288, 13824, 16384, 18432, 20736, 24576, 27648, 32768, 36864, 41472, 49152, 55296, 65536

Offset: 1

Dario Ch, Sep 01 2018

Dario Ch, Table of n, a(n) for n = 1..10000

Cf. A025612, A003596, A107326, A003597, A107364, A025616, A108201, A108238.

With[{max = 10^5}, Flatten[Table[2^i*12^j, {i, 0, Log2[max]}, {j, 0, Log[12, max/2^i]}]] // Sort] (* Amiram Eldar, Mar 29 2025 *)

from heapq import heappush, heappop
def sequence():
    pq = [1]
    seen = set(pq)
    while True:
        value = heappop(pq)
        yield value
        seen.remove(value)
        for x in 2 * value, 12 * value:
            if x not in seen:
                heappush(pq, x)
                seen.add(x)
seq = sequence()
finalsequence_list = [next(seq) for i in range(100)]  # Dario Ch, Sep 01 2018

from sympy import integer_log
def A317804(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): return n+x-sum((x//12**i).bit_length() for i in range(integer_log(x,12)[0]+1))
    return bisection(f,n,n) # Chai Wah Wu, Mar 26 2025

Sum_{n>=1} 1/a(n) = 24/11. - Amiram Eldar, Mar 29 2025

A108218 Numbers of the form (11^i)*(12^j), with i, j >= 0.

Original entry on oeis.org

1, 11, 12, 121, 132, 144, 1331, 1452, 1584, 1728, 14641, 15972, 17424, 19008, 20736, 161051, 175692, 191664, 209088, 228096, 248832, 1771561, 1932612, 2108304, 2299968, 2509056, 2737152, 2985984, 19487171, 21258732, 23191344, 25299648

Offset: 1

Douglas Winston (douglas.winston(AT)srupc.com), Jun 28 2005

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

Cf. A003596, A003597, A107988, A003598, A108201, A108698, A064476, A003599, A108238, A107788, A108687, A108779, A108090, A108771.

import Data.Set (singleton, deleteFindMin, insert)
a108218 n = a108218_list !! (n-1)
a108218_list = f $ singleton (1,0,0) where
   f s = y : f (insert (11 * y, i + 1, j) $ insert (12 * y, i, j + 1) s')
         where ((y, i, j), s') = deleteFindMin s
-- Reinhard Zumkeller, May 15 2015

With[{max = 3*10^7}, Flatten[Table[11^i*12^j, {i, 0, Log[11, max]}, {j, 0, Log[12, max/11^i]}]] // Sort] (* Amiram Eldar, Mar 29 2025 *)

Sum_{n>=1} 1/a(n) = 6/5. - Amiram Eldar, Mar 29 2025

A340268 Composite numbers k>1 such that (s-1) | (d-1) for each d | k, where s = lpf(k) = A020639(k).

Original entry on oeis.org

4, 6, 8, 9, 10, 12, 14, 15, 16, 18, 20, 21, 22, 24, 25, 26, 27, 28, 30, 32, 33, 34, 36, 38, 39, 40, 42, 44, 45, 46, 48, 49, 50, 51, 52, 54, 56, 57, 58, 60, 62, 63, 64, 65, 66, 68, 69, 70, 72, 74, 75, 76, 78, 80, 81, 82, 84, 85, 86, 87, 88, 90, 91, 92, 93, 94, 96

Offset: 1

Maxim Karimov, Jan 02 2021

nonn

Not a duplicate of A340058 because the complements A335902 and A340269 differ. - R. J. Mathar, Feb 16 2021

Cf. A000010, A000961, A020639, A340058, A335902, A340269 (complement).

Contains all composite terms of at least A003586, A003591, A003592, A003593, A003596.

n=300; % gives all terms of the sequence not exceeding n
A=[];
for i=2:n
    lpf=2;
    while mod(i,lpf)~=0
        lpf=lpf+1;
    end
    for d=1:floor(i/2)
        if mod(i,d)==0 && mod(d-1,lpf-1)~=0
            break
        elseif d==floor(i/2)
            A=[A i];
        end
    end
end

with(numtheory):
q:= n-> (f-> andmap(d-> irem(d-1, f)=0, divisors(n)))(min(factorset(n))-1):
select(not isprime and q, [$2..96])[];  # Alois P. Heinz, Feb 12 2021

Select[Range[2, 96], Function[{n, s}, And[! PrimeQ@ n, AllTrue[Divisors[n] - 1, Mod[#, s] == 0 &]]] @@ {#, FactorInteger[#][[1, 1]] - 1} &] (* Michael De Vlieger, Feb 12 2021 *)

isok(c) = if ((c>1) && !isprime(c), my(f=factor(c)[,1]); for (k=1, #f~, if ((f[k]-1) % (f[1]-1), return(0))); return(1)); \\ Michel Marcus, Jan 03 2021

cp's OEIS Frontend

A107788 Numbers of the form (8^i)*(11^j), with i, j >= 0.

Views

Author

Keywords

Links

Crossrefs

Programs

Haskell

Mathematica

Formula

A108687 Numbers of the form (9^i)*(11^j), with i, j >= 0.

Views

Author

Keywords

Links

Crossrefs

Programs

Haskell

Mathematica

Python

Formula

A025635 Numbers of form 9^i*10^j, with i, j >= 0.

Views

Author

Keywords

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Python

Formula

A107466 Numbers of the form (5^i)*(13^j).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A108698 Numbers of the form (6^i)*(11^j), with i, j >= 0.

Views

Author

Keywords

Links

Crossrefs

Programs

Haskell

Mathematica

Formula

A108056 Numbers of the form (7^i)*(13^j).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

Extensions

A108779 Numbers of the form (10^i)*(11^j), with i, j >= 0.

Views

Author

Keywords

Links

Crossrefs

Programs

Haskell

Mathematica

Formula

A317804 Numbers of form 2^i*12^j, with i, j >= 0.

Views

Author