A143202 -id:A143202 - cp's OEIS frontend

A003593 Numbers of the form 3^i*5^j with i, j >= 0.

1, 3, 5, 9, 15, 25, 27, 45, 75, 81, 125, 135, 225, 243, 375, 405, 625, 675, 729, 1125, 1215, 1875, 2025, 2187, 3125, 3375, 3645, 5625, 6075, 6561, 9375, 10125, 10935, 15625, 16875, 18225, 19683, 28125, 30375, 32805, 46875, 50625, 54675, 59049

Offset: 1

N. J. A. Sloane

nonn

Odd 5-smooth numbers (A051037). - Reinhard Zumkeller, Sep 18 2005

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

Cf. A033849, A112751-A112756, A143202, A022337 (list of j), A022336(list of i).

Cf. A264997 (partitions into), see also A264998. Cf. A108347 (odd 7-smooth).

Filtered([1..60000],n->PowerMod(15,n,n)=0); # Muniru A Asiru, Mar 19 2019

import Data.Set (singleton, deleteFindMin, insert)
a003593 n = a003593_list !! (n-1)
a003593_list = f (singleton 1) where
   f s = m : f (insert (3*m) $ insert (5*m) s') where
     (m,s') = deleteFindMin s
-- Reinhard Zumkeller, Sep 13 2011

[n: n in [1..60000] | PrimeDivisors(n) subset [3,5]]; // Bruno Berselli, Sep 24 2012

isA003593 := proc(n)
    if n = 1 then
        true;
    else
        return (numtheory[factorset](n) minus {3, 5} = {} );
    end if;
end proc:
A003593 := proc(n)
    option remember;
    if n = 1 then
        1;
    else
        for a from procname(n-1)+1 do
            if isA003593(a) then
                return a;
            end if;
        end do:
    end if;
end proc:
seq(A003593(n),n=1..30) ; # R. J. Mathar, Aug 04 2016

fQ[n_] := PowerMod[15, n, n] == 0; Select[Range[60000], fQ] (* Bruno Berselli, Sep 24 2012 *)

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

is(n)=n==3^valuation(n,3)*5^valuation(n,5) \\ Charles R Greathouse IV, Apr 23 2013

from sympy import integer_log
def A003593(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//5**i,3)[0]+1 for i in range(integer_log(x,5)[0]+1))
    return bisection(f,n,n) # Chai Wah Wu, Oct 22 2024

a(n) ~ 1/sqrt(15)*exp(sqrt(2*log(3)*log(5)*n)) asymptotically. - Benoit Cloitre, Jan 22 2002
The characteristic function of this sequence is given by Sum_{n >= 1} x^a(n) = Sum_{n >= 1} mu(15*n)*x^n/(1 - x^n), where mu(n) is the Möbius function A008683. Cf. with the formula of Hanna in A051037. - Peter Bala, Mar 18 2019
Sum_{n>=1} 1/a(n) = (3*5)/((3-1)*(5-1)) = 15/8. - Amiram Eldar, Sep 22 2020

A033849 Numbers whose prime factors are 3 and 5.

Original entry on oeis.org

15, 45, 75, 135, 225, 375, 405, 675, 1125, 1215, 1875, 2025, 3375, 3645, 5625, 6075, 9375, 10125, 10935, 16875, 18225, 28125, 30375, 32805, 46875, 50625, 54675, 84375, 91125, 98415, 140625, 151875, 164025, 234375, 253125, 273375, 295245

Offset: 1

Jeff Burch

nonn

Numbers k such that phi(k) = (8/15)*k. - Benoit Cloitre, Apr 19 2002

Subsequence of A143202. - Reinhard Zumkeller, Sep 13 2011

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

Cf. A033845, A033846, A033847, A033848, A033850, A033851, A143201, A143202.

Subsequence of A256617.

import Data.Set (singleton, deleteFindMin, insert)
a033849 n = a033849_list !! (n-1)
a033849_list = f (singleton (3*5)) where
   f s = m : f (insert (3*m) $ insert (5*m) s') where
     (m,s') = deleteFindMin s
-- Reinhard Zumkeller, Sep 13 2011

Sort[Flatten[Table[Table[3^j*5^k, {j, 1, 10}], {k, 1, 10}]]] (* Geoffrey Critzer, Dec 07 2014 *)
Select[Range[300000],FactorInteger[#][[All,1]]=={3,5}&] (* Harvey P. Dale, Oct 19 2022 *)

from sympy import integer_log
def A033849(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//5**i,3)[0]+1 for i in range(integer_log(x,5)[0]+1))
    return 15*bisection(f,n,n) # Chai Wah Wu, Oct 22 2024

From Reinhard Zumkeller, Sep 13 2011: (Start)
A143201(a(n)) = 3.
a(n) = 15*A003593(n). (End)
Sum_{n>=1} 1/a(n) = 1/8. - Amiram Eldar, Dec 22 2020

Offset and typo in data fixed by Reinhard Zumkeller, Sep 13 2011

A143201 Product of distances between prime factors in factorization of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 4, 1, 2, 1, 6, 3, 1, 1, 2, 1, 4, 5, 10, 1, 2, 1, 12, 1, 6, 1, 6, 1, 1, 9, 16, 3, 2, 1, 18, 11, 4, 1, 10, 1, 10, 3, 22, 1, 2, 1, 4, 15, 12, 1, 2, 7, 6, 17, 28, 1, 6, 1, 30, 5, 1, 9, 18, 1, 16, 21, 12, 1, 2, 1, 36, 3, 18, 5, 22, 1, 4, 1, 40, 1, 10, 13, 42, 27, 10, 1, 6, 7, 22

Offset: 1

Reinhard Zumkeller, Aug 12 2008

nonn

a(n) is the product of the sum of 1 and first differences of prime factors of n with multiplicity, with a(n) = 1 for n = 1 or prime n. - Michael De Vlieger, Nov 12 2023.
a(A007947(n)) = a(n);
A006093 and A001747 give record values and where they occur:
A006093(n)=a(A001747(n+1)) for n>1.
a(n) = 1 iff n is a prime power: a(A000961(n))=1;
a(n) = 2 iff n has exactly 2 and 3 as prime factors:
a(A033845(n))=2;
a(n) = 3 iff n is in A143202;
a(n) = 4 iff n has exactly 2 and 5 as prime factors:
a(A033846(n))=4;
a(n) = 5 iff n is in A143203;
a(n) = 6 iff n is in A143204;
a(n) = 7 iff n is in A143205;
a(n) <> A006512(k)+1 for k>1.
a(A033849(n))=3; a(A033851(n))=3; a(A033850(n))=5; a(A033847(n))=6; a(A033848(n))=10. [Reinhard Zumkeller, Sep 19 2011]

			a(86) = a(43*2) = 43-2+1 = 42;
a(138) = a(23*3*2) = (23-3+1)*(3-2+1) = 42;
a(172) = a(43*2*2) = (43-2+1)*(2-2+1) = 42;
a(182) = a(13*7*2) = (13-7+1)*(7-2+1) = 42;
a(276) = a(23*3*2*2) = (23-3+1)*(3-2+1)*(2-2+1) = 42;
a(330) = a(11*5*3*2) = (11-5+1)*(5-3+1)*(3-2+1) = 42.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Index entries for primes, gaps between

Cf. A000961, A001747, A006093, A020639, A033845, A033846.
Cf. A033847, A033848, A033849, A033850.
Cf. A143203, A143204, A143205.
Cf. A027746.

a143201 1 = 1
a143201 n = product $ map (+ 1) $ zipWith (-) (tail pfs) pfs
   where pfs = a027748_row n
-- Reinhard Zumkeller, Sep 13 2011

Table[Times@@(Differences[Flatten[Table[First[#],{Last[#]}]&/@ FactorInteger[ n]]]+1),{n,100}] (* Harvey P. Dale, Dec 07 2011 *)

a(n) = f(n,1,1) where f(n,q,y) = if n=1 then y else if q=1 then f(n/p,p,1)) else f(n/p,p,y*(p-q+1)) with p = A020639(n) = smallest prime factor of n.

A033851 Numbers whose prime factors are 5 and 7.

Original entry on oeis.org

35, 175, 245, 875, 1225, 1715, 4375, 6125, 8575, 12005, 21875, 30625, 42875, 60025, 84035, 109375, 153125, 214375, 300125, 420175, 546875, 588245, 765625, 1071875, 1500625, 2100875, 2734375, 2941225, 3828125, 4117715, 5359375, 7503125

Offset: 1

Jeff Burch

nonn

Numbers k such that phi(k)/k == 24/35. - Artur Jasinski, Nov 09 2008

Subsequence of A143202. - Reinhard Zumkeller, Sep 13 2011

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

Cf. A003595, A147571-A147575, A147576-A147580. [Artur Jasinski, Nov 09 2008]
Cf. A033845, A033846, A033847, A033848, A033849, A033850, A143201, A143202.
Subsequence of A256617.

import Data.Set (singleton, deleteFindMin, insert)
a033851 n = a033851_list !! (n-1)
a033851_list = f (singleton (5*7)) where
   f s = m : f (insert (5*m) $ insert (7*m) s') where
     (m,s') = deleteFindMin s
-- Reinhard Zumkeller, Sep 13 2011

a = {}; Do[If[EulerPhi[x]/x == 24/35, AppendTo[a, x]], {x, 1, 10000}]; a (* Artur Jasinski, Nov 09 2008 *)
Take[With[{nn=10},Sort[Flatten[Table[5^i 7^j,{i,nn},{j,nn}]]]],40] (* Harvey P. Dale, Feb 09 2013 *)

a(n) = 35*A003595(n). - Artur Jasinski, Nov 09 2008
A143201(a(n)) = 3. - Reinhard Zumkeller, Sep 13 2011
Sum_{n>=1} 1/a(n) = 1/24. - Amiram Eldar, Dec 22 2020

Offset fixed by Reinhard Zumkeller, Sep 13 2011

cp's OEIS Frontend

A003593 Numbers of the form 3^i*5^j with i, j >= 0.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

PARI

PARI

Python

Formula

A033849 Numbers whose prime factors are 3 and 5.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Mathematica

Python

Formula

Extensions

A143201 Product of distances between prime factors in factorization of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

Formula

A033851 Numbers whose prime factors are 5 and 7.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Mathematica

Formula

Extensions