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

A375512 a(n) is the number of distinct compositions of four positive integers x, u, v, y (x < u <= v < y) for which x + u + v + y = n and u*v = x*y.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 2, 1, 0, 2, 0, 2, 3, 0, 0, 4, 3, 0, 4, 3, 0, 6, 0, 3, 5, 0, 6, 9, 0, 0, 6, 8, 0, 9, 0, 5, 13, 0, 0, 13, 6, 7, 8, 6, 0, 11, 10, 12, 9, 0, 0, 23, 0, 0, 19, 10, 12, 15, 0, 8, 11, 18, 0, 27, 0, 0, 23, 9, 15, 18, 0, 25, 19
Offset: 0

Views

Author

Felix Huber, Aug 19 2024

Keywords

Comments

(bin(4,0) + bin(4,2) + bin(4,4))*a(n) = 8*a(n) is the number of distinct compositions of four integers x, u, v, y (abs(x) < abs(u) <= abs(v) < abs(y)) for which abs(x) + abs(u) + abs(v) + abs(y) = n and u*v = x*y.
a(n) is also the number of 2X2 matrices having the determinant 0 whose elements [x,u;v,y] are positive integers with x < u <= v < y and x + u + v + y = n.
a(n) is also the number of distinct linear 2X2 equation systems that do not have exactly one solution and whose coefficients [x,u;v,y] are positive integers with x < u <= v < y and x + u + v + y = n.

Examples

			a(9) = 1 because only (1, 2, 2, 4) satisfies the conditions: 1 + 2 + 2 + 4 = 9 and 2*2 = 1*4.
a(24) = 4 because (1, 2, 7, 14), (1, 3, 5, 15), (2, 4, 6, 12), (3, 5, 6, 10) satisfy the conditions: 1 + 2 + 7 + 14 = 24 and 2*7 = 1*14, 1 + 3 + 5 + 15 = 24 and 3*5 = 1*15, 2 + 4 + 6 + 12 = 24 and 4*6 = 2*12, 3 + 5 + 6 + 10 = 24 and 5*6 = 3*10.
See also linked Maple code.
		

Crossrefs

Programs

  • Maple
    See Huber link.
  • Python
    def A375512(n): return sum(1 for x in range(1,(n>>2)+1) for y in range(x+1,(n-x)//3+1) for z in range(y,(n-y>>1)+1) if xChai Wah Wu, Aug 23 2024

Formula

Conjecture: a(p) = 0 for primes p.
From Robert Israel, Aug 23 2024: (Start)
The conjecture is true, in fact for any x,y,u,v as in the definition, n has proper divisor gcd(x,u) + gcd(v,y).
Proof: Suppose x,y,u,v are positive integers with x + y + u + v = n and x*y = u*v = m. Let g = gcd(x,u). Then x = g*X and u = g*U where X and U are coprime. Since X*y = U*v = m/g, we must have y = h*U and v = h*X where h = gcd(v,y). Then n = g*X + h*U + g*U + h*X = (g+h)*(U+X).
(End)

A357260 a(n) is the number of 2 X 2 Euclid-reduced matrices having coprime elements and determinant n.

Original entry on oeis.org

1, 2, 3, 4, 5, 8, 7, 9, 9, 14, 11, 16, 13, 20, 18, 19, 17, 28, 19, 26, 26, 32, 23, 36, 25, 38, 31, 38, 29, 54, 31, 41, 42, 50, 38, 56, 37, 56, 50, 56, 41, 76, 43, 62, 58, 68, 47, 78, 49, 78, 66, 74, 53, 92, 62, 76, 74, 86, 59, 114, 61, 92, 78, 85, 74, 124, 67, 98, 90, 118
Offset: 1

Views

Author

Michel Marcus, Sep 21 2022

Keywords

Comments

See Bacher link for the definition of Euclid-reduced.

Crossrefs

Cf. A357259.

Programs

  • Mathematica
    f[n_] := DivisorSum[n, # + 1 - n/# &, #^2 >= n &]; a[n_] := DivisorSum[n, MoebiusMu[Sqrt[#]] * f[n/#] &, IntegerQ[Sqrt[#]] &]; Array[a, 100] (* Amiram Eldar, Sep 21 2022 *)
  • PARI
    f(n) = sumdiv(n, d, if (d^2 >= n, d + 1 -n/d)); \\ A357259
    a(n) = sumdiv(n, d, if (issquare(d), moebius(sqrtint(d))*f(n/d)));

Formula

a(n) = Sum_{d^2|n} moebius(d)*A357259(n/d^2).
Showing 1-2 of 2 results.