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.
%I A172143 #17 Sep 01 2024 09:37:19 %S A172143 0,1,3,5,13,21,28,53,85,113,213,227,341,453,853,909,1365,1813,1820, %T A172143 3413,3637,5461,7253,7281,13653,14549,14563,21845,29013,29125,54613, %U A172143 58197,58253,87381,116053,116501,116508,218453,232789,233013,349525,464213,466005,466033 %N A172143 a(n) = (A172126(n) - 1)/3. %C A172143 Conjecture: sequence consists of an infinite number of subsequences S(m,0) = A172241(n) = (1/18)[8^n-(-1)^n-9], m>0, S(m,n+1) = 4*S(m,n)+1. The first subsequences %C A172143 S(1,n) = A002450(n) = (4^n-1)/3 = 0, 1, 5, 21, 85, ..., %C A172143 S(2,n) = A072197(n) = (10*4^n-1)/3 = 3, 13, 53, 213, ..., %C A172143 S(3,n) = (85*4^n-1)/3 = 28, 113, 453, ..., %C A172143 S(4,n) = (682*4^n-1)/3 = 227, 909, 3637, ..., and generally, %C A172143 S(m,n) = [(3*A172241(m) + 1) * 4^n - 1]/3. %H A172143 Amiram Eldar, <a href="/A172143/b172143.txt">Table of n, a(n) for n = 1..10000</a> %t A172143 seq[max_] := Module[{kmax = Floor[Log[4, 3*max+1]], s = {}, s1, odd},Do[odd = (4^k-1)/3; s1 = 2^Range[0, Floor[Log2[max/odd]]] * odd; s = Join[s, s1], {k, 1, kmax}]; Select[(Union[s] - 1)/3, IntegerQ]]; seq[10^7] (* _Amiram Eldar_, Sep 01 2024 *) %o A172143 (PARI) for(n=1, 300000, o=3*n/2^valuation(n, 2)+1; b=ispower(o); if(b&&b%2==0&&round(sqrtn(o, b/2))==4&&(n-1)%3==0, print1((n-1)/3, ", "))) %Y A172143 Cf. A002450, A072197, A172126, A172241. %K A172143 nonn %O A172143 1,3 %A A172143 _Ralf Stephan_, Nov 19 2010 %E A172143 More terms from _Amiram Eldar_, Sep 01 2024