I’m looking for help with calculating the dice probabilities for the Heresy dice system, as featured in the Victoriana (3rd ed, in case it matters) RPG.

It’s a d6 dice pool system. You roll attribute+skill d6s. Each 1 or 6 is a success, and for each 6 you roll, you roll an additional die. Lather, rinse, repeat.

More-difficult tasks mean adding more dice into the mix, which the game calls “black dice”. Each black die subtracts a success if it rolls a 1 or 6. Black dice do **not** roll extra dice on a 6.

So, the black die are easy. That part is just

`- [count {1,6} in Bd6] `

Where B is the variable I’ve stuffed the number of black dice into.

It seems to me that the overall logic (algorithm/pseudocode) should be something like:

`Input N skill + B black initiate empty variable S loop: roll N skill dice { for each die 1: increase S by 1 6: increase S by 1, increase N by 1 } loop: roll B black dice { for each die 1,6: decrease S by 1 } `

But I’m failing to translate that into Anydice syntax. Repeatedly, in many different ways. Any help?