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.

A057897 Numbers which can be written as m^k-k, with m, k > 1.

Original entry on oeis.org

2, 5, 7, 12, 14, 23, 24, 27, 34, 47, 58, 61, 62, 77, 79, 98, 119, 121, 122, 142, 167, 194, 213, 223, 238, 248, 252, 254, 287, 322, 340, 359, 398, 439, 482, 503, 509, 527, 574, 621, 623, 674, 723, 726, 727, 782, 839, 898, 959, 997, 1014, 1019, 1022, 1087, 1154
Offset: 1

Views

Author

Henry Bottomley, Sep 26 2000

Keywords

Comments

It may be that positive integers can be written as m^k-k (with m and k > 1) in at most one way [checked up to 10000].
All numbers < 10^16 of this form have a unique representation. The uniqueness question leads to a Pillai-like exponential Diophantine equation a^x-b^y = x-y for x > y > 1 and b > a > 1, which appears to have no solutions. - T. D. Noe, Oct 06 2004

Crossrefs

Cf. A099225 (numbers of the form m^k+k, with m and k > 1), A074981 (n such that there is no solution to Pillai's equation), A099226 (numbers that can be represented as both a^x+x and b^y-y, for some a, b, x, y > 1).

Programs

  • Mathematica
    nLim=1000; lst={}; Do[k=2; While[n=m^k-k; n<=nLim, AppendTo[lst, n]; k++ ], {m, 2, Sqrt[nLim]}]; Union[lst] (* T. D. Noe, Oct 06 2004 *)
  • PARI
    ok(n)={my(e=2); while(2^e <= n+e, if(ispower(n+e, e), return(1)); e++); 0} \\ Andrew Howroyd, Oct 20 2020
    
  • PARI
    upto(lim)={my(p=logint(lim,2)); while(logint(lim+p+1,2)>p, p++); Vec(Set(concat(vector(p-1, e, e++; vector(sqrtnint(lim+e,e)-1, m, (m+1)^e-e)))))} \\ Andrew Howroyd, Oct 20 2020

A057896 Nonnegative numbers that can be written as m^k - m (with m and k nonnegative) in more than one way.

Original entry on oeis.org

0, 6, 30, 210, 240, 2184, 8190, 78120, 24299970
Offset: 1

Views

Author

Henry Bottomley, Sep 26 2000

Keywords

Comments

The next term, if it exists, is at least 2*10^17. - David Wasserman, May 01 2002
a(10) > 10^24, if it exists. The only numbers below 10^24 that can be written as m^k+m in more than one way are 30 = 5^2+5 = 3^3+3 and 130 = 5^3+5 = 2^7+2. - Giovanni Resta, Jun 21 2018
Conjectured to be finite and complete by Bennett (2001).

Examples

			30 is in the sequence since 30 = 2^5 - 2 = 6^2 - 6;
2184 is in the sequence since 2184 = 3^7 - 3 = 13^3 - 13.
		

Crossrefs

Programs

  • Maple
    res:= {0}:
    for k from 3 to 60 do
      for m from 2 while m^k-m < 2^60 do
         x:= m^k-m;
         if assigned(R[x]) or issqr(4*x+1) then res:= res union {x}
         else R[x]:= [m,k]
         fi
      od
    od:
    res; # Robert Israel, Oct 07 2015

Extensions

More terms from Jud McCranie, Oct 01 2000
Offset corrected by Joerg Arndt, Oct 07 2015

A308324 Numbers which can be written in the form m^k - m with m an odd prime and k a positive integer.

Original entry on oeis.org

0, 6, 20, 24, 42, 78, 110, 120, 156, 240, 272, 336, 342, 506, 620, 726, 812, 930, 1320, 1332, 1640, 1806, 2162, 2184, 2394, 2756, 3120, 3422, 3660, 4422, 4896, 4970, 5256, 6162, 6558, 6806, 6840, 7832, 9312
Offset: 1

Views

Author

Craig J. Beisel, May 20 2019

Keywords

Comments

Besides the trivial example a(1)=0, the only known term which has two representations is a(24) = 2184 = 3^7 - 3 = 13^3 - 13. It is conjectured by Bennett to be the only term with this property.

Examples

			a(3) = 5^2 - 5 = 20.
		

Crossrefs

Programs

  • PARI
    x=List([]);lim=10000;forprime(m=3,lim,for(k=1,100,y=(m^k-m);if(y>lim,break,i=setsearch(x,y,1);if(i>0,listinsert(x, y, i)))));for(i=1, #x,print(x[i]));

A308394 Numbers which can be written in the form m^k - m with m prime and k a positive integer.

Original entry on oeis.org

0, 2, 6, 14, 20, 24, 30, 42, 62, 78, 110, 120, 126, 156, 240, 254, 272, 336, 342, 506, 510, 620, 726, 812, 930, 1022, 1320, 1332, 1640, 1806, 2046, 2162, 2184, 2394, 2756, 3120, 3422, 3660, 4094, 4422, 4896, 4970, 5256, 6162, 6558, 6806, 6840, 7832, 8190, 9312
Offset: 1

Views

Author

Craig J. Beisel, May 24 2019

Keywords

Comments

The only known terms which have two representations where m is prime are 6 and 2184. It is conjectured by Bennett these are the only terms with this property.

Examples

			a(9) = 2^6 - 2 = 62.
For the two terms known to have two representations we have a(3) = 6 = 2^3 - 2 = 3^2 - 3 and a(33)= 2184 = 3^7 - 3 = 13^3 - 13.
		

Crossrefs

Subsequences: A000918 (2^n - 2), A036689 (p^2 - p), A058809 (3^n - 3), A178671 (5^n - 5).

Programs

  • Maple
    N:= 10^6; # to get all terms <= N
    P:= select(isprime,[2,seq(i,i=3..floor((1+sqrt(1+4*N))/2),2)]):
    S:= {0,seq(seq(m^k-m,k=2..floor(log[m](N+m))),m=P)}:
    sort(convert(S,list)); # Robert Israel, Aug 11 2019
  • PARI
    x=List([]); lim=10000; forprime(m=2, lim, for(k=1, 100, y=(m^k-m); if(y>lim, break, i=setsearch(x, y, 1); if(i>0, listinsert(x, y, i))))); for(i=1, #x, print(x[i]));
    
  • PARI
    isok(n) = {forprime(p=2, oo, my(keepk = 0); for (k=1, oo, if ((x=p^k - p) == n, return(1)); if (x > n, keepk = k; break);); if (keepk == 2, break););} \\ Michel Marcus, Aug 06 2019

A356639 Number of integer sequences b with b(1) = 1, b(m) > 0 and b(m+1) - b(m) > 0, of length n which transform under the map S into a nonnegative integer sequence. The transform c = S(b) is defined by c(m) = Product_{k=1..m} b(k) / Product_{k=2..m} (b(k) - b(k-1)).

Original entry on oeis.org

1, 1, 3, 17, 155, 2677, 73327, 3578339, 329652351
Offset: 1

Views

Author

Thomas Scheuerle, Aug 19 2022

Keywords

Comments

This sequence can be calculated by a recursive algorithm:
Let B1 be an array of finite length, the "1" denotes that it is the first generation. Let B1' be the reversed version of B1. Let C be the element-wise product C = B1 * B1'. Then B2 is a concatenation of taking each element of B1 and add all divisors of the corresponding element in C. If we start with B1 = {1} then we get this sequence of arrays: B2 = {2}, B3 = {3, 4, 6}, ... . a(n) is the length of the array Bn. In short the length of Bn+1 and so a(n+1) is the sum over A000005(Bn * Bn').
The transform used in the definition of this sequence is its own inverse, so if c = S(b) then b = S(c). The eigensequence is 2^n = S(2^n).
There exist some transformation pairs of infinite sequences in the database:
A026549 <--> A038754; A100071 <--> A001405; A058295 <--> A------;
A111286 <--> A098011; A093968 <--> A205825; A166447 <--> A------;
A079352 <--> A------; A082458 <--> A------; A008233 <--> A264635;
A138278 <--> A------; A006501 <--> A264557; A336496 <--> A------;
A019464 <--> A------; A062112 <--> A------; A171647 <--> A359039;
A279312 <--> A------; A031923 <--> A------.
These transformation pairs are conjectured:
A137326 <--> A------; A066332 <--> A300902; A208147 <--> A308546;
A057895 <--> A------; A349080 <--> A------; A019442 <--> A------;
A349079 <--> A------.
("A------" means not yet in the database.)
Some sequences in the lists above may need offset adjustment to force a beginning with 1,2,... in the transformation.
If we allowed signed rational numbers, further interesting transformation pairs could be observed. For example, 1/n will transform into factorials with alternating sign. 2^(-n) transforms into ones with alternating sign and 1/A000045(n) into A000045 with alternating sign.

Examples

			a(4) = 17. The 17 transformation pairs of length 4 are:
  {1, 2, 3, 4}  = S({1, 2, 6, 24}).
  {1, 2, 3, 5}  = S({1, 2, 6, 15}).
  {1, 2, 3, 6}  = S({1, 2, 6, 12}).
  {1, 2, 3, 9}  = S({1, 2, 6, 9}).
  {1, 2, 3, 12} = S({1, 2, 6, 8}).
  {1, 2, 3, 21} = S({1, 2, 6, 7}).
  {1, 2, 4, 5}  = S({1, 2, 4, 20}).
  {1, 2, 4, 6}  = S({1, 2, 4, 12}).
  {1, 2, 4, 8}  = S({1, 2, 4, 8}).
  {1, 2, 4, 12} = S({1, 2, 4, 6}).
  {1, 2, 4, 20} = S({1, 2, 4, 5}).
  {1, 2, 6, 7}  = S({1, 2, 3, 21}).
  {1, 2, 6, 8}  = S({1, 2, 3, 12}).
  {1, 2, 6, 9}  = S({1, 2, 3, 9}).
  {1, 2, 6, 12} = S({1, 2, 3, 6}).
  {1, 2, 6, 15} = S({1, 2, 3, 5}).
  {1, 2, 6, 24} = S({1, 2, 3, 4}).
b(1) = 1 by definition, b(2) = 1+1 as 1 has only 1 as divisor.
a(3) = A000005(b(2)*b(2)) = 3.
The divisors of b(2) are 1,2,4. So b(3) can be b(2)+1, b(2)+2 and b(2)+4.
a(4) = A000005((b(2)+1)*(b(2)+4)) + A000005((b(2)+2)*(b(2)+2)) + A000005((b(2)+4)*(b(2)+1)) = 17.
		

Crossrefs

Showing 1-5 of 5 results.