A086299 -id:A086299 - cp's OEIS frontend

A068191 Numbers n such that A067734(n)=0; complement of A002473; at least one prime-factor of n is larger than 7, it has 2 decimal digits.

11, 13, 17, 19, 22, 23, 26, 29, 31, 33, 34, 37, 38, 39, 41, 43, 44, 46, 47, 51, 52, 53, 55, 57, 58, 59, 61, 62, 65, 66, 67, 68, 69, 71, 73, 74, 76, 77, 78, 79, 82, 83, 85, 86, 87, 88, 89, 91, 92, 93, 94, 95, 97, 99, 101, 102, 103, 104, 106, 107, 109, 110, 111, 113, 114

Offset: 1

Labos Elemer, Feb 19 2002

Also numbers n such that A198487(n) = 0 and A107698(n) = 0. - Jaroslav Krizek, Nov 04 2011
A086299(a(n)) = 0. - Reinhard Zumkeller, Apr 01 2012
A262401(a(n)) < a(n). - Reinhard Zumkeller, Sep 25 2015
Numbers not in A007954. - Mohammed Yaseen, Sep 13 2022

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

Cf. A001222, A002473, A067734.
Cf. A068183, A068184, A068185, A068186, A068187, A068189, A068190.
Cf. A262401.

import Data.List (elemIndices)
a068191 n = a068191_list !! (n-1)
a068191_list = map (+ 1) $ elemIndices 0 a086299_list
-- Reinhard Zumkeller, Apr 01 2012

Select[Range@120, Last@Map[First, FactorInteger@#] > 7 &] (* Vincenzo Librandi, Sep 19 2016 *)

from sympy import integer_log
def A068191(n):
    def f(x):
        c = n
        for i in range(integer_log(x,7)[0]+1):
            i7 = 7**i
            m = x//i7
            for j in range(integer_log(m,5)[0]+1):
                j5 = 5**j
                r = m//j5
                for k in range(integer_log(r,3)[0]+1):
                    c += (r//3**k).bit_length()
        return c
    m, k = n, f(n)
    while m != k: m, k = k, f(k)
    return m # Chai Wah Wu, Sep 16 2024

A238985 Zeroless 7-smooth numbers.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 12, 14, 15, 16, 18, 21, 24, 25, 27, 28, 32, 35, 36, 42, 45, 48, 49, 54, 56, 63, 64, 72, 75, 81, 84, 96, 98, 112, 125, 126, 128, 135, 144, 147, 162, 168, 175, 189, 192, 196, 216, 224, 225, 243, 245, 252, 256, 288, 294, 315, 324, 336

Offset: 1

Charles R Greathouse IV and Reinhard Zumkeller, Mar 07 2014

A001221(a(n)) <= 3 since 10 cannot divide a(n).

It seems that this sequence is finite and contains 12615 terms. - Daniel Mondot, May 03 2022 and Jianing Song, Jan 28 2023

			a(12615) = 2^25 * 3^227 * 7^28.

Daniel Mondot, Table of n, a(n) for n = 1..12615 (terms 1..10000 from Charles R Greathouse IV)

Cf. A168046, intersection of A002473 and A052382.
A238938, A238939, A238940, A195948, A238936, A195908 are proper subsequences.
Cf. A059405 (subsequence), A350180 through A350187.

import Data.Set (singleton, deleteFindMin, fromList, union)
a238985 n = a238985_list !! (n-1)
a238985_list = filter ((== 1) . a168046) $ f $ singleton 1 where
   f s = x : f (s' `union` fromList
               (filter ((> 0) . (`mod` 10)) $ map (* x) [2,3,5,7]))
               where (x, s') = deleteFindMin s

zf(n)=vecmin(digits(n))
list(lim)=my(v=List(),t,t1); for(e=0,log(lim+1)\log(7), t1=7^e; for(f=0,log(lim\t1+1)\log(3), t=t1*3^f; while(t<=lim, if(zf(t), listput(v, t)); t<<=1)); for(f=0,log(lim\t1+1)\log(5), t=t1*5^f; while(t<=lim, if(zf(t), listput(v, t)); t*=3))); Set(v)

A086299(a(n)) * A168046(a(n)) = 1.

Keyword:fini and keyword:full removed by Jianing Song, Jan 28 2023 as finiteness is only conjectured.

A071604 a(n) is the number of 7-smooth numbers <= n.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 10, 11, 11, 12, 13, 14, 14, 15, 15, 16, 17, 17, 17, 18, 19, 19, 20, 21, 21, 22, 22, 23, 23, 23, 24, 25, 25, 25, 25, 26, 26, 27, 27, 27, 28, 28, 28, 29, 30, 31, 31, 31, 31, 32, 32, 33, 33, 33, 33, 34, 34, 34, 35, 36, 36, 36, 36, 36, 36, 37, 37, 38

Offset: 1

Benoit Cloitre, Jun 02 2002

A 7-smooth number is a number of the form 2^x*3^y*5^z*7^u, (x,y,z,u) >= 0.

In other words, a 7-smooth number is a number with no prime factor greater than 7. - Peter Munn, Nov 20 2021

			a(11) = 10 as there are 10 7-smooth numbers <= 11. Namely 1, 2, 3, 4, 5, 6, 7, 8, 9, 10. - _David A. Corneth_, Apr 19 2021

David A. Corneth, Table of n, a(n) for n = 1..10000

Partial sums of A086299.
Column 7 of A080786.
Cf. A002473, A106600, A343598.
Equivalent sequences with other limits on greatest prime factor: A070939 (2), A071521 (3), A071520 (5), A071523 (11), A080684 (13), A080685 (17), A080686 (19), A096300 (log n).

for(n=1,100,print1(sum(k=1,n,if(sum(i=5,n,if(k%prime(i),0,1)),0,1)),","))

from sympy import integer_log
def A071604(n):
    c = 0
    for i in range(integer_log(n,7)[0]+1):
        i7 = 7**i
        m = n//i7
        for j in range(integer_log(m,5)[0]+1):
            j5 = 5**j
            r = m//j5
            for k in range(integer_log(r,3)[0]+1):
                c += (r//3**k).bit_length()
    return c # Chai Wah Wu, Sep 16 2024

a(n) = Card{ k | A002473 (k) <= n }.

Name corrected by David A. Corneth, Apr 19 2021

cp's OEIS Frontend

A068191 Numbers n such that A067734(n)=0; complement of A002473; at least one prime-factor of n is larger than 7, it has 2 decimal digits.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Mathematica

Python

A238985 Zeroless 7-smooth numbers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

PARI

Formula

Extensions

A071604 a(n) is the number of 7-smooth numbers <= n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

Python

Formula

Extensions