OFFSET
1,1
COMMENTS
Numbers k such that A376301(k) = -1.
Numbers k such that every number ending in k is divisible by a prime >= 11.
If i*k is a term where 2 <= i <= 9 and i*k has the same number of digits as k, then k is a term.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
The multiplicative order of 3 mod 100 is 20, and that of 7 is 4. Of the 36 numbers from 10 to 99 that are coprime to 10, the only ones congruent to some 3^i * 7^j (mod 100) are those where i = 0 to 19 and j = 0 to 3: 21, 23, 27, 29, 41, 43, 47, 49, 61, 63, 67, 69, 81, 83, 87, and 89. The remaining 20: 11, 13, 17, 19, 31, 33, 37, 39, 51, 53, 57, 59, 71, 73, 77, 79, 91, 93, 97, 99, are members of this sequence.
MAPLE
dmax:= 3: # for terms < 10^dmax
N:= 10^dmax-1: R:= NULL:
with(priqueue):
V[0]:= 10:
for d from 1 to dmax do
A:= Array(0..10^d-1, -1);
initialize(pq);
insert([-1, 0, 0, 0, 0], pq);
while pq[0] <> 0 do
t:= extract(pq);
x:= -t[1];
xd:= x mod 10^d;
if A[xd] = -1 then
A[xd]:= x;
insert([-7*x, t[2], t[3], t[4], t[5]+1], pq);
if t[5] = 0 then
insert([-5*x, t[2], t[3], t[4]+1, 0], pq);
if t[4] = 0 then
insert([-3*x, t[2], t[3]+1, 0, 0], pq);
if t[3] = 0 then
insert([-2*x, t[2]+1, 0, 0, 0], pq);
fi fi fi fi od;
R:= R, op(select(t -> A[t] = -1, {$10^(d-1)..10^d-1}))):
od:
R;
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Robert Israel, Sep 19 2024
STATUS
approved