OFFSET
1,2
COMMENTS
The conventional a(1) = 0 represents the empty concatenation.
Due to simple concatenation, this encoding of the positive integers becomes ambiguous from n = 613 = prime(112)^1 on, which has the same encoding a(n) = 1121 as 6 = prime(1)^1*prime(2)^1. To get a unique encoding, one could use, e.g., the digit 9 as delimiter to separate indices and exponents, written in base 9 as to use only digits 0..8, as soon as a term would be the duplicate of an earlier term (or for all n >= 613). Then one would have, e.g., a(613) = prime(134_9)^1 = 13491.
Sequence A067599 is based on the same idea, but uses the primes instead of their indices. In A037276 the prime factors are repeated, instead of giving the exponent. In A080670 exponents 1 are omitted. In A124010 only the prime signature is given. In A054841 the sum e_i*10^(i-1) is given, i.e., exponents are used as digits in base 10, while they are listed individually in the rows of A067255.
EXAMPLE
2 = prime(1)^1 => a(2) = 11,
3 = prime(2)^1 => a(3) = 21,
4 = prime(1)^2 => a(4) = 12,
5 = prime(3)^1 => a(5) = 31,
6 = prime(1)^1*prime(2)^1 => a(1) = 1121,
7 = prime(3)^1 => a(7) = 41,
8 = prime(1)^3 => a(8) = 13, and so on.
MAPLE
a:= n-> `if`(n=1, 0, parse(cat(seq([numtheory[pi]
(i[1]), i[2]][], i=sort(ifactors(n)[2]))))):
seq(a(n), n=1..60); # Alois P. Heinz, Mar 16 2018
MATHEMATICA
Array[FromDigits@ Flatten@ Map[{PrimePi@ #1, #2} & @@ # &, FactorInteger@ #] &, 51] (* Michael De Vlieger, Mar 16 2018 *)
PROG
(PARI) A299400(n)=if(n=factor(n), eval(concat(apply(f->Str(primepi(f[1]), f[2]), Col(n)~))))
CROSSREFS
KEYWORD
nonn,base
AUTHOR
M. F. Hasler, Mar 08 2018
STATUS
approved