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-10 of 34 results. Next

A259986 This sequence and A259987 are base-6 analogs of A007185 and A016090, written in base 10.

Original entry on oeis.org

3, 9, 81, 81, 6561, 29889, 76545, 636417, 3995649, 24151041, 326481921, 689278977, 11573190657, 76876660737, 155240824833, 2035980763137, 4857090670593, 55637069004801, 157197025673217, 1375916505694209, 19656708706009089, 129341461907898369
Offset: 1

Views

Author

N. J. A. Sloane, Jul 13 2015

Keywords

Comments

See Schut (1991) for precise definition.
Ignoring repetitions, the subsequence of A237583 of terms ending in 3 in base 6. - Eric M. Schmidt, Jul 18 2015

References

  • C. P. Schut, Idempotents. Report AM-R9101, Centrum voor Wiskunde en Informatica, Amsterdam, 1991.

Crossrefs

Programs

  • Sage
    def a(n) : return crt(1, 0, 2^n, 3^n) # Eric M. Schmidt, Jul 18 2015

Extensions

More terms from Eric M. Schmidt, Jul 18 2015

A259991 This sequence and A259990 are base-14 analogs of A007185 and A016090, written in base 10.

Original entry on oeis.org

8, 148, 344, 36016, 151264, 1764736, 46941952, 1101076992, 15858967552, 139825248256, 2453862488064, 14602557997056, 127990382747648, 921705156001792, 56481739283791872, 523186025957228544, 15768859390622826496, 198716939766610001920, 3186868919241067200512
Offset: 1

Views

Author

N. J. A. Sloane, Jul 13 2015

Keywords

Comments

See Schut (1991) for precise definition.
Ignoring repetitions, the subsequence of A201919 of terms ending in 8 in base 14. - Eric M. Schmidt, Jul 18 2015

References

  • C. P. Schut, Idempotents. Report AM-R9101, Centrum voor Wiskunde en Informatica, Amsterdam, 1991.

Crossrefs

Programs

  • Sage
    def a(n) : return crt(0, 1, 2^n, 7^n) # Eric M. Schmidt, Jul 18 2015

Extensions

More terms from Eric M. Schmidt, Jul 18 2015

A214880 5-adic valuation of A007185.

Original entry on oeis.org

1, 2, 4, 4, 5, 6, 7, 8, 9, 10, 11, 12, 17, 14, 15, 18, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 29, 29, 30, 33, 33, 33, 35, 35, 36, 37, 38, 40, 40, 41, 42, 43, 44, 45, 46, 47, 48, 50, 51, 54, 54, 54, 54, 55, 56, 57, 58, 59, 61, 61, 63, 63, 66, 66, 66, 67
Offset: 1

Views

Author

Eric M. Schmidt, Jul 29 2012

Keywords

Comments

a(n) >= n.
It appears that the density of n such that a(n) = n is 4/5.

Crossrefs

Programs

  • GAP
    A214880 := function(max) local result, i; result := [5]; for i in [2..max] do Add(result, result[i-1]^2 mod 10^i); od; result := List(result, x->Valuation(x, 5)); return result; end;

A214882 A007185(n)/5^n.

Original entry on oeis.org

1, 1, 5, 1, 29, 57, 37, 33, 109, 841, 373, 3761, 8125, 9817, 8517, 41025, 73741, 67177, 118293, 967377, 822621, 3100537, 6492133, 16397921, 33478573, 13406601, 83211957, 177703665, 35540733, 114482329, 881889925, 176377985, 35275597, 15468937385, 3093787477
Offset: 1

Views

Author

Eric M. Schmidt, Jul 31 2012

Keywords

Comments

Conjecture: For any odd m and for any k, the density of n such that a(n) == k (mod m) is 1/m.

Crossrefs

Programs

  • GAP
    A214882 := function(max) local result, i; result := [5]; for i in [2..max] do Add(result, result[i-1]^2 mod 10^i); od; result := List([1..max], n->result[n]/5^n); return result; end;

A259987 This sequence and A259986 are base 6 analogs of A007185 and A016090, written in base 10.

Original entry on oeis.org

4, 28, 136, 1216, 1216, 16768, 203392, 1043200, 6082048, 36315136, 36315136, 1487503360, 1487503360, 1487503360, 314944159744, 785129144320, 12069568774144, 45922887663616, 452162714337280, 2280241934368768, 2280241934368768, 2280241934368768, 2280241934368768
Offset: 1

