频道栏目
首页 > 程序开发 > 软件开发 > C++ > 正文
HDU 3037 Saving Beans (Lucas定理)
2014-06-18 10:59:30      个评论    来源:HDU 3037 Saving Beans (Lucas定理)  
收藏   我要投稿

题目链接:https://acm.hdu.edu.cn/showproblem.php?pid=3037

推出公式为C(n + m, m) % p, 用Lucas定理求大组合数取模的值

代码:

#include 
#include 
#include 
using namespace std;

int t;
long long n, m, p;

long long pow(long long n, long long k) {
    if (k == 0) return 1;
    if (k == 1) return n;
    long long ans = pow(n * n % p, k>>1);
    if (k&1) ans = ans * n % p;
    return ans;
}

long long C(long long n, long long m) {
    if (m > n) return 0;
    m = min(m, n - m);
    long long zi = 1, mu = 1;
    for (long long i = 0; i < m; i++) {
	zi = zi * (n - i) % p;
	mu = mu * (i + 1) % p;
    }
    return zi * pow(mu, p - 2) % p;
}

long long Lucas(long long n, long long m, long long p) {
    if (m == 0) return 1;
    return C(n % p, m % p) * Lucas(n / p, m / p, p) % p;
}

int main() {
    scanf("%d", &t);
    while (t--) {
	scanf("%lld%lld%lld", &n, &m, &p);
	printf("%lld\n", Lucas(n + m, m, p));
    }
    return 0;
}


点击复制链接 与好友分享!回本站首页
相关TAG标签 定理
上一篇:基本排序系列之基数排序
下一篇:POJ 2506 -TILING
相关文章
图文推荐
点击排行

关于我们 | 联系我们 | 广告服务 | 投资合作 | 版权申明 | 在线帮助 | 网站地图 | 作品发布 | Vip技术培训 | 举报中心

版权所有: 红黑联盟--致力于做实用的IT技术学习网站