%I #49 Apr 04 2024 10:04:25
%S 1,7,57,586,7465,114381,2054353,42374116,987654321,25678050355,
%T 736867805641,23136292864686,789018236134297,29043982525261081,
%U 1147797409030816545,48471109094902544776,2178347851919531492065,103805969587115219182431
%N Digits 1..n in strict descending order n..1 interpreted in base n+1.
%C All odd-indexed (2n+1) terms are divisible by (2n+1). See A051847.
%C All even-indexed (2n) terms are divisible by n. - _Alexander R. Povolotsky_, Oct 20 2022
%H Alois P. Heinz, <a href="/A051846/b051846.txt">Table of n, a(n) for n = 1..200</a>
%H Chai Wah Wu, <a href="https://arxiv.org/abs/2403.20304">Pandigital and penholodigital numbers</a>, arXiv:2403.20304 [math.GM], 2024. See p. 1.
%F a(n) = Sum_{i=1..n} i*(n+1)^(i-1).
%F a(n) = ((n+1)^(n+1)*(n-1) + 1)/n^2 = A062806(n+1)/(n+1) - (n+1)^(n+1). - _Benoit Cloitre_, Sep 28 2002
%F a(n) = A028310(n-1) * A023811(n+1) + A199969(n+1). - _M. F. Hasler_, Jan 22 2013
%F a(n) = (n-1) * A058128(n+1) + 1. - _Seiichi Manyama_, Apr 10 2022
%e a(1) = 1,
%e a(2) = 2*3 + 1 = 7,
%e a(3) = 3*(4^2) + 2*4 + 1 = 57,
%e a(4) = 4*(5^3) + 3*(5^2) + 2*5 + 1 = 586.
%p a(n) := proc(n) local i; add(i*((n+1)^(i-1)),i=1..n); end;
%t Array[Sum[i*(# + 1)^(i - 1), {i, #}] &, 18] (* _Michael De Vlieger_, Apr 04 2024 *)
%o (PARI) a(n)=((n+1)^(n+1)*(n-1)+1)/n^2
%o (Maxima) makelist(((n+1)^(n+1)*(n-1) + 1)/n^2,n,1,20); /* _Martin Ettl_, Jan 25 2013 */
%o (Python)
%o def a(n): return sum((i+1)*(n+1)**i for i in range(n))
%o print([a(n) for n in range(1, 20)]) # _Michael S. Branicky_, Apr 10 2022
%Y The right edge of A051845.
%Y Cf. A023811, A051847, A058128, A062806, A199969.
%K easy,base,nonn
%O 1,2
%A _Antti Karttunen_, Dec 13 1999
%E Minor edits in formulas by _M. F. Hasler_, Oct 11 2019