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.

A065774 Numbers k such that the product of the digits of k is equal to the sum of the prime factors of k, counted with multiplicity.

This page as a plain text file.
%I A065774 #22 Apr 13 2024 14:57:30
%S A065774 2,3,4,5,7,18,25,154,329,418,442,532,1519,1826,2354,2472,2781,3383,
%T A065774 4343,4712,5282,5561,6171,6623,7922,9331,9911,11248,12328,12773,13125,
%U A065774 14193,15194,16512,17267,19852,21479,24516,26522,31816,32661,36512
%N A065774 Numbers k such that the product of the digits of k is equal to the sum of the prime factors of k, counted with multiplicity.
%H A065774 Giovanni Resta, <a href="/A065774/b065774.txt">Table of n, a(n) for n = 1..10000</a> (first 400 terms from Harry J. Smith)
%t A065774 Select[Range[2, 10^5], Times @@ IntegerDigits @ # == Plus @@ Times @@@ FactorInteger @ # &] (* _Giovanni Resta_, Apr 23 2017 *)
%o A065774 (PARI) ProdD(x)= { local(p=1); while (x>9 && p>0, p*=(x-10*(x\10)); x\=10); return(p*x) }
%o A065774 sopfr(n)= { local(f,s=0); f=factor(n); for(i=1, matsize(f)[1], s+=f[i, 1]*f[i, 2]); return(s) }
%o A065774 { n=0; for (m=1, 10^9, if (ProdD(m) == sopfr(m), write("b065774.txt", n++, " ", m); if (n==400, return)) ) } \\ _Harry J. Smith_, Oct 30 2009
%Y A065774 Cf. A067173, A001414, A007954.
%K A065774 base,easy,nonn
%O A065774 1,1
%A A065774 _Jason Earls_, Dec 04 2001
%E A065774 Offset changed from 0 to 1 by _Harry J. Smith_, Oct 30 2009