%I #7 Aug 18 2014 00:52:06
%S 1,2,5,3,4,9,11,7,6,18,10,59,20,25,16,8,50,15,32,31,12,13,38,21,41,
%T 125,85,43,17,45,52,35,22,19,103,105,33,24,14,190,68,27,66,28,161,29,
%U 80,26,54,46,177,84,258,34,180,64,90,70,507,37,196,96,39,110,430,92,78,75,600,48,40,82,213,218,71,23,87,72,51,132,30
%N Permutation of natural numbers induced when A091205 is restricted to {1} and binary codes for polynomials reducible over GF(2): a(1) = 1, a(n) = A062298(A091205(A091242(n-1))).
%H Antti Karttunen, <a href="/A245813/b245813.txt">Table of n, a(n) for n = 1..10001</a>
%H <a href="/index/Per#IntegerPermutation">Index entries for sequences that are permutations of the natural numbers</a>
%F a(1) = 1, and for n > 1, a(n) = A062298(A091205(A091242(n-1))).
%F As a composition of related permutations:
%F a(n) = A245815(A245820(n)).
%o (PARI)
%o allocatemem(234567890);
%o v091226 = vector(2^22);
%o v091242 = vector(2^22);
%o isA014580(n)=polisirreducible(Pol(binary(n))*Mod(1, 2)); \\ This function from _Charles R Greathouse IV_
%o i=0; j=0; n=2; while((n < 2^22), if(isA014580(n), i++; v091226[n] = v091226[n-1]+1, j++; v091242[j] = n; v091226[n] = v091226[n-1]); n++);
%o A062298(n) = n-primepi(n);
%o A091226(n) = v091226[n];
%o A091242(n) = v091242[n];
%o A091205(n) = if(n<=1, n, if(isA014580(n), prime(A091205(A091226(n))), {my(irfs, t); irfs=subst(lift(factor(Mod(1, 2)*Pol(binary(n)))), x, 2); irfs[,1]=apply(t->A091205(t), irfs[,1]); factorback(irfs)}));
%o A245813(n) = if(n<=1, n, A062298(A091205(A091242(n-1))));
%o for(n=1, 10001, write("b245813.txt", n, " ", A245813(n)));
%o (Scheme) (define (A245813 n) (if (<= n 1) n (A062298 (A091205 (A091242 (- n 1))))))
%Y Inverse: A245814.
%Y Cf. A062298, A091242.
%Y Related permutations: A091205, A245815, A245820.
%K nonn
%O 1,2
%A _Antti Karttunen_, Aug 16 2014