今天我刷题了吗
  • 2U's GitBook
  • 我的模板
    • I/O 及其他
      • 快速读写
      • 关闭同步流
      • 文件重定向
      • 随机数
      • 运行时间
      • Lambda
      • 万能头文件
    • 字符串
      • 字符串分割
      • KMP 算法
      • Manacher 算法
      • AC 自动机
      • 高精度
    • 数学
      • 快速幂
      • 矩阵快速幂
      • 筛法
      • 欧拉函数
      • Polya 定理
      • 逆元
      • 组合数
    • 图论
      • Dijkstra 算法
      • SPFA 算法
      • Kruskal 算法
      • Prim 算法
      • 倍增
      • 离线 Tarjan 算法
      • Tarjan 算法
      • Hungary 算法
      • A* 算法
      • EK 算法
      • Dinic 算法
    • 数据结构
      • 线段树
      • 树状数组
      • 可持久化数组
      • 主席树
      • 树堆
      • 无旋树堆
      • 伸展树
      • 树套树
      • 树链剖分
      • 点分治
      • 动态树
  • 我的题单
由 GitBook 提供支持
在本页

这有帮助吗?

  1. 我的模板
  2. 数学

组合数

https://www.zybuluo.com/ArrowLLL/note/713749

Pascal 公式打表组合数取模

const int maxn = 1e3 + 5;
const int mod = 1e5 + 3;
int tC[maxn][maxn];

inline int C(int n, int k)
{
    if (k > n) return 0;
    return tC[n][k];
}

void calcC(int n)
{
    for (int i = 0; i < n; i++) {
        tC[i][0] = 1;
        for (int j = 1; j < i; j++) {
            tC[i][j] = (C(i - 1, j - 1) + C(i - 1, j)) % mod;
        }
        tC[i][i] = 1;
    }
}
const int maxn = 1e3 + 5;
const int mod = 1e5 + 3;
int tC[maxn * maxn];

inline int loc(int n, int k)
{
    int pos = (1 + (n >> 1)) * (n >> 1);
    pos += k;
    pos += (n & 1) ? (n + 1) >> 1 : 0;
    return pos;
}

inline int C(int n, int k)
{
    if (k > n) return 0;
    k = min(n - k, k);
    return tC[loc(n, k)];
}

void calcC(int n)
{
    for (int i = 0; i < n; i++) {
        tC[loc(i, 0)] = 1;
        for (int j = 1, e = i >> 1; j <= e; j++) {
            tC[loc(i, j)] = C(i - 1, j) + C(i - 1, j - 1);
        }
    }
}

逆元求组合数取模

上一页逆元下一页图论

最后更新于4年前

这有帮助吗?