A227653
a(1) = least k such that 1/2 + 1/3 < H(k) - H(3); a(2) = least k such that H(a(1)) - H(3) < H(k) -H(a(1)), and for n > 2, a(n) = least k such that H(a(n-1)) - H(a(n-2)) > H(k) - H(a(n-1)), where H = harmonic number.
Original entry on oeis.org
8, 21, 54, 138, 352, 897, 2285, 5820, 14823, 37752, 96148, 244872, 623645, 1588311, 4045140, 10302237, 26237926, 66823230, 170186624, 433434405, 1103878665, 2811378360, 7160069791, 18235396608, 46442241368, 118279949136, 301237536249, 767197263003
Offset: 1
The first two values (a(1),a(2)) = (8,21) match the beginning of the following inequality chain: 1/2 + 1/3 < 1/4 + ... + 1/8 < 1/9 + ... + 1/21 < ...
-
z = 100; h[n_] := h[n] = HarmonicNumber[N[n, 500]]; x = 2; y = 3;
a[1] = Ceiling[w /. FindRoot[h[w] == 2 h[y] - h[x - 1], {w, 1}, WorkingPrecision -> 400]]; a[2] = Ceiling[w /. FindRoot[h[w] == 2 h[a[1]] - h[y], {w, a[1]}, WorkingPrecision -> 400]]; Do[s = 0; a[t] = Ceiling[w /. FindRoot[h[w] == 2 h[a[t - 1]] - h[a[t - 2]], {w, a[t - 1]}, WorkingPrecision -> 400]], {t, 3, z}];
m = Map[a, Range[z]] (* A227653, Peter J. C. Moses, Jul 12 2013 *)
A228016
a(1) = least k such that 1/1+1/2+1/3+1/4+1/5 < H(k) - H(5); a(2) = least k such that H(a(1)) - H(5) < H(k) -H(a(1)), and for n > 2, a(n) = least k such that H(a(n-1)) - H(a(n-2)) > H(k) - H(a(n-1)), where H = harmonic number.
Original entry on oeis.org
54, 539, 5340, 52865, 523314, 5180279, 51279480, 507614525, 5024865774, 49741043219, 492385566420, 4874114620985, 48248760643434, 477613491813359, 4727886157490160, 46801248083088245, 463284594673392294, 4586044698650834699, 45397162391834954700
Offset: 1
The first two values (a(1),a(2)) = (54,539) match the beginning of the following inequality chain: 1/1+1/2+1/3+1/4+1/5 < 1/6+...+1/54 < 1/55+...+1/539 < ...
-
z = 100; h[n_] := h[n] = HarmonicNumber[N[n, 500]]; x = 1; y = 5;
a[1] = Ceiling[w /. FindRoot[h[w] == 2 h[y] - h[x - 1], {w, 1}, WorkingPrecision -> 400]]; a[2] = Ceiling[w /. FindRoot[h[w] == 2 h[a[1]] - h[y], {w, a[1]}, WorkingPrecision -> 400]]; Do[s = 0; a[t] = Ceiling[w /. FindRoot[h[w] == 2 h[a[t - 1]] - h[a[t - 2]], {w, a[t - 1]}, WorkingPrecision -> 400]], {t, 3, z}];
m = Map[a, Range[z]] (* A227653, Peter J. C. Moses, Jul 12 2013 *)
A225605
a(1) = least k such that 1/3 < H(k) - 1/3; a(2) = least k such that H(a(1)) - H(3) < H(k) - H(a(1)), and for n > 2, a(n) = least k such that H(a(n-1)) - H(a(n-2)) > H(k) - H(a(n-1)), where H = harmonic number.
Original entry on oeis.org
5, 9, 16, 29, 53, 97, 178, 327, 601, 1105, 2032, 3737, 6873, 12641, 23250, 42763, 78653, 144665, 266080, 489397, 900141, 1655617, 3045154, 5600911, 10301681, 18947745, 34850336, 64099761, 117897841, 216847937, 398845538, 733591315, 1349284789, 2481721641
Offset: 1
The first two values (a(1),a(2)) = (5,9) match the beginning of the following inequality chain:
1/3 < 1/4 + 1/5 < 1/6 + 1/7 + 1/8 + 1/9 < ...
-
z = 100; h[n_] := h[n] = HarmonicNumber[N[n, 500]]; x = 3; y = 3;
a[1] = Ceiling[w /. FindRoot[h[w] == 2 h[y] - h[x - 1], {w, 1}, WorkingPrecision -> 400]]; a[2] = Ceiling[w /. FindRoot[h[w] == 2 h[a[1]] - h[y], {w, a[1]}, WorkingPrecision -> 400]]; Do[s = 0; a[t] = Ceiling[w /. FindRoot[h[w] == 2 h[a[t - 1]] - h[a[t - 2]], {w, a[t - 1]}, WorkingPrecision -> 400]], {t, 3, z}];
m = Map[a, Range[z]] (* A225605, Peter J. C. Moses, Jul 12 2013 *)
A228025
a(1) = least k such that 1/2+1/3+1/4+1/5 < H(k) - H(5); a(2) = least k such that H(a(1)) - H(5) < H(k) -H(a(1)), and for n > 2, a(n) = least k such that H(a(n-1)) - H(a(n-2)) > H(k) - H(a(n-1)), where H = harmonic number.
Original entry on oeis.org
20, 76, 285, 1065, 3976, 14840, 55385, 206701, 771420, 2878980, 10744501, 40099025, 149651600, 558507376, 2084377905, 7779004245, 29031639076, 108347552060, 404358569165, 1509086724601, 5631988329240, 21018866592360, 78443478040201, 292755045568445
Offset: 1
The first two values (a(1),a(2)) = (20,76) match the beginning of the following inequality chain: 1/2+1/3+1/4+1/5 < 1/6+...+1/20 < 1/21+...+1/76 < ...
-
z = 100; h[n_] := h[n] = HarmonicNumber[N[n, 500]]; x = 2; y = 5;
a[1] = Ceiling[w /. FindRoot[h[w] == 2 h[y] - h[x - 1], {w, 1}, WorkingPrecision -> 400]]; a[2] = Ceiling[w /. FindRoot[h[w] == 2 h[a[1]] - h[y], {w, a[1]}, WorkingPrecision -> 400]]; Do[s = 0; a[t] = Ceiling[w /. FindRoot[h[w] == 2 h[a[t - 1]] - h[a[t - 2]], {w, a[t - 1]}, WorkingPrecision -> 400]], {t, 3, z}];
m = Map[a, Range[z]] (* A227653, Peter J. C. Moses, Jul 12 2013 *)
Showing 1-4 of 4 results.
Comments