OFFSET
1,1
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..37 (terms < 3*10^12)
EXAMPLE
25 is a term as p = 2*5 = 10 and (25-10)*(25+10) = 525 which contains '25' as a substring.
8926 is a term as p = 8*9*2*6 = 864 and (8926-864)*(8926+864) = 78926980 which contains '8926' as a substring.
PROG
(PARI) isokp(dx, d) = {if (!#setintersect(Set(dx), Set(d)), return (0)); for (i=1, #dx - #d + 1, if (vector(#d, k, dx[k+i-1]) == d, return(1)); ); }
isokd(x, d, n) = {if (x==n, return (1)); my(dx = digits(x)); if (#dx < #d, return (0)); isokp(dx, d); }
isok(n) = {my(d = digits(n), p = vecprod(d)); if (p>0, isokd((n-p)*(n+p), d, n)); } \\ Michel Marcus, May 07 2020
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Scott R. Shannon, May 05 2020
EXTENSIONS
More terms from Giovanni Resta, May 07 2020
STATUS
approved