A050376 -id:A050376 - cp's OEIS frontend

A050377 Number of ways to factor n into "Fermi-Dirac primes" (members of A050376).

1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 2, 1, 1, 1, 4, 1, 2, 1, 2, 1, 1, 1, 2, 2, 1, 2, 2, 1, 1, 1, 4, 1, 1, 1, 4, 1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 4, 2, 2, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 1, 1, 2, 6, 1, 1, 1, 2, 1, 1, 1, 4, 1, 1, 2, 2, 1, 1, 1, 4, 4, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 4, 1, 2, 2, 4, 1, 1, 1, 2, 1

Offset: 1

Christian G. Bower, Nov 15 1999

a(n) depends only on prime signature of n (cf. A025487). So a(24) = a(375) since 24 = 2^3 * 3 and 375 = 3 * 5^3 both have prime signature (3,1).

Antti Karttunen, Table of n, a(n) for n = 1..100000 (the first 10000 terms from Reinhard Zumkeller)
Index entries for sequences computed from exponents in factorization of n.

Cf. A000123, A001055, A002110, A005117, A018819, A050376-A050380, A025487.

Cf. A108951, A330687 (positions of records), A330688 (record values), A330689, A330690, A382295.

A018819:= proc(n) option remember;
  if n::odd then procname(n-1)
  else procname(n-1) + procname(n/2)
  fi
end proc:
A018819(0):= 1:
f:= n -> mul(A018819(s[2]),s=ifactors(n)[2]):
seq(f(n),n=1..100); # Robert Israel, Jan 14 2016

b[0] = 1; b[n_] := b[n] = b[n - 1] + If[EvenQ[n], b[n/2], 0];
a[n_] := Times @@ (b /@ FactorInteger[n][[All, 2]]);
Array[a, 102] (* Jean-François Alcover, Jan 27 2018 *)

A018819(n) = if( n<1, n==0, if( n%2, A018819(n-1), A018819(n/2)+A018819(n-1))); \\ From A018819
A050377(n) = factorback(apply(e -> A018819(e), factor(n)[, 2])); \\ Antti Karttunen, Dec 28 2019

Dirichlet g.f.: Product_{n in A050376} (1/(1-1/n^s)).
a(p^k) = A000123([k/2]) for all primes p.
a(A002110(n)) = 1.
Multiplicative with a(p^e) = A018819(e). - Christian G. Bower and David W. Wilson, May 22 2005
a(n) = Sum{a(d): d^2 divides n}, a(1) = 1. - Reinhard Zumkeller, Jul 12 2007
a(A108951(n)) = A330690(n). - Antti Karttunen, Dec 28 2019
a(n) = 1 for all squarefree values of n (A005117). - Eric Fox, Feb 03 2020
G.f.: Sum_{k>=1} a(k) * x^(k^2) / (1 - x^(k^2)). - Ilya Gutkovskiy, Nov 25 2020
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Product_{p prime} f(1/p) = 1.7876368001694456669... (A382295), where f(x) = (1-x) / Product_{k>=0} (1 - x^(2^k)). - Amiram Eldar, Oct 03 2023

More terms from Antti Karttunen, Dec 28 2019

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

Antti Karttunen and Peter Munn, Nov 02 2019

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.

			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).

Antti Karttunen, Table of n, a(n) for n = 0..77; the first 12 antidiagonals of array
Wikipedia, Polynomial ring

Transpose: A329049.
Permutation of A050376.
Rows 1-4: A001146, A011764, A176594, A165425 (after the two initial terms).
Columns 1-5: A000040, A001248, A030514, A179645, A030635.
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.

