pub fn bell_number<const P: u32, const PR: u32>(n: usize) -> Vec<ConstModInt<P>>
ベル数$B_0, \dots, B_n$を列挙する。