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.

Previous Showing 51-60 of 83 results. Next

A302033 a(n) = A019565(A003188(n)).

Original entry on oeis.org

1, 2, 6, 3, 15, 30, 10, 5, 35, 70, 210, 105, 21, 42, 14, 7, 77, 154, 462, 231, 1155, 2310, 770, 385, 55, 110, 330, 165, 33, 66, 22, 11, 143, 286, 858, 429, 2145, 4290, 1430, 715, 5005, 10010, 30030, 15015, 3003, 6006, 2002, 1001, 91, 182, 546, 273, 1365, 2730, 910, 455, 65, 130, 390, 195, 39, 78, 26, 13, 221, 442, 1326, 663, 3315, 6630, 2210, 1105
Offset: 0

Views

Author

Antti Karttunen & Peter Munn, Apr 16 2018

Keywords

Comments

A squarefree analog of A207901 (and the subsequence consisting of its squarefree terms): Each term is either a divisor or a multiple of the next one, and the terms differ by a single prime factor. Compare also to A284003.
For all n >= 0, max(a(n + 1), a(n)) / min(a(n + 1), a(n)) = A094290(n + 1) = prime(valuation(n + 1, 2) + 1) = A000040(A001511(n + 1)). [See Russ Cox's Dec 04 2010 comment in A007814.] - David A. Corneth & Antti Karttunen, Apr 16 2018

Crossrefs

A permutation of A005117. Subsequence of A207901.
Cf. A302054 (gives the sum of prime divisors).
Cf. also A277811, A283475, A284003.

Programs

  • Mathematica
    Array[Times @@ Prime@ Flatten@ Position[#, 1] &@ Reverse@ IntegerDigits[BitXor[#, Floor[#/2]], 2] &, 72, 0] (* Michael De Vlieger, Apr 27 2018 *)
  • PARI
    A003188(n) = bitxor(n, n>>1);
    A019565(n) = {my(j); factorback(Mat(vector(if(n, #n=vecextract(binary(n), "-1..1")), j, [prime(j), n[j]])~))}; \\ From A019565
    A302033(n) = A019565(A003188(n));
    
  • PARI
    first(n) = {my(pr = primes(1 + logint(n, 2)), ex = vector(#pr, i, 1), res = vector(n)); res[1] = 1; for(i = 1, n-1, v = valuation(i, 2); res[i + 1] = res[i] * pr[v++] ^ ex[v]; ex[v]*=-1); res}

Formula

a(n) = A019565(A003188(n)).
a(n) = A284003(A064706(n)).
a(n+1) = A059897(a(n), A094290(n+1)). - Peter Munn, Apr 01 2019

A089913 Table T(n,k) = lcm(n,k)/gcd(n,k) = n*k/gcd(n,k)^2 read by antidiagonals (n >= 1, k >= 1).

Original entry on oeis.org

1, 2, 2, 3, 1, 3, 4, 6, 6, 4, 5, 2, 1, 2, 5, 6, 10, 12, 12, 10, 6, 7, 3, 15, 1, 15, 3, 7, 8, 14, 2, 20, 20, 2, 14, 8, 9, 4, 21, 6, 1, 6, 21, 4, 9, 10, 18, 24, 28, 30, 30, 28, 24, 18, 10, 11, 5, 3, 2, 35, 1, 35, 2, 3, 5, 11, 12, 22, 30, 36, 40, 42, 42, 40, 36, 30, 22, 12, 13, 6, 33, 10, 45
Offset: 1

Views

Author

Marc LeBrun, Nov 14 2003

Keywords

Comments

A multiplicative analog of absolute difference A049581. Exponents in prime factorization of T(n,k) are absolute differences of those of n and k. Commutative non-associative operator with identity 1. T(nx,kx)=T(n,k), T(n^x,k^x)=T(n,k)^x, etc.
The bivariate function log(T(., .)) is a distance (or metric) function. It is a weighted analog of A130836, in the sense that if e_i (resp. f_i) denotes the exponent of prime p_i in the factorization of m (resp. of n), then both log(T(m, n)) and A130836(m, n) are writable as Sum_{i} w_i * abs(e_i - f_i). For A130836, w_i = 1 for all i, whereas for log(T(., .)), w_i = log(p_i). - Luc Rousseau, Sep 17 2018
If the analog of absolute difference, as described in the first comment, is determined by factorization into distinct terms of A050376 instead of by prime factorization, the equivalent operation is defined by A059897 and is associative. The positive integers form a group under A059897. The two factorization methods give the same factorization for squarefree numbers (A005117), so that T(.,.) restricted to A005117 is associative. Thus the squarefree numbers likewise form a group under the operation defined by this sequence. - Peter Munn, Apr 04 2019

Examples

			T(6,10) = lcm(6,10)/gcd(6,10) = 30/2 = 15.
  1,  2,  3,  4,  5, ...
  2,  1,  6,  2, 10, ...
  3,  6,  1, 12, 15, ...
  4,  2, 12,  1, 20, ...
  5, 10, 15, 20,  1, ...
  ...
		

Crossrefs

Programs

  • GAP
    T:=Flat(List([1..13],n->List([1..n-1],k->Lcm(k,n-k)/Gcd(k,n-k)))); # Muniru A Asiru, Oct 24 2018
    
  • Mathematica
    Flatten[Table[LCM[i, m - i]/GCD[i, m - i], {m, 15}, {i, m - 1}]] (* Ivan Neretin, Apr 27 2015 *)
  • PARI
    A089913(n,k)=n*k/gcd(n,k)^2 \\ M. F. Hasler, Dec 06 2019

Formula

A130836(n, k) = A001222(T(n, k)). - Luc Rousseau, Sep 17 2018

A240521 a(n) = A050376(n)*A050376(n+1) where A050376(n) is the n-th number of the form p^(2^k) with p is prime and k >= 0.

Original entry on oeis.org

6, 12, 20, 35, 63, 99, 143, 208, 272, 323, 437, 575, 725, 899, 1147, 1517, 1763, 2021, 2303, 2597, 3127, 3599, 4087, 4757, 5183, 5767, 6399, 6723, 7387, 8633, 9797, 10403, 11021, 11663, 12317, 13673, 15367, 16637, 17947, 19043, 20711, 22499, 23707, 25591
Offset: 1

Views

Author

Vladimir Shevelev, Apr 07 2014

Keywords

Comments

Let m be an odd positive number. Let S_m denote the sequence {Product_{i=1..r} q_(n+t_i)}A050376%20and%20Sum">{n>=1}, where {q_i} is sequence A050376 and Sum{i=1..r} 2^(t_1 - t_i) is the binary representation of m, such that t_1 > t_2 > ... > t_r = 0. Note that {S_1, S_3, S_5, ...} is a partition of all integers > 1. Then S_1=A050376, which is obtained when we set r=1, t_1 = 0. [Formula made compatible with A240535 data by Peter Munn, Aug 10 2021]
This present sequence is S_3 in this partition. It is obtained when we set r=2, t_1=1, t_2=0.
S_m(n) = A052330(A030101(m)*2^(n-1)) = A329330(A050376(n), A052330(A030101(m))). - Peter Munn, Aug 10 2021
A minimal set of generators for A000379 as a group under A059897(.,.). - Peter Munn, Aug 11 2019

Crossrefs

Positions of 3's in A240535.
Sequences for other parts of the partition described in the first comment: A050376 (S_1), A240522 (S_5), A240524 (S_7), A240536 (S_9), A241024 (S_11), A241025 (S_13).

Programs

  • Python
    from sympy import primepi, integer_nthroot
    def A240521(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(primepi(integer_nthroot(x,1<Chai Wah Wu, Feb 18-19 2025

Formula

a(n) = A052330(3*2^(n-1)) = A329330(A050376(n), 6). - Peter Munn, Aug 10 2021

Extensions

More terms from Peter J. C. Moses, Apr 18 2014

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.
		

Crossrefs

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

A120230 Split-floor-multiplier sequence (SFMS) using multipliers 1/4 and 4. (SFMS is defined at A120229.)

Original entry on oeis.org

4, 8, 12, 1, 20, 24, 28, 2, 36, 40, 44, 3, 52, 56, 60, 64, 68, 72, 76, 5, 84, 88, 92, 6, 100, 104, 108, 7, 116, 120, 124, 128, 132, 136, 140, 9, 148, 152, 156, 10, 164, 168, 172, 11, 180, 184, 188, 192, 196, 200, 204, 13, 212, 216, 220, 14, 228, 232, 236, 15, 244, 248
Offset: 1

Views

Author

Clark Kimberling, Jun 11 2006

Keywords

Comments

Self-inverse permutation of the natural numbers.

Examples

			a(1)=1*4 because [1/4] is not positive.
a(2)=2*4 because [2/4] is not positive.
a(3)=3*4 because [3/4] is not positive.
a(4)=1=[4*(1/4)].
a(5)=5*4 because [5/4]=a(4), not new.
		

Crossrefs

Row 4 and column 4 of A059897.

Formula

a(n)=[n/4] if this is positive and new, else a(n)=4n.

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

Original entry on oeis.org

1, 2, 3, 4, 6, 7, 8, 9, 12, 13, 14, 16, 18, 19, 21, 24, 25, 26, 27, 28, 31, 32, 36, 37, 38, 39, 42, 43, 48, 49, 50, 52, 54, 55, 56, 57, 61, 62, 63, 64, 67, 72, 73, 74, 75, 76, 78, 79, 81, 84, 85, 86, 91, 93, 96, 97, 98, 100, 103, 104, 108, 109, 110, 111, 112
Offset: 1

Views

Author

Ralf Stephan, May 16 2013

Keywords

Comments

The asymptotic density of this sequence is 1/2. - Amiram Eldar, Apr 03 2022
From Peter Munn, Nov 16 2023: (Start)
Contains all nonzero squares.
Dividing by 5 the terms that are multiples of 5 gives its complement, A225838.
(A352272, 2*A352272, 3*A352272, 6*A352272) is a partition of the terms.
The terms form a subgroup of the positive integers under the operation A059897(.,.) and are the positive integers in an index 2 multiplicative subgroup of rationals that is generated by 2, 3 and integers congruent to 1 modulo 6. See A225857 and A352272 for further information about such subgroups.
(End)

Crossrefs

Complement of A225838.
Subsequences: A003136\{0}, A083854\{0}, A260488\{0}, A352272.
Symmetric difference of A026225 and A036554; of A036668 and A189716.

Programs

  • Magma
    [n: n in [1..200] | IsOne(d mod 6) where d is n div (2^Valuation(n,2)*3^Valuation(n,3))]; // Bruno Berselli, May 16 2013
    
  • Mathematica
    mx = 122; t = {}; Do[n = 2^i*3^j (6 k + 1); 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==1),print1(n,",")))
    
  • Python
    from sympy import integer_log
    def A225837(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//3**i>>j)+5)//6 for i in range(integer_log(x,3)[0]+1) for j in range((x//3**i).bit_length()))
        return bisection(f,n,n) # Chai Wah Wu, Feb 02 2025

A334110 The squares of squarefree numbers (A062503), ordered lexicographically according to their prime factors. a(n) = Product_{k in I} prime(k+1)^2, where I are the indices of nonzero binary digits in n = Sum_{k in I} 2^k.

Original entry on oeis.org

1, 4, 9, 36, 25, 100, 225, 900, 49, 196, 441, 1764, 1225, 4900, 11025, 44100, 121, 484, 1089, 4356, 3025, 12100, 27225, 108900, 5929, 23716, 53361, 213444, 148225, 592900, 1334025, 5336100, 169, 676, 1521, 6084, 4225, 16900, 38025, 152100, 8281, 33124, 74529, 298116, 207025, 828100, 1863225, 7452900, 20449, 81796, 184041
Offset: 0

Views

Author

Antti Karttunen and Peter Munn, May 01 2020

Keywords

Comments

For the lexicographic ordering, the prime factors must be written in nonincreasing order. If we write the factors in nondecreasing order, we get a lexicographically ordered set with an order type that is greater than a natural number index - the resulting sequence does not include all qualifying numbers. (Note also that the symbols used for the lexicographic order are the prime numbers, not their digits.)
a(n) is the n-th power of 4 in the monoid defined in A331590.
Conjecture: a(n) is the position of the first occurrence of n in A334109.

Examples

			The initial terms are shown below, equated with the product of their prime factors to exhibit the lexicographic ordering. The list starts with 1, since 1 is factored as the empty product and the empty list is first in lexicographic order.
    1 = .
    4 = 2*2.
    9 = 3*3.
   36 = 3*3*2*2.
   25 = 5*5.
  100 = 5*5*2*2.
  225 = 5*5*3*3.
  900 = 5*5*3*3*2*2.
   49 = 7*7.
  196 = 7*7*2*2.
  441 = 7*7*3*3.
		

Crossrefs

Cf. A000079, A019565 (square roots), A048675, A097248, A225546, A267116, A332382, A334109 (a left inverse).
Column 2 of A329332. Permutation of A062503.
After 1, the right children of the leftmost edge of A334860, or respectively, the left children of the rightmost edge of A334866.
Subsequences: A001248, A061742, A166329.
Subsequence of A052330.
A003961, A003987, A059897, A331590 are used to express relationship between terms of this sequence.

Programs

  • Mathematica
    Array[If[# == 0, 1, Times @@ Flatten@ Map[Function[{p, e}, Map[Prime[Log2@ # + 1]^(2^(PrimePi@ p - 1)) &, DeleteCases[NumberExpand[e, 2], 0]]] @@ # &, FactorInteger[3^#]]] &, 51, 0] (* Michael De Vlieger, May 26 2020 *)
  • PARI
    A334110(n) = { my(p=2,m=1); while(n, if(n%2, m *= p^2); n >>= 1; p = nextprime(1+p)); (m); };

Formula

a(n) = A019565(n)^2.
For n >= 1, a(A000079(n-1)) = A001248(n).
For all n >= 0, A334109(a(n)) = n.
a(n+k) = A331590(a(n), a(k)).
a(n XOR k) = A059897(a(n), a(k)), where XOR denotes bitwise exclusive-or, A003987.
a(n) = A225546(3^n).
a(2n) = A003961(a(n)).
a(2n+1) = 4 * a(2n).
a(2^k-1) = A061742(k).
A267116(a(n)) = 2.
A048675(a(n)) = 2n.
A097248(a(n)) = A332382(n) = A019565(2n).

A352272 Numbers whose squarefree part is congruent to 1 modulo 6.

Original entry on oeis.org

1, 4, 7, 9, 13, 16, 19, 25, 28, 31, 36, 37, 43, 49, 52, 55, 61, 63, 64, 67, 73, 76, 79, 81, 85, 91, 97, 100, 103, 109, 112, 115, 117, 121, 124, 127, 133, 139, 144, 145, 148, 151, 157, 163, 169, 171, 172, 175, 181, 187, 193, 196, 199, 205, 208, 211, 217, 220, 223, 225, 229
Offset: 1

Views

Author

Peter Munn, Mar 10 2022

Keywords

Comments

Numbers of the form 4^i * 9^j * (6k+1), i, j, k >= 0.
Closed under multiplication.
The sequence forms a subgroup of the positive integers under the commutative operation A059897(.,.), one of 8 subgroups of the form {k : A007913(k) == 1 (mod m)} - in each case m is a divisor of 24. A059897 has a relevance to squarefree parts that arises from the identity A007913(k*n) = A059897(A007913(k), A007913(n)), where A007913(n) is the squarefree part of n.
The subgroup has 8 cosets, which partition the positive integers as follows. For each i in {1, 5}, j in {1, 2, 3, 6} there is a coset {m^2 * (6k+i) * j : m >= 1, k >= 0}. See the table in the examples.
None of the 8 cosets have been entered into the database previously, but many subgroups of the quotient group (which are formed of certain combinations of cosets) are represented among earlier OEIS sequences, including 6 of the 7 subgroups of index 2 (which combine 4 cosets). This sequence can therefore be defined as the intersection of pairs or triples of these sequences in many combinations (see the cross-references). See also the table in the example section of A352273 (the coset that includes 5).

Examples

			The squarefree part of 9 is 1, which is congruent to 1 (mod 6), so 9 is in the sequence.
The squarefree part of 14 is 14, which is congruent to 2 (mod 6), so 14 is not in the sequence.
The squarefree part of 52 = 2^2 * 13 is 13, which is congruent to 1 (mod 6), so 52 is in the sequence.
The 8 cosets described in the initial comments (forming a partition of the positive integers) are shown as rows of the following table. The first half of the table corresponds to (6k+i) with i=1; the second half to i=5, with row 5 being A352273.
   1,  4,   7,   9,  13,  16,  19,  25,  28,  31,  36, ...
   2,  8,  14,  18,  26,  32,  38,  50,  56,  62,  72, ...
   3, 12,  21,  27,  39,  48,  57,  75,  84,  93, 108, ...
   6, 24,  42,  54,  78,  96, 114, 150, 168, 186, 216, ...
   5, 11,  17,  20,  23,  29,  35,  41,  44,  45,  47, ...
  10, 22,  34,  40,  46,  58,  70,  82,  88,  90,  94, ...
  15, 33,  51,  60,  69,  87, 105, 123, 132, 135, 141, ...
  30, 66, 102, 120, 138, 174, 210, 246, 264, 270, 282, ...
The product of two positive integers is in this sequence if and only if they are in the same coset. The asymptotic density of cosets (containing) 1 and 5 is 1/4; of cosets 2 and 10 is 1/8; of cosets 3 and 15 is 1/12; of cosets 6 and 30 is 1/24.
		

Crossrefs

Intersection of any 2 of A055047, A339690 and A352274.
Intersection of any 4 sets chosen from A003159, A007417, A026225, A036668, A189715 and A225837 (in most cases, only 3 sets are needed - specifically if the pairwise intersections of the 3 sets differ from each other).
Closure of A084089 under multiplication by 9.
Other subsequences: A000290\{0}, A016921, A229848 (apparently, with 55 the first difference).
A334832 lists equivalent sequences modulo other divisors of 24.

Programs

  • PARI
    isok(m) = core(m) % 6 == 1;
    
  • Python
    from itertools import count
    def A352272(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-1)//6+1
            return c
        return bisection(f,n,n) # Chai Wah Wu, Feb 14 2025

Formula

{a(n) : n >= 1} = {m >= 1 : A007913(m) == 1 (mod 6)}.
{a(n) : n >= 1} = A334832 U A334832/7 U A334832/13 U A334832/19 where A334832/k denotes {A334832(m)/k : m >= 1, k divides A334832(m)}.
Using the same denotation, {a(n) : n >= 1} = A352273/5 = {A307151(A352273(m)) : m >= 1}.

A082522 Numbers of the form p^(2^k) with p prime and k>0.

Original entry on oeis.org

4, 9, 16, 25, 49, 81, 121, 169, 256, 289, 361, 529, 625, 841, 961, 1369, 1681, 1849, 2209, 2401, 2809, 3481, 3721, 4489, 5041, 5329, 6241, 6561, 6889, 7921, 9409, 10201, 10609, 11449, 11881, 12769, 14641, 16129, 17161, 18769, 19321, 22201, 22801
Offset: 1

Views

Author

Reinhard Zumkeller, May 11 2003

Keywords

Comments

Every positive square (A000290 without 0) is the product of a unique subset of these numbers. The lexicographically earliest (when ordered) minimal set of generators for the positive squares as a group under A059897(.,.); the intersection of A050376 and A000290. - Peter Munn, Aug 25 2019

Examples

			3^(2^2) = 81, therefore 81 is a term.
		

Crossrefs

Programs

  • PARI
    lst(lim)=my(v=List(apply(n->n^2,primes(primepi(sqrtint(lim))))), t);forprime(p=2,(lim+.5)^(1/4),t=p^2;while((t=t^2)<=lim,listput(v,t)));vecsort(Vec(v)) \\ Charles R Greathouse IV, Apr 10 2012
    
  • Python
    from sympy import primepi, integer_nthroot
    def A082522(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(primepi(integer_nthroot(x,1<Chai Wah Wu, Feb 19 2025

Formula

a(n) = A050376(A181970(n)) = A050376(n)^2. - Vladimir Shevelev, Apr 05 2013
a(n) ~ n^2 log^2 n. - Charles R Greathouse IV, Oct 19 2015
Sum_{n>=1} 1/a(n) = Sum_{k>=1} P(2^k) = 0.53331724743088069672..., where P is the prime zeta function. - Amiram Eldar, Nov 26 2020

A319241 Heinz numbers of strict integer partitions of even numbers. Squarefree numbers whose prime indices sum to an even number.

Original entry on oeis.org

1, 3, 7, 10, 13, 19, 21, 22, 29, 30, 34, 37, 39, 43, 46, 53, 55, 57, 61, 62, 66, 70, 71, 79, 82, 85, 87, 89, 91, 94, 101, 102, 107, 111, 113, 115, 118, 129, 130, 131, 133, 134, 138, 139, 146, 151, 154, 155, 159, 163, 165, 166, 173, 181, 183, 186, 187, 190, 193
Offset: 1

Views

Author

Gus Wiseman, Sep 15 2018

Keywords

Comments

The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).
From Peter Munn, Feb 04 2022: (Start)
For every odd squarefree number, s, exactly one of s and 2s is a term.
Closed under the commutative operation A350066(.,.).
Closed under the commutative operation A059897(.,.) forming a subgroup of the positive integers considered as a group under A059897. As subgroups, this sequence and A028982 are each a transversal of the other.
(End)

Examples

			30 is the Heinz number of (3,2,1), which is strict and has even weight, so 30 belongs to the sequence.
The sequence of all even-weight strict partitions begins: (), (2), (4), (3,1), (6), (8), (4,2), (5,1), (10), (3,2,1), (7,1), (12), (6,2), (14), (9,1), (16), (5,3), (8,2), (18), (11,1), (5,2,1), (4,3,1).
		

Crossrefs

Complement of the union of A319242 and A013929.
Intersection of A005117 and A300061.

Programs

  • Mathematica
    Select[Range[100],And[SquareFreeQ[#],EvenQ[Total[Cases[FactorInteger[#],{p_,k_}:>k*PrimePi[p]]]]]&]
  • PARI
    isok(m) = issquarefree(m) && !(vecsum(apply(primepi, factor(m)[,1])) % 2); \\ Michel Marcus, Feb 08 2022

Formula

{a(n) : n >= 1} = {A019565(A158704(n)) : n >= 1} = {A073675(A319242(n)) : n >= 1}. - Peter Munn, Feb 04 2022
Previous Showing 51-60 of 83 results. Next