joj 2736: 二叉树计数卡特兰数取模 1<=n<=10^9,1<=m<=32768

最新推荐文章于 2020-11-16 01:45:50 发布

kongming_acm

最新推荐文章于 2020-11-16 01:45:50 发布

阅读量1.6k

点赞数

分类专栏： joj acm_数学问题文章标签： pair output input c

本文链接：https://blog.csdn.net/kongming_acm/article/details/6361176

版权

该博客探讨了如何计算含有特定数量节点的不同二叉树的总数，并将其结果对一个大素数取模。输入包括节点数量n和模数k，输出为二叉树的数量对k取模后的值。

摘要由CSDN通过智能技术生成

给出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

//师哥模板

#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 %