a=1:100 b=outer(a,a,">=") for(i in a)if(i>2)for(j in 2:(i-1))b[i,j]=b[i-j,j]+b[i-1,j-1] write(paste("p(",a,")=",rowSums(b),sep=""),T)

Note that non-ascii characters in the above source code will be escaped (such as \x9f).

download

return to the top page