1### R code from vignette source 'partitionspaper.Rnw' 2 3################################################### 4### code chunk number 1: partitionspaper.Rnw:122-122 5################################################### 6 7 8 9################################################### 10### code chunk number 2: partitionspaper.Rnw:123-124 11################################################### 12require(partitions) 13 14 15################################################### 16### code chunk number 3: partitionspaper.Rnw:132-133 17################################################### 18P(100) 19 20 21################################################### 22### code chunk number 4: partitionspaper.Rnw:137-138 23################################################### 24diffparts(10) 25 26 27################################################### 28### code chunk number 5: partitionspaper.Rnw:170-176 29################################################### 30f <- function(x){prod(factorial(x),factorial(tabulate(x)))} 31prob <- function(a,n){ 32 jj <- restrictedparts(n,a,include.zero=FALSE) 33 N <- factorial(a)*factorial(n)*sum(1/apply(jj,2,f)) 34 return(N/a^n) 35} 36 37 38