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.
%I A138407 #24 Nov 22 2022 02:33:44 %S A138407 16,162,2500,14406,146410,342732,1336336,2345778,6156502,19803868, %T A138407 27705630,67469796,113030440,143589642,224465326,410305012,702806938, %U A138407 830750460,1329973986,1778817670,2044673352,3038106318,3891582322 %N A138407 a(n) = p^4*(p-1), where p = prime(n). %H A138407 Vincenzo Librandi, <a href="/A138407/b138407.txt">Table of n, a(n) for n = 1..200</a> %H A138407 <a href="/index/Pri#prime_powers">Index to sequences related to prime powers</a>. %F A138407 a(n) = A000010(prime(n)^5). - _R. J. Mathar_, Oct 15 2017 %F A138407 From _Amiram Eldar_, Nov 22 2022: (Start) %F A138407 a(n) = prime(n)^5 - prime(n)^4 = A050997(n) - A030514(n). %F A138407 Product_{n>=1} (1 - 1/a(n)) = A065416. (End) %t A138407 a = {}; Do[p = Prime[n]; AppendTo[a, p^5 - p^4], {n, 1, 50}]; a %t A138407 f54[n_]:=Module[{c=Prime[n]},c^5-c^4]; Array[f54,30] (* _Harvey P. Dale_, Mar 29 2015 *) %o A138407 (PARI) forprime(p=2,1e3,print1(p^5-p^4", ")) \\ _Charles R Greathouse IV_, Jun 16 2011 %o A138407 (Magma) [NthPrime((n))^5 - NthPrime((n))^4: n in [1..30] ]; // _Vincenzo Librandi_, Jun 17 2011 %Y A138407 Cf. A000010, A030514, A050997, A065416. %K A138407 nonn,easy %O A138407 1,1 %A A138407 _Artur Jasinski_, Mar 19 2008 %E A138407 First Mathematica program corrected by _Harvey P. Dale_, Mar 29 2015