Views

Author

N. J. A. Sloane, Jul 13 2015

Keywords

Comments

See Schut (1991) for precise definition.
Ignoring repetitions, the subsequence of A237583 of terms ending in 4 in base 6. - Eric M. Schmidt, Jul 18 2015

References

  • C. P. Schut, Idempotents. Report AM-R9101, Centrum voor Wiskunde en Informatica, Amsterdam, 1991.

Crossrefs

Programs

  • Sage
    def a(n) : return crt(0, 1, 2^n, 3^n) # Eric M. Schmidt, Jul 18 2015

Extensions

More terms from Eric M. Schmidt, Jul 18 2015

A259988 This sequence and A259989 are base-6 analogs of A007185 and A016090, written in base 6.

Original entry on oeis.org

3, 13, 213, 213, 50213, 350213, 1350213, 21350213, 221350213, 2221350213, 52221350213, 152221350213, 5152221350213, 55152221350213, 155152221350213, 4155152221350213, 14155152221350213, 314155152221350213, 1314155152221350213, 21314155152221350213
Offset: 1

Views

Author

N. J. A. Sloane, Jul 13 2015

Keywords

Comments

See Schut (1991) for precise definition.
Ignoring repetitions, the subsequence of A201821 of terms ending in 3. - Eric M. Schmidt, Jul 18 2015

References

  • C. P. Schut, Idempotents. Report AM-R9101, Centrum voor Wiskunde en Informatica, Amsterdam, 1991.

Crossrefs

Programs

  • Sage
    def a(n) : return Integer(crt(1, 0, 2^n, 3^n).str(6)) # Eric M. Schmidt, Jul 18 2015

Extensions

More terms from Eric M. Schmidt, Jul 18 2015

A259989 This sequence and A259988 are base-6 analogs of A007185 and A016090, written in base 6.

Original entry on oeis.org

4, 44, 344, 5344, 5344, 205344, 4205344, 34205344, 334205344, 3334205344, 3334205344, 403334205344, 403334205344, 403334205344, 400403334205344, 1400403334205344, 41400403334205344, 241400403334205344, 4241400403334205344, 34241400403334205344
Offset: 1

Views

Author

N. J. A. Sloane, Jul 13 2015

Keywords

Comments

See Schut (1991) for precise definition.
Ignoring repetitions, the subsequence of A201821 of terms ending in 4. - Eric M. Schmidt, Jul 18 2015

References

  • C. P. Schut, Idempotents. Report AM-R9101, Centrum voor Wiskunde en Informatica, Amsterdam, 1991.

Crossrefs

Programs

  • Sage
    def a(n) : return Integer(crt(0, 1, 2^n, 3^n).str(6)) # Eric M. Schmidt, Jul 18 2015

Extensions

Corrected and extended by Eric M. Schmidt, Jul 18 2015

A259990 This sequence and A259991 are base-14 analogs of A007185 and A016090, written in base 10.

Original entry on oeis.org

7, 49, 2401, 2401, 386561, 5764801, 58471553, 374712065, 4802079233, 149429406721, 1595702681601, 42091354378241, 665724390506497, 10190301669556225, 99086356274020353, 1654767311852142593, 14722487338708369409, 228161914444026740737, 2789435039707847196673
Offset: 1

Views

Author

N. J. A. Sloane, Jul 13 2015

Keywords

Comments

See Schut (1991) for precise definition.
Ignoring repetitions, the subsequence of A201919 of terms ending in 7 in base 14. - Eric M. Schmidt, Jul 18 2015

References

  • C. P. Schut, Idempotents. Report AM-R9101, Centrum voor Wiskunde en Informatica, Amsterdam, 1991.

Crossrefs

Programs

  • Sage
    def a(n) : return crt(1, 0, 2^n, 7^n) # Eric M. Schmidt, Jul 18 2015

Extensions

More terms from Eric M. Schmidt, Jul 18 2015

A003226 Automorphic numbers: m^2 ends with m.

Original entry on oeis.org

