A102928 -id:A102928 - cp's OEIS frontend

A132767 a(n) = n*(n + 25).

0, 26, 54, 84, 116, 150, 186, 224, 264, 306, 350, 396, 444, 494, 546, 600, 656, 714, 774, 836, 900, 966, 1034, 1104, 1176, 1250, 1326, 1404, 1484, 1566, 1650, 1736, 1824, 1914, 2006, 2100, 2196, 2294, 2394, 2496, 2600, 2706, 2814, 2924, 3036, 3150, 3266, 3384

Offset: 0

Omar E. Pol, Aug 28 2007

a(n) is the Zagreb 1 index of the Mycielskian of the cycle graph C[n]. See p. 205 of the D. B. West reference. - Emeric Deutsch, Nov 04 2016

Douglas B. West, Introduction to Graph Theory, 2nd ed., Prentice-Hall, NJ, 2001.

G. C. Greubel, Table of n, a(n) for n = 0..5000
Felix P. Muga II, Extending the Golden Ratio and the Binet-de Moivre Formula, Preprint on ResearchGate, March 2014.
Wikipedia, Mycielskian.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A001008, A002378, A005563, A028347, A028552, A028557, A028560, A028563, A028566, A028569, A098603, A098847, A098848, A098849, A098850, A102928, A120071, A132759, A132760, A132761, A132762, A132763, A132764, A132765, A132766.

Table[n (n + 25), {n, 0, 50}] (* Bruno Berselli, Aug 22 2018 *)
LinearRecurrence[{3,-3,1},{0,26,54},60] (* Harvey P. Dale, Feb 20 2023 *)

a(n)=n*(n+25) \\ Charles R Greathouse IV, Jun 17 2017

[n*(n+25) for n in (0..50)] # G. C. Greubel, Mar 13 2022

a(n) = 2*n + a(n-1) + 24 (with a(0)=0). - Vincenzo Librandi, Aug 03 2010
a(n) = n^2 + 25*n. - Omar E. Pol, Nov 04 2016
From Chai Wah Wu, Dec 17 2016: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n > 2.
G.f.: 2*x*(13 - 12*x)/(1-x)^3. (End)
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(25)/25 = A001008(25)/A102928(25) = 34052522467/223092870000, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*log(2)/25 - 19081066231/669278610000. (End)
E.g.f.: x*(26 + x)*exp(x). - G. C. Greubel, Mar 13 2022

A175441 Denominator of the harmonic mean of the first n positive integers.

Original entry on oeis.org

1, 3, 11, 25, 137, 49, 363, 761, 7129, 7381, 83711, 86021, 1145993, 1171733, 1195757, 2436559, 42142223, 14274301, 275295799, 11167027, 18858053, 19093197, 444316699, 1347822955, 34052522467, 34395742267, 312536252003, 315404588903, 9227046511387, 9304682830147

Offset: 1

Jaroslav Krizek, May 16 2010

See A102928 - numerators of the harmonic means of the first n positive integers.
a(n) = A001008(n) for n = 1 - 19 and other n.
a(n) is also the numerator of H(n)/(n+1)+1/(n+1)^2 = -int(x^n*log(1-x), x=0..1) with H(n) = A001008(x)/A002805(n) harmonic number of order n. - Groux Roland, Jan 08 2011
a(n) coincides with A001008(n) iff n is not in the sequence A256102. For the quotient A001008(n) / a(n) if n is from A256102 see the corresponding entry of A256103. - Wolfdieter Lang, Apr 23 2015

			H(n) = 1, 4/3, 18/11, 48/25, 300/137, 120/49, 980/363, 2240/761, ...
Comparison with A001008: the first 19 entries coincide because 20 is the first entry of A256102; indeed, A001008(20) = 55835135 and a(2) = 11167027. The quotient is 5 = A256103(1). - _Wolfdieter Lang_, Apr 23 2015

Andrew Howroyd, Table of n, a(n) for n = 1..200

Cf. A102928 (numerators), A001008, A256102, A256103.

Table[Denominator[HarmonicMean[Range[n]]],{n,30}] (* Harvey P. Dale, May 21 2021 *)

a(n)={denominator(n/sum(k=1, n, 1/k))} \\ Andrew Howroyd, Jan 08 2020

a(n) = denominator(n/(Sum_{k=1..n} 1/k)). - Andrew Howroyd, Jan 08 2020

