给出n,k.求含有n个节点的不同二叉树数目%k
Input
MULTI TEST CASE!!
每行两个整数n,k(1<=n<=10^9,1<=k<=32768)
Output
一行二叉树的数目%k
Sample Input
3 100
Sample Output
5
//师哥模板
#include <stdio.h> #include <utility> #include <vector> #include <algorithm>
int Pow(int x, int b) { int ret = 1; for(int s = x; b; b >>= 1) { if(b & 1) ret *= s; s *= s; } return ret; } int PowMod(int x, int b, int p) { int ret = 1 % p; for(int s = x % p; b; b >>= 1) { if(b & 1) { ret *= s; ret %= p; } s *= s; s %= p; } return ret; }
std::pair<int, int> ExGcd(int a, int b) { if(a == 0) return std::make_pair(0, 1); std::pair<int, int> s = ExGcd(b % a, a); return std::make_pair(s.second - b / a * s.first, s.first); }
int Inv(int a, int m) // ax == 1 mod m { // assert gcd(a, m) == 1 std::pair<int, int> s = ExGcd(a, m); // s.first * a + s.second * m == 1 return s.first < 0 ? s.first + m : s.first; }
std::pair<int, int> FacMod(int n, int p, int k) { // assert p > 1, k > 0 int pk = Pow(p, k); int S[pk]; S[0] = 1 % pk; for(int i = 1; i < pk; ++i) { S[i] = S[i - 1]; if(i % p != 0) { S[i] *= i; S[i] %= pk; } } int ret = 1 % pk, et = 0, ep = 0; // S[pk - 1]^et, p^ep while(n) { et += n / pk; ret *= S[n % pk]; ret %= pk; ep += n / p; n /= p; } ret *= PowMod(S[pk - 1], et, pk); ret %