proposed
approved
proposed
approved
editing
proposed
This is GF(2)[X] analog of A005940, but note the indexing: here the domain starts from 0, even though although the range excludes zero.
Antti Karttunen, <a href="/A305417/b305417.txt">Table of n, a(n) for n = 0..16383</a>
This is GF(2)[X] analog of A005940, but note the indexing: here the domain starts from 0, while the range starts from 1.
This sequence can be represented as a binary tree. Each child to the left is obtained by applying A305421 to the parent, and each child to the right is obtained by doubling the parent:
1
|
...................2...................
3 4
7......../ \........6 5......../ \........8
/ \ / \ / \ / \
/ \ / \ / \ / \
/ \ / \ / \ / \
11 14 9 12 21 10 15 16
13 22 29 28 49 18 27 24 69 42 63 20 107 30 17 32
Sequence A305427 is obtained by scanning the same tree level by level from right to left.
nonn,tabf,changed
allocated for Antti KarttunenPermutation of natural numbers: a(0) = 1, a(2n) = A305421(a(n)), a(2n+1) = 2*a(n).
1, 2, 3, 4, 7, 6, 5, 8, 11, 14, 9, 12, 21, 10, 15, 16, 13, 22, 29, 28, 49, 18, 27, 24, 69, 42, 63, 20, 107, 30, 17, 32, 19, 26, 23, 44, 35, 58, 39, 56, 127, 98, 83, 36, 151, 54, 45, 48, 81, 138, 207, 84, 475, 126, 65, 40, 743, 214, 189, 60, 273, 34, 51, 64, 25, 38, 53, 52, 121, 46, 57, 88, 173, 70, 101, 116, 233, 78, 105, 112, 199, 254, 129
0,2
This is GF(2)[X] analog of A005940, but note the indexing: the domain starts from 0, while the range starts from 1.
<a href="/index/Ge#GF2X">Index entries for sequences operating on GF(2)[X]-polynomials</a>
<a href="/index/Per#IntegerPermutation">Index entries for sequences that are permutations of the natural numbers</a>
a(0) = 1, a(2n) = A305421(a(n)), a(2n+1) = 2*a(n).
(PARI)
A091225(n) = polisirreducible(Pol(binary(n))*Mod(1, 2));
A305420(n) = { my(k=1+n); while(!A091225(k), k++); (k); };
A305421(n) = { my(f = subst(lift(factor(Pol(binary(n))*Mod(1, 2))), x, 2)); for(i=1, #f~, f[i, 1] = Pol(binary(A305420(f[i, 1])))); fromdigits(Vec(factorback(f))%2, 2); };
A305417(n) = if(0==n, (1+n), if(!(n%2), A305421(A305417(n/2)), 2*(A305417((n-1)/2))));
allocated
nonn
Antti Karttunen, Jun 10 2018
approved
editing