a(n) = numerator(Sum_{k>0} 1/(k*(k+n))). - Mohammed Yaseen, Jun 23 2024

Terms a(25) and beyond from Andrew Howroyd, Jan 08 2020

A132768 a(n) = n*(n + 26).

Original entry on oeis.org

0, 27, 56, 87, 120, 155, 192, 231, 272, 315, 360, 407, 456, 507, 560, 615, 672, 731, 792, 855, 920, 987, 1056, 1127, 1200, 1275, 1352, 1431, 1512, 1595, 1680, 1767, 1856, 1947, 2040, 2135, 2232, 2331, 2432, 2535, 2640, 2747, 2856, 2967, 3080, 3195, 3312, 3431

Offset: 0

Omar E. Pol, Aug 28 2007

G. C. Greubel, Table of n, a(n) for n = 0..5000
Felix P. Muga II, Extending the Golden Ratio and the Binet-de Moivre Formula, Preprint on ResearchGate, March 2014.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A001008, A002378, A005563, A028347, A028552, A028557, A028560, A028563, A028566, A028569.
Cf. A098849, A098850, A098603, A098847, A098848, A102928, A120071, A132759, A132760, A132761.
Cf. A132762, A132763, A132764, A132765, A132766, A132767.

Table[n(n+26),{n,0,50}] (* or *) LinearRecurrence[{3,-3,1},{0,27,56},50] (* Harvey P. Dale, Dec 15 2018 *)

a(n)=n*(n+26) \\ Charles R Greathouse IV, Jun 17 2017

[n*(n+26) for n in (0..50)] # G. C. Greubel, Mar 13 2022

a(n) = n*(n + 26).
a(n) = 2*n + a(n-1) + 25, with a(0)=0. - Vincenzo Librandi, Aug 03 2010
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(26)/26 = A001008(26)/A102928(26) = 34395742267/232016584800, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 18051406831/696049754400. (End)
From G. C. Greubel, Mar 13 2022: (Start)
G.f.: x*(27 - 25*x)/(1-x)^3.
E.g.f.: x*(27 + x)*exp(x). (End)

A132769 a(n) = n*(n + 27).

Original entry on oeis.org

0, 28, 58, 90, 124, 160, 198, 238, 280, 324, 370, 418, 468, 520, 574, 630, 688, 748, 810, 874, 940, 1008, 1078, 1150, 1224, 1300, 1378, 1458, 1540, 1624, 1710, 1798, 1888, 1980, 2074, 2170, 2268, 2368, 2470, 2574, 2680, 2788, 2898, 3010, 3124, 3240, 3358, 3478

Offset: 0

Omar E. Pol, Aug 28 2007

Harvey P. Dale, Table of n, a(n) for n = 0..1000
Felix P. Muga II, Extending the Golden Ratio and the Binet-de Moivre Formula, Preprint on ResearchGate, March 2014.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A001008, A002378, A005563, A028347, A028552, A028557, A028560, A028563, A028566.
Cf. A028569, A098849, A098850, A098603, A098847, A098848, A102928, A120071, A132759.
Cf. A132760, A132761, A132762, A132763, A132764, A132765, A132766, A132767, A132768.
Cf. A132756.

Table[n(n+27),{n,0,50}] (* or *) LinearRecurrence[{3,-3,1},{0,28,58},50] (* Harvey P. Dale, Oct 14 2012 *)

a(n)=n*(n+27) \\ Charles R Greathouse IV, Oct 07 2015

a(n) = 2*n + a(n-1) + 26, with a(0)=0. - Vincenzo Librandi, Aug 03 2010
a(0)=0, a(1)=28, a(2)=58; for n > 2, a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Harvey P. Dale, Oct 14 2012
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(27)/27 = A001008(27)/A102928(27) = 312536252003/2168462696400, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*log(2)/27 - 57128792093/2168462696400. (End)
From Elmo R. Oliveira, Nov 29 2024: (Start)
G.f.: 2*x*(14 - 13*x)/(1 - x)^3.
E.g.f.: exp(x)*x*(28 + x).
a(n) = 2*A132756(n). (End)

A256102 Numbers m such that gcd(A001008(m), m) > 1, in increasing order.

Original entry on oeis.org

