A147576 -id:A147576 - cp's OEIS frontend

A087207 A binary representation of the primes that divide a number, shown in decimal.

0, 1, 2, 1, 4, 3, 8, 1, 2, 5, 16, 3, 32, 9, 6, 1, 64, 3, 128, 5, 10, 17, 256, 3, 4, 33, 2, 9, 512, 7, 1024, 1, 18, 65, 12, 3, 2048, 129, 34, 5, 4096, 11, 8192, 17, 6, 257, 16384, 3, 8, 5, 66, 33, 32768, 3, 20, 9, 130, 513, 65536, 7, 131072, 1025, 10, 1, 36, 19, 262144, 65, 258

Offset: 1

Mitch Cervinka (puritan(AT)planetkc.com), Oct 26 2003

The binary representation of a(n) shows which prime numbers divide n, but not the multiplicities. a(2)=1, a(3)=10, a(4)=1, a(5)=100, a(6)=11, a(10)=101, a(30)=111, etc.
For n > 1, a(n) gives the (one-based) index of the column where n is located in array A285321. A008479 gives the other index. - Antti Karttunen, Apr 17 2017
From Antti Karttunen, Jun 18 & 20 2017: (Start)
A268335 gives all n such that a(n) = A248663(n); the squarefree numbers (A005117) are all the n such that a(n) = A285330(n) = A048675(n).
For all n > 1 for which the value of A285331(n) is well-defined, we have A285331(a(n)) <= floor(A285331(n)/2), because then n is included in the binary tree A285332 and a(n) is one of its ancestors (in that tree), and thus must be at least one step nearer to its root than n itself.
Conjecture: Starting at any n and iterating the map n -> a(n), we will always reach 0 (see A288569). This conjecture is equivalent to the conjecture that at any n that is neither a prime nor a power of two, we will eventually hit a prime number (which then becomes a power of two in the next iteration). If this conjecture is false then sequence A285332 cannot be a permutation of natural numbers. On the other hand, if the conjecture is true, then A285332 must be a permutation of natural numbers, because all primes and powers of 2 occur in definite positions in that tree. This conjecture also implies the conjectures made in A019565 and A285320 that essentially claim that there are neither finite nor infinite cycles in A019565.
If there are any 2-cycles in this sequence, then both terms of the cycle should be present in A286611 and the larger one should be present in A286612.
(End)
Binary rank of the distinct prime indices of n, where the binary rank of an integer partition y is given by Sum_i 2^(y_i-1). For all prime indices (with multiplicity) we have A048675. - Gus Wiseman, May 25 2024

			a(38) = 129 because 38 = 2*19 = prime(1)*prime(8) and 129 = 2^0 + 2^7 (in binary 10000001).
a(140) = 13, binary 1101 because 140 is divisible by the first, third and fourth primes and 2^(1-1) + 2^(3-1) + 2^(4-1) = 13.

N. J. A. Sloane, Table of n, a(n) for n = 1..10000 [First 1000 terms from _T. D. Noe_]
Index entries for sequences related to binary expansion of n
Index entries for sequences computed from indices in prime factorization

For partial sums see A288566.
Cf. A000040, A000120, A001221, A005117, A008479, A019565, A055396, A285320, A285321, A285329, A285330, A285332.
Sequences with related definitions: A007947, A008472, A027748, A048675, A248663, A276379 (same sequence shown in base 2), A288569, A289271, A297404.
Cf. A286608 (numbers n for which a(n) < n), A286609 (n for which a(n) > n), and also A286611, A286612.
A003986, A003961, A059896 are used to express relationship between terms of this sequence.
Related to A267116 via A225546.
Positions of particular values are: A000079\{1} (1), A000244\{1} (2), A033845 (3), A000351\{1} (4), A033846 (5), A033849 (6), A143207 (7), A000420\{1} (8), A033847 (9), A033850 (10), A033851 (12), A147576 (14), A147571 (15), A001020\{1} (16), A033848 (17).
A048675 gives binary rank of prime indices.
A061395 gives greatest prime index, least A055396.
A112798 lists prime indices, length A001222, reverse A296150, sum A056239.
Binary indices (listed A048793):
- length A000120, complement A023416
- min A001511, opposite A000012
- sum A029931, product A096111
- max A029837 or A070939, opposite A070940
- complement A368494, sum A359400
- opposite complement A371571, sum A359359
- opposite A371572, sum A230877
Cf. A000720, A005940, A018819, A023506, A071814, A225620, A277319, A277905, A304818, A372689, A372890.

