A366260 - OEIS

%I #10 Oct 06 2023 08:34:48

%S 1,2,4,3,9,5,6,8,16,7,10,12,15,27,25,18,54,11,14,20,21,45,35,30,24,32,

%T 49,50,36,75,81,125,625,13,22,28,33,63,55,42,40,48,77,70,60,105,135,

%U 175,90,162,121,98,100,147,225,245,150,72,64,343,250,108,375,243,729,17,26,44,39,99,65,66,56,80,91,110

%N Doudna sequence permuted by May code: a(n) = A005940(1+A303767(n)).

%H Antti Karttunen, <a href="/A366260/b366260.txt">Table of n, a(n) for n = 0..16384</a>

%H <a href="/index/Bi#binary">Index entries for sequences related to binary expansion of n</a>

%H <a href="/index/Per#IntegerPermutation">Index entries for sequences that are permutations of the natural numbers</a>

%o (PARI)

%o A005940(n) = { my(p=2, t=1); n--; until(!n\=2, if((n%2), (t*=p), p=nextprime(p+1))); t };

%o A209229(n) = (n && !bitand(n,n-1));

%o A053644(n) = { my(k=1); while(k<=n, k<<=1); (k>>1); }; \\ From A053644

%o A303767(n) = if(!n,n,if(A209229(n),n+A303767(n-1),A053644(n)+A303767(n-A053644(n)-1)));

%o A366260(n) = A005940(1+A303767(n));

%Y Cf. A005940, A303767, A366261, A366262.

%Y Cf. also A243353, A366263.

%K nonn

%O 0,2

%A _Antti Karttunen_, Oct 05 2023