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 40 results. Next

A075693 Difference between 10-adic numbers defined in A018248 & A018247.

Original entry on oeis.org

1, 5, -3, 9, -9, -7, 5, 7, 5, -7, 7, -9, -9, -1, 5, -3, -1, 5, 5, -9, 3, -5, -5, -9, -9, 7, -3, -3, 9, 7, 1, 9, 9, -9, 9, 7, -3, -9, -7, 9, 3, 5, 3, 5, 1, 3, 5, 1, -5, -1, -1, 9, 9, 9, 7, 7, -7, 3, -3, -7, 9, -7, -1, -9, 9, -1, -3, -3, 7, 5, -3, 9, 9, -9, -7, -9, 9, -1, -7, 3, -9, 5, 9, -7
Offset: 0

Views

Author

Robert G. Wilson v, Sep 26 2002

Keywords

Comments

Numbers in A018247 and A018248 are known as automorphic numbers in base 10, meaning that the infinite integers a=(...256259918212890625) or b=(...743740081787109376) provides a nontrivial solution to x*x == x (mod any power of 10).
Read backwards so as to match their counterparts (A007185 & A016090), A018247(0)+A018248(0) = 11 & A018247(n)+A018248(n) = 9 for all n's > 0 and their product is A076308.
All entries must be odd.
Is the accumulative sum equally positive and negative, i.e. does the sum equal 0 infinitely often?

Crossrefs

Programs

  • Mathematica
    (* execute the programming in both A018247 & A018248 *) Reverse[b - a]
    aa[n_] := For[t = 5; k = 1, True, k++, t = Mod[t^2, 10^k]; If[k == n, Return[ Quotient[t, 10^(n-1)]]]]; bb[n_] := Reap[ For[t = 6; k = 1, k <= n , k++, t = Mod[t^5, 10^k]; Sow[ Quotient[10*t, 10^k]]]][[2, 1, n]]; a[n_] := bb[n] - aa[n]; Table[a[n], {n, 1, 84}](* Jean-François Alcover, May 25 2012, after Paul D. Hanna *)

Formula

a(n) = A018248(n) - A018247(n). - Seiichi Manyama, Jul 26 2017

A076308 Product of 10-adic numbers defined in A018247 and A018248.

Original entry on oeis.org

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

Views

Author

Robert G. Wilson v, Oct 05 2002

Keywords

Crossrefs

Programs

  • Mathematica
    (* execute the programming in both A018247 & A018248 *) IntegerDigits[ FromDigits[ Reverse[a]]*FromDigits[ Reverse[b]]]

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

A317905 Convergence speed of m^^m, where m = A067251(n) and n >= 2. a(n) = f(m, m) - f(m, m - 1), where f(x, y) corresponds to the maximum value of k, such that x^^y == x^^(y + 1) (mod 10^k).

Original entry on oeis.org

0, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 4, 1, 1, 2, 1, 1, 1, 1, 2, 3, 2, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 2, 2, 1, 1, 1, 3, 1, 3, 1, 1, 1, 1, 1, 1, 6, 1, 1, 3, 1, 1, 1, 1, 2, 2, 2, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 1, 5, 1, 1
Offset: 2

Views

Author

Marco Ripà, Aug 10 2018

Keywords

Comments

