A233231
a(n) = 10*a(n-3) - a(n-6) + 4 for n>5, a(0)=2, a(1)=3, a(2)=5, a(3)=12, a(4)=29, a(5)=51.
Original entry on oeis.org
2, 3, 5, 12, 29, 51, 122, 291, 509, 1212, 2885, 5043, 12002, 28563, 49925, 118812, 282749, 494211, 1176122, 2798931, 4892189, 11642412, 27706565, 48427683, 115248002, 274266723, 479384645, 1140837612, 2714960669, 4745418771, 11293128122, 26875339971
Offset: 0
a(5)=29 as the triangle with sides (5,29,30) has integer area 72.
-
seq[n_] := seq[n]=Which[n==0, 2, n==1, 3, n==2, 5, n==3, 12, n==4, 29, n==5, 51, True, 10seq[n-3]-seq[n-6]+4]; Table[seq[m], {m, 0, 100}]
LinearRecurrence[{1, 0, 10, -10, 0, -1, 1}, {2, 3, 5, 12, 29, 51, 122}, 30] (* T. D. Noe, Dec 09 2013 *)
A233232
Primes p such that a Heronian triangle with a fixed side length of 5 contains p as another side length.
Original entry on oeis.org
3, 5, 13, 29, 509, 1213, 4892189, 111790443613, 8585304626467575518951161931678989213, 196181078582773936644856635510156388051229, 18087581947968179558090719299773079036323829315869
Offset: 1
a(4)=29 because the triangle with sides (5, 29, 30) is Heronian, 29 is prime and is the 4th occurrence of such a prime.
-
seq[n_] := seq[n]=Which[n==0, 2, n==1, 3, n==2, 5, n==3, 12, n==4, 29, n==5, 51, True, 10seq[n-3]-seq[n-6]+4]; lst={}; Do[Which[PrimeQ[seq[m]], AppendTo[lst, seq[m]], PrimeQ[seq[m]+1], AppendTo[lst, seq[m]+1], True, Null], {m, 1, 400}]; lst
t = LinearRecurrence[{1, 0, 10, -10, 0, -1, 1}, {2, 3, 5, 12, 29, 51, 122}, 400]; Select[Union[t, t + 1], PrimeQ[#] &] (* T. D. Noe, Dec 09 2013 *)
A237518
Least primes that together with prime(n) forms a Heronian triangle, starting at n = 2.
Original entry on oeis.org
5, 3, 4729, 13, 5, 17, 37, 5280071830550089, 5, 97, 13, 17, 61, 1824001, 53, 109, 11, 3301, 1009, 19, 241, 241, 17, 11, 29, 409, 6841, 11, 17, 3169, 181, 41, 157, 3, 457, 13, 10369, 231781748893580717709514473745694370721, 173, 277, 19, 7297, 31, 53, 3049, 373
Offset: 2
a(18)=11 as prime(18)=61, the triple (11, 60, 61) is Heronian and right angled with area=330 and 61 is the least such prime. prime(18)=61==1 mod 4 and a(18)=11==3 mod 4 and prime(18)>a(18).
- Noam Elkies, Heronian triangles with two sides that are prime, Mathoverflow, 2013.
- Eugen J. Ionascu, Florian Luca, Pantelimon Stanica, Heron triangles with two fixed sides, arXiv:math/0608185 [math.NT], 2006.
- Eugen J. Ionascu, Florian Luca, Pantelimon Stanica, Heron triangles with two fixed sides, Journal of Number Theory, Volume 126, Issue 1, September 2007, Pages 52-67.
- Pantelimon Stanica, Santanu Sarkar, Sourav Sen Gupta, Subhamoy Maitra, and Nirupam Kar, Counting Heron triangles with Constraints, Integers, Vol. 13, 2013.
-
maxn = 150000; nn=Prime[Range[maxn]]; lst={}; nn1=Prime[Range[2, 100]]; Do[Do[s=(a+b+c)/2; If[IntegerQ[s], area2=s(s-a)(s-b)(s-c); If[area2>0 && IntegerQ[Sqrt[area2]], (AppendTo[lst, b]; Break[])]]; If[b==Prime[maxn], AppendTo[lst, 0]; Break[]], {b, nn}, {a, b-c+2, b+c-2, 2}], {c, nn1}]; lst
(* 1st Program *)
q=23; d=1; nextpair[{y0_, x0_}] := (y=23; x=4; y1=y*y0+x*x0*33; x1=x0*y+y0*x; {y1, x1}); pair=nextpair[{0, q}]; While[!PrimeQ[(pair[[2]]-d)/2] && !PrimeQ[(pair[[2]]-d)/2+d], pair=nextpair[pair]]; primepair={(pair[[2]]-d)/2, (pair[[2]]-d)/2+d}; primepair(* 2nd Program *)
q=167; d=25; y=88751; x=2150; nextpair[{y0_, x0_}] := (If[IntegerQ[(q^2-d^2)/16], k=(q^2-d^2)/16, k=(q^2-d^2)/4]; y1=y*y0+x*x0*k; x1=x0*y+y0*x; {y1, x1}); pair=nextpair[{0, q}]; While[!PrimeQ[(pair[[2]]-d)/2] && !PrimeQ[(pair[[2]]+d)/2], pair=nextpair[pair]]; primepair={(pair[[2]]-d)/2, (pair[[2]]+d)/2}; primepair(* 3rd Program *)
Showing 1-3 of 3 results.
Comments