%I #12 Sep 18 2015 23:35:08
%S 1,1,3,1,7,1,9,7,3,7,9,11,3,11,17,3,13,1,27,1,37,3,17,9,19,1,39,7,19,
%T 3,31,19,7,27,7,33,7,39,11,23,3,47,9,29,3,49,1,49,9,37,9,41,9,47,21,1,
%U 51,13,19,9,53,23,9,67,3,53,33,13,21,11,29,17,21,13,27,11,51,19,13,61
%N a(1) = 1, a(n) = smallest number such that the concatenation of two successive terms gives a prime which has not occurred earlier.
%C Conjecture: Every odd prime occurs in this sequence infinitely many times.
%H Alois P. Heinz, <a href="/A077202/b077202.txt">Table of n, a(n) for n = 1..10000</a>
%p b:= proc() true end:
%p a:= proc(n) option remember; local h, k, p;
%p if n=1 then 1
%p else h:= a(n-1);
%p for k do p:=parse(cat(h,k));
%p if b(p) and isprime(p) then break fi
%p od; b(p):= false; k
%p fi
%p end:
%p seq(a(n), n=1..100); # _Alois P. Heinz_, Sep 18 2015
%o (PARI) A077202(nmax)= { local(a,tst,hadp,hSet) ; a=[1] ; hadp=[1] ; for(n=2,nmax, for(new=1,10000, tst=Str(eval(a[n-1]) eval(new)) ; tst=eval(tst) ; if(isprime(tst), hSet=Set(hadp) ; if( setsearch(hSet,tst)==0, hadp=concat(hadp,tst) ; a=concat(a,new) ; break ; ) ; ) ; ) ; ) ; return(a) ; } { print(A077202(80)) ; } - _R. J. Mathar_, May 19 2006
%Y Cf. A082238.
%K base,nonn
%O 1,3
%A _Amarnath Murthy_, Nov 02 2002
%E Corrected and extended by _R. J. Mathar_, May 19 2006
%E Offset corrected by _Alois P. Heinz_, Sep 18 2015
|