OFFSET
1,2
LINKS
FORMULA
PROG
(PARI)
allocatemem(234567890);
vecsize = (2^24)-4;
uplim = 2^25;
v091226 = vector(vecsize); A091226 = n -> v091226[n];
v206074 = vector(vecsize); A206074 = n -> v206074[n];
v205783 = vector(vecsize); A205783 = n -> v205783[n];
isA014580(n) = polisirreducible(Pol(binary(n))*Mod(1, 2)); \\ This function from Charles R Greathouse IV
isA206074(n) = polisirreducible(Pol(binary(n)));
v091226[1] = 0; v205783[1] = 1; i=0; j=1; n=2; while((n < uplim), if(!(n%65536), print1(n, ", ")); if(isA206074(n), i++; v206074[i] = n, j++; v205783[j] = n); if(isA014580(n), v091226[n] = v091226[n-1]+1, v091226[n] = v091226[n-1]); n++); print(n);
A260425(n) = if(1==n, 1, if(isA014580(n), A206074(A260425(A091226(n))), A205783(1+A260425(A091245(n)))));
for(n=1, 4244, write("b260425.txt", n, " ", A260425(n)));
CROSSREFS
Inverse: A260426.
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jul 26 2015
STATUS
approved