0, 1, 5, 6, 25, 76, 376, 625, 9376, 90625, 109376, 890625, 2890625, 7109376, 12890625, 87109376, 212890625, 787109376, 1787109376, 8212890625, 18212890625, 81787109376, 918212890625, 9918212890625, 40081787109376, 59918212890625, 259918212890625, 740081787109376
Offset: 1

Views

Author

Keywords

Comments

Also called curious numbers.
For entries after the second, two successive terms sum up to a total having the form 10^n + 1. - Lekraj Beedassy, Apr 29 2005 [This comment is clearly wrong as stated. The sums of two consecutive terms are 1, 6, 11, 31, 101, 452, 1001, 10001, 100001, 200001, 1000001, 3781250, .... - T. D. Noe, Nov 14 2010]
If a d-digit number n is in the sequence, then so is 10^d+1-n. However, the same number can be 10^d+1-n for different n in the sequence (e.g., 10^3+1-376 = 10^4+1-9376 = 625), which spoils Beedassy's comment. - Robert Israel, Jun 19 2015
Substring of both its square and its cube not congruent to 0 (mod 10). See A029943. - Robert G. Wilson v, Jul 16 2005
a(n)^k ends with a(n) for k > 0; see also A029943. - Reinhard Zumkeller, Nov 26 2011
Apart from initial term, a subsequence of A046831. - M. F. Hasler, Dec 05 2012
This is also the sequence of numbers such that the n-th m-gonal number ends in n for any m == 0,4,8,16 (mod 20). - Robert Dawson, Jul 09 2018
Apart from 6, a subsequence of A301912. - Robert Dawson, Aug 01 2018

