OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
5823 is a term as 5823, 5824 and 5825 are divisible by 9, 16 and 25 respectively.
MAPLE
with(numtheory):for n from 1 to 200000 do: k:=0:q:=floor(sqrt(n)):for m from
2 to q do: p1:=m^2:p2:=(m+1)^2:p3:=(m+2)^2:if irem(n, p1)=0 and irem(n+1, p2)=0
and irem(n+2, p3)=0 and k=0 then k:=1:printf(`%d, `, n):else fi:od:od:
MATHEMATICA
f[p_, e_] := p^Floor[e/2]; maxsq[n_] := Times@@ (f @@@ FactorInteger[n]); aQ[n_] := (s = maxsq[n]) > 1 && AnyTrue[Rest @ Divisors[s], Divisible[n+1, (#+1)^2] && Divisible[n+2, (#+2)^2] &]; Select[Range[50000], aQ] (* Amiram Eldar, Sep 09 2019 *)
PROG
(Sage) is_A178919 = lambda n: any(all(((k+x)**2).divides(n+x) for x in range(3)) for k in divisors(n) if k > 1) # D. S. McNeil, Dec 29 2010
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Dec 29 2010
STATUS
approved