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.

Showing 1-7 of 7 results.

A177334 Largest factor in the factorization of n! over distinct terms of A050376.

Original entry on oeis.org

2, 3, 4, 5, 16, 16, 16, 81, 256, 256, 256, 256, 256, 256, 256, 256, 65536, 65536, 65536, 65536, 65536, 65536, 65536, 65536, 65536, 65536, 65536, 65536, 65536, 65536, 65536, 65536, 4294967296, 4294967296, 4294967296, 4294967296, 4294967296, 4294967296
Offset: 2

Views

Author

Vladimir Shevelev, May 06 2010

Keywords

Comments

Each number >=2 has a unique factorization over distinct terms of A050376.
This is obtained from the standard prime factor representation by splitting the exponents into a sum of powers of 2, and further factorization according to the nonzero term of this base-2 representation.
The largest factor of this representation of A000142(n) defines this sequence.

References

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

Crossrefs

Programs

  • Maple
    A177334 := proc(n) local a,p,pow2 ; a := 1 ; 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 := max(a,op(1,p)^(2^(j-1))) ; end if; end do: end do: return a ; end proc:
    seq(A177334(n),n=2..60) ; # R. J. Mathar, Jun 16 2010
  • Mathematica
    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}]; Max[(Prime/@Range[np])^(b/@v) // Flatten]]; Array[a, 38, 2]  (* Amiram Eldar, Sep 17 2019 *)

Extensions

a(18) and a(19) corrected and sequence extended by R. J. Mathar, Jun 16 2010

A241139 Number of nonprimes in factorization of n! over distinct terms of A050376.

Original entry on oeis.org

0, 0, 1, 1, 2, 2, 3, 3, 3, 3, 4, 4, 5, 6, 7, 7, 4, 4, 5, 5, 6, 6, 8, 9, 10, 10, 9, 9, 11, 11, 12, 12, 10, 9, 8, 8, 9, 10, 11, 11, 12, 12, 11, 12, 14, 14, 16, 15, 15, 15, 13, 13, 14, 14, 14, 14, 16, 16, 16, 16, 17, 19, 21, 21, 18, 18, 19, 16, 14, 14, 16, 16, 17
Offset: 2

Views

Author

Vladimir Shevelev, Apr 16 2014

Keywords

Examples

			Factorization of 4! over distinct terms of A050376 is 4! = 2*3*4. This factorization contains only one A050376-nonprime. So a(4)=1.
		

References

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

Crossrefs

Programs

  • Mathematica
    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}]; Length[Select[(b /@ v) // Flatten, # > 1 &]]]; Array[a, 73, 2]  (* Amiram Eldar, Sep 17 2019 *)
  • PARI
    a(n)={my(f=factor(n!)[,2]); sum(i=1, #f~, hammingweight(f[i]>>1))} \\ Andrew Howroyd, Sep 17 2019

Formula

a(n) = A177329(n) - A055460(n).

Extensions

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

A241148 Number of factorials k!, 0<=k<=n, relatively prime to n! in Fermi-Dirac arithmetic.

Original entry on oeis.org

1, 2, 2, 2, 2, 2, 5, 5, 2, 2, 7, 7, 4, 4, 4, 2, 2, 2, 5, 5, 7, 4, 3, 3, 4, 4, 2, 2, 4, 4, 4, 4, 2, 2, 3, 3, 4, 4, 3, 2, 4, 4, 3, 3, 2, 4, 5, 5, 4, 4, 2, 2, 2, 2, 6, 5, 2, 2, 3, 3, 7, 7, 3, 2, 2, 2, 3, 3, 3, 4, 3, 3, 4, 4, 2, 2, 2, 2, 6, 6, 4, 4, 2, 2, 2, 3, 4
Offset: 0

Views

Author

Vladimir Shevelev, Apr 16 2014

Keywords

Comments

Or, equivalently, the number of factorials k!, 0<=k<=n, for which k! and n! have no common A050376-factors in their factorizations over distinct terms of A050376.
Note that 1 (=0!=1!) corresponds to an empty subset of A050376.

Examples

			0!=1, 1!=1; further we have the following factorizations of k! over distinct terms of A050376 for k = 2,3,4,5,6:
2!=2, 3!=2*3, 4!=2*3*4, 5!=2*3*4*5, 6!=5*9*16. Thus, in the sense of the factorizations being considered, 6! is relatively prime to 0!,1!,2!,3!, and 4!. So a(6)=5.
		

References

  • V. S. Shevelev, Multiplicative functions in the Fermi-Dirac arithmetic, Izvestia Vuzov of the North-Caucasus region, Nature sciences 4 (1996), 28-43 (Russian; MR 2000f: 11097, pp. 3912-3913).

Crossrefs

Programs

  • Mathematica
    b[n_] := 2^(-1 + Position[Reverse@IntegerDigits[n, 2], ?(# == 1 &)]) // Flatten; infp[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}]; (Prime /@ Range[np])^(b /@ v) // Flatten]; infCoprimeQ[x_, y_] := Intersection[infp[x], infp[y]] == {}; a[n_] := Length @ Select[Range[0, n], infCoprimeQ[n, #] & ]; Array[a, 87, 0] (* Amiram Eldar, Sep 17 2019 *)

Extensions

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

A241289 Numbers n for which in the factorization of n! over distinct terms of A050376, the numbers of prime and nonprime terms are equal.

Original entry on oeis.org

7, 8, 9, 13, 18, 22, 37, 57, 71
Offset: 1

Views

Author

Vladimir Shevelev, Apr 18 2014

Keywords

Comments

a(10), if it exists, should be more than 5000. Is a(9)=71 the last term of sequence? - Peter J. C. Moses, Apr 19 2014
One can prove that a(9)=71 indeed is the last term of this sequence. - Vladimir Shevelev, Apr 19 2014.

Examples

			7 is in the sequence, since 7! in the considered factorization is 5*7*9*16, and here we have 2 primes and 2 nonprimes.
		

References

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

Crossrefs

Extensions

Terms a(7) - a(9) from Peter J. C. Moses, Apr 19 2014

A235623 Numbers n for which in the prime power factorization of n!, the numbers of exponents 1 and >1 are equal.

Original entry on oeis.org

0, 1, 4, 7, 8, 9, 13, 19, 20, 21
Offset: 1

Views

Author

Vladimir Shevelev, Apr 20 2014

Keywords

Comments

Number n is in the sequence, if and only if pi(n) = 2*pi(n/2), where pi(x) is the number of primes<=x. Indeed, all primes from interval (n/2, n] appear in prime power factorization of n! with exponent 1, while all primes from interval (0, n/2] appear in n! with exponents >1. However, it follows from Ehrhart's link that, for n>=22, pi(n) < 2*pi(n/2). Therefore, a(9)=21 is the last term of the sequence.
m is in this sequence if and only if the number of prime divisors of [m/2]! equals the number of unitary prime divisors of m! - Peter Luschny, Apr 29 2014

Examples

			21! = 2^20*3^9*5^4*7^3*11*13*17*19. Here 4 primes with exponent 1 and 4 primes with exponents >1, so 21 is in the sequence.
		

Crossrefs

Programs

  • Maple
    with(numtheory): a := proc(n) factorset(n!); factorset(iquo(n,2)!);
    `if`(nops(%% minus %) = nops(%), n, NULL) end: seq(a(n), n=0..30); # Peter Luschny, Apr 28 2014
  • PARI
    isok(n) = {f = factor(n!); sum(i=1, #f~, f[i,2] == 1) == sum(i=1, #f~, f[i,2] > 1);} \\ Michel Marcus, Apr 20 2014

A241270 Numbers with the property that in their factorization over distinct terms of A050376, the sums of prime and nonprime terms of A050376 are equal.

Original entry on oeis.org

126, 468, 624, 792, 880, 1056, 1150, 2900, 3264, 4606, 5824, 6375, 6624, 8320, 9856, 10388, 11375, 12798, 13650, 16400, 16704, 19250, 20925, 30135, 32625, 36720, 39150, 39900, 53784, 56446, 56925, 57000, 59500, 63455, 65520, 71400, 71500, 72471
Offset: 1

Views

Author

Vladimir Shevelev, Apr 18 2014

Keywords

Comments

The corresponding sequence of the sum over the primes, which equals the sum over the nonprimes, is 9, 13, 16, 13, 16, 16, 25, 29, 20, 49, 20, 25, 25, 20, 20, 53, 25, 81, 25, 41, 29, 25, 34, 49, 34, 25, 34, 29, 85, 169, 34, 29, 29, 49, 25, 29, 29, 49, ... - Wolfdieter Lang, Apr 25 2014

Examples

			126 and 468 are in the sequence since the factorizations are 2*7*9 and 4*9*13 respectively, and 2+7=9, 4+9=13.
		

References

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

Crossrefs

Extensions

More terms from Peter J. C. Moses, Apr 18 2014
New extension from Wolfdieter Lang, Apr 25 2014

A265626 Let S be the set of factorizations of n! where the largest factor is the largest prime <= n, and let f(s) be the least factor in the factorization s. Then a(n) = max f(S).

Original entry on oeis.org

2, 2, 2, 2, 3, 4, 3, 3, 4, 7, 6, 7, 7, 7, 7, 10, 11, 11, 11, 11, 11, 13, 13, 13, 13, 13, 13, 17, 17, 17, 17, 17, 17, 17, 17, 19, 19, 19, 19, 23, 23, 23, 23, 23, 23, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 31, 31, 31, 31, 31, 31, 31, 31, 37, 37
Offset: 2

Views

Author

Keywords

Examples

			2! = 2
3! = 2 * 3
4! = 2^3 * 3
5! = 2 * 3 * 4 * 5
6! = 3^2 * 4^2 * 5
7! = 4 * 5 * 6^2 * 7
8! = 3 * 4^3 * 5 * 6 * 7
9! = 3 * 4^2 * 5 * 6^3 * 7
10! = 4^2 * 5^2 * 6^4 * 7
11! = 7 * 8^2 * 9^2 * 10^2 * 11
12! = 6^5 * 7 * 8 * 10^2 * 11
13! = 7 * 8^2 * 9^2 * 10^2 * 11 * 12 * 13
14! = 7^2 * 8 * 9 * 10^2 * 11 * 12^3 * 13
15! = 7^2 * 9 * 10^3 * 11 * 12^4 * 13
16! = 7^2 * 10^3 * 11 * 12^6 * 13
17! = 10 * 11 * 12^4 * 13 * 14^2 * 15^2 * 16 * 17
		

Crossrefs

Programs

  • PARI
    f(n,mn,mx)=if(n%mn,return(0)); n/=mn; if(n==1, return(1)); for(k=mn,mx, if(f(n,k,mx), return(1))); 0
    a(n)=if(n<6,return(2)); my(p=precprime(n),q=nextprime(p/2),N=n!); forprime(r=q+1,p-1, N/=r^valuation(N,r)); forstep(k=q,3,-1, if(f(N,k,p), return(k)))

Formula

a(n) > 5 for n > 10. a(n) < A007918(ceiling(A007917(n)/2)).
Showing 1-7 of 7 results.