20, 42, 77, 110, 156, 272, 342, 506, 812, 930, 1247, 1332, 1640, 1806, 2162, 2756, 3422, 3660, 4422, 4970, 5256, 6162, 6806, 7832, 9312, 9328, 10100, 10506, 11342, 11772, 12656, 16002, 17030, 18632, 19182, 22052, 22650, 24492, 26406, 27722, 29756, 31862, 32580, 36290, 37056, 38612, 39402

Offset: 1

Wolfdieter Lang, Apr 16 2015

For the corresponding values of GCD(A001008(a(n)), a(n)) see A256103(n).
A001008(a(n))/A175441(a(n)) = A256103(n), n >= 1.
This means that for all values n not in the present sequence the numerator of the harmonic sum (HS) of the first n positive integers coincides with the denominator of the harmonic mean (HM) of the first n positive integers. That is, n divides the HM(n) numerator A102928(n) for n not in the present sequence.
Of course, HS(n)*HM(n) = n, n >= 1, where HS(n) = A001008(n)/A002805(n) and HM(n) = A102928(n)/A175441(n).
All terms are composite. Sequences contains all numbers of the form p*(p - 1), where p is a prime >= 5. This is because p^2 divides numerator(Sum_{i=1..p-1} 1/(k*p + i)), and p divides numerator(Sum_{i=1..p-1} 1/(i*p)), so p divides numerator(Sum_{i=1..p*(p-1)} 1/i). - Jianing Song, Dec 24 2018

			n = 1: gcd(A001008(20), 20) = gcd(55835135, 20) = 5 = A256103(1) > 1.
A001008(20)/A175441(20) = 55835135/11167027 = 5 = A256103(1).
Because 19 is not in this sequence 1 = gcd(A001008(19), 19) = gcd(275295799, 19).

Amiram Eldar, Table of n, a(n) for n = 1..500

Cf. A256103, A001008, A002805, A102928, A175441.

