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.

A371631 Primes whose product of nonzero digits divided by the sum of its digits is also prime.

This page as a plain text file.
%I A371631 #43 Sep 25 2024 13:16:16
%S A371631 167,257,523,541,617,761,1447,1607,1861,2053,2251,2503,2521,2851,4051,
%T A371631 5023,5281,5821,6701,8161,8521,10067,10607,10861,11273,11471,12713,
%U A371631 13127,13217,13721,14407,16007,17123,17231,17321,18061,20507,20521,21247,21317,21713,22051
%N A371631 Primes whose product of nonzero digits divided by the sum of its digits is also prime.
%C A371631 No term N can have a "9" digit. [Proof: The sum of the digits of N is not a multiple of 3, but the numerator would be a multiple of 9, and so the number would be a multiple of 9, so not a prime.]
%H A371631 Harvey P. Dale, <a href="/A371631/b371631.txt">Table of n, a(n) for n = 1..1000</a>
%e A371631 167 (prime) is a term because 1*6*7/(1+6+7)=42/14=3 (prime).
%t A371631 pQ[n_] := Block[{idp = DeleteCases[IntegerDigits[n], 0]}, PrimeQ[Times @@ idp/Total@ idp]]; Cases[Prime@ Range@ PrimePi[10^5], _?pQ]
%t A371631 Select[Prime[Range[2500]],PrimeQ[Times@@(IntegerDigits[#]/.(0->1))/Total[ IntegerDigits[ #]]]&] (* _Harvey P. Dale_, Sep 24 2024 *)
%Y A371631 Subsequence of A038367.
%Y A371631 Equals prime terms of A138566.
%K A371631 nonn,base,easy
%O A371631 1,1
%A A371631 _Mikk Heidemaa_, May 24 2024