|
%I #12 Oct 28 2017 21:39:30
%S 1,1,1,1,1,1,1,1,1,1,2,4,8,12,18,27,36,48,64,80,100,125,175,245,343,
%T 490,700,1000,1400,1960,2744,3724,5054,6859,9386,12844,17576,24336,
%U 33696,46656,64800,90000,125000,172500,238050,328509,452295,622725,857375,1182275,1630295,2248091,3106141,4291691,5929741,8190250,11312500,15625000,21562500
%N Number of subsets of {1,2,...,n-9} without differences equal to 3, 6 or 9.
%C a(n) is the number of permutations (p(1), p(2), ..., p(n)) satisfying -k <= p(i)-i <= r and p(i)-i in the set I, i=1..n, with k=3, r=9, I={-3,0,9}.
%H G. C. Greubel, <a href="/A224814/b224814.txt">Table of n, a(n) for n = 0..1000</a>
%H Vladimir Baltic, <a href="http://pefmath.etf.rs/vol4num1/AADM-Vol4-No1-119-135.pdf">On the number of certain types of strongly restricted permutations</a>, Applicable Analysis and Discrete Mathematics Vol. 4, No 1 (April, 2010), 119-13
%F a(n) = a(n-1) +a(n-4) -a(n-6) +a(n-7) -3*a(n-9) +4*a(n-10) +5*a(n-12) -2*a(n-13) +3*a(n-15) -8*a(n-16) +a(n-18) -4*a(n-19) +3*a(n-21) -4*a(n-22) -3*a(n-24) -5*a(n-27) +8*a(n-28) +7*a(n-30) -2*a(n-31) -9*a(n-33) +2*a(n-34) +5*a(n-36) +4*a(n-37) +a(n-39) -6*a(n-40) -3*a(n-42) +2*a(n-43) +2*a(n-45) +a(n-46) -a(n-48) -a(n-49) +a(n-52).
%F G.f.: (1 -x^4 -x^5 -x^7 -x^8 +2*x^9 -x^10 -3*x^12 -x^13 -2*x^15 +3*x^16 +3*x^17 +2*x^18 -x^20 -4*x^21 +x^23 +3*x^24 +3*x^25 +x^27 -4*x^28 -x^29 -2*x^30 +x^31 +2*x^33 +x^34 -x^36 -x^37 +x^40 )/((1-x-x^4)*(1-x^9-x^12)*(1 +x^6 +4*x^9 -4*x^12 -2*x^15 +4*x^18 -3*x^21 -3*x^24 +7*x^27 -6*x^30 +3*x^33 -x^36)).
%F a(3*k) = (A003269(k))^3,
%F a(3*k+1) = (A003269(k))^2 * A003269(k+1),
%F a(3*k+2) = A003269(k) * (A003269(k+1))^2.
%t CoefficientList[Series[(1 - x^4 - x^5 - x^7 - x^8 + 2*x^9 - x^10 - 3*x^12 - x^13 - 2*x^15 + 3*x^16 + 3*x^17 + 2*x^18 - x^20 - 4*x^21 + x^23 + 3*x^24 + 3*x^25 + x^27 - 4*x^28 - x^29 - 2*x^30 + x^31 + 2*x^33 +x^34 - x^36 - x^37 + x^40)/((1 - x - x^4)*(1 - x^9 - x^12)*(1 + x^6 + 4*x^9 - 4*x^12 - 2*x^15 + 4*x^18 - 3*x^21 - 3*x^24 + 7*x^27 - 6*x^30 + 3*x^33 - x^36)), {x, 0, 50}], x] (* _G. C. Greubel_, Oct 28 2017 *)
%o (PARI) x='x+O('x^50); Vec((1 -x^4 -x^5 -x^7 -x^8 +2*x^9 -x^10 -3*x^12 -x^13 -2*x^15 +3*x^16 +3*x^17 +2*x^18 -x^20 -4*x^21 +x^23 +3*x^24 +3*x^25 +x^27 -4*x^28 -x^29 -2*x^30 +x^31 +2*x^33 +x^34 -x^36 -x^37 +x^40 )/((1-x-x^4)*(1-x^9-x^12)*(1 +x^6 +4*x^9 -4*x^12 -2*x^15 +4*x^18 -3*x^21 -3*x^24 +7*x^27 -6*x^30 +3*x^33 -x^36))) \\ _G. C. Greubel_, Oct 28 2017
%Y Cf. A002524-A002529, A072827, A072850-A072856, A079955-A080014, A217694, A224808-A224815.
%K nonn,easy
%O 0,11
%A _Vladimir Baltic_, May 18 2013
|