Select[Range[10^4], !CoprimeQ[#, Numerator @ HarmonicNumber[#]] &] (* Amiram Eldar, Feb 24 2020 *)

a(n) is the n-th smallest element of the set M:= {m positive inter | gcd(A001008(m), m) > 1}, n >= 1.

A132770 a(n) = n*(n + 28).

Original entry on oeis.org

0, 29, 60, 93, 128, 165, 204, 245, 288, 333, 380, 429, 480, 533, 588, 645, 704, 765, 828, 893, 960, 1029, 1100, 1173, 1248, 1325, 1404, 1485, 1568, 1653, 1740, 1829, 1920, 2013, 2108, 2205, 2304, 2405, 2508, 2613, 2720, 2829, 2940, 3053, 3168, 3285, 3404, 3525

Offset: 0

Omar E. Pol, Aug 28 2007

G. C. Greubel, Table of n, a(n) for n = 0..5000
Felix P. Muga II, Extending the Golden Ratio and the Binet-de Moivre Formula, Preprint on ResearchGate, March 2014.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A001008, A002378, A005563, A028347, A028552, A028557, A028560, A028563, A028566, A028569.
Cf. A098603, A098847, A098848, A098849, A098850, A102928, A120071, A132759, A132760, A132761.
Cf. A132762, A132763, A132764, A132765, A132766, A132767, A132768, A132769.

#(#+28)&/@Range[0,50] (* or *) LinearRecurrence[{3,-3,1},{0,29,60},50] (* Harvey P. Dale, Apr 30 2018 *)

a(n)=n*(n+28) \\ Charles R Greathouse IV, Jun 17 2017

[n*(n+28) for n in (0..50)] # G. C. Greubel, Mar 13 2022

a(n) = 2*n + a(n-1) + 27, with a(0)=0. - Vincenzo Librandi, Aug 03 2010
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(28)/28 = A001008(28)/A102928(28) = 315404588903/2248776129600, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 7751493599/321253732800. (End)
G.f.: x*(29 - 27*x)/(1-x)^3. - Harvey P. Dale, Aug 03 2021
E.g.f.: x*(29 + x)*exp(x). - G. C. Greubel, Mar 13 2022

A132771 a(n) = n*(n + 29).

Original entry on oeis.org

0, 30, 62, 96, 132, 170, 210, 252, 296, 342, 390, 440, 492, 546, 602, 660, 720, 782, 846, 912, 980, 1050, 1122, 1196, 1272, 1350, 1430, 1512, 1596, 1682, 1770, 1860, 1952, 2046, 2142, 2240, 2340, 2442, 2546, 2652, 2760, 2870, 2982, 3096, 3212, 3330, 3450, 3572

Offset: 0

Omar E. Pol, Aug 28 2007

G. C. Greubel, Table of n, a(n) for n = 0..5000
Felix P. Muga II, Extending the Golden Ratio and the Binet-de Moivre Formula, Preprint on ResearchGate, March 2014.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A001008, A002378, A005563, A028347, A028552, A028557, A028560, A028563, A028566, A028569.
Cf. A098849, A098850, A098603, A098847, A098848, A102928, A120071, A132759, A132760, A132761.
Cf. A132762, A132763, A132764, A132765, A132766, A132767, A132768, A132769, A132770.

Table[n(n+29), {n,0,50}] (* Bruno Berselli, Apr 03 2015 *)
LinearRecurrence[{3,-3,1},{0,30,62},50] (* Harvey P. Dale, Oct 18 2024 *)

a(n)=n*(n+29) \\ Charles R Greathouse IV, Jun 17 2017

[n*(n+29) for n in (0..50)] # G. C. Greubel, Mar 13 2022

a(n) = 2*n + a(n-1) + 28 (with a(0)=0). - Vincenzo Librandi, Aug 03 2010
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(29)/29 = A001008(29)/A102928(29) = 9227046511387/67543597321200, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*log(2)/29 - 236266661971/9649085331600. (End)
From G. C. Greubel, Mar 13 2022: (Start)
G.f.: 2*(15*x - 14*x^2)/(1-x)^3.
E.g.f.: x*(30 + x)*exp(x). (End)

A232193 Numerators of the expected value of the length of a random cycle in a random n-permutation.

Original entry on oeis.org

1, 3, 23, 55, 1901, 4277, 198721, 16083, 14097247, 4325321, 2132509567, 4527766399, 13064406523627, 905730205, 13325653738373, 362555126427073, 14845854129333883, 57424625956493833, 333374427829017307697, 922050973293317, 236387355420350878139797

Offset: 1

Geoffrey Critzer, Nov 20 2013

In this experiment we randomly select (uniform distribution) an n-permutation and then randomly select one of the cycles from that permutation. Cf. A102928/A175441 which gives the expected cycle length when we simply randomly select a cycle.

			Expectations for n=1,... are 1/1, 3/2, 23/12, 55/24, 1901/720, 4277/1440, 198721/60480, 16083/4480, ... = A232193/A232248
For n=3 there are 6 permutations.  We have probability 1/6 of selecting (1)(2)(3) and the cycle size is 1.  We have probability 3/6 of selecting a permutation with cycle type (1)(23) and (on average) the cycle length is 3/2.  We have probability 2/6 of selecting a permutation of the form (123) and the cycle size is 3.  1/6*1 + 3/6*3/2 + 2/6*3 = 23/12.

Alois P. Heinz, Table of n, a(n) for n = 1..250

Denominators are A232248.
Cf. A028417(n)/n! the expected value of the length of the shortest cycle in a random n-permutation.
Cf. A028418(n)/n! the expected value of the length of the longest cycle in a random n-permutation.

with(combinat):
b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
      expand(add(multinomial(n, n-i*j, i$j)/j!*(i-1)!^j
      *b(n-i*j, i-1) *x^j, j=0..n/i))))
    end:
a:= n->numer((p->add(coeff(p, x, i)/i, i=1..n))(b(n$2))/(n-1)!):
seq(a(n), n=1..30);  # Alois P. Heinz, Nov 21 2013
# second Maple program:
a:= n-> numer(add(abs(combinat[stirling1](n, i))/i, i=1..n)/(n-1)!):
seq(a(n), n=1..30);  # Alois P. Heinz, Nov 23 2013

