       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