The OEIS mourns the passing of Jim Simons and is grateful to the Simons Foundation for its support of research in many branches of science, including the OEIS.
login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A367451 Squares in A353729. 2
1, 4, 9, 324, 1444, 26244, 114244, 116964, 191844, 213444, 311364, 412164, 777924, 2196324, 3392964, 6441444, 11142244, 13264164, 16499844, 21734244, 23213124, 23444964, 24621444, 26193924, 27436644, 36216324, 44116164, 47361924, 61121124, 61496964, 63393444 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
The sequence is infinite because the squares of numbers of the form (10^n + 14)/3, n >= 2, 38, 338, 3338, ... are 1444, 114244, 11142244, ... and have only the digits 1, 2, 4. Deleting any digit results in a number divisible by 4.
a(1) = 1 and a(3) = 9 are the only terms with an odd last digit. There are no other terms with the last digit 1, 5, 9 because the squares are odd, the tens digit must be even, and by deleting it the remaining number is odd.
The last digit of the terms cannot be 6 because the last two digits of the squares would be 16, 36, 56, 76, 96, so the penultimate digit is odd, and removing the 6 would leave an odd number. Thus, the terms a(n) with n >= 4, have the last digit 4.
LINKS
EXAMPLE
324 = 18^2 is a term, since 3 divides 24, 2 divides 34 and 4 divides 32.
MATHEMATICA
A353729Q[n_] := Block[{d = IntegerDigits[n], i=0}, FreeQ[d, 0] && (While[++i <= Length[d] && Divisible[FromDigits[Drop[d, {i}]], d[[i]]]]; i > Length[d])];
Select[Range[15000]^2, A353729Q] (* Paolo Xausa, Feb 27 2024 *)
PROG
(Magma) ints:=func<n| n eq 0 select[0] else Intseq(n)>; f:=func<n, i|Seqint([ints(n)[j]: j in [1..#ints(n)]|j ne (#ints(n)-i+1)])>; [p:p in [s*s:s in [1..8000]] |not 0 in ints(p) and forall{i:i in [1..#ints(p)]|IsIntegral(f(p, i)/ints(p)[(#Intseq(p)-i+1)])}];
(PARI) isok(k) = if (issquare(k), my(d=digits(k)); if (vecmin(d), for(i=1, #d, my(list=List(d)); listpop(list, i); if (fromdigits(Vec(list)) % d[i], return(0)); ); return(1); ); ); \\ Michel Marcus, Dec 19 2023
CROSSREFS
Intersection of A000290 and A353729.
Sequence in context: A085065 A068708 A370610 * A063248 A027522 A330580
KEYWORD
nonn,base
AUTHOR
Marius A. Burtea, Dec 16 2023
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified May 26 10:56 EDT 2024. Contains 372824 sequences. (Running on oeis4.)