The OEIS is supported by the many generous donors to the OEIS Foundation. Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A067206 Numbers n such that the digits of n end in phi(n). 4
 1, 1320, 1640, 1768, 1996, 2640, 3960, 13200, 16400, 19984, 19996, 26400, 39600, 132000, 164000, 199996, 264000, 396000, 1320000, 1640000, 1999936, 2640000, 3960000, 13200000, 16400000, 16666240, 17999488, 18515584, 19999984, 19999996 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,2 COMMENTS Comments from Farideh Firoozbakht, Dec 30 2006: (Start) "(1). If n is in the sequence and 10 divides n then for each natural number k, n*10^k is in the sequence. So since 1320, 1640, 2640, 3960 & 16666240 are in the sequence, for each natural number k, 132*10^k, 164*10^k, 264*10^k, 396*10^k & 1666624*10^k are in the sequence. Hence the sequence is infinite. "(2). If 5*10^k-1 is prime then 4*(5*10^k-1) is in the sequence. So 4*A093945 is a subsequence of this sequence. "(3). If p=125*10^k-1 is prime then 16*p is in the sequence. For k = 1, 4, 5, 8, 13, 19, 25, 26, 76, 88, 167, 290, 389, ... p is prime. "(4). If p=3125*10^k-1 is prime then 64*p is in the sequence. For k = 1, 3, 9, 33, 121, 223, 357, 363, 447, ... p is prime." (End) REFERENCES Pickover, C. "Wonders of Numbers". Oxford Univ. Press, 2001. LINKS Giovanni Resta, Table of n, a(n) for n = 1..66 (terms < 10^12, first 33 terms from Farideh Firoozbakht) C. A. Pickover, "Wonders of Numbers, Adventures in Mathematics, Mind and Meaning," Zentralblatt review EXAMPLE The digits of 1768 end in phi(1768) = 768, so 1768 is a term of the sequence. MAPLE with(numtheory); P:=proc(q) local a, b, n; for n from 1 to q do a:=phi(n); b:=0; while a>0 do b:=b+1; a:=trunc(a/10); od; if phi(n)=(n mod 10^b) then print(n); fi; od; end: P(10^9); # Paolo P. Lava, Mar 13 2014 MATHEMATICA (*returns true if a ends in b, false o.w.*) f[a_, b_] := Module[{c, d, e, g, h, i, r}, r = False; c = ToString[a]; d = ToString[b]; e = StringLength[c]; g = StringPosition[c, d]; h = Length[g]; If[h > 0, i = g[[h]]; If[i[] == e, r = True]]; r]; Select[Range[10^5], f[ #, EulerPhi[ # ]] &] CROSSREFS Cf. A093945, A091439. Cf. A066663. - R. J. Mathar, Sep 30 2008 Sequence in context: A255799 A023318 A066663 * A248857 A260839 A157266 Adjacent sequences:  A067203 A067204 A067205 * A067207 A067208 A067209 KEYWORD base,nonn AUTHOR Joseph L. Pe, Feb 19 2002 EXTENSIONS More terms from Farideh Firoozbakht, Dec 30 2006 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified January 18 10:47 EST 2022. Contains 350454 sequences. (Running on oeis4.)