for(k=1,10^10,if(issquare(51+163*k),print(k))) 10 91 255 498 826 1231 1723 2290 2946 3675 4495 5386 6370 7423 8571 9786 11098 12475 13951 15490 17130 18831 20635 22498 24466 26491 28623 30810 33106 ... ...