3n+1小于2009，5n+1可以表示为完全平方数，求正整数n的解有多少个？

报告～！
22个～！

具体算法
delphi实现
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
procedure TForm1.Button1Click(Sender: TObject);
var
n:integer;
i:integer;
begin
n:=1;
i:=0;
while((3*n+1)<2009) do
begin
if sqrt(5*n+1)= trunc(sqrt(5*n+1)) then
begin
i:=i+1;
end;
n:=n+1;
end;
self.Edit1.Text:=intTostr(i);
end;

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