A200510 Least m>0 such that n = 10^x-y^2 (mod m) has no solution, or 0 if no such m exists.

0, 0, 3, 5, 9, 3, 0, 9, 3, 0, 0, 3, 35, 5, 3, 11, 9, 3, 5, 0, 3, 16, 9, 3, 11, 9, 3, 20, 5, 3, 16, 9, 3, 5, 9, 3, 0, 11, 3, 0, 9, 3, 20, 5, 3, 32, 11, 3, 5, 9, 3, 0, 9, 3, 28, 37, 3, 11, 5, 3, 200, 9, 3, 5, 0, 3, 16, 9, 3, 16, 9, 3, 35, 5, 3, 0, 9, 3, 5, 9

Offset

0,3

Comments

If such an m>0 exists, this proves that n is not in A051221, i.e., not of the form 10^x-y^2. On the other hand, if n is in A051221, i.e., there are integers x, y such that n=10^x-y^2, then we know that a(n)=0.

Links

M. F. Hasler, Table of n, a(n) for n = 0..398

Formula

a(111)=11111.

Example

See A200507.

Prog

(PARI) A200510(n, b=10, p=3)={ my( x=0, qr, bx, seen ); for( m=3, 9e9, while( x^p < m, issquare(b^x-n) & return(0); x++); qr=vecsort(vector(m, i, i^2+n)%m, , 8); seen=0; bx=1; until( bittest(seen+=1<<bx, bx=bx*b%m), for(i=1, #qr, qr[i]<bx & next; qr[i]>bx & break; next(3))); return(m))}

Crossrefs

Cf. A051204-A051221, A200505-A200524.

Sequence in context: A210946 A319984 A303772 * A138055 A236556 A189044

Adjacent sequences: A200507 A200508 A200509 * A200511 A200512 A200513

Keyword

nonn

Author

M. F. Hasler, Nov 18 2011

Status

approved