A088111 -id:A088111 - cp's OEIS frontend

A046080 a(n) is the number of integer-sided right triangles with hypotenuse n.

0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 2, 1, 0, 0, 1, 1, 0, 0, 0, 1, 1, 0, 1, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 2, 1, 1, 1, 0, 1, 0, 0, 1, 0, 1, 1, 0, 0, 0, 4, 0, 0, 1, 0, 1, 0, 0, 1, 1, 2, 0, 0, 1, 0, 1, 0, 1, 0, 0, 4, 0, 1, 0, 1, 1, 1, 0, 0, 0, 1, 0, 1, 0, 0

Offset: 1

Eric W. Weisstein

nonn

Or number of ways n^2 can be written as the sum of two positive squares: a(5) = 1: 3^2 + 4^2 = 5^2; a(25) = 2: 7^2 + 24^2 = 15^2 + 20^2 = 25^2. - Alois P. Heinz, Aug 01 2019

A. H. Beiler, Recreations in the Theory of Numbers, New York: Dover, pp. 116-117, 1966.

Stanislav Sykora, Table of n, a(n) for n = 1..20000
Ron Knott, Pythagorean Triples and Online Calculators
F. Richman, Pythagorean Triples
A. Tripathi, On Pythagorean triples containing a fixed integer, Fib. Q., 46/47 (2008/2009), 331-340. See Theorem 7.
Eric Weisstein's World of Mathematics, Pythagorean Triple

First differs from A083025 at n=65.
Cf. A000290, A006339, A024362, A046079, A046081, A009000.
A088111 gives records; A088959 gives where records occur.
Cf. A046109, A063725.
Cf. A004144, A084647, A084648, A084649.
Partial sums: A224921.

f:= proc(n) local F,t;
  F:= select(t -> t[1] mod 4 = 1, ifactors(n)[2]);
  1/2*(mul(2*t[2]+1, t=F)-1)
end proc:
map(f, [$1..100]); # Robert Israel, Jul 18 2016

