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.

A082059 Numbers k such that k, sigma(k) and phi(k) have the same decimal digits.

This page as a plain text file.
%I A082059 #13 Jun 18 2022 19:29:11
%S A082059 1,1021,1201,2011,7687,7867,8677,10211,10271,10321,12011,12041,12071,
%T A082059 12101,12301,12343,12401,12433,12601,14323,16787,17021,20011,20101,
%U A082059 20161,20411,20611,21001,21011,21031,21061,21101,21401,21433,21601,21701,22343,22433
%N A082059 Numbers k such that k, sigma(k) and phi(k) have the same decimal digits.
%H A082059 Donovan Johnson, <a href="/A082059/b082059.txt">Table of n, a(n) for n = 1..1000</a>
%e A082059 n = 89767, sigma(n) = 89768, phi(n) = 89766.
%t A082059 Select[Range[25000],Union[IntegerDigits[#]] == Union[ IntegerDigits[ DivisorSigma[ 1,#]]] == Union[IntegerDigits[EulerPhi[#]]]&] (* _Harvey P. Dale_, Jun 12 2014 *)
%o A082059 (PARI) c=0; for(n=1, 324773, m=n; nd=vector(10); for(j=1, 6, r=m%10; nd[r+1]=1; m=m\10; if(m==0, j=6)); s=sigma(n); m=s; sd=vector(10); for(j=1, 6, r=m%10; sd[r+1]=1; m=m\10; if(m==0, j=6)); for(j=1, 10, if(nd[j]<>sd[j], next(2))); p=eulerphi(n); m=p; pd=vector(10); for(j=1, 6, r=m%10; pd[r+1]=1; m=m\10; if(m==0, j=6)); for(j=1, 10, if(nd[j]<>pd[j], next(2))); c++; write("b082059.txt", c " " n)) \\ _Donovan Johnson_, Feb 28 2013
%Y A082059 Cf. A000203 (sigma), A000010 (phi).
%K A082059 nonn,base
%O A082059 1,2
%A A082059 _Labos Elemer_, Apr 04 2003