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.

A060291 Number of twin Harshads, including overlaps, whose sum is prime and where the 2nd Harshad is <= 10^n.

This page as a plain text file.
%I A060291 #13 Jan 19 2021 03:58:34
%S A060291 7,8,21,60,278,1520,9583,61835,419705,3023007,22597508
%N A060291 Number of twin Harshads, including overlaps, whose sum is prime and where the 2nd Harshad is <= 10^n.
%F A060291 Generate the twin Harshads whose sum is prime. Count how many there are where the 2nd Harshad in the pair is <= to a consecutive power of 10.
%e A060291 a(1)=7 because there are 7 pairs of Harshads whose sum is prime and the 2nd Harshad in the pair is <= 10; these are 1+2=3, 2+3=5, 3+4=7, 5+6=11, 6+7=13, 8+9=17, 9+10=19. (Another sequence does not permit overlapping Harshad numbers.)
%t A060291 harshadQ[n_] := Divisible[n, Plus @@ IntegerDigits[n]]; c = 0; p = 10; s = {}; n = 0; k = 2; q1 = True; While[n < 6, q2 = harshadQ[k]; If[q1 && q2 && PrimeQ[2*k-1], c++; If[k > p, n++; AppendTo[s, c-1]; p *= 10]]; q1 = q2; k++]; s (* _Amiram Eldar_, Jan 19 2021 *)
%Y A060291 Cf. A060288, A060289, A060290, A005349, A060159.
%K A060291 nonn,base,more
%O A060291 1,1
%A A060291 _Enoch Haga_, Mar 24 2001
%E A060291 a(8)-a(11) from _Amiram Eldar_, Jan 19 2021