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

A371699 The smallest composite number which divides the concatenation of its descending ordered prime factors, with repetition, when written in base n.

Original entry on oeis.org

42, 4, 42, 4, 374, 4, 9, 4, 378, 4, 609, 4, 9, 4, 3525, 4, 343, 4, 9, 4, 70, 4, 25, 4, 9, 4, 195, 4, 343, 4, 9, 4, 25, 4, 130, 4, 9, 4, 366, 4, 3562, 4, 9, 4, 42, 4, 49, 4, 9, 4, 474, 4, 25, 4, 9, 4, 238, 4, 1131, 4, 9, 4, 25, 4, 555, 4, 9, 4, 14405, 4, 12207
Offset: 2

Views

Author

Chai Wah Wu, Apr 12 2024

Keywords

Comments

Conjecture: a(n) <= A371641(n) for n >= 2.

Examples

			a(2) = 42 since 42 = 7*3*2 = 111_2 * 11_2 * 10 _2 and 42 divides 126 = 1111110_2.
a(10) = 378 since 278 = 7*3*3*3*2 and 278 divides 73332.
		

Crossrefs

Programs

  • Python
    from itertools import count
    from sympy import factorint, integer_log
    def A371699(n):
        for m in count(4):
            f = factorint(m)
            if sum(f.values()) > 1:
                c = 0
                for p in sorted(f,reverse=True):
                    a = pow(n,integer_log(p,n)[0]+1,m)
                    for _ in range(f[p]):
                        c = (c*a+p)%m
                if not c:
                    return m

Formula

If p is prime, then a(p*(m+2)-1) <= p^2 for all m >= 0.
If n+1 is composite, then a(n) <= q^2, where q = A020639(n+1) is the smallest prime factor of n+1. This implies that if n > 2 is odd, then a(n) = A371641(n) = 4.
The first few terms n where n+1 is composite and a(n) < A020639(n+1)^2 are a(288) = 70, a(298) = 42, a(340) = 42, a(360) = 182, ...
If n is even, then a(n) >= 9. This is true as it is easy to verify that a(n) cannot be equal to 4, 6 or 8 in this case.
Suppose n>2 is even. 2 concatenated twice in base n is 2(n+1) which is not divisible by 4.
Next, 3 concatenated by 2 is 3*n+2 which is not divisible by 6. Finally 2 concatenated 3 times is 2(n^3-1)/(n-1) which is not divisible by 8 since n^3-1 is odd.
This implies that if n = 6*k+2 for some k > 0, then a(n) = A371641(n) = 9.

A371666 Composite numbers which equal the reverse of the concatenation of their ascending ordered prime factors, with repetition, when written in binary.

Original entry on oeis.org

623, 78205, 101239, 80073085, 1473273719
Offset: 1

Views

Author

Scott R. Shannon, Apr 02 2024

Keywords

Comments

A subsequence of A355973. The first five numbers each have only two prime factors - do terms exist with three or more?

Examples

			101239 is a term as 101239_10 = 29_10 * 3491_10 = 11101_2 * 110110100011_2 = "11101110110100011"_2 which when reversed is "11000101101110111"_2 = 101239_10.
		

Crossrefs

Programs

  • Python
    from itertools import count, islice
    from sympy import factorint
    def A371666_gen(startvalue=4): # generator of terms >= startvalue
        for n in count(max(startvalue,4)):
            f = factorint(n)
            if sum(f.values()) > 1:
                c = 0
                for p in sorted(f,reverse=True):
                    a, q = p.bit_length(), int(bin(p)[:1:-1],2)
                    for _ in range(f[p]):
                        c = (c<A371666_list = list(islice(A371666_gen(),3)) # Chai Wah Wu, Apr 13 2024

A372277 Composite numbers that divide the concatenation of the reverse of their ascending order prime factors, with repetition, when written in binary.

Original entry on oeis.org

87339, 332403, 9813039
Offset: 1

Views

Author

Scott R. Shannon, Apr 25 2024

Keywords

Comments

The base-2 version of A372046.
a(4) > 120000000000. - Robert P. P. McKone, May 20 2024

Examples

			332403 is a term as 332403 = 3 * 179 * 619 = 11_2 * 10110011_2 * 1001101011_2 = "11"_2 * "11001101"_2 * "1101011001"_2 when each prime factor is reversed. This gives "11110011011101011001"_2 when concatenated, and 11110011011101011001_2 = 997209 which is divisible by 332403.
		

Crossrefs

Programs

  • Mathematica
    a[n_Integer] := Module[{f}, f = Flatten[ConstantArray @@@ FactorInteger[n]]; If[Length[f] < 2, Return[False]]; Mod[FromDigits[StringJoin[StringReverse[IntegerString[#, 2]] & /@ f], 2], n] == 0];
    Select[Range[2, 10^5], a] (* Robert P. P. McKone, May 02 2024 *)
  • Python
    from itertools import count, islice
    from sympy import factorint
    def A372277_gen(startvalue=4): # generator of terms >= startvalue
        for n in count(max(startvalue,4)):
            f = factorint(n)
            if sum(f.values()) > 1:
                c = 0
                for p in sorted(f):
                    a = pow(2,len(s:=bin(p)[2:]),n)
                    q = int(s[::-1],2)
                    for _ in range(f[p]):
                        c = (c*a+q)%n
                if not c:
                    yield n
    A372277_list = list(islice(A372277_gen(),3)) # Chai Wah Wu, Apr 25 2024
Showing 1-3 of 3 results.