MathGroup Archive 2008

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Inclusion-Exclusion Principle in Mathematica

  • To: mathgroup at smc.vnet.net
  • Subject: [mg91666] Re: Inclusion-Exclusion Principle in Mathematica
  • From: dh <dh at metrohm.ch>
  • Date: Wed, 3 Sep 2008 06:45:31 -0400 (EDT)
  • References: <g88v23$grc$1@smc.vnet.net>

Hi,

you want to calculate telescopic sums: Sum(a[i1,i2,..]) where 0<i1<n, 

i1<i2<=n, ...

Let's start small with two netsed sums:

nsum=5;

Sum[

a[i,j]

,{i,1,nsum},{j,i+1,nsum}]

this gives: 

a[1,2]+a[1,3]+a[1,4]+a[1,5]+a[2,3]+a[2,4]+a[2,5]+a[3,4]+a[3,5]+a[4,5]

Obviously, we need to create the iterators automatically. We will do 

this outside sum and then "insert" the iterators into Sum using 

"Sequence". Therefore, what we need to create looks like:

{ {x1,1,n}, {x2,x1+1,n},.. }

It is easier to use unique identifiers (using Unique[]) than x1,x2... 

The following creates a list "iters" of 3 iterators and a list "vars" of 

variables:

nsum=5;

niters=3;

last1=0;

vars={};

iters=Table[{last2=Unique[],last1+1,last1=last2;AppendTo[vars,last1];nsum},{niters}]

the variables are e.g.: {$7,$8,$9} and the corresponding ieterators: 

{{$7,1,5},{$8,1+$7,5},{$9,1+$8,5}}. Now we insert this into the sum (use 

lowercase sum to prevent evaluation in order to see what we get. If this 

is o.k. change to uppercase and reevaluate):

sum[

a[Sequence@@ vars]

,Evaluate[Sequence@@ bounds]]

this gives e.g.:

sum[a[$7,$8,$9],{$7,1,5},{$8,1+$7,5},{$9,1+$8,5}]

this looks fine and we replace sum by Sum and get:

a[1,2,3]+a[1,2,4]+a[1,2,5]+a[1,3,4]+a[1,3,5]+a[1,4,5]+a[2,3,4]+a[2,3,5]+a[2,4,5]+a[3,4,5]

hope this helps, Daniel



steaxauce wrote:

> Hey guys, first post!

> I believe I've just solved a problem, and I'm trying to use mathematica to test my solution for known results, but I haven't been able to write my expression into mathematica. It looks a lot like the inclusion-exclusion principle, which is equation (3) on the following mathworld page:

> 

> http://mathworld.wolfram.com/Inclusion-ExclusionPrinciple.html

> 

> My expression is similar in that it contains alternating series of multiple series, where the multiple series are of increasing "depth." I hope that's clear. I can enter it in the way it's written on that page, but that wouldn't be practical for me. I can't have it in the form sum1[] - sum2[] + sum3[] - sum4[], etc. It needs to be in the form sum[-1^(j-1)sumj[]], because I'm working with a very large and varying number of sets and I can't type in all of the series manually. To do this, I need to know how to specify the "depth," or number of series in a multiple series, with a number. I've gone through the documentation, contacted mathematica tech support and posted in the student support forum, but I haven't gotten a solution yet. Can anyone who's more familiar with mathematica that me help me out? I just downloaded the program a few days ago!

> 

> Thanks guys!

> 





-- 



Daniel Huber

Metrohm Ltd.

Oberdorfstr. 68

CH-9100 Herisau

Tel. +41 71 353 8585, Fax +41 71 353 8907

E-Mail:<mailto:dh at metrohm.com>

Internet:<http://www.metrohm.com>




  • Prev by Date: Programmatic headers and footers
  • Next by Date: Controlling the size of the integral sign
  • Previous by thread: Programmatic headers and footers
  • Next by thread: Controlling the size of the integral sign