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.

A127629 Numbers m such that a divisor, together with its quotient and remainder, are consecutive terms (in that order) in a geometric sequence.

Original entry on oeis.org

9, 28, 34, 58, 65, 75, 110, 126, 132, 201, 205, 217, 224, 246, 254, 258, 294, 344, 384, 399, 436, 498, 502, 513, 516, 520, 579, 657, 680, 690, 730, 786, 810, 866, 880, 978, 979, 1001, 1008, 1028, 1038, 1105, 1128, 1164, 1330, 1332, 1365, 1370, 1374, 1388
Offset: 1

Views

Author

Nick Hobson, Jan 20 2007

Keywords

Comments

The sequence misses the primes.
When m is a term, then m = d*q + r and rBernard Schott, May 15 2020

Examples

			58 is in the sequence because 58 = 9*6 + 4, where 9, 6 and 4 are consecutive terms in a geometric sequence.
For a(4) = 58 with noninteger ratio = 3/2:
     58 | 9          58 | 6
        ------          ------
      4 | 6           4 | 9
For a(16) = 258 with integer ratio = 4:
    258 | 32         258 |  8
        ------           -------
      2 |  8           2 | 32
		

Crossrefs

Disjoint union of A334185 and A334186.
Subsequence: A001093 \ {0, 1, 2} (for remainder = 1).

