A090570 -id:A090570 - cp's OEIS frontend

A218864 Numbers of the form 9k^2 + 8k, k an integer.

0, 1, 17, 20, 52, 57, 105, 112, 176, 185, 265, 276, 372, 385, 497, 512, 640, 657, 801, 820, 980, 1001, 1177, 1200, 1392, 1417, 1625, 1652, 1876, 1905, 2145, 2176, 2432, 2465, 2737, 2772, 3060, 3097, 3401, 3440, 3760, 3801, 4137, 4180, 4532, 4577, 4945, 4992

Offset: 1

Jason Kimberley, Nov 08 2012

Numbers m such that 9*m + 16 is a square. - Vincenzo Librandi, Apr 07 2013
Equivalently, integers of the form h*(h + 8)/9 (nonnegative values of h are listed in A090570). - Bruno Berselli, Jul 15 2016
Generalized 20-gonal (or icosagonal) numbers: r*(9*r - 8) with r = 0, +1, -1, +2, -2, +3, -3, ... - Omar E. Pol, Jun 06 2018
Partial sums of A317316. - Omar E. Pol, Jul 28 2018
Exponents in expansion of Product_{n >= 1} (1 + x^(18*n-17))*(1 + x^(18*n-1))*(1 - x^(18*n)) = 1 + x + x^17 + x^20 + x^52 + .... - Peter Bala, Dec 10 2020

Jason Kimberley, Table of n, a(n) for n = 1..2000
S. Cooper and M. D. Hirschhorn, Results of Hurwitz type for three squares. Discrete Math., Vol. 274, No. 1-3 (2004), pp. 9-24. See C(q).
Index entries for linear recurrences with constant coefficients, signature (1,2,-2,-1,1).

Characteristic function is A205987.
Numbers of the form 9*m^2+k*m, for integer n: A016766 (k=0), A132355 (k=2), A185039 (k=4), A057780 (k=6), this sequence (k=8).
Cf. A074377 (numbers m such that 16*m+9 is a square).
Cf. A317316.
For similar sequences of numbers m such that 9*m+i is a square, see list in A266956.
Cf. sequences of the form m*(m+i)/(i+1) listed in A274978. [Bruno Berselli, Jul 25 2016]
Sequences of generalized k-gonal numbers: A001318 (k=5), A000217 (k=6), A085787 (k=7), A001082 (k=8), A118277 (k=9), A074377 (k=10), A195160 (k=11), A195162 (k=12), A195313 (k=13), A195818 (k=14), A277082 (k=15), A274978 (k=16), A303305 (k=17), A274979 (k=18), A303813 (k=19), this sequence (k=20), A303298 (k=21), A303299 (k=22), A303303 (k=23), A303814 (k=24), A303304 (k=25), A316724 (k=26), A316725 (k=27), A303812 (k=28), A303815 (k=29), A316729 (k=30).

a:=func; [0]cat[a(n*m): m in [-1,1], n in [1..20]];

