b = 2; s = {{1}};
Do[NestWhile[# + 1 &, 0, ! (PrimeQ[FromDigits[tmp = Join[Last[s], (nn = #; IntegerDigits[nn - Sum[b^n, {n, l = NestWhile[# + 1 &, 1, ! (nn - (Sum[b^n, {n, #}]) < 0) &] - 1}], b, l + 1])], b]]) &];
AppendTo[s, tmp], {30}]; Map[FromDigits, s]
Map[FromDigits[#, b] &, s] (*
A261201 *)