It is possible to anticipate the convergence speed of m^^m, where ^^ indicates tetration or hyper-4 (e.g., 3^^4 = 3^(3^(3^3))), simply looking at the congruence (mod 25) of m. In fact, assuming m > 2, a(n) = 1 for any m == 2, 3, 4, 6, 8, 9, 11, 12, 13, 14, 16, 17, 19, 21, 22, 23 (mod 25), and a(n) >= 2 otherwise.
It follows that 32/45 = 71.11% of the a(n) assume unitary value.
You can also obtain an arbitrary high convergence speed, such as taking the beautiful base b = 999...99 (9_9_9... n times), which gives a(n) = len(b), for any len(b) > 1. Thus, 99...9^^m == 99...9^^(m + 1) (mod m*10^len(b)), as proved by Ripà in "La strana coda della serie n^n^...^n", pages 25-26. In fact, m = 99...9 == 24 (mod 25) and a(m=24) > 1.
From Marco Ripà, Dec 19 2021: (Start)
Knowing the "constant congruence speed" of a given base (a.k.a. the convergence speed of the base m, assuming m > 2) is very useful in order to calculate the exact number of stable digits of all its tetrations of height b > 1. As an example, let us consider all the a(n) such that n is congruent to 4 (mod 9) (i.e., all the tetration bases belonging to the congruence class 5 (mod 10)). Then, the exact number of stable digits (#S(m, b)) of any tetration m^^b (i.e., the number of its last "frozen" digits) such that m is congruent to 5 (mod 10), for any b >= 3, can automatically be calculated by simply knowing that (under the stated constraint) the congruence speed of the m corresponds to the 2-adic valuation of (m^2 - 1) minus 1. Thus, let k = 1, 2, 3, ..., and we have that
If m = 20*k - 5, then #S(m, b > 2) = b*(v_2(m^2 - 1) - 1) + 1 = b*(v_2(m + 1) + 1);
If m = 20*k + 5, then #S(m, b > 2) = (b + 1)*(v_2(m^2 - 1) - 1) = (b + 1)*(v_2(m - 1));
If m = 5, then #S(m, 1) = 1, #S(m, 2) = 4, #S(m, b > 2) = 8 + 2*(b - 3).
(End)
For any n > 2, the value of a(n) depends on the congruence modulo 18 of n, since the constant congruence speed of m arises from the 14 nontrivial solutions of the fundamental equation y^5 = y in the (commutative) ring of decadic integers (e.g., y = -1 = ...9999 is a solution of y^5 = y, so it originates the law a(n) = min(v_2(m + 1), v_5(m + 1)) concerning every n belonging to the congruence class 0 modulo 18, as stated in the "Formula" section of the present sequence). - Marco Ripà, Feb 17 2022
a(n) satisfies the following multiplicative constraint: for each pair (m_1, m_2) of terms of A067251, a(m_1*m_2) is necessarily greater than or equal to the minimum between a(m_1) and a(m_2) (see Equation 2.4 and Appendix of "A Compact Notation for Peculiar Properties Characterizing Integer Tetration" in Links). - Gabriele Di Pietro, Apr 29 2025

Examples

			For m = 25, a(23) = 3 implies that 25^^(25 + i) freezes 3*i "new" rightmost digits (i >= 0).
		

References

  • Marco Ripà, La strana coda della serie n^n^...^n, Trento, UNI Service, Nov 2011. ISBN 978-88-6178-789-6

Crossrefs

Programs

  • PARI
    \\ uses reducetower.gp from links
    f2(x,y) = my(k=0); while(reducetower(x, 10^k, y) == reducetower(x, 10^k, y+1), k++); k;
    f1(n) = polcoef(x*(x+1)*(x^4-x^3+x^2-x+1)*(x^4+x^3+x^2+x+1) / ((x-1)^2*(x^2+x+1)*(x^6+x^3+1)) + O(x^(n+1)), n, x); \\ A067251
    a(n) =  my(m=f1(n)); f2(m, m) - f2(m, m-1);
    lista(nn) = {for (n=2, nn, print1(a(n), ", "););} \\ Michel Marcus, Jan 27 2021

Formula

Let n > 2. For any integer c >= 0, if n is an element of the set {5, 7, 14, 17, 22, 23, 24, 29, 32, 39, 41, 45, 46}, then a(n + 45*c) >= 2; whereas a(n) = 1 otherwise. - Marco Ripà, Sep 28 2018
If n == 5 (mod 9), then a(n) = v_2(a(n)^2 - 1) - 1, where v_2(x) indicates the 2-adic valuation of x. - Marco Ripà, Dec 19 2021
If n == 1 (mod 18) and n <> 1, then a(n) = min(v_2(m - 1), v_5(m - 1)) (i.e., 1 plus the number of trailing zeros, if any, next to the rightmost digit of m);
if n == 10 (mod 18), then a(n) = min(v_2(m + 1), v_5(m - 1));
if n == {2,8}(mod 9) and n <> 2, then a(n) = v_5(m^2 + 1);
if n == {3,7}(mod 18), then a(n) = min(v_2(m + 1), v_5(n^2 + 1));
if n == {12,16}(mod 18), then a(n) = min(v_2(m - 1), v_5(n^2 + 1));
if n == 4 (mod 9), then a(n) = v_5(m + 1);
if n == 5 (mod 18), then a(n) = v_2(m - 1);
if n == 14 (mod 18), then a(n) = v_2(m + 1);
if n == 6 (mod 9), then a(n) = v_5(m - 1);
if n == 9 (mod 18), then a(n) = min(v_2(m - 1), v_5(m + 1));
if n == 0 (mod 18), then a(n) = min(v_2(m + 1), v_5(m + 1)) (i.e., number of digits of the rightmost repunit "9's" of m); where v_2(x) and v_5(x) indicates the 2-adic valuation of (x) and the 5-adic valuation of (x), respectively. - Marco Ripà, Feb 17 2022

Extensions

Edited by Jinyuan Wang, Aug 30 2020

A007185 Automorphic numbers ending in digit 5: a(n) = 5^(2^n) mod 10^n.

Original entry on oeis.org

5, 25, 625, 625, 90625, 890625, 2890625, 12890625, 212890625, 8212890625, 18212890625, 918212890625, 9918212890625, 59918212890625, 259918212890625, 6259918212890625, 56259918212890625, 256259918212890625, 2256259918212890625, 92256259918212890625
Offset: 1

Views

Author

Keywords

Comments

Conjecture: For any m coprime to 10 and for any k, the density of n such that a(n) == k (mod m) is 1/m. - Eric M. Schmidt, Aug 01 2012
a(n) is the unique positive integer less than 10^n such that a(n) is divisible by 5^n and a(n) - 1 is divisible by 2^n. - Eric M. Schmidt, Aug 18 2012

Examples

			625 is in the sequence because 625^2 = 390625, which ends in 625.
90625 is in the sequence because 90625^2 = 8212890625, which ends in 90625.
90635 is not in the sequence because 90635^2 = 8214703225, which does not end in 90635.
		

References

  • 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-4.
  • Ya. I. Perelman, Algebra can be fun, pp. 97-98.
  • 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

A018247 gives the associated 10-adic number.
A003226 = {0, 1} union (this sequence) union A016090.

Programs

Formula

a(n) = 5^(2^n) mod 10^n.
a(n)^2 == a(n) (mod 10^n), that is, a(n) is an idempotent in Z[10^n].
a(n+1) = a(n)^2 mod 10^(n+1). - Eric M. Schmidt, Jul 28 2012
a(2n) = (3*a(n)^2 - 2*a(n)^3) mod 10^(2n). - Sylvie Gaudel, Mar 10 2018

Extensions

More terms from David W. Wilson
Edited by David W. Wilson, Sep 26 2002
Further edited by N. J. A. Sloane, Jul 21 2010
Comment moved to name by Alonso del Arte, Mar 10 2018

A018248 The 10-adic integer x = ...1787109376 satisfies x^2 = x.

Original entry on oeis.org

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

Views

Author

Yoshihide Tamori (yo(AT)salk.edu)

Keywords

Comments

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

Examples

			x equals the limit of the (n+1) trailing digits of 6^(5^n):
6^(5^0)=(6), 6^(5^1)=77(76), 6^(5^2)=28430288029929701(376), ...
x = ...9442576576769103890995893380022607743740081787109376.
From _Peter Bala_, Nov 05 2022: (Start)
Trailing digits of 2^(10^n), 4^(10^n) and 6^(10^n) for n = 5:
2^(10^5) = ...9883(109376);
4^(10^5) = ...7979(109376);
6^(10^5) = ...4155(109376). (End)
		

References

  • W. W. R. Ball, Mathematical Recreations & Essays, N.Y. Macmillan Co, 1947.
  • R. Cuculière, Jeux Mathématiques, in Pour la Science, No. 6 (1986), 10-15.
  • V. deGuerre and R. A. Fairbairn, Automorphic numbers, J. Rec. Math., 1 (No. 3, 1968), 173-179.
  • M. Kraitchik, Sphinx, 1935, p. 1.
  • A. M. Robert, A Course in p-adic Analysis, Springer, 2000; see pp. 63, 419.

Crossrefs

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

Programs

  • Maple
    a := proc (n) option remember; if n = 1 then 2 else irem(a(n-1)^10, 10^n) end if; end proc:
    # display the digits of a(100) from right to left
    S := convert(a(100), string):
    with(ListTools):
    the_List := [seq(parse(S[i]), i = 1..length(S))]:
    Reverse(the_List); # Peter Bala, Nov 04 2022
  • Mathematica
    b = {6}; g[n_] := Block[{k = 0, c}, While[c = FromDigits[Prepend[b, k]]; Mod[c^2, 10^n] != c, k++ ]; b = Prepend[b, k]]; Do[ g[n], {n, 2, 105}]; Reverse[b]
    With[{n = 150}, Reverse[IntegerDigits[PowerMod[16, 5^n, 10^n]]]] (* IWABUCHI Yu(u)ki, Feb 16 2024 *)
  • PARI
    {a(n)=local(b=6,v=[]);for(k=1,n+1,b=b^5%10^k;v=concat(v,(10*b\10^k)));v[n+1]} \\ Paul D. Hanna, Jul 06 2006
    
  • PARI
    Vecrev(digits(lift(chinese(Mod(0, 2^100), Mod(1, 5^100))))) \\ Seiichi Manyama, Aug 07 2019

Formula

x = r^4 where r=...1441224165530407839804103263499879186432 (A120817). x = 10-adic limit_{n->oo} 6^(5^n). - Paul D. Hanna, Jul 06 2006
For n >= 2, the final n+1 digits of either 2^(10^n), 4^(10^n) or 6^(10^n), when read from right to left, give the first n+1 entries in the sequence. - Peter Bala, Nov 05 2022

Extensions

More terms from David W. Wilson
Edited by David W. Wilson, Sep 26 2002

A091664 10-adic integer x=.....06619977392256259918212890624 satisfying x^3 = x.

Original entry on oeis.org

4, 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
Offset: 0

Views

Author

Edoardo Gueglio (egueglio(AT)yahoo.it), Jan 28 2004

Keywords

Comments

Let a,b be integers defined in A018247, A018248 satisfying a^2=a, b^2=b, obviously a^3=a, b^3=b; let c,d,e,f be integers defined in A091661, A063006, A091663, A091664; then c^3=c, d^3=d, e^3=e, f^3=f, c+d=1, a+e=1, b+f=1, b+c=a, d+f=e, a+f=c, a=f+1, b=e+1, cd=-1, af=-1, gh=-1 where -1=.....999999999.

Examples

			x equals the limit of the (n+1) trailing digits of 4^(5^n):
4^(5^0) = (4), 4^(5^1) = 10(24), 4^(5^2) = 1125899906842(624), ...
x = ...0557423423230896109004106619977392256259918212890624.
		

Crossrefs

Programs

  • Mathematica
    To calculate c, d, e, f use Mathematica algorithms for a, b and equations: c=a-b, d=1-c, e=b-1, f=a-1.
  • PARI
    {a(n)=local(b=4,v=[]); for(k=1, n+1, b=b^5%10^k; v=concat(v,(10*b\10^k))); v[n+1]} \\ Paul D. Hanna, Jul 06 2006
    
  • PARI
    (A091664_vec(n)=Vecrev(digits(lift(chinese(Mod(0,2^n),Mod(-1,5^n))))))(99) \\ M. F. Hasler, Jan 26 2020

Formula

x = r^2 where r = ...1441224165530407839804103263499879186432 (A120817). x = 10-adic lim_{n->oo} 4^(5^n). - Paul D. Hanna, Jul 06 2006
For n > 0, a(n) = 9 - A018248(n) = A018247(n). - Seiichi Manyama, Jul 28 2017

A091663 10-adic integer x=.....93380022607743740081787109375 satisfying x^3 = x.

Original entry on oeis.org

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

Views

Author

Edoardo Gueglio (egueglio(AT)yahoo.it), Jan 28 2004

Keywords

Comments

Let a,b be integers defined in A018247, A018248 satisfying a^2=a, b^2=b, obviously a^3=a, b^3=b; let c,d,e,f be integers defined in A091661, A063006, A091663, A091664 then c^3=c, d^3=d, e^3=e, f^3=f, c+d=1, a+e=1, b+f=1, b+c=a, d+f=e, a+f=c, a=f+1, b=e+1, cd=-1, af=-1, gh=-1 where -1=.....999999999.

Crossrefs

Programs

  • Mathematica
    To calculate c, d, e, f use Mathematica algorithms for a, b and equations: c=a-b, d=1-c, e=b-1, f=a-1.

Formula

For n > 0, a(n) = 9 - A018247(n) = A018248(n). - Seiichi Manyama, Jul 28 2017

A373387 Constant congruence speed of the tetration base n (in radix-10), or -1 if n is a multiple of 10.

Original entry on oeis.org

0, 0, 1, 1, 1, 2, 1, 2, 1, 1, -1, 1, 1, 1, 1, 4, 1, 1, 2, 1, -1, 1, 1, 1, 2, 3, 2, 1, 1, 1, -1, 1, 2, 1, 1, 2, 1, 1, 1, 1, -1, 1, 1, 2, 1, 2, 1, 1, 1, 2, -1, 2, 1, 1, 1, 3, 1, 3, 1, 1, -1, 1, 1, 1, 1, 6, 1, 1, 3, 1, -1, 1, 1, 1, 2, 2, 2, 1, 1, 1, -1, 1, 2, 1
Offset: 0

Views

Author

Marco Ripà, Jun 02 2024

Keywords

Comments

It has been proved that this sequence contains arbitrarily large entries, while a(0) = a(1) = 0 by definition (given the fact that 0^0 = 1 is a reasonable choice and then 0^^b is 1 if b is even, whereas 0^^b is 0 if b is even). For any nonnegative integer n which is not a multiple of 10, a(n) is given by Equation (16) of the paper "Number of stable digits of any integer tetration" (see Links).
Moreover, a sufficient condition for having a constant congruence speed of any tetration base n, greater than 1 and not a multiple of 10, is that b >= 2 + v(n), where v(n) is equal to
u_5(n - 1) iff n == 1 (mod 5),
u_5(n^2 + 1) iff n == 2,3 (mod 5),
u_5(n + 1) iff n == 4 (mod 5),
u_2(n^2 - 1) - 1 iff n == 5 (mod 10)
(u_5 and u_2 indicate the 5-adic and the 2-adic valuation of the argument, respectively).
Therefore b >= n + 1 is always a sufficient condition for the constancy of the congruence speed (as long as n > 1 and n <> 0 (mod 10)).
As a trivial application of this property, we note that the constant congruence speed of the tetration 3^^b is 1 for any b > 1, while 3^3 is not congruent to 3 modulo 10. Thus, we can easily calculate the exact number of the rightmost digits of Graham’s number, G(64) (see A133613), that are the same of the homologous rightmost digits of 3^3^3^... since 3^3 is not congruent to 3 modulo 10, while the congruence speed of n = 3 is constant from height 2 (see A372490). This means that the last slog_3(G(64))-1 digits of G(64) are the same slog_3(G(64))-1 final digits of 3^3^3^..., whereas the difference between the slog_3(G(64))-th digit of G(64) and the slog_3(G(64))-th digit of 3^3^3^... is congruent to 6 modulo 10.
The constant congruence speed of tetration satisfies the following multiplicative constraint: for each pair (n_1, n_2) of nonnegative integers whose product is not divisible by 10, a(n_1*n_2) is necessarily greater than or equal to the minimum between a(n_1) and a(n_2) (see Equation 2.4 and Appendix of "A Compact Notation for Peculiar Properties Characterizing Integer Tetration" in Links). - Marco Ripà, Apr 26 2025

Examples

			a(3) = 1 since 3^^b := 3^3^3^... freezes 1 more rightmost digit for each unit increment of b, starting from b = 2.
		

References

  • Marco Ripà, La strana coda della serie n^n^...^n, Trento, UNI Service, Nov 2011. ISBN 978-88-6178-789-6.

Crossrefs

Programs

  • Python
    def v2(n):
        count = 0
        while n % 2 == 0 and n > 0:
            n //= 2
            count += 1
        return count
    def v5(n):
        count = 0
        while n % 5 == 0 and n > 0:
            n //= 5
            count += 1
        return count
    def V(a):
        mod_20 = a % 20
        mod_10 = a % 10
        if mod_20 == 1:
            return min(v2(a - 1), v5(a - 1))
        elif mod_20 == 11:
            return min(v2(a + 1), v5(a - 1))
        elif mod_10 in {2, 8}:
            return v5(a ** 2 + 1)
        elif mod_20 in {3, 7}:
            return min(v2(a + 1), v5(a ** 2 + 1))
        elif mod_20 in {13, 17}:
            return min(v2(a - 1), v5(a ** 2 + 1))
        elif mod_10 == 4:
            return v5(a + 1)
        elif mod_20 == 5:
            return v2(a - 1)
        elif mod_20 == 15:
            return v2(a + 1)
        elif mod_10 == 6:
            return v5(a - 1)
        elif mod_20 == 9:
            return min(v2(a - 1), v5(a + 1))
        elif mod_20 == 19:
            return min(v2(a + 1), v5(a + 1))
    def generate_sequence():
        sequence = []
        for a in range(1026):
            if a == 0 or a == 1:
                sequence.append(0)
            elif a % 10 == 0:
                sequence.append(-1)
            else:
                sequence.append(V(a))
        return sequence
    sequence = generate_sequence()
    print("a(0), a(1), a(2), ..., a(1025) =", ", ".join(map(str, sequence)))

Formula

a(n) = -1 iff n == 0 (mod 10), a(n) = 0 iff n = 1 or 2. Otherwise, a(n) >= 1 and it is given by Equation (16) from Ripà and Onnis.

A091661 Coefficients in a 10-adic square root of 1.

Original entry on oeis.org

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

Views

Author

Edoardo Gueglio (egueglio(AT)yahoo.it), Jan 28 2004

Keywords

Comments

10-adic integer x=.....239954784512519836425781249 satisfying x^3 = x.
Let a,b be integers defined in A018247, A018248 satisfying a^2=a, b^2=b, obviously a^3=a, b^3=b; let c,d,e,f be integers defined in A091661, A063006, A091663, A091664 then c^3=c, d^3=d, e^3=e, f^3=f, c+d=1, a+e=1, b+f=1, b+c=a, d+f=e, a+f=c, a=f+1, b=e+1, cd=-1, af=-1, gh=-1 where -1=.....999999999.

Crossrefs

Another 10-adic root of 1 is given by A063006.

Programs

  • Mathematica
    To calculate c, d, e, f use Mathematica algorithms for a, b and equations: c=a-b, d=1-c, e=b-1, f=a-1.
  • Ruby
    def A(s, n)
      n.times{|i|
        m = 10 ** (i + 1)
        (0..9).each{|j|
          k = j * m + s
          if (k ** 2 - k) % (m * 10) == 0
            s = k
            break
          end
        }
      }
      s
    end
    def A091661(n)
      str = (10 ** (n + 1) + A(5, n) - A(6, n)).to_s.reverse
      (0..n).map{|i| str[i].to_i}
    end
    p A091661(100) # Seiichi Manyama, Jul 31 2017

Formula

For n>0, a(n) = 9 - A063006(n).
Showing 1-10 of 40 results. Next