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.

A071970 List the positive rationals in the order in which they are produced by the Stern sequence A002487 and apply the Sagher map to turn them into integers.

This page as a plain text file.
%I A071970 #12 Oct 03 2017 03:32:50
%S A071970 1,2,4,3,18,12,9,8,48,45,50,20,75,72,16,5,200,112,147,288,320,175,98,
%T A071970 28,245,800,192,63,392,80,25,6,180,675,648,176,847,490,300,99,3872,
%U A071970 832,845,600,1008,1323,162,108,567,1176,720,325,5408,704,363,90,700,539
%N A071970 List the positive rationals in the order in which they are produced by the Stern sequence A002487 and apply the Sagher map to turn them into integers.
%C A071970 The Sagher map sends Product p_i^e_i / Product q_i^f_i (p_i and q_i being distinct primes) to Product p_i^(2e_i) * Product q_i^(2f_i-1). This is multiplicative.
%H A071970 Y. Sagher, <a href="http://www.jstor.org/stable/2324846">Counting the rationals</a>, Amer. Math. Monthly, 96 (1989), p. 823.
%e A071970 The first few rationals and their images are 1/1 -> 1, 1/2 -> 2, 2/1 -> 4, 1/3 -> 3, 3/2 -> 18, 2/3 -> 12, 3/1 -> 9, 1/4 -> 8, ...
%t A071970 nmax = 58; s[0] = 0; s[1] = 1; s[n_?EvenQ] := s[n/2]; s[n_] := s[(n-1)/2] + s[(n+1)/2]; v = Table[ FactorInteger /@ {s[n] , s[n+1]}, {n, 1, nmax}]; a[n_] := Times @@ (#[[1]]^(2*#[[2]])&) /@ v[[n, 1]]*Times @@ (#[[1]]^(2*#[[2]]-1)&) /@ v[[n, 2]]; Table[a[n], {n, 1, nmax}] (* _Jean-François Alcover_, Nov 25 2011, after Pari *)
%o A071970 (PARI) s(n)=if(n<2,n>0,if(n%2,s((n+1)/2)+s((n-1)/2),s(n/2))) /* A002487(n) */
%o A071970 (PARI) a(n)=local(v); if(n,v=factor(s(n)/s(n+1))~; prod(k=1,length(v),v[1,k]^if(v[2,k]<0,-1-2*v[2,k],2*v[2,k])),0)
%Y A071970 Cf. A002487, A060837.
%K A071970 nonn,nice,easy
%O A071970 1,2
%A A071970 _N. J. A. Sloane_, Jun 19 2002
%E A071970 More terms from _Michael Somos_, Jul 19 2002