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

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

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