|
%I #26 Feb 27 2020 03:03:01
%S 9,81,18,81,9,702,9,171,27,72,18,693,18,72,27,171,9,702,9,81,18,81,9,
%T 5913,9,81,18,81,9,1602,9,261,36,63,27,594,18,162,36,162,18,603,9,171,
%U 27,72,18,5814,9,171,27,72,18,603,9,261,36,63,27,1584,27,63,36,261,9
%N First difference of A050289 = numbers whose digits are a permutation of (1,...,9).
%C This sequence is the natural extension of A107346 (and others, see below) from 5!-1 to 9!-1 terms, which is the natural (since maximal) length, given that OEIS sequence data are stored as decimal numbers. On the other hand, it is quite different from A219664 in many aspects, not only for the reason that the other sequence is infinite and therefore differs from this one in all terms beyond n = 9!-1.
%C The sequence is finite, with 9!-1 terms, and symmetric: a(n)=a(9!-n).
%C All terms are multiples of 9, cf. formula.
%C The subsequence of the first n!-1 terms (n=2,...,9) yields the first differences of the sequence of numbers whose digits are a permutation of (1,...,n):
%C The first 8!-1 terms yield the first differences of A178478: numbers whose digits are a permutation of 12345678.
%C The first 7!-1 terms yield the first differences of A178477: numbers whose digits are a permutation of 1234567.
%C The first 6!-1 terms yield the first differences of A178476: numbers whose digits are a permutation of 123456.
%C The first 5!-1 terms yield A107346, the first differences of A178475: numbers whose digits are a permutation of 12345.
%F a(n) = A219664(n) = 9*A217626(n) (for n < 9!). - _M. F. Hasler_, Jan 12 2013
%F a(n) = a(m!-n) for any m < 10 such that n < m!.
%e The same initial terms are obtained for the permutations of any set of the form {1,...,m}, e.g., {1,2,3} or {1,...,9}: In the first case we have P = (123,132,213,231,312,321) and P(4)-P(3) = 231 - 213 = 18 = a(3), and in the latter case P(4)-P(3) = 123456897 - 123456879 = 18, again. - _M. F. Hasler_, Jan 12 2013
%t Take[Differences[Sort[FromDigits/@Permutations[Range[9]]]],70] (* _Harvey P. Dale_, Mar 31 2018 *)
%o (PARI) A209280_list(N=5)={my(v=vector(N,i,10^(N-i))~); v=vecsort(vector(N!,k,numtoperm(N,k)*v)); vecextract(v,"^1")-vecextract(v,"^-1")} \\ return the N!-1 first terms as a vector
%o (PARI) A209280(n)={if(a209280=='a209280 || #a209280<n, a209280=A209280_list(A090529(n+1)));a209280[n]}
%Y Cf. A030299, A219664.
%K easy,nonn,base,fini
%O 1,1
%A _M. F. Hasler_, Jan 12 2013
|
