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-5 of 5 results.

A114258 Numbers k such that k^2 contains exactly 2 copies of each digit of k.

Original entry on oeis.org

72576, 406512, 415278, 494462, 603297, 725760, 3279015, 4065120, 4152780, 4651328, 4915278, 4927203, 4944620, 4972826, 4974032, 4985523, 4989323, 5002245, 5016125, 6032970, 6214358, 6415002, 6524235, 7257600, 9883667
Offset: 1

Views

Author

Giovanni Resta, Nov 18 2005

Keywords

Comments

From Chai Wah Wu, Feb 27 2024: (Start)
If k is a term, then k == 0 (mod 9) or k == 2 (mod 9) (see A370676).
First decimal digit of each term is 3 or larger. (End)

Examples

			72576 is in the sequence since its square 5267275776 contains four 7's, two 2's, two 5's and two 6's.
		

Crossrefs

Programs

  • Python
    from math import isqrt
    from itertools import count, islice
    def A114258_gen(): # generator of terms
        for l in count(1):
            a = isqrt(10**((l<<1)-1))
            if (a9:=a%9):
                a -= a9
            for b in range(a,10**l,9):
                for c in (0,2):
                    k = b+c
                    if sorted(str(k)*2)==sorted(str(k**2)):
                        yield k
    A114258_list = list(islice(A114258_gen(),20)) # Chai Wah Wu, Feb 27 2024

A199631 Numbers having each digit once and whose cube has each digit three times.

Original entry on oeis.org

4680215379, 4752360918, 4765380219, 4915280637, 5063248197, 5164738920, 5382417906, 5426370189, 5429013678, 5628130974, 5679321048, 5697841320, 5762831940, 5783610492, 5786430129, 5903467821, 6019285734, 6053147982, 6095721483, 6143720958, 6158723094
Offset: 1

Views

Author

T. D. Noe, Nov 09 2011

Keywords

Examples

			4680215379^3 = 102517384602327906545167884939.
		

Crossrefs

Cf. A050278 (pandigital numbers), A199630, A365144, A199632, A199633. Subsequence of A114259.

