# must first enter n = number of dice, N = number of repetitions 
#
score = rep(NA, times = N)
die = (1:6)
for (i in 1:N) {
   dice = sample(die, n, replace=TRUE)
   isitone = (dice == 1)
   anyones = (sum(isitone) > 0)
   if(anyones == TRUE) {score[i]=0} else {score[i] = sum(dice)}
}
table(score)
hist(score)
mean(score)