a087207 = sum . map ((2 ^) . (subtract 1) . a049084) . a027748_row
-- Reinhard Zumkeller, Jul 16 2013

a[n_] := Total[ 2^(PrimePi /@ FactorInteger[n][[All, 1]] - 1)]; a[1] = 0; Table[a[n], {n, 1, 69}] (* Jean-François Alcover, Dec 12 2011 *)

a(n) = {if (n==1, 0, my(f=factor(n), v = []); forprime(p=2, vecmax(f[,1]), v = concat(v, vecsearch(f[,1], p)!=0);); fromdigits(Vecrev(v), 2));} \\ Michel Marcus, Jun 05 2017

A087207(n)=vecsum(apply(p->1<M. F. Hasler, Jun 23 2017

from sympy import factorint, primepi
def a(n):
    return sum(2**primepi(i - 1) for i in factorint(n))
print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Jun 06 2017

(definec (A087207 n) (if (= 1 n) 0 (+ (A000079 (+ -1 (A055396 n))) (A087207 (A028234 n))))) ;; This uses memoization-macro definec
(define (A087207 n) (A048675 (A007947 n))) ;; Needs code from A007947 and A048675. - Antti Karttunen, Jun 19 2017

Additive with a(p^e) = 2^(i-1) where p is the i-th prime. - Vladeta Jovovic, Oct 29 2003
a(n) gives the m such that A019565(m) = A007947(n). - Naohiro Nomoto, Oct 30 2003
A000120(a(n)) = A001221(n); a(n) = Sum(2^(A049084(p)-1): p prime-factor of n). - Reinhard Zumkeller, Nov 30 2003
G.f.: Sum_{k>=1} 2^(k-1)*x^prime(k)/(1-x^prime(k)). - Franklin T. Adams-Watters, Sep 01 2009
From Antti Karttunen, Apr 17 2017, Jun 19 2017 & Dec 06 2018: (Start)
a(n) = A048675(A007947(n)).
a(1) = 0; for n > 1, a(n) = 2^(A055396(n)-1) + a(A028234(n)).
A000035(a(n)) = 1 - A000035(n). [a(n) and n are of opposite parity.]
A248663(n) <= a(n) <= A048675(n). [XOR-, OR- and +-variants.]
a(A293214(n)) = A218403(n).
a(A293442(n)) = A267116(n).
A069010(a(n)) = A287170(n).
A007088(a(n)) = A276379(n).
A038374(a(n)) = A300820(n) for n >= 1.
(End)
From Peter Munn, Jan 08 2020: (Start)
a(A059896(n,k)) = a(n) OR a(k) = A003986(a(n), a(k)).
a(A003961(n)) = 2*a(n).
a(n^2) = a(n).
a(n) = A267116(A225546(n)).
a(A225546(n)) = A267116(n).
(End)

More terms from Don Reble, Ray Chandler and Naohiro Nomoto, Oct 28 2003

Name clarified by Antti Karttunen, Jun 18 2017

A147571 Numbers with exactly 4 distinct prime divisors {2,3,5,7}.

Original entry on oeis.org

210, 420, 630, 840, 1050, 1260, 1470, 1680, 1890, 2100, 2520, 2940, 3150, 3360, 3780, 4200, 4410, 5040, 5250, 5670, 5880, 6300, 6720, 7350, 7560, 8400, 8820, 9450, 10080, 10290, 10500, 11340, 11760, 12600, 13230, 13440, 14700, 15120, 15750, 16800

Offset: 1

Artur Jasinski, Nov 07 2008

nonn

Successive numbers k such that EulerPhi(x)/x = m:
( Family of sequences for successive n primes )
m=1/2 numbers with exactly 1 distinct prime divisor {2} see A000079
m=1/3 numbers with exactly 2 distinct prime divisors {2,3} see A033845
m=4/15 numbers with exactly 3 distinct prime divisors {2,3,5} see A143207
m=8/35 numbers with exactly 4 distinct prime divisors {2,3,5,7} see A147571
m=16/77 numbers with exactly 5 distinct prime divisors {2,3,5,7,11} see A147572
m=192/1001 numbers with exactly 6 distinct prime divisors {2,3,5,7,11,13} see A147573
m=3072/17017 numbers with exactly 7 distinct prime divisors {2,3,5,7,11,13,17} see A147574
m=55296/323323 numbers with exactly 8 distinct prime divisors {2,3,5,7,11,13,17,19} see A147575

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

Cf. A002473, A086780, A143207, A147572, A147573, A147574, A147575, A147576, A147577, A147578, A147579, A147580.

[n: n in [1..2*10^4] | PrimeDivisors(n) eq [2,3,5,7]]; // Vincenzo Librandi, Sep 15 2015

a = {}; Do[If[EulerPhi[x]/x == 8/35, AppendTo[a, x]], {x, 1, 100000}]; a
Select[Range[20000],PrimeNu[#]==4&&Max[FactorInteger[#][[;;,1]]]<11&] (* Harvey P. Dale, Nov 05 2024 *)

is(n)=n%210==0 && n==2^valuation(n,2) * 3^valuation(n,3) * 5^valuation(n,5) * 7^valuation(n,7) \\ Charles R Greathouse IV, Jun 19 2016

a(n) = 210 * A002473(n). - David A. Corneth, May 14 2019

Sum_{n>=1} 1/a(n) = 1/48. - Amiram Eldar, Nov 12 2020

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

A147580 Numbers with exactly 7 distinct odd prime divisors {3,5,7,11,13,17,19}.

Original entry on oeis.org

4849845, 14549535, 24249225, 33948915, 43648605, 53348295, 63047985, 72747675, 82447365, 92147055, 101846745, 121246125, 130945815, 160044885, 169744575, 189143955, 218243025, 237642405, 247342095, 266741475, 276441165, 305540235, 315239925, 363738375, 373438065

Offset: 1

Artur Jasinski, Nov 07 2008

nonn

Numbers k such that phi(k)/k = m
( Family of sequences for successive n odd primes )
m=2/3 numbers with exactly 1 distinct prime divisor {3} see A000244
m=8/15 numbers with exactly 2 distinct prime divisors {3,5} see A033849
m=16/35 numbers with exactly 3 distinct prime divisors {3,5,7} see A147576
m=32/77 numbers with exactly 4 distinct prime divisors {3,5,7,11} see A147577
m=384/1001 numbers with exactly 5 distinct prime divisors {3,5,7,11,13} see A147578
m=6144/17017 numbers with exactly 6 distinct prime divisors {3,5,7,11,13,17} see A147579
m=3072/323323 numbers with exactly 7 distinct prime divisors {3,5,7,11,13,17,19} see A147580
m=110592/323323 numbers with exactly 8 distinct prime divisors {3,5,7,11,13,17,19,23} see A147581

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

Cf. A060735, A143207, A147571-A147575, A147576-A147580.

a = {}; Do[If[EulerPhi[4849845 x] == 1658880 x, AppendTo[a, 4849845 x]], {x, 1, 100}]; a

Sum_{n>=1} 1/a(n) = 1/1658880. - Amiram Eldar, Dec 22 2020

More terms from Amiram Eldar, Mar 11 2020

A147575 Numbers with exactly 8 distinct prime divisors {2,3,5,7,11,13,17,19}.

Original entry on oeis.org

9699690, 19399380, 29099070, 38798760, 48498450, 58198140, 67897830, 77597520, 87297210, 96996900, 106696590, 116396280, 126095970, 135795660, 145495350, 155195040, 164894730, 174594420, 184294110, 193993800, 203693490, 213393180, 232792560, 242492250, 252191940

Offset: 1

Artur Jasinski, Nov 07 2008

nonn

Successive numbers k such that EulerPhi(x)/x = m:
( Family of sequences for successive n primes )
m=1/2 numbers with exactly 1 distinct prime divisor {2} see A000079
m=1/3 numbers with exactly 2 distinct prime divisors {2,3} see A033845
m=4/15 numbers with exactly 3 distinct prime divisors {2,3,5} see A143207
m=8/35 numbers with exactly 4 distinct prime divisors {2,3,5,7} see A147571
m=16/77 numbers with exactly 5 distinct prime divisors {2,3,5,7,11} see A147572
m=192/1001 numbers with exactly 6 distinct prime divisors {2,3,5,7,11,13} see A147573
m=3072/17017 numbers with exactly 7 distinct prime divisors {2,3,5,7,11,13,17} see A147574
m=55296/323323 numbers with exactly 8 distinct prime divisors {2,3,5,7,11,13,17,19} see A147575

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

Cf. A060735, A080682, A143207, A147571-A147575, A147576-A147580.

a = {}; Do[If[EulerPhi[9699690 x] == 1658880 x, AppendTo[a, 9699690 x]], {x, 1, 100}]; a

a(n) = 9699690 * A080682(n). - Amiram Eldar, Mar 10 2020

Sum_{n>=1} 1/a(n) = 1/1658880. - Amiram Eldar, Nov 12 2020

More terms from Amiram Eldar, Mar 10 2020

A147572 Numbers with exactly 5 distinct prime divisors {2,3,5,7,11}.

Original entry on oeis.org

2310, 4620, 6930, 9240, 11550, 13860, 16170, 18480, 20790, 23100, 25410, 27720, 32340, 34650, 36960, 41580, 46200, 48510, 50820, 55440, 57750, 62370, 64680, 69300, 73920, 76230, 80850, 83160, 92400, 97020, 101640, 103950, 110880, 113190, 115500, 124740, 127050

Offset: 1

Artur Jasinski, Nov 07 2008

nonn

Successive numbers k such that EulerPhi(x)/x = m:
( Family of sequences for successive n primes )
m=1/2 numbers with exactly 1 distinct prime divisor {2} see A000079
m=1/3 numbers with exactly 2 distinct prime divisors {2,3} see A033845
m=4/15 numbers with exactly 3 distinct prime divisors {2,3,5} see A143207
m=8/35 numbers with exactly 4 distinct prime divisors {2,3,5,7} see A147571
m=16/77 numbers with exactly 5 distinct prime divisors {2,3,5,7,11} see A147572
m=192/1001 numbers with exactly 6 distinct prime divisors {2,3,5,7,11,13} see A147573
m=3072/17017 numbers with exactly 7 distinct prime divisors {2,3,5,7,11,13,17} see A147574
m=55296/323323 numbers with exactly 8 distinct prime divisors {2,3,5,7,11,13,17,19} see A147575

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

Cf. A051038, A060735, A143207, A147571-A147575, A147576-A147580.

a = {}; Do[If[EulerPhi[x]/x == 16/77, AppendTo[a, x]], {x, 1, 100000}]; a
Select[Range[130000],FactorInteger[#][[All,1]]=={2,3,5,7,11}&] (* Harvey P. Dale, Oct 04 2020 *)

from sympy import integer_log, prevprime
def A147572(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 g(x,m): return sum((x//3**i).bit_length() for i in range(integer_log(x,3)[0]+1)) if m==3 else sum(g(x//(m**i),prevprime(m))for i in range(integer_log(x,m)[0]+1))
    def f(x): return n+x-g(x,11)
    return 2310*bisection(f,n,n) # Chai Wah Wu, Sep 16 2024

a(n) = 2310 * A051038(n). - Amiram Eldar, Mar 10 2020

Sum_{n>=1} 1/a(n) = 1/480. - Amiram Eldar, Nov 12 2020

More terms from Amiram Eldar, Mar 10 2020

A147577 Numbers with exactly 4 distinct odd prime divisors {3,5,7,11}.

Original entry on oeis.org

1155, 3465, 5775, 8085, 10395, 12705, 17325, 24255, 28875, 31185, 38115, 40425, 51975, 56595, 63525, 72765, 86625, 88935, 93555, 114345, 121275, 139755, 144375, 155925, 169785, 190575, 202125, 218295, 259875, 266805, 280665, 282975, 317625

Offset: 1

Artur Jasinski, Nov 07 2008

nonn

Numbers k such that phi(k)/k = m
( Family of sequences for successive n odd primes )
m=2/3 numbers with exactly 1 distinct prime divisor {3} see A000244
m=8/15 numbers with exactly 2 distinct prime divisors {3,5} see A033849
m=16/35 numbers with exactly 3 distinct prime divisors {3,5,7} see A147576
m=32/77 numbers with exactly 4 distinct prime divisors {3,5,7,11} see A147577
m=384/1001 numbers with exactly 5 distinct prime divisors {3,5,7,11,13} see A147578
m=6144/17017 numbers with exactly 6 distinct prime divisors {3,5,7,11,13,17} see A147579
m=3072/323323 numbers with exactly 7 distinct prime divisors {3,5,7,11,13,17,19} see A147580
m=110592/323323 numbers with exactly 8 distinct prime divisors {3,5,7,11,13,17,19,23} see A147581

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

Cf. A060735, A143207, A147571-A147575, A147576-A147580.

a = {}; Do[If[EulerPhi[x]/x == 32/77, AppendTo[a, x]], {x, 1, 1000000}]; a
Select[Range[350000],EulerPhi[#]/#==32/77&] (* Harvey P. Dale, Mar 25 2016 *)

from sympy import integer_log
def A147577(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):
        c = n+x
        for i11 in range(integer_log(x,11)[0]+1):
            for i7 in range(integer_log(x11:=x//11**i11,7)[0]+1):
                for i5 in range(integer_log(x7:=x11//7**i7,5)[0]+1):
                    c -= integer_log(x7//5**i5,3)[0]+1
        return c
    return 1155*bisection(f,n,n) # Chai Wah Wu, Oct 22 2024

Sum_{n>=1} 1/a(n) = 1/480. - Amiram Eldar, Dec 22 2020

A147578 Numbers with exactly 5 distinct odd prime divisors {3,5,7,11,13}.

Original entry on oeis.org

15015, 45045, 75075, 105105, 135135, 165165, 195195, 225225, 315315, 375375, 405405, 495495, 525525, 585585, 675675, 735735, 825825, 945945, 975975, 1126125, 1156155, 1216215, 1366365, 1486485, 1576575, 1756755, 1816815, 1876875, 2027025, 2147145, 2207205, 2477475

Offset: 1

Artur Jasinski, Nov 07 2008

nonn

Numbers k such that phi(k)/k = m
( Family of sequences for successive n odd primes )
m=2/3 numbers with exactly 1 distinct prime divisor {3} see A000244;
m=8/15 numbers with exactly 2 distinct prime divisors {3,5} see A033849;
m=16/35 numbers with exactly 3 distinct prime divisors {3,5,7} see A147576;
m=32/77 numbers with exactly 4 distinct prime divisors {3,5,7,11} see A147577;
m=384/1001 numbers with exactly 5 distinct prime divisors {3,5,7,11,13} see A147578;
m=6144/17017 numbers with exactly 6 distinct prime divisors {3,5,7,11,13,17} see A147579;
m=3072/323323 numbers with exactly 7 distinct prime divisors {3,5,7,11,13,17,19} see A147580;
m=110592/323323 numbers with exactly 8 distinct prime divisors {3,5,7,11,13,17,19,23} see A147581.

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

Cf. A000010 (EulerPhi), A060735, A143207, A147571-A147575, A147576-A147580.

a = {}; Do[If[EulerPhi[x]/x == 384/1001, AppendTo[a, x]], {x, 1, 1000000}]; a

Sum_{n>=1} 1/a(n) = 1/5760. - Amiram Eldar, Dec 22 2020

More terms from Amiram Eldar, Mar 11 2020

A147579 Numbers with exactly 6 distinct odd prime divisors {3,5,7,11,13,17}.

Original entry on oeis.org

255255, 765765, 1276275, 1786785, 2297295, 2807805, 3318315, 3828825, 4339335, 5360355, 6381375, 6891885, 8423415, 8933925, 9954945, 11486475, 12507495, 13018005, 14039025, 16081065, 16591575, 19144125, 19654635, 20675655, 21696675, 23228205, 25270245, 26801775

Offset: 1

Artur Jasinski, Nov 07 2008

nonn

Numbers k such that phi(k)/k = m
( Family of sequences for successive n odd primes )
m=2/3 numbers with exactly 1 distinct prime divisor {3} see A000244
m=8/15 numbers with exactly 2 distinct prime divisors {3,5} see A033849
m=16/35 numbers with exactly 3 distinct prime divisors {3,5,7} see A147576
m=32/77 numbers with exactly 4 distinct prime divisors {3,5,7,11} see A147577
m=384/1001 numbers with exactly 5 distinct prime divisors {3,5,7,11,13} see A147578
m=6144/17017 numbers with exactly 6 distinct prime divisors {3,5,7,11,13,17} see A147579
m=3072/323323 numbers with exactly 7 distinct prime divisors {3,5,7,11,13,17,19} see A147580
m=110592/323323 numbers with exactly 8 distinct prime divisors {3,5,7,11,13,17,19,23} see A147581

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

Cf. A060735, A143207, A147571-A147575, A147576-A147580.

a = {}; Do[If[EulerPhi[255255 x] == 92160 x, AppendTo[a, 255255 x]], {x, 1, 100}]; a

Sum_{n>=1} 1/a(n) = 1/92160. - Amiram Eldar, Dec 22 2020

More terms from Amiram Eldar, Mar 11 2020

A147581 Numbers with exactly 8 distinct odd prime divisors {3,5,7,11,13,17,19,23}.

Original entry on oeis.org

111546435, 334639305, 557732175, 780825045, 1003917915, 1227010785, 1450103655, 1673196525, 1896289395, 2119382265, 2342475135, 2565568005, 2788660875, 3011753745, 3681032355, 3904125225, 4350310965, 5019589575, 5465775315, 5688868185, 6135053925, 6358146795

Offset: 1

Artur Jasinski, Nov 07 2008

nonn

Numbers k such that phi(k)/k = m
( Family of sequences for successive n odd primes )
m=2/3 numbers with exactly 1 distinct prime divisor {3} see A000244
m=8/15 numbers with exactly 2 distinct prime divisors {3,5} see A033849
m=16/35 numbers with exactly 3 distinct prime divisors {3,5,7} see A147576
m=32/77 numbers with exactly 4 distinct prime divisors {3,5,7,11} see A147577
m=384/1001 numbers with exactly 5 distinct prime divisors {3,5,7,11,13} see A147578
m=6144/17017 numbers with exactly 6 distinct prime divisors {3,5,7,11,13,17} see A147579
m=3072/323323 numbers with exactly 7 distinct prime divisors {3,5,7,11,13,17,19} see A147580
m=110592/323323 numbers with exactly 8 distinct prime divisors {3,5,7,11,13,17,19,23} see A147581

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

Cf. A060735, A143207, A147571-A147575, A147576-A147580.

a = {}; Do[If[EulerPhi[111546435 x] == 36495360 x, AppendTo[a, 111546435 x]], {x, 1, 100}]; a

from sympy import integer_log
def A147581(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):
        c = n+x
        for i23 in range(integer_log(x,23)[0]+1):
            for i19 in range(integer_log(x23:=x//23**i23,19)[0]+1):
                for i17 in range(integer_log(x19:=x23//19**i19,17)[0]+1):
                    for i13 in range(integer_log(x17:=x19//17**i17,13)[0]+1):
                        for i11 in range(integer_log(x13:=x17//13**i13,11)[0]+1):
                            for i7 in range(integer_log(x11:=x13//11**i11,7)[0]+1):
                                for i5 in range(integer_log(x7:=x11//7**i7,5)[0]+1):
                                    c -= integer_log(x7//5**i5,3)[0]+1
        return c
    return 111546435*bisection(f,n,n) # Chai Wah Wu, Oct 22 2024

Sum_{n>=1} 1/a(n) = 1/36495360. - Amiram Eldar, Dec 22 2020

More terms from Amiram Eldar, Mar 11 2020

cp's OEIS Frontend

A087207 A binary representation of the primes that divide a number, shown in decimal.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

PARI

Python

Scheme

Formula

Extensions

A147571 Numbers with exactly 4 distinct prime divisors {2,3,5,7}.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

A033851 Numbers whose prime factors are 5 and 7.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Mathematica

Formula

Extensions

A147580 Numbers with exactly 7 distinct odd prime divisors {3,5,7,11,13,17,19}.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A147575 Numbers with exactly 8 distinct prime divisors {2,3,5,7,11,13,17,19}.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A147572 Numbers with exactly 5 distinct prime divisors {2,3,5,7,11}.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Python

Formula

Extensions

A147577 Numbers with exactly 4 distinct odd prime divisors {3,5,7,11}.

Views

Author

Keywords