cp's OEIS Frontend

A247714 Position of A036561(n) in sequence A003586.

%I A247714 #12 Sep 23 2014 05:07:37
%S A247714 1,2,3,4,5,7,6,8,10,12,9,11,14,16,19,13,15,18,21,24,27,17,20,23,26,30,
%T A247714 33,37,22,25,29,32,36,40,44,49,28,31,35,39,43,47,52,57,62,34,38,42,46,
%U A247714 51,55,60,66,71,77,41,45,50,54,59,64,69,75,81,87,93,48
%N A247714 Position of A036561(n) in sequence A003586.
%C A247714 Motivated by _L. Edson Jeffery_ comment in A036561 that says it is a permutation of A003586.
%H A247714 Reinhard Zumkeller, <a href="/A247714/b247714.txt">Table of n, a(n) for n = 0..10000</a>
%p A247714 N0:= 10: # to get the first (N0+1)*(N0+2)/2 terms
%p A247714 V:= 3^N0:
%p A247714 S:= {seq(seq(2^i*3^j, i=0..ilog2(V/3^j)),j=0..N0)}:
%p A247714 # in Maple 11 or earlier, uncomment the next line and comment out the previous one
%p A247714 # S:= sort([seq(seq(2^i*3^j, i=0..ilog2(V/3^j)),j=0..N0)]):
%p A247714 for k from 1 to nops(S) do
%p A247714   r:= S[k];
%p A247714   jr:= padic[ordp](r,3);
%p A247714   ir:= jr + padic[ordp](r,2);
%p A247714   A[1+jr+ir*(ir+1)/2] := k;
%p A247714 od:
%p A247714 seq(A[k],k=1..(N0+1)*(N0+2)/2); # _Robert Israel_, Sep 22 2014
%o A247714 (PARI) lista(nn) = {w = readvec("b036561.txt"); v = readvec("b003586.txt"); for (i=1, nn, print1(setsearch(v, w[i], 0), ", "););}
%o A247714 (Haskell)
%o A247714 import Data.List (elemIndex); import Data.Maybe (fromJust)
%o A247714 a247714 = (+ 1) . fromJust .
%o A247714                   (`elemIndex` a003586_list) . (a036561_list !!)
%o A247714 -- _Reinhard Zumkeller_, Sep 23 2014
%Y A247714 Cf. A003586, A036561.
%K A247714 nonn
%O A247714 0,2
%A A247714 _Michel Marcus_, Sep 22 2014