OFFSET
1,1
COMMENTS
Though the first 27 terms match those of A089392, the next term of A089392 (2221) is not a member of this sequence. Conjecture: sequence is finite.
No more terms < 10^8. - David Wasserman, Oct 04 2005
EXAMPLE
863 is a member 863, 8 + 63, 86 + 3, 8 + 6 + 3 are all prime.
MAPLE
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=choose([seq(j, j=2..d)]))]: for n from 10^(d-1) to 10^d-1 do sn:=convert(n, base, 10): fl:=0: for s in sch do m:=add(j, j=[seq(ds(sn[s[i]..s[i+1]-1]), i=1..nops(s)-1)]): if not isprime(m) then fl:=1: break fi od: if fl=0 then printf("%d, ", n) fi od od: # C. Ronaldo
# second Maple program:
b:= proc(s) option remember; (n-> {s, seq(seq(seq(""||x||"+"||y,
y=b(s[i+1..n])), x=b(s[1..i])), i=1..n-1)})(length(s))
end:
q:= n-> andmap(isprime, map(parse, b(""||n))):
select(q, [$1..10000])[]; # Alois P. Heinz, Oct 29 2021
MATHEMATICA
Select[Prime@Range@1010, And@@PrimeQ[n=#; Total/@(FromDigits/@#&/@Union[DeleteCases[SplitBy[#, #==-1&], {-1}]&/@(Insert[IntegerDigits@n, -1, #]&/@(List/@#&/@Rest@Subsets[Range@IntegerLength@n]))])]&] (* Giorgos Kalogeropoulos, Oct 29 2021 *)
CROSSREFS
KEYWORD
base,nonn,more
AUTHOR
Amarnath Murthy, Nov 10 2003
EXTENSIONS
Corrected and extended by C. Ronaldo (aga_new_ac(AT)hotmail.com), Dec 25 2004
STATUS
approved