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

A173426 a(n) is obtained by starting with 1, sequentially concatenating all decimal numbers up to n, and then, starting from n-1, sequentially concatenating all decimal numbers down to 1.

Original entry on oeis.org

1, 121, 12321, 1234321, 123454321, 12345654321, 1234567654321, 123456787654321, 12345678987654321, 12345678910987654321, 123456789101110987654321, 1234567891011121110987654321, 12345678910111213121110987654321, 123456789101112131413121110987654321
Offset: 1

Views

Author

Umut Uludag, Feb 18 2010

Keywords

Comments

The first prime in this sequence is the 20-digit number a(10) = 12345678910987654321. On Jul 20 2015, Shyam Sunder Gupta reported on the Number Theory Mailing List that he has found what is probably the second prime in the sequence. This is the 2446th term, namely the 17350-digit probable prime 1234567..244524462445..7654321. See A359148. - N. J. A. Sloane, Jul 29 2015 - Aug 03 2015
There are no other (PR)prime members in this sequence for n<60000. - Serge Batalov, Jul 29 2015
David Broadhurst gives heuristic arguments which suggest that this sequence contains infinitely many primes.
See A075023 and A075024 for the smallest and largest prime factor of the terms. - M. F. Hasler, Jul 29 2015
Using summation in decimal length clades, one can obtain analytical expressions for the sequence:
a(n) = A002275(n)^2, for 1 <= n < 10;
a(n) = (120999998998*10^(4*n-28) - 2*10^(2*n-9) + 8790000000121)/99^2, for 10 <= n < 10^2;
a(n) = (120999998998*10^(6*n-227) - (1099022*10^(6*n-406) + 242*10^(3*n-108) - 1087789*10^191)/111^2 + 8790000000121)/99^2, for 10^2 <= n < 10^3; etc. - Serge Batalov, Jul 29 2015
Curiously, 1234567891010987654321 is also a prime (see A259937). - N. J. A. Sloane, Nov 30 2021

References

  • D. Broadhurst, Primes from concatenation: results and heuristics, Number Theory List, Aug 01 2015 and later postings.

Crossrefs

This sequence and A002477 (Wonderful Demlo numbers) agree up to the 9th term.

