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-5 of 5 results.

A034836 Number of ways to write n as n = x*y*z with 1 <= x <= y <= z.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 3, 2, 2, 1, 4, 1, 2, 2, 4, 1, 4, 1, 4, 2, 2, 1, 6, 2, 2, 3, 4, 1, 5, 1, 5, 2, 2, 2, 8, 1, 2, 2, 6, 1, 5, 1, 4, 4, 2, 1, 9, 2, 4, 2, 4, 1, 6, 2, 6, 2, 2, 1, 10, 1, 2, 4, 7, 2, 5, 1, 4, 2, 5, 1, 12, 1, 2, 4, 4, 2, 5, 1, 9, 4, 2, 1, 10, 2, 2, 2, 6, 1, 10, 2, 4, 2, 2, 2, 12, 1, 4, 4, 8
Offset: 1

Views

Author

Erich Friedman

Keywords

Comments

Number of boxes with integer edge lengths and volume n.
Starts the same as, but is different from, A033273. First values of n such that a(n) differs from A033273(n) are 36,48,60,64,72,80,84,90,96,100. - Benoit Cloitre, Nov 25 2002
a(n) depends only on the signature of n; the sorted exponents of n. For instance, a(12) and a(18) are the same because both 12 and 18 have signature (1,2). - T. D. Noe, Nov 02 2011
Number of 3D grids of n congruent cubes, in a box, modulo rotation (cf. A007425 and A140773 for boxes instead of cubes; cf. A038548 for the 2D case). - Manfred Boergens, Apr 06 2021

Examples

			a(12) = 4 because we can write 12 = 1*1*12 = 1*2*6 = 1*3*4 = 2*2*3.
a(36) = 8 because we can write 36 = 1*1*36 = 1*2*18 = 1*3*12 = 1*4*9 = 1*6*6 = 2*2*9 = 2*3*6 = 3*3*4.
For n = p*q, p < q primes: a(n) = 2 because we can write n = 1*1*pq = 1*p*q.
For n = p^2, p prime: a(n) = 2 because we can write n = 1*1*p^2 = 1*p*p.

Links

Crossrefs

See also: writing n = x*y (A038548, unordered, A000005, ordered), n = x*y*z (this sequence, unordered, A007425, ordered), n = w*x*y*z (A007426, ordered)
Cf. A001358, A008578, A010057, A046951, A088432, A088433, A088434.
Differs from A033273 and A226378 for the first time at n=36.

Programs

  • Maple
    f:=proc(n) local t1,i,j,k; t1:=0; for i from 1 to n do for j from i to n do for k from j to n do if i*j*k = n then t1:=t1+1; fi; od: od: od: t1; end;
# second Maple program:
A034836:=proc(n)
   local a,b,i;
   a:=0;
   b:=(l,x,h)->l<=x and x<=h;
   for i in select(`<=`,NumberTheory:-Divisors(n),iroot(n,3)) do
      a:=a+nops(select[2](b,i,NumberTheory:-Divisors(n/i),isqrt(n/i)))
   od;
   return a
end proc;
seq(A034836(n),n=1..100); # Felix Huber, Oct 02 2024
  • Mathematica
    Table[c=0; Do[If[i<=j<=k && i*j*k==n,c++],{i,t=Divisors[n]},{j,t},{k,t}]; c,{n,100}] (* Jayanta Basu, May 23 2013 *)
(* Similar to the first Mathematica code but with fewer steps in Do[..] *)
b=0; d=Divisors[n]; r=Length[d];
Do[If[d[[h]] d[[i]] d[[j]]==n, b++], {h, r}, {i, h, r}, {j, i, r}]; b (* Manfred Boergens, Apr 06 2021 *)
a[1] = 1; a[n_] := Module[{e = FactorInteger[n][[;; , 2]]}, If[IntegerQ[Surd[n, 3]], 1/3, 0] + (Times @@ ((e + 1)*(e + 2)/2))/6 + (Times @@ (Floor[e/2] + 1))/2]; Array[a, 100] (* Amiram Eldar, Apr 19 2024 *)
  • PARI
    A038548(n)=sumdiv(n, d, d*d<=n) /* <== rhs from A038548 (Michael Somos) */
