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.

A306676 Number of distinct acute triangles with prime sides and largest side = prime(n).

Original entry on oeis.org

1, 2, 3, 5, 5, 8, 9, 11, 13, 12, 18, 17, 21, 27, 30, 28, 30, 38, 38, 43, 56, 53, 59, 59, 56, 64, 79, 85, 100, 106, 79, 90, 96, 115, 102, 123, 124, 130, 144, 147, 152, 177, 161, 188, 199, 225, 193, 175, 195, 228, 248, 247, 280, 259, 267, 277, 288, 324
Offset: 1

Views

Author

César Eliud Lozada, Mar 04 2019

Keywords

Examples

			For n=3, prime(n)=5. Acute triangles: {2,5,5}, {3,5,5}, {5,5,5} (Total=3=a(3)).
For n=4, prime(n)=7. Acute triangles: {2,7,7}, {3,7,7}, {5,5,7}, {5, 7, 7}, {7, 7, 7} (Total=5=a(4)).
		

Crossrefs

Programs

  • Maple
    #nType=1 for acute triangles, nType=2 for obtuse triangles
    #nType=0 for both triangles
    CountPrimeTriangles := proc (n, nType := 1)
      local aa, oo, j, k, sg, a, b, c, tt, lAcute;
      aa := {}; oo := {};
      a := ithprime(n);
      for j from n by -1 to 1 do
        b := ithprime(j);
        for k from j by -1 to 1 do
          c := ithprime(k);
          if a < b+c and abs(b-c) < a and b < c+a and abs(c-a) < b and c < a+b and abs(a-b) < c then
            lAcute := evalb(0 < b^2+c^2-a^2);
            tt := sort([a, b, c]);
            if lAcute then aa := {op(aa), tt} else oo := {op(oo), tt} end if
          end if
        end do
      end do;
      return sort(`if`(nType = 1, aa, `if`(nType = 2, oo, `union`(aa, oo))))
    end proc: