login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A058433 Numbers k such that k^2 contains only digits {0,3,9}, not ending with zero. 2
3, 969071253 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

No more terms up to 10^23. - Charles R Greathouse IV, Jul 27 2009

LINKS

Table of n, a(n) for n=1..2.

P. De Geest, Index to related sequences

Hisanori Mishima, Sporadic tridigital solutions

PROG

(PARI) admissibleMod(M=1000, t=[3, 9], debug=0)={ local( tt=1 ); while( M > tt*=10, t=concat([t, t+vector(#t, i, t[1])*tt, t+vector(#t, i, t[2])*tt])); if( debug, print("t="t)); t=Set(t); tt=[];

for(i=1, M, if( setsearch(t, i^2%M), tt=concat(tt, i))); concat(tt, M+tt[1])} A058433(Nmax=1e10, N=1, addMod=100000, debug=1)={ local( a=[], addNext=admissibleMod(addMod=10^round(log(addMod)/log(10)), [3, 9]), d=1, add=vector(addMod, i, if(i-1>addNext[d], d++); addNext[d]-i+1), pow10 = vector( d=#Str((Nmax=round(Nmax))^2), i, 10^(i-1)), nextOK=[0, 2, 1, 0, 5, 4, 3, 2, 1, 0], Nnext, N2, numDigits, place );

nextOK = vector( #nextOK, i, if( nextOK[i], nextOK[i]*pow10)); N=round(N); while( Nmax >= N, numDigits = #Str(N2=N^2); if( place = nextOK[1+ d=N2 \ pow10[numDigits] ], N = max(sqrtint( place[numDigits]+d*pow10[numDigits] ), N+1); next); Nnext = min( Nmax, sqrtint( round((d+1)*pow10[numDigits]) ));

if( debug, print( "checking from "N" to "Nnext": <= ", 1+max(0, Nnext-N)*(#addNext-1)\ addMod, " candidates.")); N += add[1+ N%addMod]; place=1; while( Nnext >= N, dr = divrem( N2=N^2, pow10[ place=numDigits ] ); while( place-- & !d=nextOK[1+ (dr = divrem( dr[2], pow10[ place ] ))[1]], );

if( !place, break); N = sqrtint( N2 - dr[2] + d[ place ])+1; N+=add[1+N%addMod]; ); if( !place, if( debug, print( N, "^2 = ", N^2 )); a=concat(a, N)); N = Nnext*3\2+1 ); a} - M. F. Hasler, May 14 2007

CROSSREFS

Cf. A058434.

Sequence in context: A058469 A081508 A067481 * A154998 A036236 A235357

Adjacent sequences:  A058430 A058431 A058432 * A058434 A058435 A058436

KEYWORD

nonn,base,hard,more,bref

AUTHOR

Patrick De Geest, Nov 15 2000

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified September 22 06:01 EDT 2021. Contains 347605 sequences. (Running on oeis4.)