1,1

with(combinat): ds:=proc(s) local j: RETURN(add(s[j]*10^(j-1), j=1..nops(s))):end: for d from 1 to 6 do sch:=[seq([1, op(i), d+1], i=[[], seq([j], j=2..d)])]: for n from 10^d-1 by -1 to 10^(d-1) do sn:=convert(n, base, 10): fl:=0: for s in sch do m:=mul(j, j=[seq(ds(sn[s[i]..s[i+1]-1]), i=1..nops(s)-1)])+1: if not isprime(m) then fl:=1: break fi od: if fl=0 then printf("%d, ", n):break fi od od: (C. Ronaldo)

Cf. A089392, A089393, A089394, A089395, A089396.

Sequence in context: A121790 A032491 A135624 this_sequence A089701 A068253 A155811

Adjacent sequences: A089394 A089395 A089396 this_sequence A089398 A089399 A089400

base,nonn

Amarnath Murthy (amarnath_murthy(AT)yahoo.com), Nov 10 2003

More terms from C. Ronaldo (aga_new_ac(AT)hotmail.com), Dec 25 2004