Programs

  • Mathematica
    mx = 1388; m = Ceiling @ Sqrt[mx]; s={}; Do[r = Select[Divisors[k^2], #Amiram Eldar, Aug 28 2019 *)
  • PARI
    is(n)={for(d=1, n, if((n\d)*(n%d)==d^2, return(1))); return(0)}

A334185 Let m = d*q + r be the Euclidean division of m by d. The terms m of this sequence satisfy that r, q, d are consecutive positive integer terms in a geometric progression with a common integer ratio.

Original entry on oeis.org

9, 28, 34, 65, 75, 110, 126, 132, 205, 217, 246, 258, 294, 344, 399, 436, 502, 513, 520, 579, 657, 680, 730, 810, 866, 978, 979, 1001, 1028, 1128, 1164, 1330, 1332, 1365, 1374, 1582, 1605, 1729, 1736, 1815, 1947, 2004, 2050, 2064, 2196, 2198, 2310, 2329, 2610, 2710
Offset: 1

Views

Author

Bernard Schott, Apr 18 2020

Keywords

Comments

Inspired by the problem 141 of Project Euler (see the link).
If b is the common ratio, then b is an integer >= 2.
So, when b >= 2 and r >= 1, q=r*b, d=r*b^2, then every m = r * (1+r*b^3) is a term, and the division becomes: r*(1+r*b^3) = (r*b^2) * (r*b) + r. The integers (r, r*b, r*b^2) are in geometric progression.
When (r < q < d) is solution with m = d * q + r, then, with d' = q and q' = d, m = d' * q' + r and (r < d' < q') is also a solution with another order between remainder, divisor and quotient (see last example).
m is a term if m = r * (1+r*b^3) with r >= 1 and b >= 2; so, when r = 1, A001093(n) for n > 1 are terms (see 1st example).

Examples

			a(2) = 28 = 9*3 + 1 with (1,3,9) and ratio = 3;
a(5) = 75 = 12*6 + 3 with (3,6,12) and ratio = 2;
a(12) = 258 = 32*8 + 2 with (2,8,32) and ratio = 4;
a(42) = 2004 = 100*20 + 4 with (r=4, q=20, d=100) but also 2004 = 20*100 + 4 with (r=4, d'=20, q'=100) both with ratio = 5:
2004 | 100             2004 |  20
     +-----                 +-----
   4 |  20                4 | 100
		

Crossrefs

Cf. A334186 (similar, with b is an irreducible fraction).
Subsequence: A001093 \ {0, 1, 2} (for r = 1).

Programs

  • Mathematica
    Select[Range[2000], Length @ Solve[r * (1 + r*b^3) == # && r >=1 && b >= 2, {r, b}, Integers] > 0 &] (* Amiram Eldar, Apr 18 2020 *)
  • PARI
    isok(m) = {for (d=1, m, if (m % d, q = m\d; r = m % d; if (!(d % q) && (d/q == q/r), return (1));););} \\ Michel Marcus, Apr 19 2020

Extensions

Name improved by Michel Marcus, Apr 19 2020
More terms from Michel Marcus, Apr 19 2020

A335064 Let m = d*q + r be the Euclidean division of m by d. The terms m of this sequence satisfy that q, r, d are consecutive positive integer terms in a geometric progression with a noninteger common ratio > 1.

Original entry on oeis.org

42, 110, 156, 210, 240, 342, 420, 462, 506, 600, 702, 812, 930, 1122, 1190, 1260, 1332, 1482, 1560, 1640, 1806, 1980, 2070, 2162, 2352, 2550, 2652, 2756, 2970, 3080, 3192, 3306, 3422, 3660, 3906, 4032, 4290, 4422, 4692, 4830, 4970, 5256, 5550, 5700, 5852, 6006, 6162
Offset: 1

Views

Author

Bernard Schott, May 22 2020

Keywords

Comments

Inspired by the problem 141 of Project Euler (see the link).
The terms of this sequence are oblong numbers m = k*(k+1) with k in A024619.
When q < r < d are consecutive terms of a geometric progression of constant b = p/s noninteger, with b>1, s>=2, p>s, it is necessary that q is a multiple of s^2, so q = q' * s^2 with q' >= 1; the Euclidean division of a term m by q becomes
p*s*q' * (1+p*s*q') = (p^2*q') * (s^2*q') + p*s*q' with k = p*s*q',
so (q, r, d) = (s^2*q', p*s*q', p^2*q') is solution. (see examples).
But, as these terms are oblong, there exists also another division where the constant ratio is the integer psq' and (q,r,d) = (1, p*s*q', (p*s*q')^2) are in geometric progression.

Examples

			Examples for 42, 110 and 156 with consecutive ratios 3/2, 5/2, 4/3:
   42 | 9         110 | 25         156 | 16
      -----           -----            -----
    6 | 4    ,     10 |  4     ,    12 |  9 ,
then with consecutive ratios 2, 10 and 12:
   42 | 12        110 | 100        156 | 144
      -----           -----            ------
    6 |  3   ,     10 |   1    ,    12 |   1.
		

Crossrefs

Subsequence of A002378 and of A335065.

Programs

  • Mathematica
    Table[n*(n + 1), {n, Select[Range[80], PrimeNu[#] > 1 &]}] (* Amiram Eldar, May 23 2020 *)
  • PARI
    apply(x->x*(x+1), select(x->!isprimepower(x), [2..80])) \\ Michel Marcus, May 23 2020

Formula

a(n) = A024619(n) * (1+A024619(n)).
a(n) = A002378(A024619(n)). - Michel Marcus, May 23 2020

A335065 Let m = d*q + r be the Euclidean division of m by d. The terms m of this sequence satisfy that d, q, r are consecutive positive integer terms in a geometric progression but not necessarily in that order.

Original entry on oeis.org

6, 9, 12, 20, 28, 30, 34, 42, 56, 58, 65, 72, 75, 90, 110, 126, 132, 156, 182, 201, 205, 210, 217, 224, 240, 246, 254, 258, 272, 294, 306, 342, 344, 380, 384, 399, 420, 436, 462, 498, 502, 506, 513, 516, 520, 552, 579, 600, 650, 657, 680, 690, 702, 730, 756, 786
Offset: 1

Views

Author

Bernard Schott, May 23 2020

Keywords

Comments

Inspired by the problem 141 of Project Euler (see link).
There exist 3 possibilities to get such terms m that satisfy that d, q, r are consecutive positive integer terms in a geometric progression but not necessarily in that order:
-> the geometric progression is r < q < d (A127629).
-> the geometric progression is r < d < q (same terms of A127629).
-> the geometric progression is q < r < d (A002378 \ {0,2} = oblong numbers >= 6).
Some numbers have a geometric progression solution in the 3 cases (132, 1332, 6162, ...) [see examples].

Examples

			Examples with r < q < d, r < d < q, q < r <d:
   34 | 8        75 |  6           42 | 12
      ----          -----             -----
    2 | 4     ,   3 | 12     ,      6 |  3
The 3 possible divisions by 132:
  132 | 16      132 |  8          132 | 121
      -----         ------            ------
    4 |  8   ,    4 | 16     ,     11 |   1.
		

Crossrefs

Equals A127629 Union A002378 \ {0,2}.
Subsequences: A334185, A334186, A335064.

Programs

  • Mathematica
    mx = 800; Union@ Reap[ Do[y = x+1; While[(z = y^2/x) < mx, If[ IntegerQ@ z, If[(m = z y + x) <= mx, Sow@ m]; If[(m = z x + y) <= mx, Sow@ m]]; y++], {x, mx}]][[2, 1]] (* Giovanni Resta, May 24 2020 *)
  • PARI
    isok(n) = {my(r, d); for (q=2, n-1, if (r=(n % q), d = n\q; if ((r*d == q^2) || (r*q == d^2) || (q*d == r^2), return (1));););} \\ Michel Marcus, May 25 2020

A335272 For m to be a term there must exist three Euclidean divisions of m by d, d', and d", m = d*q + r = d'*q' + r' = d"*q" + r", such that (r, q, d), (r', d', q'), and (q", r", d") are three geometric progressions.

Original entry on oeis.org

110, 132, 1332, 6162, 10712, 12210, 35156, 60762, 67340, 152490, 296480, 352242, 354620, 357006, 648830, 771762, 932190, 1197930, 2057790, 2803950, 3241800, 3310580, 4458432, 6454140, 7865220, 9613100, 10814232, 13976382, 16382256, 19267710, 53824232, 55138050
Offset: 1

Views

Author

Bernard Schott, May 30 2020

Keywords

Comments

Inspired by Project Euler, Problem 141 (see link).
The terms are necessary oblong numbers >= 6.

Examples

			For 110:
   110 | 18             110 |  6               110 | 100
       -----               ------               ---------
     2 |  6       ,       2 | 18       ,        10 |   1
For 132, see A335065.
For 1332:
  1332 | 121           1332 |  11            1332 | 1296
       ------              -------                -------
     1 |  11      ,       1 | 121      ,       36 |   1   .
		

Crossrefs

Intersection of A127629 and A002378.

Programs

  • Mathematica
    Select[(#^2 + #) & /@ Range[2000], (n = #; AnyTrue[ Range[1 + Sqrt@ n], #^2 == Mod[n, #] Floor[n/#] &]) &] (* Giovanni Resta, Jun 03 2020 *)
  • PARI
    isob(n) = my(m=sqrtint(n)); m*(m+1)==n; \\ A002378
    isgd(n) = {for(d=1, n, if((n\d)*(n%d)==d^2, return(1))); return(0)}; \\ A127629
    isok(n) = isob(n) && isgd(n); \\ Michel Marcus, May 30 2020

Extensions

More terms from Michel Marcus, May 30 2020
a(18)-a(26) from Jinyuan Wang, May 30 2020
Terms a(27) and beyond from Giovanni Resta, Jun 03 2020
Showing 1-5 of 5 results.