login
(0, 1, 2, 3, 2^2, 5, 2*3, 7, 2^3, 3^2, 2*5, 11, 2^2*3, 13, ...) becomes (0^(1+2), 3^(2+2), 5^(2+3), 7^(2+3), 3^(2+2), 5^(11+2), 2^(3+13), ...).
1

%I #7 Oct 05 2015 02:59:24

%S 0,81,3125,16807,81,1220703125,65536,1024,15625,1419857,2097152,256,

%T 96889010407,6436343,2187,65536,81,157775382034845806615042743,

%U 150094635296999121,256,61159090448414546291,1953125,32

%N (0, 1, 2, 3, 2^2, 5, 2*3, 7, 2^3, 3^2, 2*5, 11, 2^2*3, 13, ...) becomes (0^(1+2), 3^(2+2), 5^(2+3), 7^(2+3), 3^(2+2), 5^(11+2), 2^(3+13), ...).

%e 0^(1 + 2) = 0^3 = 0 = a(1).

%e 3^(2 + 2) = 3^4 = 81 = a(2).

%e 5^(2 + 3) = 5^5 = 3125 = a(3).

%e 7^(2 + 3) = 7^5 = 16807 = a(4).

%e 3^(2 + 2) = 3^4 = 81 = a(5).

%e 5^(11 + 2) = 5^13 = 1220703125 = a(6).

%e 2^(3 + 13) = 2^16 = 65536 = a(7).

%e 2^(7 + 3) = 2^10 = 1024 = a(8), etc.

%p pflat2 := proc(nmax) local a, ifs, n, p, c ; a := [0,1] ; for n from 2 to nmax do ifs := ifactors(n)[2] ; for p in ifs do a := [op(a),op(1,p)] ; if op(2,p) > 1 then a := [op(a),op(2,p)] ; fi; od: od: a ; end: pL := pflat2(120) : for n from 1 to nops(pL)-4 by 3 do printf("%d,", op(n, pL)^(op(n+1, pL)+op(n+2,pL)) ) ; od: # _R. J. Mathar_, Nov 06 2008

%Y Cf. A141269, A141270.

%K nonn

%O 1,2

%A _Juri-Stepan Gerasimov_, Nov 01 2008

%E Extended by _R. J. Mathar_, Nov 06 2008