Programs

  • Mathematica
    t = Select[Permutations[Range[0, 9]], #[[1]] > 0 &]; t2 = Select[t, Union[DigitCount[FromDigits[#]^3]] == {3} &]; FromDigits /@ t2

A114260 Numbers k such that the 4th power of k contains exactly 4 copies of each digit of k.

Original entry on oeis.org

5702631489, 7264103985, 7602314895, 7824061395, 8105793624, 8174035962, 8304269175, 8904623175, 8923670541, 9451360827, 9785261403, 9804753612, 9846032571, 57026314890, 59730829461, 60947591328, 64017823995, 65190218436, 67024081935, 70645192839, 72641039850, 74991208356
Offset: 1

Views

Author

Giovanni Resta, Nov 18 2005

Keywords

Comments

First 13 terms of the sequence are also pandigital, i.e., they contain all the 10 digits at least once. This is probably accidental, but quite curious.
If a(n) is in the sequence, then 10*a(n) is also in the sequence. So 10^k*a(n) is also in the sequence for positive integers k. Thus this sequence differs from A365144. - Ray Chandler, Aug 23 2023
From David A. Corneth, Aug 30 2023: (Start)
Not all terms are pandigital, for example 65190218436, 75932056341 and 83581076421 are not.
For any k ends in, say, 425742 (which has six digits) the last six digits of k^4 are fixed. So if k ends in 425742 then k^4 ends in 425742^2 mod 10^6 = 318096 and so it must have a 3, 1, 8, 0, 9 and 6 none of which are contained in 425742.
Therefore if k ends in 425742 then it must have at least 12 digits. In a search for terms <= 10^11 this eases the search, in a search for terms <= 10^12 this leaves only 6! * 5 / 6 = 600 cases to check instead of 10^6.
An additional idea one might use is that the number of digits of k^4 must be a multiple of 4. I.e. 10^(4*m + 3) <= k^4 <= 10^(4*m + 4) so 10^m * 10^0.75 < m < 10^(m + 1) (all strict inequalities as 10^0.75 (cf. A210522) is irrational) which tells us a bunch about the leading digits of k.
Checking pandigital numbers separately might ease the search. That way if the union of some k and last q digits of k^4 is all decimal digits one could end the search there. This goes for example for 100426. If a term ends in 100426 then it has all decimal digits. (End)
All terms are divisible by 9. First decimal digit of a term is 5 or larger. - Chai Wah Wu, Feb 27 2024

Examples

			5702631489 is a term since its 4th power 1057550783692741389295697108242363408641 contains four 5's, four 7's, four 0's and so on.
		

Crossrefs

Programs

  • Python
    from itertools import count, islice
    from sympy import integer_nthroot
    def A114260_gen(): # generator of terms
        for l in count(1):
            a = integer_nthroot(10**(4*l-1),4)[0]
            if (a9:=a%9):
                a += 9-a9
            for b in range(a,10**l,9):
                if sorted(str(b)*4)==sorted(str(b**4)):
                    yield b
    A114260_list = list(islice(A114260_gen(), 5)) # Chai Wah Wu, Feb 27 2024

Extensions

a(14) from Ray Chandler, Aug 24 2023
More terms from David A. Corneth, Aug 30 2023

A114261 Numbers k such that the 5th power of k contains exactly 5 copies of each digit of k.

Original entry on oeis.org

961527834, 7351062489, 8105632794, 8401253976, 8731945026, 9164072385, 9238750614, 9615278340, 9847103256, 72308154699, 73510624890, 81056327940, 83170652949, 83792140506, 84012539760, 87319450260, 91602408573, 91640723850, 92387506140, 96152783400, 98471032560
Offset: 1

Views

Author

Giovanni Resta, Nov 18 2005

Keywords

Comments

Some of the early terms of the sequence are also pandigital, i.e. they contain all the 10 digits once. This is probably accidental, but quite curious!
All terms are divisible by 9. First decimal digit of a term is 6 or larger. - Chai Wah Wu, Feb 27 2024

Examples

			E.g. 961527834 is in the sequence since its 5th power 821881685441327565743977956591832631269739424 contains five 9's, five 6's, five 1's and so on.
		

Crossrefs

Programs

  • Python
    from itertools import count, islice
    from sympy import integer_nthroot
    def A114261_gen(): # generator of terms
        for l in count(1):
            a = integer_nthroot(10**(5*l-1),5)[0]
            if (a9:=a%9):
                a += 9-a9
            for b in range(a,10**l,9):
                if sorted(str(b)*5)==sorted(str(b**5)):
                    yield b
    A114261_list = list(islice(A114261_gen(),5)) # Chai Wah Wu, Feb 27 2024

Extensions

a(8)-a(9) from Ray Chandler, Aug 23 2023
a(10)-a(21) from Chai Wah Wu, Feb 28 2024

A363160 Smallest positive integer m with all digits distinct such that m^n contains each digit of m exactly n times, or -1 if no such m exists.

Original entry on oeis.org

1, 406512, 516473892, 5702631489, 961527834, 7025869314, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1
Offset: 1

Views

Author

Jean-Marc Rebert, Sep 07 2023

Keywords

Comments

For 7 <= n <= 185, I tried all possibilities with at most 10 distinct digits and I found no solution.
9876543210^186 has only 1859 < 186 * 10 = 1860 digits, so a(n) = -1 for n = 186.
So 9876543210^n has fewer than 10*n digits for n >= 186, so a(n) = -1 for n >= 186.

Examples

			a(1) = 1, because 1^1 = 1 has each digit of 1, 1 time, and no lesser number > 0 satisfies this.
a(2) = 406512, because 406512 has distinct digits, 406512^2 = 165252006144 has each digit of 406512, 2 times, and no lesser number satisfies this.
n a(n)        a(n)^n
1 1           1
2 406512      165252006144
3 516473892   137766973511455269432948288
4 5702631489  1057550783692741389295697108242363408641
5 961527834   821881685441327565743977956591832631269739424
6 7025869314  120281934463386157260042215510596389732740014997586987548736
		

Crossrefs

Programs

  • Mathematica
    hasDistinctDigitsQ[m_Integer?NonNegative]:=Length@IntegerDigits@m==Length@DeleteDuplicates@IntegerDigits@m;validNumberQ[n_Integer?NonNegative,m_Integer?NonNegative]:=AllTrue[Tally@IntegerDigits@m,Function[{digitFreq},MemberQ[Tally@IntegerDigits[m^n],{digitFreq[[1]],n*digitFreq[[2]]}]]];a[n_Integer?Positive,ex_Integer?Positive]:=Module[{m=1},Monitor[While[True,If[hasDistinctDigitsQ[m]&&validNumberQ[n,m],Return[m]];m++;If[m>10^(ex*n),Return[-1]];];m,m]];Table[a[n,7],{n,1,7}] (* Robert P. P. McKone, Sep 09 2023 *)
Showing 1-5 of 5 results.