Table[Prime[#]^(2^k) &[m - k + 1], {m, 0, 7}, {k, m, 0, -1}] // Flatten (* Michael De Vlieger, Dec 28 2019 *)

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

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.

Example annotated for clarity by Peter Munn, Feb 12 2020

A177329 Number of factors in the representation of n! as a product of distinct terms of A050376.

Original entry on oeis.org

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

Offset: 2

Vladimir Shevelev, May 06 2010

nonn

Vladimir S. Shevelev, Multiplicative functions in the Fermi-Dirac arithmetic, Izvestia Vuzov of the North-Caucasus region, Nature sciences 4 (1996), 28-43 [Russian].

Chai Wah Wu, Table of n, a(n) for n = 2..10000 (terms 2..1000 from Amiram Eldar)
Simon Litsyn and Vladimir Shevelev, On factorization of integers with restrictions on the exponent, INTEGERS: Electronic Journal of Combinatorial Number Theory, 7 (2007), #A33, 1-36.

Cf. A000120, A001358, A050292, A050376, A064380, A064547, A176472, A176509, A176525.

read("transforms") ; A064547 := proc(n) f := ifactors(n)[2] ; a := 0 ; for p in f do a := a+wt(op(2,p)) ; end do: a ; end proc:
A177329 := proc(n) A064547(n!) ; end proc: seq(A177329(n),n=2..80) ; # R. J. Mathar, May 28 2010

f[p_, e_] := DigitCount[e, 2, 1]; a[1] = 0; a[n_] := Plus @@ f @@@ FactorInteger[n!]; Array[a, 100, 2] (* Amiram Eldar, Aug 24 2024 *)

a(n) = vecsum(apply(x -> hammingweight(x), factor(n!)[,2])); \\ Amiram Eldar, Aug 24 2024

from collections import Counter
from sympy import factorint
def A177329(n): return sum(map(int.bit_count,sum((Counter(factorint(i)) for i in range(2,n+1)),start=Counter()).values())) # Chai Wah Wu, Jul 18 2024

a(n) = Sum_{i} A000120(e_i), where n! = Product_{i} p_i^e_i is the prime factorization of n!.

a(n) = A064547(n!). - R. J. Mathar, May 28 2010

a(20)=10 inserted by Vladimir Shevelev, May 08 2010

Terms from a(14) onwards replaced according to the formula - R. J. Mathar, May 28 2010

A240695 a(n) is the smallest k such that a unique product of distinct terms of A050376 which is equal to k! contains at least the first n terms of A050376.

Original entry on oeis.org

2, 3, 4, 5, 125, 125, 138, 220, 220, 1766, 5526, 10351, 12365, 65653, 65653, 202738, 490333, 808762, 1478432, 1971352, 1971352, 1971352, 14798206, 14798206, 14798206, 14798206, 161974053, 547880880, 1619543840, 1619543840, 1619543840, 2103844465, 6435961044

Offset: 1

Vladimir Shevelev, Apr 10 2014

nonn

By the definition, the representation of a(n)! as a product of distinct terms of A050376 should contain the first n terms of A050376 and there is no restriction on the distribution of other factors of this product.

a(38) > 2 * 10^11. - Hiroaki Yamanouchi, Oct 01 2014

			5! = 2*3*4*5. We have the first 4 terms of A050376, so a(4) = 5.

Hiroaki Yamanouchi, Table of n, a(n) for n = 1..37

Cf. A240537, A240606, A240619, A240620, A240668, A240669, A240670, A240672.

bad[n_, pp_, mo_] := Catch[Do[If[ Mod[(n - Total@ IntegerDigits[n, pp[[i]]]) /(pp[[i]] - 1), mo[[i]] + 1] != mo[[i]], Throw@ True], {i, Length@ pp}]; False]; a[n_]:= Block[{fa, mo, pp, k},fa = FactorInteger[ Times @@ Select[Range[2, Prime[n]], (f = FactorInteger@# ; Length[f] == 1 && IntegerQ[Log[2, f[[1, 2]]]]) &, n]]; pp = First /@ fa; mo = Last /@ fa; k = fa[[-1, 1]]; While[ bad[k, pp, mo], k++]; k]; Array[a,15] (* Giovanni Resta, Apr 11 2014 *)

a(5)-a(23) from Giovanni Resta, Apr 11 2014

a(24)-a(33) from Hiroaki Yamanouchi, Oct 01 2014

A302781 Divisor-or-multiple permutation of natural numbers constructed from two-dimensional Hilbert curve (A163357) and Fermi-Dirac primes (A050376).

Original entry on oeis.org

1, 2, 6, 3, 15, 5, 10, 30, 120, 40, 20, 60, 12, 24, 8, 4, 28, 84, 168, 56, 14, 7, 21, 42, 210, 105, 35, 70, 280, 840, 420, 140, 1260, 3780, 7560, 2520, 630, 315, 945, 1890, 378, 189, 63, 126, 504, 1512, 756, 252, 36, 72, 216, 108, 540, 180, 360, 1080, 270, 90, 45, 135, 27, 54, 18, 9, 117, 351, 702, 234, 936, 468

Offset: 0

Antti Karttunen, Apr 14 2018

nonn

Note that the starting offset is 0, to align with A052330 and A207901.
Shares with A064736, A207901, A298480, A302350, A302783, A303771, etc. the property that a(n) is either a divisor or a multiple of a(n+1). Permutations satisfying such property are called "divisor-or-multiple permutations" in the OEIS, although Mazet & Saias call them "chain permutations" in their paper. [Edited by Antti Karttunen, Aug 26 2018]
One way to construct such permutations is by composing A052330 from the right with any such permutation like A003188 or A302846 where the binary expansions of a(n) and a(n+1) always differ by just a single bit-position.
Further permutations satisfying the same condition could be constructed from higher-dimensional versions (i.e., greater than 2) of Hilbert's space-filling curves, where the coordinates of each dimension would be Gray coded separately and then interleaved together. Permutation A207901 is essentially a one-dimensional variant of the same idea, while this is constructed from the 2-dimensional curve A163357, which is a Hamiltonian path on N X N grid.
See Peter Munn's A300012 for another idea for constructing such a permutation. - Antti Karttunen, Aug 26 2018

Antti Karttunen, Table of n, a(n) for n = 0..16383
Pierre Mazet and Eric Saias, Etude du graphe divisoriel 4, arXiv:1803.10073 [math.NT], 2018.
Various, Discussion on SeqFan-list, April 2018.
Index entries for sequences that are permutations of the natural numbers

Cf. A302782 (inverse).

Cf. A003188, A050376, A052330, A059252, A059253, A064706, A163356, A163357, A207901, A298480, A300012, A302844, A302846, A302783, A303771.

up_to_e = 14;
v050376 = vector(up_to_e);
A050376(n) = v050376[n];
ispow2(n) = (n && !bitand(n,n-1));
i = 0; for(n=1,oo,if(ispow2(isprimepower(n)), i++; v050376[i] = n); if(i == up_to_e,break));
A052330(n) = { my(p=1,i=1); while(n>0, if(n%2, p *= A050376(i)); i++; n >>= 1); (p); };
A064706(n) = bitxor(n, n>>2);
A057300(n) = { my(t=1,s=0); while(n>0, if(1==(n%4),n++,if(2==(n%4),n--)); s += (n%4)*t; n >>= 2; t <<= 2); (s); };
A163356(n) = if(!n,n,my(i = (#binary(n)-1)\2, f = 4^i, d = (n\f)%4, r = (n%f)); (((((2+(i%2))^d)%5)-1)*f) + if(3==d,f-1-A163356(r),A057300(A163356(r))));
A302781(n) = A052330(A064706(A163356(n)));

a(n) = A052330(A302846(n)), where A302846(n) = A000695(A003188(A059253(n))) + 2*A000695(A003188(A059252(n))).

a(n) = A207901(A302844(n)) = A052330(A064706(A163356(n))).

Name edited by Antti Karttunen, Aug 26 2018

A305829 Factor n into distinct Fermi-Dirac primes (A050376), normalize by replacing every instance of the k-th Fermi-Dirac prime with k, then multiply everything together.

Original entry on oeis.org

1, 1, 2, 3, 4, 2, 5, 3, 6, 4, 7, 6, 8, 5, 8, 9, 10, 6, 11, 12, 10, 7, 12, 6, 13, 8, 12, 15, 14, 8, 15, 9, 14, 10, 20, 18, 16, 11, 16, 12, 17, 10, 18, 21, 24, 12, 19, 18, 20, 13, 20, 24, 21, 12, 28, 15, 22, 14, 22, 24, 23, 15, 30, 27, 32, 14, 24, 30, 24, 20, 25

Offset: 1

Gus Wiseman, Jun 10 2018

Let f(n) = A050376(n) be the n-th Fermi-Dirac prime. Every positive integer n has a unique factorization of the form n = f(s_1)*...*f(s_k) where the s_i are strictly increasing positive integers. Then a(n) = s_1 * ... * s_k.

Multiplicative because for coprime m and n the Fermi-Dirac factorizations of m and n are disjoint and their union is the Fermi-Dirac factorization of m * n. - Andrew Howroyd, Aug 02 2018

Andrew Howroyd, Table of n, a(n) for n = 1..1000
Gus Wiseman, Tree of x -> a(x) for n = 1...75

Cf. A003963, A050376, A064547, A213925, A279065, A279614, A299755, A299756, A299757, A305830, A305831, A305832.

nn=100;
FDfactor[n_]:=If[n===1,{},Sort[Join@@Cases[FactorInteger[n],{p_,k_}:>Power[p,Cases[Position[IntegerDigits[k,2]//Reverse,1],{m_}->2^(m-1)]]]]];
FDprimeList=Array[FDfactor,nn,1,Union];FDrules=MapIndexed[(#1->#2[[1]])&,FDprimeList];
Table[Times@@(FDfactor[n]/.FDrules),{n,nn}]

\\ here isfd is membership test for A050376.
isfd(n)={my(e=isprimepower(n)); e && e == 1<Andrew Howroyd, Aug 02 2018

A366242 Numbers that are products of "Fermi-Dirac primes" (A050376) that are powers of primes with exponents that are powers of 4.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 10, 11, 13, 14, 15, 16, 17, 19, 21, 22, 23, 26, 29, 30, 31, 32, 33, 34, 35, 37, 38, 39, 41, 42, 43, 46, 47, 48, 51, 53, 55, 57, 58, 59, 61, 62, 65, 66, 67, 69, 70, 71, 73, 74, 77, 78, 79, 80, 81, 82, 83, 85, 86, 87, 89, 91, 93, 94, 95, 96, 97

Offset: 1

Amiram Eldar, Oct 05 2023

A subsequence of A252895, and first differs from it at n = 172. A252895(172) = 256 = 2^(2^3) is not a term of this sequence.
Equivalently, numbers that are products of "Fermi-Dirac primes" that are powers of primes with exponents that are powers of 2 with even exponents.
Products of distinct numbers of the form p^(2^(2*k)), where p is prime and k >= 0.
Numbers whose prime factorization has exponents that are positive terms of the Moser-de Bruijn sequence (A000695).
Every integer k has a unique representation as a product of 2 numbers: one is in this sequence and the other is in A366243: k = A366244(k) * A366245(k).
The asymptotic density of this sequence is 1/c = 0.65531174251481086750..., where c is given in the formula section.

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

Cf. A000695, A050376, A366243, A366244, A366245.

Subsequence of A252895.

mdQ[n_] := AllTrue[IntegerDigits[n, 4], # < 2 &]; Select[Range[100], AllTrue[FactorInteger[#][[;; , 2]], mdQ] &]

ismd(n) = {my(d = digits(n, 4)); for(i = 1, #d, if(d[i] > 1, return(0))); 1;}
is(n) = {my(e = factor(n)[ ,2]); for(i = 1, #e, if(!ismd(e[i]), return(0))); 1;}

a(n) ~ c * n, where c = Product_{k>=0} zeta(2^(2*k+1))/zeta(2^(2*k+2)) = 1.52599127273749217982... .

A366243 Numbers that are products of "Fermi-Dirac primes" (A050376) that are powers of primes with exponents that are not powers of 4.

Original entry on oeis.org

1, 4, 9, 25, 36, 49, 100, 121, 169, 196, 225, 256, 289, 361, 441, 484, 529, 676, 841, 900, 961, 1024, 1089, 1156, 1225, 1369, 1444, 1521, 1681, 1764, 1849, 2116, 2209, 2304, 2601, 2809, 3025, 3249, 3364, 3481, 3721, 3844, 4225, 4356, 4489, 4761, 4900, 5041, 5329

Offset: 1

Amiram Eldar, Oct 05 2023

Equivalently, numbers that are products of "Fermi-Dirac primes" that are powers of primes with exponents that are powers of 2 with odd exponents.
Products of distinct numbers of the form p^(2^(2*k-1)), where p is prime and k >= 1.
Numbers whose prime factorization has exponents that are positive terms of A062880.
Every integer k has a unique representation as a product of 2 numbers: one is in this sequence and the other is in A366242: k = A366245(k) * A366244(k).

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

Cf. A062880, A050376, A366242, A366244, A366245.

A062503 is a subsequence.

mdQ[n_] := AllTrue[IntegerDigits[n, 4], # < 2 &]; q[e_] := EvenQ[e] && mdQ[e/2]; Select[Range[6000], # == 1 || AllTrue[FactorInteger[#][[;; , 2]], q] &]

ismd(n) = {my(d = digits(n, 4)); for(i = 1, #d, if(d[i] > 1, return(0))); 1;}
is(n) = {my(e = factor(n)[ ,2]); for(i = 1, #e, if(e[i]%2 || !ismd(e[i]/2), return(0))); 1;}

Sum_{n>=1} 1/a(n) = Product_{k>=0} zeta(2^(2*k+1))/zeta(2^(2*k+2)) = 1.52599127273749217982... (this is the constant c in A366242).

A177333 Smallest factor in the factorization of n! over distinct terms of A050376.

Original entry on oeis.org

2, 2, 2, 2, 5, 5, 2, 2, 7, 7, 3, 3, 2, 2, 2, 2, 5, 5, 4, 3, 2, 2, 4, 4, 2, 2, 2, 2, 4, 4, 2, 2, 3, 3, 3, 3, 2, 2, 4, 4, 2, 2, 2, 2, 3, 3, 4, 4, 2, 2, 2, 2, 4, 4, 2, 2, 3, 3, 7, 7, 2, 2, 2, 2, 3, 3, 3, 4, 2, 2, 4, 4, 2, 2, 2, 2, 4, 4, 4, 4, 2, 2, 2, 2, 3, 4, 2, 2, 4, 4, 5, 3, 2, 2, 4, 4, 2, 2, 2, 2, 3, 3, 2, 2, 4

Offset: 2

Vladimir Shevelev, May 06 2010

nonn

			The factorization of 10! = 3628800 is 2^8*3^4*5^2*7^1, where 2^8 > 3^4 > 5^2 > 7, so a(10)=7 is the smallest of these 4 factors.

V. S. Shevelev, Multiplicative functions in the Fermi-Dirac arithmetic, Izvestia Vuzov of the North-Caucasus region, Nature sciences 4 (1996), 28-43 [Russian].

Amiram Eldar, Table of n, a(n) for n = 2..1000
S. Litsyn and V. S. Shevelev, On factorization of integers with restrictions on the exponent, INTEGERS: Electronic Journal of Combinatorial Number Theory, 7 (2007), #A33, 1-36.

Cf. A050376, A177329, A055460, A176525, A001358, A176472, A176509, A064380, A050292.

A177333 := proc(n) local a,p,pow2 ; a := n! ; for p in ifactors(n!)[2] do pow2 := convert( op(2,p),base,2) ; for j from 1 to nops(pow2) do if op(j,pow2) <> 0 then a := min(a,op(1,p)^(2^(j-1))) ; end if; end do: end do: return a ; end proc:
seq(A177333(n),n=2..120) ; # R. J. Mathar, Jun 16 2010

b[n_] :=2^(-1+Position[ Reverse@IntegerDigits[n, 2],?(#==1&)])//Flatten; a[n] := Module[{np = PrimePi[n]}, v=Table[0,{np}]; Do[p = Prime[k]; Do[v[[k]] += IntegerExponent[j, p], {j,2,n}],  {k,1,np}]; Min[(Prime/@Range[np])^(b/@v) // Flatten]]; Array[a, 105, 2] (* Amiram Eldar, Sep 17 2019 *)

Corrected from a(10) on and extended beyond a(30) by R. J. Mathar, Jun 16 2010

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

Vladimir Shevelev, Apr 07 2014

nonn

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

Eric Weisstein's World of Mathematics, Group.
Wikipedia, Generating set of a group.

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).
Cf. A000379, A030101, A052330, A059897, A329330.

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

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

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

cp's OEIS Frontend

A050377 Number of ways to factor n into "Fermi-Dirac primes" (members of A050376).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

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.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A177329 Number of factors in the representation of n! as a product of distinct terms of A050376.

Views

Author

Keywords

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Formula

Extensions

A240695 a(n) is the smallest k such that a unique product of distinct terms of A050376 which is equal to k! contains at least the first n terms of A050376.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Extensions

A302781 Divisor-or-multiple permutation of natural numbers constructed from two-dimensional Hilbert curve (A163357) and Fermi-Dirac primes (A050376).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Formula

Extensions

A305829 Factor n into distinct Fermi-Dirac primes (A050376), normalize by replacing every instance of the k-th Fermi-Dirac prime with k, then multiply everything together.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

A366242 Numbers that are products of "Fermi-Dirac primes" (A050376) that are powers of primes with exponents that are powers of 4.

Views

Author

Keywords

Comments

Links

Crossrefs