OFFSET
1,1
COMMENTS
There exists a subsequence of squares such that 16, 36, 49, 64, 81, 100, 121, 144, 169, 196, 225, 256, 289, 324, ...
There exists a subsequence of primes such that 11, 13, 19, 23, 31, 41, 47, 53, 59, 61, 67, 71, 73, 83, 89, 97, 109, 113, 131, 137, 139, 149,... but the subsequence of primes 17, 29, 37, 43, 101, 317, 433, 439, 487, 569,... is not included in the sequence.
LINKS
Michel Lagneau, Table of n, a(n) for n = 1..10000
MAPLE
with(numtheory):for n from 1 to 100 do:x:=convert(n, base, 10):n1:=nops(x): s:=sum('x[i]*10^(n1-i)', 'i'=1..n1):y:=n*10^n1+s:z:=factorset(y):n2:=nops(z):s1:=sum('z[j]', 'j'=1..n2):printf(`%d, `, s1):od:
MATHEMATICA
Join[{11}, d[n_]:=IntegerDigits[n]; Rest[Total[Transpose[FactorInteger[Plus[FromDigits[Join[x=d[#], Reverse[x]]]]]][[1]]]&/@Range[100]]]
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Michel Lagneau, Apr 06 2014
STATUS
approved