OFFSET
1,2
COMMENTS
Seems to be a subsequence of the positive numbers primitively represented by the binary quadratic form (1, 6, -3) with discriminant 48 (see A244291, A243168). - Peter Luschny, Jun 25 2014
LINKS
N. J. A. Sloane et al., Binary Quadratic Forms and OEIS (Index to related sequences, programs, references)
EXAMPLE
MATHEMATICA
(* b = A083752 *) b[n_] := b[n] = For[k = n+1, True, k++, If[IntegerQ[Sqrt[(4k+3n)(4n+3k)]], Return[k]]]; Reap[For[n = 1, n < 600, n++, mm = Most @ Divisors[n]; If[NoneTrue[mm, b[n] == (n/#) b[#] &], Print[n]; Sow[n]]]][[2, 1]] (* Jean-François Alcover, Oct 31 2016 *)
PROG
(Sage)
def is_A085018(n):
for d in divisors(n):
if d < n:
return false
return true
filter(is_A085018, (1..577)) # Peter Luschny, Jun 25 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Jun 18 2003
EXTENSIONS
Edited and extended by Stefan Steinerberger, Jul 30 2007
More terms from Peter Luschny, Jun 25 2014
STATUS
approved