Array[(18 # (# - 1) - 7 (-1)^#*(2 # - 1) - 7)/8 &, 48] (* or *)
CoefficientList[Series[x (1 + 16 x + x^2)/((1 + x)^2*(1 - x)^3), {x, 0, 47}], x] (* Michael De Vlieger, Jun 06 2018 *)

a(n) = (18*n*(n - 1) - 7*(-1)^n*(2*n - 1) - 7)/8. - Bruno Berselli, Nov 13 2012
G.f.: x*(1 + 16*x + x^2)/((1 + x)^2*(1 - x)^3). - Bruno Berselli, Nov 14 2012
Sum_{n>=2} 1/a(n) = (9 + 8*Pi*cot(Pi/9))/64. - Amiram Eldar, Feb 28 2022

A058369 Numbers k such that k and k^2 have same digit sum.

Original entry on oeis.org

0, 1, 9, 10, 18, 19, 45, 46, 55, 90, 99, 100, 145, 180, 189, 190, 198, 199, 289, 351, 361, 369, 379, 388, 450, 451, 459, 460, 468, 495, 496, 550, 558, 559, 568, 585, 595, 639, 729, 739, 775, 838, 855, 900, 954, 955, 990, 999, 1000, 1098, 1099, 1179, 1188, 1189

Offset: 1

G. L. Honaker, Jr., Dec 17 2000

It is interesting that the graph of this sequence appears almost identical as the maximum value of n increases by factors of 10. Compare the graph of the b-file (having numbers up to 10^6) with the plot of the terms up to 10^8. - T. D. Noe, Apr 28 2012
If iterated digit sum (A010888, A056992) is used instead of just digit sum (A007953, A004159), we get A090570 of which this sequence is a subset. - Jeppe Stig Nielsen, Feb 18 2015
Hare, Laishram, & Stoll show that this sequence (indeed, even its subsequence A254066) is infinite. In particular for each k in {846, 847, 855, 856, 864, 865, 873, ...} there are infinitely many terms in this sequence not divisible by 10 that have digit sum k. - Charles R Greathouse IV, Aug 25 2015
There are infinitely many n such that both n and n+1 are in the sequence.  This includes A002283. - Robert Israel, Aug 26 2015

			Digit sum of 9 = 9 9^2 = 81, 8+1 = 9 digit sum of 145 = 1+4+5 = 10 145^2 = 21025, 2+1+0+2+5 = 10 digit sum of 954 = 9+5+4 = 18 954^2 = 910116, 9+1+0+1+1+6 = 18. - Florian Roeseler (hazz_dollazz(AT)web.de), May 03 2010

Zak Seidov, Table of n, a(n) for n = 1..8354 (to a(n) = 10^6)
Code Golf StackExchange, Equality in the sum of digits, coding challenge started Mar 11 2016.
K. G. Hare, S. Laishram, and T. Stoll, The sum of digits of n and n^2, International Journal of Number Theory 7:7 (2011), pp. 1737-1752.
T. D. Noe, Plot of terms up to 10^8

Cf. A002283, A007953, A058370, A058852, A077436, A254066.
Cf. A147523 (number of numbers in each decade).
Subsequence of A090570.

import Data.List (elemIndices)
import Data.Function (on)
a058369 n = a058369_list !! (n-1)
a058369_list =
   elemIndices 0 $ zipWith ((-) `on` a007953) [0..] a000290_list
-- Reinhard Zumkeller, Aug 17 2011

[n: n in [0..1200] |(&+Intseq(n)) eq (&+Intseq(n^2))]; // Vincenzo Librandi, Aug 26 2015

sd := proc (n) options operator, arrow: add(convert(n, base, 10)[j], j = 1 .. nops(convert(n, base, 10))) end proc: a := proc (n) if sd(n) = sd(n^2) then n else end if end proc; seq(a(n), n = 0 .. 1400); # Emeric Deutsch, May 11 2010
select(t -> convert(convert(t,base,10),`+`)=convert(convert(t^2,base,10),`+`),
[seq(seq(9*i+j,j=0..1),i=0..1000)]); # Robert Israel, Aug 26 2015

Select[Range[0,1200],Total[IntegerDigits[#]]==Total[IntegerDigits[ #^2]]&] (* Harvey P. Dale, Jun 14 2011 *)

is(n)=sumdigits(n)==sumdigits(n^2) \\ Charles R Greathouse IV, Aug 25 2015

def ds(n): return sum(map(int, str(n)))
def ok(n): return ds(n) == ds(n**2)
def aupto(nn): return [m for m in range(nn+1) if ok(m)]
print(aupto(1189)) # Michael S. Branicky, Jan 09 2021

A007953(a(n)) = A004159(a(n)). - Reinhard Zumkeller, Apr 25 2009

Edited by N. J. A. Sloane, May 30 2010

A145389 Digital roots of triangular numbers.

Original entry on oeis.org

0, 1, 3, 6, 1, 6, 3, 1, 9, 9, 1, 3, 6, 1, 6, 3, 1, 9, 9, 1, 3, 6, 1, 6, 3, 1, 9, 9, 1, 3, 6, 1, 6, 3, 1, 9, 9, 1, 3, 6, 1, 6, 3, 1, 9, 9, 1, 3, 6, 1, 6, 3, 1, 9, 9, 1, 3, 6, 1, 6, 3, 1, 9, 9, 1, 3, 6, 1, 6, 3, 1, 9, 9, 1, 3, 6, 1, 6, 3, 1, 9, 9, 1, 3, 6, 1, 6

Offset: 0

Reinhard Zumkeller, Oct 10 2008

Decimal expansion of 45387733/3333333330. - Enrique Pérez Herrero, Nov 14 2021

Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,0,0,0,0,1).

Cf. A000217, A004157, A010888.

digitalRoot[n_Integer?Positive] := FixedPoint[Plus@@IntegerDigits[#]&,n]; Table[If[n==0,0,digitalRoot[n(n+1)/2]], {n,0,100}] (* Vladimir Joseph Stephan Orlovsky, May 02 2011 *)

a(n)=if(n, n=n*(n+1)/2%9; if(n, n, 9), 0) \\ Charles R Greathouse IV, Dec 19 2016

def A145389(n): return (9, 1, 3, 6, 1, 6, 3, 1, 9)[n%9] if n else 0 # Chai Wah Wu, Feb 09 2023

a(n) = A010888(A000217(n)).
Periodic sequence for n>0: a(n+9) = a(n);
a(A016777(n)) = 1; a(A007494(n)) <> 1;
a(A090570(n)) = A010888(A090570(n)).
a(n) = 1 + ((n^2 + n - 2)/2) mod 9. - Ant King, Apr 25 2009
G.f.: x(1 + 3x + 6x^2 + x^3 + 6x^4 + 3x^5 + x^6 + 9x^7 + 9x^8)/((1-x)(1 + x + x^2)(1 + x^3 + x^6)). - Ant King, Nov 16 2010

A274406 Numbers m such that 9 divides m*(m + 1).

Original entry on oeis.org

0, 8, 9, 17, 18, 26, 27, 35, 36, 44, 45, 53, 54, 62, 63, 71, 72, 80, 81, 89, 90, 98, 99, 107, 108, 116, 117, 125, 126, 134, 135, 143, 144, 152, 153, 161, 162, 170, 171, 179, 180, 188, 189, 197, 198, 206, 207, 215, 216, 224, 225, 233, 234, 242, 243, 251, 252, 260, 261, 269

Offset: 1

Bruno Berselli, Jun 20 2016

Equivalently, numbers congruent to 0 or 8 mod 9.

Terms of A007494 with indices in A047264. Also, terms of A060464 with indices in A047335.

Bruno Berselli, Table of n, a(n) for n = 1..1000
Index entries for linear recurrences with constant coefficients, signature (1,1,-1).

Cf. A008591 (first bisection), A010689 (first differences), A017257 (second bisection).
Cf. similar sequences in which m*(m+1) is divisible by k: A014601 (k=4), A047208 (k=5), A007494 (k=3 and 6), A047335 (k=7), A047521 (k=8), this sequence (k=9).
Cf. A301451: numbers congruent to {1, 7} mod 9; A193910: numbers congruent to {2, 6} mod 9.

[n: n in [0..300] | IsDivisibleBy(n*(n+1),9)];

Select[Range[0, 300], Divisible[# (# + 1), 9] &]

for(n=0, 300, if(n*(n+1)%9==0, print1(n", ")))

[n for n in range(300) if 9.divides(n*(n+1))]

G.f.: x^2*(8 + x)/((1 + x)*(1 - x)^2).
a(n) = (18*n + 7*(-1)^n - 11)/4. Therefore: a(2*m) = 9*m-1, a(2*m+1) = 9*m. It follows that a(j)+a(k) and a(j)*a(k) belong to the sequence if j and k are not both even.
a(n) = -A090570(-n+2).
a(n) = a(n-1) + a(n-2) - a(n-3).
a(2*r+1) + a(2*r+s+1) = a(4*r+s+1) and a(2*r) + a(2*r+2*s+1) = a(4*r+2*s). A particular case provided by these identities: a(n) = a(n - 2*floor(n/6)) + a(2*floor(n/6) + 1).
E.g.f.: 1 + ((9*x - 2)*cosh(x) + 9*(x - 1)*sinh(x))/2. - Stefano Spezia, Apr 24 2021

A254066 Primitive numbers n such that the sums of the digits of n and n^2 coincide.

Original entry on oeis.org

1, 9, 18, 19, 45, 46, 55, 99, 145, 189, 198, 199, 289, 351, 361, 369, 379, 388, 451, 459, 468, 495, 496, 558, 559, 568, 585, 595, 639, 729, 739, 775, 838, 855, 954, 955, 999, 1098, 1099, 1179, 1188, 1189, 1198, 1269, 1468, 1485, 1494, 1495, 1585, 1738, 1747

Offset: 1

Nikhil Mahajan, Jan 25 2015

Members of A058369 not congruent to 0 (mod 10).
This sequence is to A058369 what A114135 is to A111434.
Hare, Laishram, & Stoll show that this sequence is infinite. In particular for each k in {846, 847, 855, 856, 864, 865, 873, ...} there are infinitely many terms in this sequence with digit sum k. - Charles R Greathouse IV, Aug 25 2015

			9 is in the sequence because the digit sum of 9^2 = 81 is 9.
18 is in the sequence because the digit sum of 18^2 = 324 is 9, same as the digit sum of 18.

Nikhil Mahajan, Table of n, a(n) for n = 1..10000
K. G. Hare, S. Laishram, and T. Stoll, The sum of digits of n and n^2, International Journal of Number Theory 7:7 (2011), pp. 1737-1752.

Cf. A058369, A114135, A111434.

Subsequence of A090570.

[n: n in [1..1000] | &+Intseq(n) eq &+Intseq(n^2) and not IsZero(n mod 10)]; // Bruno Berselli, Jan 29 2015

Select[Range[1000],!Divisible[#,10]&&Total[IntegerDigits[#]] == Total[ IntegerDigits[#^2]]&] (* Harvey P. Dale, Dec 27 2015 *)

is(n)=sumdigits(n)==sumdigits(n^2) \\ Charles R Greathouse IV, Aug 25 2015

list(lim)=my(v=List()); forstep(n=1,lim,[8, 9, 1, 8, 1, 8, 1, 8, 1, 8, 1, 8, 1, 8, 1, 8, 1, 9], if(sumdigits(n)==sumdigits(n^2), listput(v,n))); Vec(v) \\ Charles R Greathouse IV, Aug 26 2015

[n for n in [0..1000] if sum(n.digits())==sum((n^2).digits()) and n%10!=0] # Tom Edgar, Jan 27 2015

More terms from Harvey P. Dale, Dec 27 2015

A151972 Numbers that are congruent to {0, 1, 6, 10} mod 15.

Original entry on oeis.org

0, 1, 6, 10, 15, 16, 21, 25, 30, 31, 36, 40, 45, 46, 51, 55, 60, 61, 66, 70, 75, 76, 81, 85, 90, 91, 96, 100, 105, 106, 111, 115, 120, 121, 126, 130, 135, 136, 141, 145, 150, 151, 156, 160, 165, 166, 171, 175, 180, 181, 186, 190, 195, 196, 201, 205, 210, 211, 216, 220, 225

Offset: 1

N. J. A. Sloane, Aug 23 2009

Also, numbers n such that n^2 - n is divisible by 15.

Also, numbers n such that n^2 - n is divisible by 30.

Index entries for linear recurrences with constant coefficients, signature (1,0,0,1,-1).

For m^2 == m (mod n), see: n=2: A001477, n=3: A032766, n=4: A042948, n=5: A008851, n=6: A032766, n=7: A047274, n=8: A047393, n=9: A090570, n=10: A008851, n=11: A112651, n=12: A112652, n=13: A112653, n=14: A047274, n=15: A151972, n=16: A151977, n=17: A151978, n=18: A090570, n=19: A151979, n=20: A151980, n=21: A151971, n=22, A112651, n=24: A151973, n=26: A112653, n=30: A151972, n=32: A151983, n=34: A151978, n=38: A151979, n=42: A151971, n=48: A151981, n=64: A151984.

Cf. A215202.

[ n : n in [0..1000] | n mod 15 in [0, 1, 6, 10]]; // Vincenzo Librandi, Apr 02 2011, simplified by Eric M. Schmidt, Aug 05 2012

[ n: n in [0..1000] | (n^2-n) mod (15) eq 0 ]; // Vincenzo Librandi, Apr 03 2011, altered by Eric M. Schmidt, Aug 05 2012

A151972:=n->(30*n-41-5*I^(2*n)+(3+3*I)*I^(-n)+(3-3*I)*I^n)/8: seq(A151972(n), n=1..100); # Wesley Ivan Hurt, Jun 07 2016

Select[Range[0, 300], Divisible[#^2-#, 15]&]  (* Harvey P. Dale, Apr 01 2011, altered by Eric M. Schmidt, Aug 05 2012 *)

G.f.: x^2*(1+5*x+4*x^2+5*x^3) / ( (1+x)*(1+x^2)*(x-1)^2 ). - R. J. Mathar, Oct 25 2011
From Wesley Ivan Hurt, Jun 07 2016: (Start)
a(n) = (30*n-41-5*i^(2*n)+(3+3*i)*i^(-n)+(3-3*i)*i^n)/8 where i=sqrt(-1).
a(n) = a(n-1) + a(n-4) - a(n-5) for n>5. (End)
E.g.f.: (20 + (15*x - 23)*cosh(x) + 3*(sin(x) + cos(x) + (5*x - 6)*sinh(x)))/4. - Ilya Gutkovskiy, Jun 07 2016

This is a merge of two identical sequences, A151972 and A151975.

A151971 Numbers n such that n^2 - n is divisible by 21.

Original entry on oeis.org

0, 1, 7, 15, 21, 22, 28, 36, 42, 43, 49, 57, 63, 64, 70, 78, 84, 85, 91, 99, 105, 106, 112, 120, 126, 127, 133, 141, 147, 148, 154, 162, 168, 169, 175, 183, 189, 190, 196, 204, 210, 211, 217, 225, 231, 232, 238, 246, 252, 253, 259, 267, 273, 274, 280, 288, 294, 295, 301, 309

Offset: 1

N. J. A. Sloane, Aug 23 2009

Equivalently, numbers that are congruent to {0, 1, 7, 15} mod 21. - Bruno Berselli, Aug 06 2012

Vincenzo Librandi, Table of n, a(n) for n = 1..1000
Index entries for linear recurrences with constant coefficients, signature (1,0,0,1,-1).

For m^2 == m (mod n), see: n=2: A001477; n=3: A032766; n=4: A042948; n=5: A008851; n=6: A032766; n=7: A047274; n=8: A047393; n=9: A090570; n=10: A008851; n=11: A112651; n=12: A112652; n=13:A112653; n=14: A047274; n=15: A151972; n=16: A151977; n=17: A151978; n=18: A090570; n=19: A151979; n=20: A151980; n=21: A151971; n=22: A112651; n=24: A151973; n=26: A112653; n=30: A151972; n=32: A151983; n=34: A151978; n=38: A151979; n=42: A151971; n=48: A151981; n=64: A151984.

Cf. A215202.

[n: n in [0..309] | IsZero((n^2-n) mod 21)]; // Bruno Berselli, Aug 06 2012

A151971:=n->(42*n+14*I^((n-1)*n)-3*I^(2*n)-3)/8-7: seq(A151971(n), n=1..100); # Wesley Ivan Hurt, Jun 07 2016

Select[Range[0,400], Divisible[#^2-#,21]&] (* Harvey P. Dale, Jun 04 2012 *)

makelist((42*n+14*%i^((n-1)*n)-3*(-1)^n-3)/8-7, n, 1, 60); /* Bruno Berselli, Aug 06 2012 */

From Bruno Berselli, Aug 06 2012: (Start)
G.f.: x^2*(1+6*x+8*x^2+6*x^3)/((1+x)*(1-x)^2*(1+x^2)).
a(n) = (42*n +14*i^((n-1)*n) -3*(-1)^n -3)/8 -7, where i=sqrt(-1). (End)
a(n) = a(n-1) + a(n-4) - a(n-5) for n>5. - Wesley Ivan Hurt, Jun 07 2016
E.g.f.: (24 + (21*x - 31)*cosh(x) + 7*(sin(x) + cos(x) + (3*x - 4)*sinh(x)))/4. - Ilya Gutkovskiy, Jun 07 2016

A215202 Irregular triangle in which n-th row gives m in 1, ..., n-1 such that m^2 == m (mod n).

Original entry on oeis.org

1, 1, 1, 1, 1, 3, 4, 1, 1, 1, 1, 5, 6, 1, 1, 4, 9, 1, 1, 7, 8, 1, 6, 10, 1, 1, 1, 9, 10, 1, 1, 5, 16, 1, 7, 15, 1, 11, 12, 1, 1, 9, 16, 1, 1, 13, 14, 1, 1, 8, 21, 1, 1, 6, 10, 15, 16, 21, 25, 1, 1, 1, 12, 22, 1, 17, 18, 1, 15, 21, 1, 9, 28, 1, 1, 19, 20, 1, 13

Offset: 2

Eric M. Schmidt, Aug 05 2012

The n-th row has length A034444(n) - 1.
If m appears in row n, then gcd(n,m) appears in the n-th row of A077610. Moreover, if m', distinct from m, also appears in row n, then gcd(n, m) does not equal gcd(n, m').
For odd n and any integer m, m^2 == m (mod n) iff m^2 == m (mod 2n).
Let P(1)={1} and for integers x > 1, let P(x) be the set of distinct prime divisors of x. We can define an equivalence relation ~ on the set of elements in the ring (Z_n, +mod n,*mod n): for all a,b in Z_n (where a,b are the least nonnegative residues modulo n) a ~ b iff P(gcd(a,n)) intersect P(n) is equal to P(gcd(b,n)) intersect P(n). If we include 0 in each row then these elements can represent the equivalence classes. They form a commutative monoid. - Geoffrey Critzer, Feb 13 2016

			Triangle begins:
1;
1;
1;
1;
1, 3, 4;
1;
1;
1;
1, 5, 6;
1;
1, 4, 9;
1;
1, 7, 8;
1, 6, 10;
1;
1;
1, 9, 10; etc.  - _Bruno Berselli_, Aug 06 2012

Eric M. Schmidt, Rows 2..2000, flattened

For m^2 == m (mod n), see: n=2: A001477; n=3: A032766; n=4: A042948; n=5: A008851; n=6: A032766; n=7: A047274; n=8: A047393; n=9: A090570; n=10: A008851; n=11: A112651; n=12: A112652; n=13: A112653; n=14: A047274; n=15: A151972; n=16: A151977; n=17: A151978; n=18: A090570; n=19: A151979; n=20: A151980; n=21: A151971; n=22: A112651; n=24: A151973; n=26: A112653; n=30: A151972; n=32: A151983; n=34: A151978; n=38: A151979; n=42: A151971; n=48: A151981; n=64: A151984; n=100: A008852; n=1000: A008853.

[m: m in [1..n-1], n in [2..40] | m^2 mod n eq m]; // Bruno Berselli, Aug 06 2012

Table[Select[Range[n], Mod[#^2, n] == # &], {n, 2, 30}] // Grid (* Geoffrey Critzer, May 26 2015 *)

def A215202(n) : return [m for m in range(1, n) if m^2 % n == m];

cp's OEIS Frontend

A218864 Numbers of the form 9*k^2 + 8*k, k an integer.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

Formula

A058369 Numbers k such that k and k^2 have same digit sum.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Python

Formula

Extensions

A145389 Digital roots of triangular numbers.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A274406 Numbers m such that 9 divides m*(m + 1).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Sage

Formula

A254066 Primitive numbers n such that the sums of the digits of n and n^2 coincide.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

Sage

Extensions

A151972 Numbers that are congruent to {0, 1, 6, 10} mod 15.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Magma

Maple

A218864 Numbers of the form 9k^2 + 8k, k an integer.