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.

Showing 1-4 of 4 results.

A320486 Keep just the digits of n that appear exactly once; write 0 if all digits disappear.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 0, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 0, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 0, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 0, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 0, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 0, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 0, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 0, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 0, 1, 0, 102, 103, 104, 105, 106, 107, 108, 109, 0, 0, 2, 3, 4, 5, 6, 7, 8, 9, 120, 2, 1, 123, 124, 125, 126, 127, 128, 129, 130, 3
Offset: 0

Views

Author

N. J. A. Sloane, Oct 24 2018

Keywords

Comments

Digits that appear more than once in n are erased. Leading zeros are erased unless the result is 0. If all digits are erased, we write 0 for the result (A320485 is another version, which uses -1 for the empty string).
More than the usual number of terms are shown in order to reach some interesting examples.
a(n) = 0 mostly. - David A. Corneth, Oct 24 2018
The number of d-digit numbers n for which a(n) > 0 is at most d*9^d, so in this sense most a(n) are 0. - Robert Israel, Oct 24 2018
The set of numbers with the property that their digits appear at least twice is of asymptotic density 1 (and the set of numbers that have a digit that occurs only once is of density 0), so in that sense it is rather exceptional for large n to have a(n) > 0. - M. F. Hasler, Oct 24 2018

Examples

			1231 becomes 23, 1123 becomes 23, 11231 becomes 23, and 11023 becomes 23 (as we don't accept leading zeros). Note that 112323 disappears immediately and we get 0.
101, 110, 11000, 10001 all become 0.
		

References

  • Eric Angelini, Posting to Sequence Fans Mailing List, Oct 24 2018

Crossrefs

See A320485 for a different version.

Programs

  • Maple
    f:= proc(n) local F,S;
      F:= convert(n,base,10);
      S:= select(t -> numboccur(t,F)>1, [$0..9]);
      if S = {} then return n fi;
      F:= subs(seq(s=NULL,s=S),F);
      add(F[i]*10^(i-1),i=1..nops(F))
    end proc:
    map(f, [$0..200]); # Robert Israel, Oct 24 2018
  • Mathematica
    Table[If[(c=Select[b=IntegerDigits[n],Count[b,#]==1&])=={},0,FromDigits@c],{n,0,131}] (* Giorgos Kalogeropoulos, May 09 2021 *)
    d1[n_]:=Module[{idn=IntegerDigits[n]},FromDigits[If[DigitCount[n,10,#]>1,Nothing,#]&/@ idn]]; Array[d1,150,0] (* Harvey P. Dale, Jun 23 2023 *)
  • PARI
    a(n) = {my(d=digits(n), v = vector(10), res = 0); for(i=1,#d, v[d[i]+1]++); for(i=1,#d,if(v[d[i]+1]==1, res=10*res+d[i]));res}
    
  • PARI
    A320486(n,D=digits(n))=fromdigits(select(d->#select(t->t==d,D)<2,D)) \\ M. F. Hasler, Oct 24 2018
    
  • Python
    def A320486(n):
        return int('0'+''.join(d if str(n).count(d) == 1 else '' for d in str(n))) # Chai Wah Wu, Nov 19 2018

A320485 Keep just the digits of n that appear exactly once; write -1 if all digits disappear.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, -1, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, -1, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, -1, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, -1, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, -1, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, -1, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, -1, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, -1, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, -1, 1, 0, 102, 103, 104, 105, 106, 107, 108, 109, 0, -1, 2, 3, 4, 5, 6, 7, 8, 9, 120, 2, 1, 123, 124, 125, 126, 127, 128, 129, 130, 3
Offset: 0

Views

Author

N. J. A. Sloane, Oct 24 2018

Keywords

Comments

Digits that appear more than once in n are erased. Leading zeros are erased unless the result is 0. If all digits are erased, we write -1 for the result.
The map n -> a(n) was invented by Eric Angelini and described in a posting to the Sequence Fans Mailing List on Oct 24 2018.
More than the usual number of terms are shown in order to reach some interesting examples.
a(n) = -1 mostly. - David A. Corneth, Oct 24 2018

Examples

			1231 becomes 23, 1123 becomes 23, 11231 becomes 23, and 11023 becomes 23 (as we don't accept leading zeros). Note that 112323 disappears immediately and we get -1.
101 and 110 become 0 while 11000 and 10001 become -1.
		

References

  • Eric Angelini, Posting to Sequence Fans Mailing List, Oct 24 2018

Crossrefs

See A320486 for another version.

Programs

  • Maple
    f:= proc(n) local F,S;
      F:= convert(n,base,10);
      S:= select(t -> numboccur(t,F)>1, [$0..9]);
      if S = {} then return n fi;
      F:= subs(seq(s=NULL,s=S),F);
      if F = [] then -1
      else add(F[i]*10^(i-1),i=1..nops(F))
      fi
    end proc:
    map(f, [$0..200]); # Robert Israel, Oct 24 2018
  • Mathematica
    Array[If[# == {}, -1, FromDigits@ #] &@ Map[If[#[[-1]] > 1, -1, #[[1]] ] /. -1 -> Nothing &, Tally@ IntegerDigits[#]] &, 131] (* Michael De Vlieger, Oct 24 2018 *)
  • PARI
    a(n) = {my(d=digits(n), v = vector(10), res = 0, t = 0); for(i=1, #d, v[d[i]+1]++); for(i=1, #d, if(v[d[i]+1]==1, t = 1; res=10 * res + d[i])); res - !t + !n} \\ David A. Corneth, Oct 24 2018
    
  • Python
    def A320485(n):
        return (lambda x: int(x) if x != '' else -1)(''.join(d if str(n).count(d) == 1 else '' for d in str(n))) # Chai Wah Wu, Nov 19 2018

Formula

From Rémy Sigrist, Oct 24 2018: (Start)
a(n) = n iff n belong to A010784.
a(n) <= 9876543210 with equality iff n = 9876543210.
(End)
If n > 9876543210, then a(n) < n. If a(n) < n, then a(n) <= 99n/1000. - Chai Wah Wu, Oct 24 2018

A321011 Trajectory of 86 under repeated application of the map k -> A320486(k^2).

Original entry on oeis.org

86, 7396, 547816, 12985, 805, 648025, 1325, 1762, 3106, 94723, 86, 7396, 547816, 12985, 805, 648025, 1325, 1762, 3106, 94723, 86, 7396, 547816, 12985, 805, 648025, 1325, 1762, 3106, 94723, 86, 7396, 547816, 12985, 805, 648025, 1325, 1762, 3106, 94723
Offset: 1

Views

Author

N. J. A. Sloane, Nov 04 2018

Keywords

Comments

k -> A320486(k) is Eric Angelini's remove-repeated-digits map.
Lars Blomberg has discovered that if we start with any positive integer and repeatedly apply the map k -> A320486(k^2) then we will eventually either:
- reach 0,
- reach one of the four fixed points 1, 1465, 4376, 89476 (see A321010)
- reach the period-10 cycle shown in A321011, or
- reach the period-9 cycle shown in A321012.
Since there are only finitely many possible starting values with all digits distinct, it should not be difficult to check that this is true (and indeed, Lars Blomberg may by now have completed the proof).

Examples

			The cycle of length 10 is (86, 7396, 547816, 12985, 805, 648025, 1325, 1762, 3106, 94723).
		

References

  • Eric Angelini, Postings to Sequence Fans Mailing List, Oct 24 2018 and Oct 26 2018.

Crossrefs

Programs

  • Mathematica
    LinearRecurrence[{0,0,0,0,0,0,0,0,0,1},{86,7396,547816,12985,805,648025,1325,1762,3106,94723},40] (* or *) PadRight[ {},40,{86,7396,547816,12985,805,648025,1325,1762,3106,94723}] (* Harvey P. Dale, Nov 05 2020 *)
  • PARI
    Vec(x*(86 + 7396*x + 547816*x^2 + 12985*x^3 + 805*x^4 + 648025*x^5 + 1325*x^6 + 1762*x^7 + 3106*x^8 + 94723*x^9) / ((1 - x)*(1 + x)*(1 - x + x^2 - x^3 + x^4)*(1 + x + x^2 + x^3 + x^4)) + O(x^40)) \\ Colin Barker, Nov 04 2018

Formula

From Colin Barker, Nov 04 2018: (Start)
G.f.: x*(86 + 7396*x + 547816*x^2 + 12985*x^3 + 805*x^4 + 648025*x^5 + 1325*x^6 + 1762*x^7 + 3106*x^8 + 94723*x^9) / ((1 - x)*(1 + x)*(1 - x + x^2 - x^3 + x^4)*(1 + x + x^2 + x^3 + x^4)).
a(n) = a(n-10) for n>10.
(End)

A321010 Numbers k such that f(k^2) = k, where f is Eric Angelini's remove-repeated-digits map x->A320486(x).

Original entry on oeis.org

0, 1, 1465, 4376, 89476
Offset: 1

Views

Author

N. J. A. Sloane, Nov 03 2018

Keywords

Comments

Lars Blomberg has discovered that if we start with any positive integer and repeatedly apply the map m -> A320486(m^2) then we will eventually either:
- reach 0,
- reach one of the four fixed points 1, 1465, 4376, 89476 (this sequence),
- reach the period-10 cycle shown in A321011, or
- reach the period-9 cycle shown in A321012.
From Lars Blomberg, Nov 17 2018: (Start)
Verified by testing all possible 8877690 start values that these are the only fixed points and cycles.
Detailed counts are:
- 561354 reach 0,
- 963738 reach one of the four fixed points 1, 1465, 4376, 89476 (counts 946109, 434, 17065, 130),
- 7271337 reach the period-10 cycle, and
- 81261 reach the period-9 cycle. (End)

References

  • Eric Angelini, Postings to Sequence Fans Mailing List, Oct 24 2018 and Oct 26 2018.

Crossrefs

Showing 1-4 of 4 results.