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 11-20 of 24 results. Next

A062838 Cubes of squarefree numbers.

Original entry on oeis.org

1, 8, 27, 125, 216, 343, 1000, 1331, 2197, 2744, 3375, 4913, 6859, 9261, 10648, 12167, 17576, 24389, 27000, 29791, 35937, 39304, 42875, 50653, 54872, 59319, 68921, 74088, 79507, 97336, 103823, 132651, 148877, 166375, 185193, 195112, 205379, 226981, 238328
Offset: 1

Views

Author

Jason Earls, Jul 21 2001

Keywords

Comments

Cubefull numbers (A036966) all of whose nonunitary divisors are not cubefull (A362147). - Amiram Eldar, May 13 2023

Crossrefs

Other powers of squarefree numbers: A005117(1), A062503(2), A113849(4), A072774(all).
A329332 column 3 in ascending order.

Programs

  • Mathematica
    Select[Range[70], SquareFreeQ]^3 (* Harvey P. Dale, Jul 20 2011 *)
  • PARI
    for(n=1,35, if(issquarefree(n),print(n*n^2)))
    
  • PARI
    a(n) = my(m, c); if(n<=1, n==1, c=1; m=1; while(cAltug Alkan, Dec 03 2015
    
  • Python
    from math import isqrt
    from sympy import mobius
    def A062838(n):
        def f(x): return n+x-sum(mobius(k)*(x//k**2) for k in range(1, isqrt(x)+1))
        m, k = n, f(n)
        while m != k: m, k = k, f(k)
        return m**3 # Chai Wah Wu, Sep 11 2024

Formula

A055229(a(n)) = A005117(n) and A055229(m) < A005117(n) for m < a(n). - Reinhard Zumkeller, Apr 09 2010
a(n) = A005117(n)^3. - R. J. Mathar, Dec 03 2015
{a(n)} = {A225546(A000400(n)) : n >= 0}, where {a(n)} denotes the set of integers in the sequence. - Peter Munn, Oct 31 2019
Sum_{n>=1} 1/a(n) = zeta(3)/zeta(6) = 945*zeta(3)/Pi^6 (A157289). - Amiram Eldar, May 22 2020

Extensions

More terms from Dean Hickerson, Jul 24 2001
More terms from Vladimir Joseph Stephan Orlovsky, Aug 15 2008

A097248 a(n) is the eventual stable point reached when iterating k -> A097246(k), starting from k = n.

Original entry on oeis.org

1, 2, 3, 3, 5, 6, 7, 6, 5, 10, 11, 5, 13, 14, 15, 5, 17, 10, 19, 15, 21, 22, 23, 10, 7, 26, 15, 21, 29, 30, 31, 10, 33, 34, 35, 15, 37, 38, 39, 30, 41, 42, 43, 33, 7, 46, 47, 15, 11, 14, 51, 39, 53, 30, 55, 42, 57, 58, 59, 7, 61, 62, 35, 15, 65, 66, 67, 51, 69, 70, 71, 30, 73, 74, 21
Offset: 1

Views

Author

Reinhard Zumkeller, Aug 03 2004

Keywords

Comments

a(n) = r(n,m) with m such that r(n,m)=r(n,m+1), where r(n,k) = A097246(r(n,k-1)), r(n,0)=n. (The original definition.)
A097248(n) = r(n,a(n)).
From Antti Karttunen, Nov 15 2016: (Start)
The above remark could be interpreted to mean that A097249(n) <= a(n).
All terms are squarefree, and the squarefree numbers are the fixed points.
These are also fixed points eventually reached when iterating A277886.
(End)

Crossrefs

Range of values is A005117.
A003961, A225546, A277885, A277886, A331590 are used to express relationship between terms of this sequence.
The formula section also details how the sequence maps the terms of A007913, A260443, A329050, A329332.
See comments/formulas in A283475, A283478, A331751 giving their relationship to this sequence.

Programs

  • Mathematica
    Table[FixedPoint[Times @@ Map[#1^#2 & @@ # &, Partition[#, 2, 2] &@ Flatten[FactorInteger[#] /. {p_, e_} /; e >= 2 :> {If[OddQ@ e, {p, 1}, {1, 1}], {NextPrime@ p, Floor[e/2]}}]] &, n], {n, 75}] (* Michael De Vlieger, Mar 18 2017 *)
  • PARI
    A097246(n) = { my(f=factor(n)); prod(i=1, #f~, (nextprime(f[i,1]+1)^(f[i,2]\2))*((f[i,1])^(f[i,2]%2))); };
    A097248(n) = { my(k=A097246(n)); while(k<>n, n = k; k = A097246(k)); k; };
    \\ Antti Karttunen, Mar 18 2017
    
  • Python
    from sympy import factorint, nextprime
    from operator import mul
    def a097246(n):
        f=factorint(n)
        return 1 if n==1 else reduce(mul, [(nextprime(i)**int(f[i]/2))*(i**(f[i]%2)) for i in f])
    def a(n):
        k=a097246(n)
        while k!=n:
            n=k
            k=a097246(k)
        return k # Indranil Ghosh, May 15 2017
  • Scheme
    ;; with memoization-macro definec
    ;; Two implementations:
    (definec (A097248 n) (if (not (zero? (A008683 n))) n (A097248 (A097246 n))))
    (definec (A097248 n) (if (zero? (A277885 n)) n (A097248 (A277886 n))))
    ;; Antti Karttunen, Nov 15 2016
    

Formula

a(A005117(n)) = A005117(n).
From Antti Karttunen, Nov 15 2016: (Start)
If A008683(n) <> 0 [when n is squarefree], a(n) = n, otherwise a(n) = a(A097246(n)).
If A277885(n) = 0, a(n) = n, otherwise a(n) = a(A277886(n)).
A007913(a(n)) = a(n).
a(A007913(n)) = A007913(n).
A048675(a(n)) = A048675(n).
a(A260443(n)) = A019565(n).
(End)
From Peter Munn, Feb 06 2020: (Start)
a(1) = 1; a(p) = p, for prime p; a(m*k) = A331590(a(m), a(k)).
a(A331590(m,k)) = A331590(a(m), a(k)).
a(n^2) = a(A003961(n)) = A003961(a(n)).
a(A225546(n)) = a(n).
a(n) = A225546(2^A048675(n)) = A019565(A048675(n)).
a(A329050(n,k)) = prime(n+k-1) = A000040(n+k-1).
a(A329332(n,k)) = A019565(n * k).
Equivalently, a(A019565(n)^k) = A019565(n * k).
(End)
From Antti Karttunen, Feb 22-25 & Mar 01 2020: (Start)
a(A019565(x)*A019565(y)) = A019565(x+y).
a(A332461(n)) = A332462(n).
a(A332824(n)) = A019565(n).
a(A277905(n,k)) = A277905(n,1) = A019565(n), for all n >= 1, and 1 <= k <= A018819(n).
(End)

Extensions

Name changed and the original definition moved to the Comments section by Antti Karttunen, Nov 15 2016

A329050 Square array A(n,k) = prime(n+1)^(2^k), read by descending antidiagonals (0,0), (0,1), (1,0), (0,2), (1,1), (2,0), ...; Fermi-Dirac primes (A050376) in matrix form, sorted into rows by their prime divisor.

Original entry on oeis.org

2, 4, 3, 16, 9, 5, 256, 81, 25, 7, 65536, 6561, 625, 49, 11, 4294967296, 43046721, 390625, 2401, 121, 13, 18446744073709551616, 1853020188851841, 152587890625, 5764801, 14641, 169, 17, 340282366920938463463374607431768211456, 3433683820292512484657849089281, 23283064365386962890625, 33232930569601, 214358881, 28561, 289, 19
Offset: 0

Views

Author

Antti Karttunen and Peter Munn, Nov 02 2019

Keywords

Comments

This sequence is a permutation of A050376, so every positive integer is the product of a unique subset, S_factors, of its terms. If we restrict S_factors to be chosen from a subset, S_0, consisting of numbers from specified rows and/or columns of this array, there are notable sequences among those that may be generated. See the examples. Other notable sequences can be generated if we restrict the intersection of S_factors with specific rows/columns to have even cardinality. In any of the foregoing cases, the numbers in the resulting sequence form a group under the binary operation A059897(.,.).
Shares with array A246278 the property that columns grow downward by iterating A003961, and indeed, this array can be obtained from A246278 by selecting its columns 1, 2, 8, 128, ..., 2^((2^k)-1), for k >= 0.
A(n,k) is the image of the lattice point with coordinates X=n and Y=k under the inverse of the bijection f defined in the first comment of A306697. This geometric relationship can be used to construct an isomorphism from the polynomial ring GF(2)[x,y] to a ring over the positive integers, using methods similar to those for constructing A297845 and A306697. See A329329, the ring's multiplicative operator, for details.

Examples

			The top left 5 X 5 corner of the array:
  n\k |   0     1       2           3                   4
  ----+-------------------------------------------------------
   0  |   2,    4,     16,        256,              65536, ...
   1  |   3,    9,     81,       6561,           43046721, ...
   2  |   5,   25,    625,     390625,       152587890625, ...
   3  |   7,   49,   2401,    5764801,     33232930569601, ...
   4  |  11,  121,  14641,  214358881,  45949729863572161, ...
Column 0 continues as a list of primes, column 1 as a list of their squares, column 2 as a list of their 4th powers, and so on.
Every nonnegative power of 2 (A000079) is a product of a unique subset of numbers from row 0; every squarefree number (A005117) is a product of a unique subset of numbers from column 0. Likewise other rows and columns generate the sets of numbers from sequences:
Row 1:                 A000244 Powers of 3.
Column 1:              A062503 Squares of squarefree numbers.
Row 2:                 A000351 Powers of 5.
Column 2:              A113849 4th powers of squarefree numbers.
Union of rows 0 and 1:     A003586 3-smooth numbers.
Union of columns 0 and 1:  A046100 Biquadratefree numbers.
Union of row 0 / column 0: A122132 Oddly squarefree numbers.
Row 0 excluding column 0:  A000302 Powers of 4.
Column 0 excluding row 0:  A056911 Squarefree odd numbers.
All rows except 0:         A005408 Odd numbers.
All columns except 0:      A000290\{0} Positive squares.
All rows except 1:         A001651 Numbers not divisible by 3.
All columns except 1:      A252895 (have odd number of square divisors).
If, instead of restrictions on choosing individual factors of the product, we restrict the product to be of an even number of terms from each row of the array, we get A262675. The equivalent restriction applied to columns gives us A268390; applied only to column 0, we get A028260 (product of an even number of primes).
		

Crossrefs

Transpose: A329049.
Permutation of A050376.
Rows 1-4: A001146, A011764, A176594, A165425 (after the two initial terms).
Antidiagonal products: A191555.
Subtable of A182944, A242378, A246278, A329332.
A000290, A003961, A225546 are used to express relationship between terms of this sequence.
Related binary operations: A059897, A306697, A329329.
See also the table in the example section.

Programs

  • Mathematica
    Table[Prime[#]^(2^k) &[m - k + 1], {m, 0, 7}, {k, m, 0, -1}] // Flatten (* Michael De Vlieger, Dec 28 2019 *)
  • PARI
    up_to = 105;
    A329050sq(n,k) = (prime(1+n)^(2^k));
    A329050list(up_to) = { my(v = vector(up_to), i=0); for(a=0, oo, for(col=0, a, i++; if(i > up_to, return(v)); v[i] = A329050sq(col, a-col))); (v); };
    v329050 = A329050list(up_to);
    A329050(n) = v329050[1+n];
    for(n=0,up_to-1,print1(A329050(n),", ")); \\ Antti Karttunen, Nov 06 2019

Formula

A(0,k) = 2^(2^k), and for n > 0, A(n,k) = A003961(A(n-1,k)).
A(n,k) = A182944(n+1,2^k).
A(n,k) = A329332(2^n,2^k).
A(k,n) = A225546(A(n,k)).
A(n,k+1) = A000290(A(n,k)) = A(n,k)^2.

Extensions

Example annotated for clarity by Peter Munn, Feb 12 2020

A329329 Multiplicative operator of a ring over the positive integers that has A059897(.,.) as additive operator and is isomorphic to GF(2)[x,y] with A329050(i,j) the image of x^i * y^j.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 4, 5, 4, 1, 1, 5, 9, 9, 5, 1, 1, 6, 7, 16, 7, 6, 1, 1, 7, 15, 25, 25, 15, 7, 1, 1, 8, 11, 36, 11, 36, 11, 8, 1, 1, 9, 27, 49, 35, 35, 49, 27, 9, 1, 1, 10, 25, 64, 13, 10, 13, 64, 25, 10, 1, 1, 11, 21, 81, 125, 77, 77, 125, 81
Offset: 1

Views

Author

Peter Munn, Nov 11 2019

Keywords

Comments

Square array A(n,k), n >= 1, k >= 1, read by descending antidiagonals.
The group defined by the binary operation A059897(.,.) over the positive integers is commutative with all elements self-inverse, and isomorphic to the additive group of GF(2) polynomial rings such as GF(2)[x,y]. There is a unique isomorphism extending each bijective mapping between respective minimal generating sets. The lexicographically earliest minimal generating set for the A059897 group is A050376, often called the Fermi-Dirac primes. This set has a natural arrangement in a square array, given by A329050(i,j) = prime(i+1)^(2^j), i >= 0, j >= 0. The most meaningful generating set for the additive group of GF(2)[x,y] is {x^i * y^j: i >= 0, j >= 0}, which similarly forms a square array. All this makes A329050(i,j) especially appropriate to be the image (under an isomorphism) of the GF(2) polynomial x^i * y^j.
Using g to denote the intended isomorphism, we specify g(x^i * y^j) = A329050(i,j). This maps minimal generating sets of the additive groups, so the definition of g is completed by specifying g(a+b) = A059897(g(a), g(b)). We then calculate the image under g of polynomial multiplication in GF(2)[x,y], giving us this sequence as the matching multiplicative operator for an isomorphic ring over the positive integers. Using f to denote the inverse of g, A[n,k] = g(f(n) * f(k)).
See the formula section for an alternative definition based on the A329050 array, independent of GF(2)[x,y].
Closely related to A306697 and A297845. If A059897 is replaced in the alternative definition by A059896 (and the definition is supplemented by the derived identity for the absorbing element, shown in the formula section), we get A306697; if A059897 is similarly replaced by A003991 (integer multiplication), we get A297845. This sequence and A306697, considered as multiplicative operators, are carryless arithmetic equivalents of A297845. A306697 uses a method analogous to binary-OR when there would be a multiplicative carry, while this sequence uses a method analogous to binary exclusive-OR. In consequence A(n,k) <> A297845(n,k) exactly when A306697(n,k) <> A297845(n,k). This relationship is not symmetric between the 3 sequences: there are n and k such that A(n,k) = A306697(n,k) <> A297845(n,k). For example A(54,72) = A306697(54,72) = 273375000 <> A297845(54,72) = 22143375000.

Examples

			Square array A(n, k) begins:
  n\k|  1   2   3    4    5    6    7     8     9    10    11    12
  ---+-------------------------------------------------------------
    1|  1   1   1    1    1    1    1     1     1     1     1     1
    2|  1   2   3    4    5    6    7     8     9    10    11    12
    3|  1   3   5    9    7   15   11    27    25    21    13    45
    4|  1   4   9   16   25   36   49    64    81   100   121   144
    5|  1   5   7   25   11   35   13   125    49    55    17   175
    6|  1   6  15   36   35   10   77   216   225   210   143   540
    7|  1   7  11   49   13   77   17   343   121    91    19   539
    8|  1   8  27   64  125  216  343    32   729  1000  1331  1728
    9|  1   9  25   81   49  225  121   729   625   441   169  2025
   10|  1  10  21  100   55  210   91  1000   441    22   187  2100
   11|  1  11  13  121   17  143   19  1331   169   187    23  1573
   12|  1  12  45  144  175  540  539  1728  2025  2100  1573    80
		

Crossrefs

A059897, A225546, A329050 are used to express relationship between terms of this sequence.
Related binary operations: A297845/A003991, A306697/A059896.

Programs

  • PARI
    \\ See Links section.

Formula

Alternative definition: (Start)
A(A329050(i_1, j_1), A329050(i_2, j_2)) = A329050(i_1+i_2, j_1+j_2).
A(A059897(n,k), m) = A059897(A(n,m), A(k,m)).
A(m, A059897(n,k)) = A059897(A(m,n), A(m,k)).
(End)
Derived identities: (Start)
A(n,1) = A(1,n) = 1 (1 is an absorbing element).
A(n,2) = A(2,n) = n.
A(n,k) = A(k,n).
A(n, A(m,k)) = A(A(n,m), k).
(End)
A(A019565(i), 2^j) = A019565(i)^j = A329332(i,j).
A(A225546(i), A225546(j)) = A225546(A(i,j)).
A(n,k) = A306697(n,k) = A297845(n,k), for n = A050376(i), k = A050376(j).
A(n,k) <= A306697(n,k) <= A297845(n,k).
A(n,k) < A297845(n,k) if and only if A306697(n,k) < A297845(n,k).

A331591 a(n) is the number of distinct prime factors of A225546(n), or equally, number of distinct prime factors of A293442(n).

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 2, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 1, 2, 2, 1, 1, 1, 1, 2, 1
Offset: 1

Views

Author

Antti Karttunen and Peter Munn, Jan 21 2020

Keywords

Comments

a(n) is the number of terms in the unique factorization of n into powers of squarefree numbers with distinct exponents that are powers of 2. See A329332 for a description of the relationship between this factorization, canonical (prime power) factorization and A225546.
The result depends only on the prime signature of n.
a(n) is the number of distinct bit-positions where there is a 1-bit in the binary representation of an exponent in the prime factorization of n. - Antti Karttunen, Feb 05 2020
The first 3 is a(128) = a(2^1 * 2^2 * 2^4) = 3 and in general each m occurs first at position 2^(2^m-1) = A058891(m+1). - Peter Munn, Mar 07 2022

Examples

			From _Peter Munn_, Jan 28 2020: (Start)
The factorization of 6 into powers of squarefree numbers with distinct exponents that are powers of 2 is 6 = 6^(2^0) = 6^1, which has 1 term. So a(6) = 1.
Similarly, 40 = 10^(2^0) * 2^(2^1) = 10^1 * 2^2 = 10 * 4, which has 2 terms. So a(40) = 2.
Similarly, 320 = 5^(2^0) * 2^(2^1) * 2^(2^2) = 5^1 * 2^2 * 2^4 = 5 * 4 * 16, which has 3 terms. So a(320) = 3.
10^100 (a googol) factorizes in this way as 10^4 * 10^32 * 10^64. So a(10^100) = 3.
(End)
		

Crossrefs

Sequences with related definitions: A001221, A331309, A331592, A331593, A331740.
Positions of records: A058891.
Positions of 1's: A340682.
Sequences used to express relationships between the terms: A007913, A008833, A059796, A331590.

Programs

  • Mathematica
    Array[PrimeNu@ If[# == 1, 1, Times @@ Flatten@ Map[Function[{p, e}, Map[Prime[Log2@ # + 1]^(2^(PrimePi@ p - 1)) &, DeleteCases[NumberExpand[e, 2], 0]]] @@ # &, FactorInteger[#]]] &, 105] (* Michael De Vlieger, Jan 24 2020 *)
    f[e_] := Position[Reverse[IntegerDigits[e, 2]], 1] // Flatten; a[n_] := CountDistinct[Flatten[f /@ FactorInteger[n][[;; , 2]]]]; a[1] = 0; Array[a, 100] (* Amiram Eldar, Dec 23 2023 *)
  • PARI
    A331591(n) = if(1==n,0,my(f=factor(n),u=#binary(vecmax(f[, 2])),xs=vector(u),m=1,e); for(i=1,u,for(k=1,#f~, if(bitand(f[k,2],m),xs[i]++)); m<<=1); #select(x -> (x>0),xs));
    
  • PARI
    A331591(n) = if(1==n, 0, hammingweight(fold(bitor, factor(n)[, 2]))); \\ Antti Karttunen, Feb 05 2020
    
  • PARI
    A331591(n) = if(n==1, 0, (core(n)>1) + A331591(core(n,1)[2])) \\ Peter Munn, Mar 08 2022

Formula

a(n) = A001221(A293442(n)) = A001221(A225546(n)).
From Peter Munn, Jan 28 2020: (Start)
a(n) = A000120(A267116(n)).
a(n) = a(A007913(n)) + a(A008833(n)).
For m >= 2, a(A005117(m)) = 1.
a(n^2) = a(n).
(End)
a(n) <= A331740(n) <= A048675(n) <= A293447(n). - Antti Karttunen, Feb 05 2020
From Peter Munn, Mar 07 2022: (Start)
a(n) <= A299090(n).
a(A337533(n)) = A299090(A337533(n)).
a(A337534(n)) < A299090(A337534(n)).
max(a(n), a(k)) <= a(A059796(n, k)) = a(A331590(n, k)) <= a(n) + a(k).
(End)

A182944 Square array A(i,j), i >= 1, j >= 1, of prime powers prime(i)^j, by descending antidiagonals.

Original entry on oeis.org

2, 4, 3, 8, 9, 5, 16, 27, 25, 7, 32, 81, 125, 49, 11, 64, 243, 625, 343, 121, 13, 128, 729, 3125, 2401, 1331, 169, 17, 256, 2187, 15625, 16807, 14641, 2197, 289, 19, 512, 6561, 78125, 117649, 161051, 28561, 4913, 361, 23
Offset: 1

Views

Author

Clark Kimberling, Dec 14 2010

Keywords

Comments

We alternatively refer to this sequence as a triangle T(.,.), with T(n,k) = A(k,n-k+1) = prime(k)^(n-k+1).
The monotonic ordering of this sequence, prefixed by 1, is A000961.
The joint-rank array of this sequence is A182869.
Main diagonal gives A062457. - Omar E. Pol, Sep 11 2018

Examples

			Square array A(i,j) begins:
  i \ j: 1      2      3      4      5  ...
  ---\-------------------------------------
  1:     2,     4,     8,    16,    32, ...
  2:     3,     9,    27,    81,   243, ...
  3:     5,    25,   125,   625,  3125, ...
  4:     7,    49,   343,  2401, 16807, ...
  ...
The triangle T(n,k) begins:
  n\k:  1     2     3     4     5     6  ...
  1:    2
  2:    4     3
  3:    8     9     5
  4:   16    27    25     7
  5:   32    81   125    49    11
  6:   64   243   625   343   121    13
  ...
		

Crossrefs

Cf. A000961, A006939 (row products of triangle), A062457, A182945, A332979 (row maxima of triangle).
Columns: A000040 (1), A001248 (2), A030078 (3), A030514 (4), A050997 (5), A030516 (6), A092759 (7), A179645 (8), A179665 (9), A030629 (10).
A319075 extends the array with 0th powers.
Subtable of A242378, A284457, A329332.

Programs

  • Mathematica
    TableForm[Table[Prime[n]^j,{n,1,14},{j,1,8}]]

Formula

From Peter Munn, Dec 29 2019: (Start)
A(i,j) = A182945(j,i) = A319075(j,i).
A(i,j) = A242378(i-1,2^j) = A329332(2^(i-1),j).
A(i,i) = A062457(i).
(End)

Extensions

Clarified in respect of alternate reading as a triangle by Peter Munn, Aug 28 2022

A191555 a(n) = Product_{k=1..n} prime(k)^(2^(n-k)).

Original entry on oeis.org

1, 2, 12, 720, 3628800, 144850083840000, 272760108249915378892800000000, 1264767303092594444142256488682840323816161280000000000000000
Offset: 0

Views

Author

Rick L. Shepherd, Jun 06 2011

Keywords

Comments

x^(2^n) - a(n) is the minimal polynomial over Q for the algebraic number sqrt(p(1)*sqrt(p(2)*...*sqrt(p(n-1)*sqrt(p(n)))...)), where p(k) is the k-th prime. Each such monic polynomial is irreducible by Eisenstein's Criterion (using p = p(n)).
A prime version of Somos's quadratic recurrence sequence A052129(n) = A052129(n-1)^2 * n = Product_{k=1..n} k^(2^(n-k)). - Jonathan Sondow, Mar 29 2014
All positive integers have unique factorizations into powers of distinct primes, and into powers of squarefree numbers with distinct exponents that are powers of 2. (See A329332 for a description of the relationship between the two.) a(n) is the least number such that both factorizations have n factors. - Peter Munn, Dec 15 2019
From Peter Munn, Jan 24 2020 to Feb 06 2020: (Start)
For n >= 0, a(n+1) is the n-th power of 12 in the monoid defined by A306697.
a(n) is the least positive integer that cannot be expressed as the product of fewer than n terms of A072774 (powers of squarefree numbers).
All terms that are less than the order of the Monster simple group (A003131) are divisors of the group's order, with a(6) exceeding its square root.
(End)
It is remarkable that 4 of the first 5 terms are factorials. - Hal M. Switkay, Jan 21 2025

Examples

			a(1) = 2^1 = 2 and x^2 - 2 is the minimal polynomial for the algebraic number sqrt(2).
a(4) = 2^8*3^4*5^2*7^1 = 3628800 and x^16 - 3628800 is the minimal polynomial for the algebraic number sqrt(2*sqrt(3*sqrt(5*sqrt(7)))).
		

Crossrefs

Sequences with related definitions: A006939, A052129, A191554, A239350 (and thence A239349), A252738, A266639.
A000290, A003961, A059896, A306697 are used to express relationship between terms of this sequence.
Subsequence of A025487, A138302, A225547, A267117 (apart from a(1) = 2), A268375, A331593.
Antidiagonal products of A329050.

Programs

  • Magma
    [n le 1 select 2 else Self(n-1)^2*NthPrime(n): n in [1..10]]; // Vincenzo Librandi, Feb 06 2016
  • Maple
    a:= proc(n) option remember;
          `if`(n=0, 1, a(n-1)^2*ithprime(n))
        end:
    seq(a(n), n=0..8);  # Alois P. Heinz, Mar 05 2020
  • Mathematica
    RecurrenceTable[{a[1] == 2, a[n] == a[n-1]^2 Prime[n]}, a, {n, 10}] (* Vincenzo Librandi, Feb 06 2016 *)
    Table[Product[Prime[k]^2^(n-k),{k,n}],{n,0,10}] (* or *) nxt[{n_,a_}]:={n+1,a^2 Prime[n+1]}; NestList[nxt,{0,1},10][[All,2]] (* Harvey P. Dale, Jan 07 2022 *)
  • PARI
    a(n) = prod(k=1, n, prime(k)^(2^(n-k)))
    
  • Scheme
    ;; Two variants, both with memoization-macro definec.
    (definec (A191555 n) (if (= 1 n) 2 (* (A000040 n) (A000290 (A191555 (- n 1)))))) ;; After the original recurrence.
    (definec (A191555 n) (if (= 1 n) 2 (* (A000079 (A000079 (- n 1))) (A003961 (A191555 (- n 1)))))) ;; After the alternative recurrence - Antti Karttunen, Feb 06 2016
    

Formula

For n > 0, a(n) = a(n-1)^2 * prime(n); a(0) = 1. [edited to extend to a(0) by Peter Munn, Feb 13 2020]
a(0) = 1; for n > 0, a(n) = 2^(2^(n-1)) * A003961(a(n-1)). - Antti Karttunen, Feb 06 2016, edited Feb 13 2020 because of the new prepended starting term.
For n > 1, a(n) = A306697(a(n-1),12) = A059896(a(n-1)^2, A003961(a(n-1))). - Peter Munn, Jan 24 2020

Extensions

a(0) added by Peter Munn, Feb 13 2020

A113850 Numbers whose prime factors are raised to the fifth power.

Original entry on oeis.org

32, 243, 3125, 7776, 16807, 100000, 161051, 371293, 537824, 759375, 1419857, 2476099, 4084101, 5153632, 6436343, 11881376, 20511149, 24300000, 28629151, 39135393, 45435424, 52521875, 69343957, 79235168, 90224199, 115856201
Offset: 1

Views

Author

Cino Hilliard, Jan 25 2006

Keywords

Examples

			7776 = 32*243 = 2^5*3^5 so the prime factors, 2 and 3, are raised to the fifth power.
		

Crossrefs

Proper subset of A000584.
Nonunit terms of A329332 column 5 in ascending order.

Programs

  • Mathematica
    Select[ Range@41^5, Union[Last /@ FactorInteger@# ] == {5} &] (* Robert G. Wilson v *)
    Rest[Select[Range[100], SquareFreeQ]^5] (* Vaclav Kotesovec, May 22 2020 *)
  • PARI
    allpwrfact(n,p) = \All prime factors are raised to the power p { local(x,j,ln,y,flag); for(x=4,n, y=Vec(factor(x)); ln = length(y[1]); flag=0; for(j=1,ln, if(y[2][j]==p,flag++); ); if(flag==ln,print1(x",")); ) }
    
  • Python
    from math import isqrt
    from sympy import mobius
    def A113850(n):
        def f(x): return int(n+x+1-sum(mobius(k)*(x//k**2) for k in range(1, isqrt(x)+1)))
        m, k = n, f(n)
        while m != k: m, k = k, f(k)
        return m**5 # Chai Wah Wu, Sep 13 2024

Formula

Sum_{k>=1} 1/a(k) = zeta(5)/zeta(10) - 1 = A157291 - 1. - Amiram Eldar, May 22 2020
a(n) = A005117(n+1)^5. - Chai Wah Wu, Sep 13 2024

Extensions

More terms from Robert G. Wilson v, Jan 26 2006
Offset corrected by Chai Wah Wu, Sep 13 2024

A331592 a(n) is the smaller of the number of terms in the factorizations of n into (1) powers of distinct primes and (2) powers of squarefree numbers with distinct exponents that are powers of 2.

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 1, 2, 2, 1, 1, 1, 1, 2, 1
Offset: 1

Views

Author

Antti Karttunen and Peter Munn, Jan 21 2020

Keywords

Comments

See A329332 for a description of the relationship between the two factorizations. From this relationship we get the formula a(n) = min(A001221(n), A001221(A225546(n))).
The result depends only on the prime signature of n.
k first appears at A191555(k).

Examples

			The factorization of 6 into powers of distinct primes is 6 = 2^1 * 3^1 = 2 * 3, which has 2 terms. Its factorization into powers of squarefree numbers with distinct exponents that are powers of 2 is 6 = 6^(2^0) = 6^1, which has 1 term. So a(6) is min(2,1) = 1.
The factorization of 40 into powers of distinct primes is 40 = 2^3 * 5^1 = 8 * 5, which has 2 terms. Its factorization into powers of squarefree numbers with distinct exponents that are powers of 2 is 40 = 10^(2^0) * 2^(2^1) = 10^1 * 2^2 = 10 * 4, which has 2 terms. So a(40) is min(2,2) = 2.
		

Crossrefs

Sequences with related definitions: A331308, A331591, A331593.
A003961, A225546 are used to express relationship between terms of this sequence.
Differs from = A071625 for the first time at n=216, where a(216) = 2, while A071625(216) = 1.

Programs

Formula

a(n) = min(A001221(n), A331591(n)) = min(A001221(n), A001221(A293442(n))).
a(A225546(n)) = a(n).
a(A003961(n)) = a(n).
a(n^2) = a(n).

A113851 Numbers whose prime factors are raised to the sixth power.

Original entry on oeis.org

64, 729, 15625, 46656, 117649, 1000000, 1771561, 4826809, 7529536, 11390625, 24137569, 47045881, 85766121, 113379904, 148035889, 308915776, 594823321, 729000000, 887503681, 1291467969, 1544804416, 1838265625, 2565726409, 3010936384, 3518743761, 4750104241
Offset: 1

Views

Author

Cino Hilliard, Jan 25 2006

Keywords

Crossrefs

Subset of A001014. Superset of A030516.
Nonunit terms of A329332 column 6 in ascending order.

Programs

  • Maple
    for n from 2 to 100 do if(numtheory[issqrfree](n))then printf("%d, ", n^6): fi: od: # Nathaniel Johnston, Jun 21 2011
  • Mathematica
    Select[ Range@37^6, Union[Last /@ FactorInteger@# ] == {6} &] (* Robert G. Wilson v *)
    Select[Range[2, 37], SquareFreeQ]^6 (* Amiram Eldar, Oct 13 2020 *)
  • Python
    from math import isqrt
    from sympy import mobius
    def A113851(n):
        def f(x): return int(n+1-sum(mobius(k)*(x//k**2) for k in range(2, isqrt(x)+1)))
        m, k = n, f(n)
        while m != k: m, k = k, f(k)
        return m**6 # Chai Wah Wu, Feb 25 2025

Formula

a(n) = A005117(n+1)^6. - Nathaniel Johnston, Jun 21 2011
Sum_{n>=1} 1/a(n) = zeta(6)/zeta(12) - 1 = A269404 - 1. - Amiram Eldar, Oct 13 2020

Extensions

More terms from Robert G. Wilson v, Jan 26 2006
Previous Showing 11-20 of 24 results. Next