[go: up one dir, main page]

login
A079416
a(n) = round(prime(n)/n).
5
2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5
OFFSET
1,1
COMMENTS
The sequence is not monotone, see example and A079417.
LINKS
EXAMPLE
a(20) = round(prime(20)/20) = round(71/20) = round(3.55) = 4;
a(21) = round(prime(21)/21) = round(73/21) = round(3.476190...) = 3;
a(22) = round(prime(22)/22) = round(79/22) = round(3.590909...) = 4.
MATHEMATICA
f[n_] := Round[ Prime[n]/n]; Array[f, 105] (* Robert G. Wilson v, Oct 23 2015 *)
PROG
(PARI) vector(100, n, round(prime(n)/n)) \\ G. C. Greubel, Jan 18 2019
(Magma) [Round(NthPrime(n)/n): n in [1..100]]; // G. C. Greubel, Jan 18 2019
(Sage) [round(nth_prime(n)/n) for n in (1..100)] # G. C. Greubel, Jan 18 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Jan 07 2003
STATUS
approved