cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A322991 a(n) = A289272(2*A289271(n)).

This page as a plain text file.
%I A322991 #7 Jan 02 2019 11:55:10
%S A322991 1,3,4,5,7,12,8,9,11,15,13,20,16,21,28,17,19,24,23,35,36,33,25,44,27,
%T A322991 39,29,40,31,60,32,37,52,48,56,45,41,51,68,63,43,84,47,55,77,57,49,76,
%U A322991 53,69,92,65,59,75,91,72,100,87,61,140,64,93,88,67,112,132,71,80,108,105,73,99,79,96,116,85,104,156,81,119,83,111,89
%N A322991 a(n) = A289272(2*A289271(n)).
%H A322991 Antti Karttunen, <a href="/A322991/b322991.txt">Table of n, a(n) for n = 1..16384</a>
%F A322991 a(n) = A289272(2*A289271(n)).
%o A322991 (PARI)
%o A322991 A289271(n) = { my(v=0,i=0,x=1); for(d=2,oo,if(n==1, return(v)); if(1==gcd(x,d)&&1==omega(d), if(!(n%d)&&1==gcd(d,n/d), v += 2^i; n /= d; x *= d); i++)); }; \\ After _Rémy Sigrist_'s program for A289271.
%o A322991 A289272(n) = { my(m=1, pp=1); while(n>0, pp++; while(!isprimepower(pp)||(gcd(pp,m)>1), pp++); if(n%2, m *= pp); n >>=1); (m); };
%o A322991 A322991(n) = A289272(2*A289271(n));
%Y A322991 Permutation of A042965 (apart from zero).
%Y A322991 Cf. A289271, A289272, A322990, A322992.
%Y A322991 Cf. also A003963, A300841.
%K A322991 nonn
%O A322991 1,2
%A A322991 _Antti Karttunen_, Jan 01 2019