%I
%S 2,3,9,101,682,2483,31774,532976
%N a(n) = smallest m such that b(m) is prime and b(m+2)/b(m) = prime(n), where b() = A098550().
%C So far a(n) is very roughly 2^(3n6). Of course it is not known at present if this sequence is infinite (see A098550). It is not even known if there are infinitely many prime terms in A098550. [It is now known that A098550 is a permutation of the natural numbers, so all of the primes appear. See Applegate et al.  _L. Edson Jeffery_, Dec 30 2014]
%C a(8) appears to deviate from the formula 2^(3n6) and is closer to 2^19 (than 2^18).  _Chai Wah Wu_, Dec 16 2014
%H David L. Applegate, Hans Havermann, Bob Selcoe, Vladimir Shevelev, N. J. A. Sloane, and Reinhard Zumkeller, <a href="http://arxiv.org/abs/1501.01669">The Yellowstone Permutation</a>, arXiv preprint arXiv:1501.01669 [math.NT], 2015 and <a href="https://cs.uwaterloo.ca/journals/JIS/VOL18/Sloane/sloane9.html">J. Int. Seq. 18 (2015) 15.6.7</a>.
%t f[lst_] := Block[{k = 4}, While[GCD[lst[[2]], k] == 1  GCD[lst[[1]], k] > 1  MemberQ[lst, k], k++]; Append[lst, k]]; A098850 = Nest[f, {1, 2, 3}, max = 3000]; b[n_] := A098850[[n]];
%t a[n_] := For[m = 1, m <= max, m++, If[PrimeQ[b[m]] && b[m+2]/b[m] == Prime[n], Return[m]]];
%t Do[Print[n, " ", a[n]], {n, 1, 6}] (* _JeanFrançois Alcover_, Dec 15 2018, after _Robert G. Wilson v_ in A098850 *)
%Y Cf. A098550, A251542.
%K nonn,more
%O 1,1
%A _David Applegate_ and _N. J. A. Sloane_, Dec 15 2014
%E a(8) from _Chai Wah Wu_, Dec 16 2014