References

  • J.-M. De Koninck, Ces nombres qui nous fascinent, Entry 76, p. 26, Ellipses, Paris 2008.
  • V. deGuerre and R. A. Fairbairn, Automorphic numbers, J. Rec. Math., 1 (No. 3, 1968), 173-179.
  • R. A. Fairbairn, More on automorphic numbers, J. Rec. Math., 2 (No. 3, 1969), 170-174.
  • Jan Gullberg, Mathematics, From the Birth of Numbers, W. W. Norton & Co., NY, page 253-254.
  • B. A. Naik, 'Automorphic numbers' in 'Science Today'(subsequently renamed '2001') May 1982 pp. 59, Times of India, Mumbai.
  • Ya. I. Perelman, Algebra can be fun, pp. 97-98.
  • Clifford A. Pickover, A Passion for Mathematics, John Wiley & Sons, Hoboken, 2005, p. 64.
  • C. P. Schut, Idempotents. Report AM-R9101, Centrum voor Wiskunde en Informatica, Amsterdam, 1991.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Haskell
    import Data.List (isSuffixOf)
    a003226 n = a003226_list !! (n-1)
    a003226_list = filter (\x -> show x `isSuffixOf` show (x^2)) a008851_list
    -- Reinhard Zumkeller, Jul 27 2011
    
  • Magma
    [n: n in [0..10^7] | Intseq(n^2)[1..#Intseq(n)] eq Intseq(n)]; // Vincenzo Librandi, Jul 03 2015
    
  • Maple
    V:= proc(m) option remember;
      select(t -> t^2 - t mod 10^m = 0, map(s -> seq(10^(m-1)*j+s, j=0..9), V(m-1)))
    end proc:
    V(0):= {0,1}:
    V(1):= {5,6}:
    sort(map(op,[V(0),seq(V(i) minus V(i-1),i=1..50)])); # Robert Israel, Jun 19 2015
  • Mathematica
    f[k_] := (r = Reduce[0 < 10^k < n < 10^(k + 1) && n^2 == m*10^(k + 1) + n, {n, m}, Integers]; If[Head[r] === And, n /. ToRules[r], n /. {ToRules[r]}]); Flatten[ Join[{0, 1}, Table[f[k], {k, 0, 13}]]] (* Jean-François Alcover, Dec 01 2011 *)
    Union@ Join[{1}, Array[PowerMod[5, 2^#, 10^#] &, 16, 0], Array[PowerMod[16, 5^#, 10^#] &, 16, 0]] (* Robert G. Wilson v, Jul 23 2018 *)
  • PARI
    is_A003226(n)={n<2 || 10^valuation(n^2-n,10)>n} \\ M. F. Hasler, Dec 05 2012
    
  • PARI
    A003226(n)={ n<3 & return(n-1); my(i=10,j=10,b=5,c=6,a=b); for( k=4,n, while(b<=a, b=b^2%i*=10); while(c<=a, c=(2-c)*c%j*=10); a=min(b,c)); a } \\ M. F. Hasler, Dec 06 2012
    
  • Python
    from itertools import count, islice
    from sympy.ntheory.modular import crt
    def A003226_gen(): # generator of terms
        a = 0
        yield from (0,1)
        for n in count(0):
            b = sorted((int(crt(m:=(1< a:
                yield from b
                a = b[1]
            elif b[1] > a:
                yield b[1]
                a = b[1]
    A003226_list = list(islice(A003226_gen(),15)) # Chai Wah Wu, Jul 25 2022
  • Sage
    def automorphic(maxdigits, pow, base=10) :
        morphs = [[0]]
        for i in range(maxdigits):
            T=[d*base^i+x for x in morphs[-1] for d in range(base)]
            morphs.append([x for x in T if x^pow % base^(i+1) == x])
        res = list(set(sum(morphs, []))); res.sort()
        return res
    # call with pow=2 for this sequence, Eric M. Schmidt, Feb 09 2014
    

Formula

Equals {0, 1} union A007185 union A016090.

Extensions

More terms from Michel ten Voorde, Apr 11 2001
Edited by David W. Wilson, Sep 26 2002
Incorrect statement removed from title by Robert Dawson, Jul 09 2018

A018247 The 10-adic integer x = ...8212890625 satisfying x^2 = x.

Original entry on oeis.org

5, 2, 6, 0, 9, 8, 2, 1, 2, 8, 1, 9, 9, 5, 2, 6, 5, 2, 2, 9, 3, 7, 7, 9, 9, 1, 6, 6, 0, 1, 4, 0, 0, 9, 0, 1, 6, 9, 8, 0, 3, 2, 3, 2, 4, 3, 2, 4, 7, 5, 5, 0, 0, 0, 1, 1, 8, 3, 6, 8, 0, 8, 5, 9, 0, 5, 6, 6, 1, 2, 6, 0, 0, 9, 8, 9, 0, 5, 8, 3, 9, 2, 0, 8, 9, 6, 1, 8, 0, 1, 9, 1, 3, 7, 0, 0, 3, 5, 9, 3, 0, 9, 3, 6, 2, 4, 6, 7
Offset: 0

Views

Author

Yoshihide Tamori (yo(AT)salk.edu)

Keywords

Comments

The 10-adic numbers a and b defined in this sequence and A018248 satisfy a^2=a, b^2=b, a+b=1, ab=0. - Michael Somos

Examples

			x = ...0863811000557423423230896109004106619977392256259918212890625.
		

References

  • W. W. R. Ball, Mathematical Recreations & Essays, N.Y. Macmillan Co, 1947.
  • V. deGuerre and R. A. Fairbairn, Jnl. Rec. Math., No. 3, (1968), 173-179.
  • M. Kraitchik, Sphinx, 1935, p. 1.

Crossrefs

A007185 gives associated automorphic numbers.
The difference between A018248 & this sequence is A075693 and their product is A075693.
The six examples given by deGuerre and Fairbairn are A055620, A054869, A018247, A018248, A259468, A259469.

Programs

  • Mathematica
    a = {5}; f[n_] := Block[{k = 0, c}, While[c = FromDigits[Prepend[a, k]]; Mod[c^2, 10^n] != c, k++ ]; a = Prepend[a, k]]; Do[ f[n], {n, 2, 105}]; Reverse[a]
    With[{n = 150}, Reverse[IntegerDigits[PowerMod[5, 2^n, 10^n]]]] (* IWABUCHI Yu(u)ki, Feb 16 2024 *)
  • PARI
    a(n)=local(t=5);for(k=1,n+1,t=t^2%10^k);t\10^n \\ Paul D. Hanna, Jul 08 2006
    
  • PARI
    Vecrev(digits(lift(chinese(Mod(1, 2^100), Mod(0, 5^100))))) \\ Seiichi Manyama, Aug 07 2019

Formula

x = 10-adic lim_{n->oo} 5^(2^n) mod 10^(n+1). - Paul D. Hanna, Jul 08 2006

Extensions

More terms from David W. Wilson
Edited by David W. Wilson, Sep 26 2002
Showing 1-10 of 34 results. Next