? for(m=2,10^5,for(n=1,m-1,p2=m*(m+5)/(n*(n+5));if(issquare(p2)&&p2==floor(p2),p=sqrtint(p2);if(isprime(p),print1("("n","m","p"), ")))))
(1,3,2), (2,9,3), (1,10,5), (7,16,2), (6,22,3), (9,121,11), (19,147,7),