OFFSET
1,1
COMMENTS
The first term ending in a 9 seems to be 1225449, and the first term ending in a 1 is 136161.
For 1 <= m <= 10^4 and 1 <= k <= m, there are 9217 numbers of the form (m*k)^2/(m^2-k^2). Of these numbers, only 679 are odd.
If a(n) is not a square, then m = 9*k or m = 7*k. If a(n) is a square, m does not appear to be a multiple of k.
Let a(n) be a square generated by m_1 and k_1. If a(n-1) is generated by m_2 and k_2, then k_1 = k_2 and m_1 < m_2.
The reciprocals of these numbers can be represented as the difference in the reciprocals of two squares (i.e., there exists two distinct integers m and k satisfying 1/a(n) = 1/m^2 - 1/k^2).
EXAMPLE
(84*12)^2/(84^2-12^2) = 84^2/48 = 147. So 147 is a member of this sequence. (Note that k=12 and m=84 and so m=7*k.)
PROG
(PARI) v=[]; for(m=1, 7500, for(n=1, m-1, if(type(s=(m*n)^2/(m^2-n^2))=="t_INT"&&(s%2), v=concat(v, s)))); vecsort(v, , 8)
CROSSREFS
KEYWORD
nonn
AUTHOR
Derek Orr, Jun 23 2015
STATUS
approved