OFFSET
1,1
COMMENTS
All such numbers may be written as n*(n+1)-k or n*(n+1)+k with k odd < n and n > 1.
The sequence may be constructed by reading the infinite square array A(n,k) by upward antidiagonals where A(n,k) = (n+k-1)*(n+k)+2*k-1 for k < n and A(n,k) = (n+k)*(n+k+1)-2*n+1 for k >= n.
LINKS
Paolo Xausa, Table of n, a(n) for n = 1..10000
EXAMPLE
As an array A(n,k) read by upward antidiagonals:
5 11 19 29 41 ...
7 17 27 39 53 ...
13 23 37 51 67 ...
21 33 47 65 83 ...
31 45 61 79 101 ...
...
MATHEMATICA
Select[Range[1, 181, 2], NoneTrue[Sqrt[{#, # + 1}], IntegerQ] &] (* or, as an array: *)
A113909[n_, k_] := If[k < n, #*(# - 1) + 2*k - 1, #*(# + 1) - 2*n + 1] & [n + k];
Table[A113909[n - k + 1, k], {n, 12}, {k, n}] (* Paolo Xausa, Aug 22 2025 *)
PROG
(PARI) isok(n)= {n%2 && !issquare(n) && !issquare(n+1)} \\ Andrew Howroyd, Aug 17 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Pierre CAMI, Jan 29 2006
EXTENSIONS
Edited by N. J. A. Sloane, Jul 18 2009
Edited by Andrew Howroyd, Aug 17 2025
STATUS
approved