a(n)=sumdiv(n, d, if(d^3<=n, A038548(n/d) - sumdiv(n/d, d0, d0Rick L. Shepherd, Aug 27 2006
  • PARI
    a(n) = {my(e = factor(n)[,2]); (2 * ispower(n, 3) + vecprod(apply(x -> (x+1)*(x+2)/2, e)) + 3 * vecprod(apply(x -> x\2 + 1, e))) / 6;} \\ Amiram Eldar, Apr 19 2024

Formula

From Ton Biegstraaten, Jan 04 2016: (Start)
Given a number n, let s(1),...,s(m) be the signature list of n, and a(n) the resulting number in the sequence.
Then np = Product_{k=1..m} binomial(2+s(k),2) is the total number of products solely based on the combination of exponents. The multiplicity of powers is not taken into account (e.g., all combinations of 1,2,4 (6 times) but (2,2,2) only once). See next formulas to compute corrections for 3rd and 2nd powers.
Let ntp = Product_{k=1..m} (floor((s(k) - s(k) mod(3))/s(k))) if the number is a 3rd power or not resulting in 1 or 0.
Let nsq = Product_{k=1..m} (floor(s(k)/2) + 1) is the number of squares.
Conjecture: a(n) = (np + 3*(nsq - ntp) + 5*ntp)/6 = (np + 3*nsq + 2*ntp)/6.
Example: n = 1728; s = [3,6]; np = 10*28 = 280; nsq = 2*4 = 8; ntp = 1 so a(1728)=51 (as in the b-file).
(End)
a(n) >= A226378(n) for all n >= 1. - Antti Karttunen, Aug 30 2017
From Bernard Schott, Dec 12 2021: (Start)
a(n) = 1 iff n = 1 or n is prime (A008578).
a(n) = 2 iff n is semiprime (A001358) (see examples). (End)
a(n) = (2 * A010057(n) + A007425(n) + 3 * A046951(n))/6 (Andrica and Ionascu, 2013, p. 19, eq. 11). - Amiram Eldar, Apr 19 2024

Extensions

Definition simplified by Jonathan Sondow, Oct 03 2013

A122179 Number of ways to write n as n = x*y*z with 1

Original entry on oeis.org

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

Views

Author

Rick L. Shepherd, Aug 23 2006

Keywords

Comments

x,y,z are proper factors of n. a(n) > 0 iff n is a term of A033942; a(n) = 0 iff n is a term of A037143.

Examples

			a(24) = 2 because 24 = 2*2*6 = 2*3*4, two products of three proper factors of 24.

Links

Crossrefs

Cf. A034836, A033942, A037143, A088432, A088433, A088434, A122180.

Programs

  • PARI
    for(n=1,105, t=0; for(x=2,n-1, for(y=x,n-1, for(z=y,n-1, if(x*y*z==n, t++)))); print1(t,", "))
  • PARI
    A122179(n) = { my(s=0); fordiv(n, x, if((x>1)&&(xAntti Karttunen, Aug 24 2017

A122180 Number of ways to write n as n = x*y*z with 1 < x < y < z < n.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 3, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 3, 0, 0, 0, 0, 0, 1, 0, 2, 0, 0, 0, 3, 0, 0, 0, 1, 0, 3, 0, 0, 0, 0, 0, 4, 0, 0, 0, 1, 0, 1, 0, 1, 1
Offset: 1

Views

Author

Rick L. Shepherd, Aug 23 2006

Keywords

Comments

x,y,z are distinct proper factors of n. See A122181 for n such that a(n) > 0.
If n has at most five divisors then a(n) = 0. - David A. Corneth, Oct 24 2024

Examples

			a(48) = 2 because 48 = 2*3*8 = 2*4*6, two products of three distinct proper factors of 48.

Links

Crossrefs

Cf. A034836, A088432, A088433, A088434, A122179, A122181, A200214.

Programs

  • PARI
    for(n=1,105, t=0; for(x=2,n-1, for(y=x+1,n-1, for(z=y+1,n-1, if(x*y*z==n, t++)))); print1(t,", "))
  • PARI
    A122180(n) = { my(s=0); fordiv(n, x, if((x>1)&&(xAntti Karttunen, Jul 08 2017
  • PARI
    a(n) = {
	my(d = divisors(n));
	if(#d <= 5, return(0));
	my(res = 0, q);
	for(i = 2, #d,
		q = d[#d + 1 - i];
		if(d[i]^2 > q,
			return(res)
		);
		for(j = i + 1, #d,
			qj = q/d[j];
			if(qj <= d[j],
				next(2)
			);
			if(denominator(qj) == 1 && n % qj == 0,
				res++
			);
		);
	);
	res
} \\ David A. Corneth, Oct 24 2024

Formula

a(n) = A200214(n)/6. - Antti Karttunen, Jul 08 2017

A088434 Number of ways to write n as n = u*v*w with 1 <= u < v < w.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 2, 0, 1, 1, 1, 0, 2, 0, 2, 1, 1, 0, 4, 0, 1, 1, 2, 0, 4, 0, 2, 1, 1, 1, 4, 0, 1, 1, 4, 0, 4, 0, 2, 2, 1, 0, 6, 0, 2, 1, 2, 0, 4, 1, 4, 1, 1, 0, 8, 0, 1, 2, 3, 1, 4, 0, 2, 1, 4, 0, 8, 0, 1, 2, 2, 1, 4, 0, 6, 1, 1, 0, 8, 1, 1, 1, 4, 0, 8, 1, 2, 1, 1, 1, 9, 0, 2, 2, 4, 0, 4, 0, 4, 4, 1, 0, 8, 0, 4, 1, 6, 0, 4, 1, 2, 2, 1, 1, 14
Offset: 1

Views

Author

Reinhard Zumkeller, Oct 01 2003

Keywords

Comments

a(n)=0 iff n=1 or n prime or n prime^2: a(A000430(n)) = 0.
The integers a(n)+1 equal A045778(n) for n < 120 and differ at all n that admit factorization into 4 or more distinct factors, the smallest ones being n = 120 = 2*3*4*5, n = 144 = 2*3*4*6, n = 168 = 2*3*4*7, n = 180 = 2*3*5*6, ..., later continuing n = 312 = 2*3*4*13, n = 320 = 2*4*5*8, n = 324 = 2*3*6*9, n = 330 = 2*3*5*11, ... Coincidentally, A068350(5) to A068350(19) start this list. - R. J. Mathar, Jul 19 2007

Examples

			n=12: (1,2,6), (1,3,4): therefore a(12)=2;
n=18: (1,2,9), (1,3,6): therefore a(18)=2.

Links

Crossrefs

Cf. A034836, A088432, A088433, A122179, A122180.

Programs

  • Mathematica
    Table[Length[Select[Cases[Subsets[Divisors[n],{3}],{x_,y_,z_}->x*y*z],#==n &]],{n,102}] (* Jayanta Basu, May 23 2013 *)
  • PARI
    A088434(n) = { my(s=0); fordiv(n, u, for(v=u+1, n-1, for(w=v+1, n, if(u*v*w==n, s++)))); (s); }; \\ Antti Karttunen, Aug 24 2017

Extensions

Data section extended to 120 terms by Antti Karttunen, Aug 24 2017

A088433 Number of ways to write n as n = u*v*w with 1<=u<=v

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 4, 1, 2, 2, 3, 1, 3, 1, 4, 2, 2, 1, 6, 1, 2, 2, 4, 1, 5, 1, 4, 2, 2, 2, 7, 1, 2, 2, 6, 1, 5, 1, 4, 4, 2, 1, 8, 1, 3, 2, 4, 1, 6, 2, 6, 2, 2, 1, 10, 1, 2, 4, 5, 2, 5, 1, 4, 2, 5, 1, 11, 1, 2, 3, 4, 2, 5, 1, 9, 3, 2, 1, 10, 2, 2, 2, 6, 1, 10, 2, 4, 2, 2, 2, 12, 1, 3, 4, 6
Offset: 1

Views

Author

Reinhard Zumkeller, Oct 01 2003

Keywords

Comments

a(n) = 1 iff n prime or n prime^2: a(A000430(n))=1.

Examples

			n=12: (1,1,12), (1,2,6), (1,3,4), (2,2,3): therefore a(12)=4;
n=18: (1,1,18), (1,2,9), (1,3,6): therefore a(18)=3.

Links

Crossrefs

Cf. A034836, A088432, A088434, A122179, A122180.

Programs

  • PARI
    A088433(n) = { my(s=0); fordiv(n, u, for(v=u, n-1, for(w=v+1, n, if(u*v*w==n, s++)))); (s); }; \\ Antti Karttunen, Aug 24 2017
Showing 1-5 of 5 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.