Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A108505
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A108505 Number of palindromic semiprimes less than 10^n. +0
2
0, 3, 7, 36, 50, 269, 367, 2181, 2816, 18391, 23617, 160773, 203733, 1429749, 1788486 (list; graph; listen)
OFFSET

0,2

COMMENT

Does the limit n-> inf. a(n+2)/a(n) =~ 8*Pi^2/9?

MATHEMATICA

NextPalindrome[n_] := Block[ {l = Floor[ Log[10, n] + 1], idn = IntegerDigits[n]}, If[ Union[ idn] == {9}, Return[n + 2], If[l < 2, Return[n + 1], If[ FromDigits[ Reverse[ Take[ idn, Ceiling[l/2]]]] > FromDigits[ Take[ idn, -Ceiling[l/2]]], FromDigits[ Join[ Take[ idn, Ceiling[l/2]], Reverse[ Take[ idn, Floor[l/2]]] ]], idfhn = FromDigits[ Take[ idn, Ceiling[l/2]]] + 1; idp = FromDigits[ Join[ IntegerDigits[ idfhn], Drop[ Reverse[ IntegerDigits[ idfhn]], Mod[l, 2]]]] ]]]]; fQ[n_] := Plus @@ Last /@ FactorInteger[n] == 2; c = np = 0; Do[ While[np < 10^n, If[ fQ[np], c++ ]; np = NextPalindrome[np]]; Print[c], {n, 0, 12}]

CROSSREFS

Cf. A046328.

Sequence in context: A081555 A027624 A063042 this_sequence A047158 A102917 A156465

Adjacent sequences: A108502 A108503 A108504 this_sequence A108506 A108507 A108508

KEYWORD

base,nonn

AUTHOR

Robert G. Wilson v (rgwv(AT)rgwv.com), Jun 06 2005

page 1

Search completed in 0.002 seconds

Lookup | Welcome | Find friends | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
More pages | Superseeker | Maintained by N. J. A. Sloane (njas@research.att.com)

Last modified November 24 23:16 EST 2009. Contains 167481 sequences.


AT&T Labs Research