Table[Numerator[Total[Map[Total[#]!/Product[#[[i]],{i,1,Length[#]}]/Apply[Times,Table[Count[#,k]!,{k,1,Max[#]}]]/(Total[#]-1)!/Length[#]&,Partitions[n]]]],{n,1,25}]

a(n) = Numerator( 1/(n-1)! * Sum_{i=1..n} A132393(n,i)/i ). - Alois P. Heinz, Nov 23 2013

a(n) = numerator(Sum_{k=0..n} A002657(k)/A091137(k)) (conjectured). - Michel Marcus, Jul 19 2019

A132772 a(n) = n*(n + 30).

Original entry on oeis.org

0, 31, 64, 99, 136, 175, 216, 259, 304, 351, 400, 451, 504, 559, 616, 675, 736, 799, 864, 931, 1000, 1071, 1144, 1219, 1296, 1375, 1456, 1539, 1624, 1711, 1800, 1891, 1984, 2079, 2176, 2275, 2376, 2479, 2584, 2691, 2800, 2911, 3024, 3139, 3256, 3375, 3496, 3619

Offset: 0

Omar E. Pol, Aug 28 2007

Harvey P. Dale, Table of n, a(n) for n = 0..1000
Felix P. Muga II, Extending the Golden Ratio and the Binet-de Moivre Formula, Preprint on ResearchGate, March 2014.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A001008, A002378, A005563, A028347, A028552, A028557, A028560, A028563, A028566.
Cf. A028569, A067079, A098849, A098850, A098603, A098847, A098848, A102928, A120071.
Cf. A132759, A132760, A132761, A132762, A132763, A132764, A132765, A132766, A132767.
Cf. A132768, A132769, A132770, A132771.

Table[n(n+30),{n,0,50}] (* or *) LinearRecurrence[{3,-3,1},{0,31,64},50] (* Harvey P. Dale, Mar 06 2015 *)

a(n)=n*(n+30) \\ Charles R Greathouse IV, Jun 17 2017

[n*(n+30) for n in (0..50)] # G. C. Greubel, Mar 13 2022

G.f.: x*(31-29*x)/(1-x)^3. - R. J. Mathar, Nov 14 2007
a(n) = 2*n + a(n-1) + 29 (with a(0)=0). - Vincenzo Librandi, Aug 03 2010
a(0)=0, a(1)=31, a(2)=64, a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Harvey P. Dale, Mar 06 2015
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(30)/30 = A001008(30)/A102928(30) = 9304682830147/69872686884000, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 225175759291/9981812412000. (End)
E.g.f.: x*(31 + x)*exp(x). - G. C. Greubel, Mar 13 2022

A368372 a(n) = numerator of AM(n)-HM(n), where AM(n) and HM(n) are the arithmetic and harmonic means of the first n positive integers.

Original entry on oeis.org

0, 1, 4, 29, 111, 103, 472, 2369, 12965, 30791, 197346, 452993, 3337271, 7485915, 4160656, 18358463, 170991927, 124184839, 1278605110, 110351535, 98802055, 211524139, 2595194516, 16562041459, 219589922071, 464651871609, 2207044831642, 4649180818987, 70862100349605, 148699793966557

Offset: 1

N. J. A. Sloane, Jan 24 2024

			0, 1/6, 4/11, 29/50, 111/137, 103/98, 472/363, 2369/1522, 12965/7129, 30791/14762, 197346/83711, 452993/172042, 3337271/1145993, 7485915/2343466, 4160656/1195757, 18358463/4873118, ...

Paolo Xausa, Table of n, a(n) for n = 1..2000

Cf. A102928/A175441, A368366, A368373.

AM:=proc(n) local i; (add(i,i=1..n)/n); end;
HM:=proc(n) local i; (add(1/i,i=1..n)/n)^(-1); end;
s1:=[seq(AM(n)-HM(n),n=1..50)];

A368372[n_] := Numerator[(n+1)/2 - n/HarmonicNumber[n]];
Array[A368372, 35] (* Paolo Xausa, Jan 29 2024 *)

a368372(n) = numerator((n+1)/2 - n/harmonic(n)) \\ Hugo Pfoertner, Jan 25 2024

from fractions import Fraction
from itertools import count, islice
def agen(): # generator of terms
    A = H = 0
    for n in count(1):
        A += n
        H += Fraction(1, n)
        yield ((A*Fraction(1, n) - n/H)).numerator
print(list(islice(agen(), 30))) # Michael S. Branicky, Jan 24 2024

from fractions import Fraction
from sympy import harmonic
def A368372(n): return (Fraction(n+1,2)-Fraction(n,harmonic(n))).numerator # Chai Wah Wu, Jan 25 2024

cp's OEIS Frontend

A132767 a(n) = n*(n + 25).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Mathematica

PARI

Sage

Formula

A175441 Denominator of the harmonic mean of the first n positive integers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A132768 a(n) = n*(n + 26).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Sage

Formula

A132769 a(n) = n*(n + 27).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A256102 Numbers m such that gcd(A001008(m), m) > 1, in increasing order.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A132770 a(n) = n*(n + 28).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Sage

Formula

A132771 a(n) = n*(n + 29).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Sage