a[1] = 0; a[n_] := With[{fi = Select[ FactorInteger[n], Mod[#[[1]], 4] == 1 & ][[All, 2]]}, (Times @@ (2*fi+1)-1)/2]; Table[a[n], {n, 1, 99}] (* Jean-François Alcover, Feb 06 2012, after first formula *)

a(n)={my(m=0,k=n,n2=n*n,k2,l2);
while(1,k=k-1;k2=k*k;l2=n2-k2;if(l2>k2,break);if(issquare(l2),m++));return(m)} \\ brute force, Stanislav Sykora, Mar 18 2015

{a(n) = if( n<1, 0, sum(k=1, sqrtint(n^2 \ 2), issquare(n^2 - k^2)))}; /* Michael Somos, Mar 29 2015 */

a(n) = {my(f = factor(n/(2^valuation(n, 2)))); (prod(k=1, #f~, if ((f[k,1] % 4) == 1, 2*f[k,2] + 1, 1)) - 1)/2;} \\ Michel Marcus, Mar 08 2016

from math import prod
from sympy import factorint
def A046080(n): return prod((e<<1)+1 for p,e in factorint(n).items() if p&3==1)>>1 # Chai Wah Wu, Sep 06 2022

Let n = 2^e_2 * product_i p_i^f_i * product_j q_j^g_j where p_i == 1 mod 4, q_j == 3 mod 4; then a(n) = (1/2)*(product_i (2*f_i + 1) - 1). - Beiler, corrected
8*a(n) + 4 = A046109(n) for n > 0. - Ralf Stephan, Mar 14 2004
a(n) = 0 for n in A004144. - Lekraj Beedassy, May 14 2004
a(A084645(k)) = 1. - Ruediger Jehn, Jan 14 2022
a(A084646(k)) = 2. - Ruediger Jehn, Jan 14 2022
a(A084647(k)) = 3. - Jean-Christophe Hervé, Dec 01 2013
a(A084648(k)) = 4. - Jean-Christophe Hervé, Dec 01 2013
a(A084649(k)) = 5. - Jean-Christophe Hervé, Dec 01 2013
a(n) = A063725(n^2) / 2. - Michael Somos, Mar 29 2015
a(n) = Sum_{k=1..n} Sum_{i=1..k} [i^2 + k^2 = n^2], where [ ] is the Iverson bracket. - Wesley Ivan Hurt, Dec 10 2021
a(A002144(k)^n) = n. - Ruediger Jehn, Jan 14 2022

A088959 Lowest numbers which are d-Pythagorean decomposable, i.e., square is expressible as sum of two positive squares in more ways than for any smaller number.

Original entry on oeis.org

1, 5, 25, 65, 325, 1105, 5525, 27625, 32045, 160225, 801125, 1185665, 5928325, 29641625, 48612265, 243061325, 1215306625, 2576450045, 12882250225, 64411251125, 157163452745, 785817263725, 3929086318625, 10215624428425, 11472932050385, 51078122142125

Offset: 1

Lekraj Beedassy, Dec 01 2003

nonn

These are also the integer radii of circles around the origin that contain record numbers of lattice points. See A071383 for radii that are not necessarily integer. - Günter Rote, Sep 14 2023

			From _Petros Hadjicostas_, Jul 21 2019: (Start)
Squares 1^2, 2^2, 3^2, and 4^2 have 0 representations as the sum of two positive squares. (Thus, A088111(1) = 0 for the number of representations of 1^2.) Thus a(1) = 1.
Square 5^2 can be written as 3^2 + 4^2 only (here A088111(2) = 1). Thus, a(2) = 5.
Looking at sequence A046080, we see that for 5 <= n <= 24, only n^2 = 5^2, 10^2, 13^2, 15^2, 17^2, 20^2 can be written as a sum of two positive squares (in a single way) because 5^2 = 3^2 + 4^2, 10^2 = 6^2 + 8^2, 13^2 = 5^2 + 12^2, 17^2 = 8^2 + 15^2, and 20^2 = 12^2 + 16^2.
Since A046080(25) = 2 and A088111(3) = 2, we have that 25^2 can be written as a sum of two positive squares in two ways. Indeed, 25^2 = 7^2 + 24^2 = 15^2 + 20^2. Thus, a(3) = 25.
For 26 <= n <= 64, we see from sequence A046080 that n^2 cannot be written in more than 2 ways as a sum of two positive squares.
Since A046080(65) = 4, we see that 65^2 can be written as the sum of two positive squares in 4 ways. Indeed, 65^2 = 16^2 + 63^2 = 25^2 + 60^2 = 33^2 + 56^2 = 39^2 + 52^2. Thus, a(4) = 65.
(End)

R. M. Sternheimer, Additional Remarks Concerning The Pythagorean Triplets, Journal of Recreational Mathematics, Vol. 30, No. 1, pp. 45-48, 1999-2000, Baywood NY.

Ray Chandler, Table of n, a(n) for n = 1..307

Cf. A052199. Subsequence of A054994. Number of ways: see A088111. Where records occur in A046080.

from math import prod
from sympy import isprime
primes_congruent_1_mod_4 = [5]
def prime_4k_plus_1(i): # the i-th prime that is congruent to 1 mod 4
    while i>=len(primes_congruent_1_mod_4): # generate primes on demand
        n = primes_congruent_1_mod_4[-1]+4
        while not isprime(n): n += 4
        primes_congruent_1_mod_4.append(n)
    return primes_congruent_1_mod_4[i]
def generate_A054994():
    TO_DO = {(1,())}
    while True:
        radius, exponents = min(TO_DO)
        yield radius, exponents
        TO_DO.remove((radius, exponents))
        TO_DO.update(successors(radius,exponents))
def successors(r,exponents):
    for i,e in enumerate(exponents):
        if i==0 or exponents[i-1]>e:
            yield (r*prime_4k_plus_1(i), exponents[:i]+(e+1,)+exponents[i+1:])
    if exponents==() or exponents[-1]>0:
        yield (r*prime_4k_plus_1(len(exponents)), exponents+(1,))
n,record=0,-1
for radius,expo in generate_A054994():
    num_pyt = (prod((2*e+1) for e in expo)-1)//2
    if num_pyt>record:
        record = num_pyt
        n += 1
        print(radius, end="") # or record, for A088111
        if n==26: break # stop after 26 entries
        print(end=", ")
print() # Günter Rote, Sep 13 2023

Corrected and extended by Ray Chandler, Jan 12 2012

Name edited by Petros Hadjicostas, Jul 21 2019

cp's OEIS Frontend

A046080 a(n) is the number of integer-sided right triangles with hypotenuse n.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

PARI

Python

Formula

A088959 Lowest numbers which are d-Pythagorean decomposable, i.e., square is expressible as sum of two positive squares in more ways than for any smaller number.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Python

Extensions

A365620 Number of integer grid points on the circle around (0,0) with radius A088959(n).

Views

Author

Keywords

Crossrefs

Formula