Programs

  • Maple
    a:= n-> parse(cat($1..n, n-i$i=1..n-1)):
    seq(a(n), n=1..14);  # Alois P. Heinz, Dec 01 2021
  • Mathematica
    Table[FromDigits[Flatten[IntegerDigits/@Join[Range[n],Reverse[Range[ n-1]]]]],{n,15}] (* Harvey P. Dale, Sep 02 2015 *)
  • PARI
    A173426(n)=eval(concat(vector(n*2-1,k,if(kM. F. Hasler, Jul 29 2015
    
  • Python
    def A173426(n): return int(''.join(str(d) for d in range(1,n+1))+''.join(str(d) for d in range(n-1,0,-1))) # Chai Wah Wu, Dec 01 2021

Formula

a(n) = concatenate(1,2,3,...,n-2,n-1,n,n-1,n-2,...,3,2,1).

Extensions

More terms from and minor edits by M. F. Hasler, Jul 29 2015

A075024 a(n) is the largest prime divisor of the number A173426(n) = concatenate(1,2,...,n-1,n,n-1,...,2,1); a(1) = 1.

Original entry on oeis.org

1, 11, 37, 101, 271, 37, 4649, 137, 333667, 12345678910987654321, 17636684157301569664903, 2799473675762179389994681, 2354041513534224607850261, 2068140300159522133, 498056174529497, 112240064764214229701, 4188353169004802474320231191377
Offset: 1

Views

Author

Amarnath Murthy, Sep 01 2002

Keywords

Comments

Also for 1 < n < 10, a(n) is the common prime divisor for all A010785(m) which consist of n digits. - Alexander R. Povolotsky, Jun 05 2014, corrected by M. F. Hasler, Jul 30 2015
According to the definition (and given terms), this is the greatest prime factor (A006530) of A173426 and not of A002477, as an earlier formula asserted and which may have been an assumption of the preceding comment. - M. F. Hasler, Jul 29 2015

Examples

			a(5) = 271 as 123454321 = 41*41*271*271.
a(25) = 12471243489559387823527232424981012432152516319410549 is the larger factor of the semiprime A173426(24) = A075023(25) * a(n).
		

Crossrefs

Programs

  • Mathematica
    Table[FactorInteger[FromDigits[Join[Flatten[IntegerDigits/@Range[ n]], Flatten[ IntegerDigits/@Range[n-1,1,-1]]]]][[-1,1]],{n,20}] (* Harvey P. Dale, May 20 2016 *)
  • PARI
    a(n) = {if (n == 1, return (1)); s = ""; for (i=1, n, s = concat(s, Str(i));); forstep (i=n-1, 1, -1, s = concat(s, Str(i));); f = factor(eval(s)); f[#f~, 1];} \\ Michel Marcus, Jun 05 2014
    
  • PARI
    A075024(n)=A006530(A173426(n)) \\ A006530 should provide efficient code and also covers the case n=1. - M. F. Hasler, Jul 29 2015

Formula

a(n) = A006530(A173426(n)). - Michel Marcus, Jun 05 2014, corrected by M. F. Hasler, Jul 29 2015

Extensions

More terms from Sascha Kurz, Jan 03 2003
a(16)-a(17) from Michel Marcus, Jun 05 2014
More terms from M. F. Hasler, Jul 29 2015

A260587 Number of distinct prime factors of A173426(n) = concatenation of (1, 2, ..., n, n-1, ..., 1).

Original entry on oeis.org

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

Views

Author

M. F. Hasler, Jul 29 2015

Keywords

Examples

			a(2) = 1 since A173426(2) = 121 = 11*11 has only one distinct prime factor, 11.
a(21) = 5 since A173426(21) = 3^2 * 7 * 828703 * 94364768151913037621 * 250591098443370396365457961250972909.
a(25) = 2 since A173426(25) = A075023(n) * A075024(n) is a semiprime.
		

Crossrefs

Cf. A001221.
See A260588 for the variant where prime factors are counted with multiplicity.
See also A075023 and A075024 for the smallest and largest prime factor of the terms.

Programs

Extensions

a(38)-a(44) added using factordb.com by Jinyuan Wang, Mar 05 2020

A260588 Number of prime factors, counted with multiplicity, of A173426(n) = concatenation of (1, 2, ..., n, n-1, ..., 1).

Original entry on oeis.org

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

Views

Author

M. F. Hasler, Jul 29 2015

Keywords

Examples

			a(2) = 2 since A173426(2) = 121 = 11*11 has twice the factor 11.
a(21) = 6 since A173426(21) = 3 * 3 * 7 * 828703 * 94364768151913037621 * 250591098443370396365457961250972909.
		

Crossrefs

See A260587 for the variant where only distinct prime factors are counted.
See also A075023 and A075024 for the smallest and largest prime factor of the terms.
Cf. A001222.

Programs

Formula

a(n) = A001222(A173426(n)).

Extensions

Terms beyond a(30) via factorization results by Serge Batalov, added by M. F. Hasler, Jul 30 2015
a(38)-a(44) added using factordb.com by Jinyuan Wang, Mar 05 2020

A260589 Irregular table read by rows: n-th row lists the prime factors of A173426(n), with repetition.

Original entry on oeis.org

11, 11, 3, 3, 37, 37, 11, 11, 101, 101, 41, 41, 271, 271, 3, 3, 7, 7, 11, 11, 13, 13, 37, 37, 239, 239, 4649, 4649, 11, 11, 73, 73, 101, 101, 137, 137, 3, 3, 3, 3, 37, 37, 333667, 333667, 12345678910987654321, 7, 17636684157301569664903, 3, 3, 7, 7, 2799473675762179389994681, 1109, 4729
Offset: 1

Views

Author

M. F. Hasler, Jul 29 2015

Keywords

Comments

Row lengths are given by A260588(n). In particular, row n = 1 would have length 0, i.e., no element, because A173426(1) = 1 has no prime factors. Therefore the sequence can be considered to start with row n = 2. (The offset refers to the k-th element of the "flattened" sequence.)
For n = 1 through n = 9, A173426(n) is the square of the repunit 1...1 of length n, therefore every prime factor appears twice. This is no longer the case for n > 9.

Crossrefs

Programs

Formula

n | A173426(n) | factors = n-th row of this table
1 | 1 | []
2 | 121 | [11, 11]
3 | 12321 | [3, 3, 37, 37]
4 | 1234321 | [11, 11, 101, 101]
5 | 123454321 | [41, 41, 271, 271]
6 | 12345654321 | [3, 3, 7, 7, 11, 11, 13, 13, 37, 37]

A261411 a(1)=1; thereafter a(n) = smallest prime factor of A261570(n).

Original entry on oeis.org

1, 11, 3, 11, 41, 3, 239, 11, 3, 11, 11, 3, 11, 11, 3, 11, 11, 3, 71, 21557, 19, 17, 31, 181, 17, 353, 19, 31, 19, 29, 17, 29, 11616377, 214141, 19, 5471, 17, 13883, 3, 7, 421219193, 3, 17, 7, 3, 7, 101, 3, 634324033999, 13, 19, 13, 83, 13, 23, 13, 19, 13, 19
Offset: 1

Views

Author

N. J. A. Sloane, Aug 24 2015, based on Robert G. Wilson v's comment in A261570

Keywords

Examples

			A261570(11) = 12345678911987654321 = (11)(59)(34631)(43117)(6373)(1999), so a(10) = 11.
Note that a(2007) = A261570(2007) is a 21233-digit (probable) prime.
		

Crossrefs

Extensions

a(41)-a(59) from Hiroaki Yamanouchi, Aug 24 2015
Showing 1-6 of 6 results.