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.

A133965 Home primes whose homeliness is greater than 5.

This page as a plain text file.
%I A133965 #5 Apr 28 2015 15:00:40
%S A133965 3389,3761,13367,23251,23761,23929,31193,31397,33797,36389,37199,
%T A133965 37463,37547,37607,37643,37717,37853,37907,37951,37967,41863,41887,
%U A133965 41941,73331,74759,74771,77269,77711,77773,77797,83383,112909,113779,114773
%N A133965 Home primes whose homeliness is greater than 5.
%C A133965 Number of terms < 10^n: 0, 0, 0, 2, 31, 341, 2695, ....
%e A133965 (118, 259, 737, 801, 1167, 3389) -> 3389, etc.
%t A133965 lst = {}; f[n_] := FromDigits@ Flatten[ IntegerDigits@ Table[ #[[1]], {#[[2]]}] & /@ FactorInteger@n, 2]; h[n_] := NestWhileList[f@# &, n, !PrimeQ@# &, 1, 28]; Do[p = h[n][[ -1]]; If[ PrimeQ@p && p < 10^7 && p != n, Print[{n, p}]; AppendTo[lst, p]], {n, 2, 1000}];
%t A133965 d = 6 - 2; lsu = {}; Do[ If[ lst[[n]] == lst[[n + d]], AppendTo[ lsu, lst[[n]] ]], {n, 188004 - d - 1}]; Take[ Union@ lsu, 35]
%Y A133965 Cf. A037274, A118756, A133957, A133959, A133961, A133963, A133967, A133969, A133971, A133973, A133975, A133977, A133979.
%K A133965 nonn,base
%O A133965 1,1
%A A133965 _Robert G. Wilson v_, Sep 30 2007