热门标签 | HotTags
当前位置:  开发笔记 > 编程语言 > 正文

求组合数

第一类时mod为素数1typedeflonglongLL;23LLn,m,p;45LLquick_mod(LLa,LLb)6{7LLans1;8a%p;9while(b

第一类时mod为素数

,,
 1 typedef long long LL;
 2 
 3 LL n,m,p;
 4 
 5 LL quick_mod(LL a, LL b)
 6 {
 7     LL ans = 1;
 8     a %= p;
 9     while(b)
10     {
11         if(b & 1)
12         {
13             ans = ans * a % p;
14             b--;
15         }
16         b >>= 1;
17         a = a * a % p;
18     }
19     return ans;
20 }
21 
22 LL C(LL n, LL m)
23 {
24     if(m > n) return 0;
25     LL ans = 1;
26     for(int i=1; i<=m; i++)
27     {
28         LL a = (n + i - m) % p;
29         LL b = i % p;
30         ans = ans * (a * quick_mod(b, p-2) % p) % p;
31     }
32     return ans;
33 }
34 
35 LL Lucas(LL n, LL m)
36 {
37     if(m == 0) return 1;
38     return C(n % p, m % p) * Lucas(n / p, m / p) % p;
39 }
p == prime
第二类时mod可以为合数
,,
 1 typedef long long LL;
 2 const int N = 200005;
 3 
 4 bool prime[N];
 5 int p[N];
 6 int cnt;
 7 
 8 void isprime()
 9 {
10     cnt = 0;
11     memset(prime,true,sizeof(prime));
12     for(int i=2; i)
13     {
14         if(prime[i])
15         {
16             p[cnt++] = i;
17             for(int j=i+i; ji)
18                 prime[j] = false;
19         }
20     }
21 }
22 
23 LL quick_mod(LL a,LL b,LL m)
24 {
25     LL ans = 1;
26     a %= m;
27     while(b)
28     {
29         if(b & 1)
30         {
31             ans = ans * a % m;
32             b--;
33         }
34         b >>= 1;
35         a = a * a % m;
36     }
37     return ans;
38 }
39 
40 LL Work(LL n,LL p)
41 {
42     LL ans = 0;
43     while(n)
44     {
45         ans += n / p;
46         n /= p;
47     }
48     return ans;
49 }
50 
51 LL Solve(LL n,LL m,LL P)
52 {
53     LL ans = 1;
54     for(int i=0; i)
55     {
56         LL x = Work(n, p[i]);
57         LL y = Work(n - m, p[i]);
58         LL z = Work(m, p[i]);
59         x -= (y + z);
60         ans *= quick_mod(p[i],x,P);
61         ans %= P;
62     }
63     return ans;
64 }
65 
66 int main()
67 {
68     int T;
69     isprime();
70     cin>>T;
71     while(T--)
72     {
73         LL n,m,P;
74         cin>>n>>m>>P;
75         n += m - 2;
76         m--;
77         cout<endl;
78     }
79     return 0;
80 }
View Code

求组合数


推荐阅读
  • 本文内容为asp.net微信公众平台开发的目录汇总,包括数据库设计、多层架构框架搭建和入口实现、微信消息封装及反射赋值、关注事件、用户记录、回复文本消息、图文消息、服务搭建(接入)、自定义菜单等。同时提供了示例代码和相关的后台管理功能。内容涵盖了多个方面,适合综合运用。 ... [详细]
  • 本文介绍了lua语言中闭包的特性及其在模式匹配、日期处理、编译和模块化等方面的应用。lua中的闭包是严格遵循词法定界的第一类值,函数可以作为变量自由传递,也可以作为参数传递给其他函数。这些特性使得lua语言具有极大的灵活性,为程序开发带来了便利。 ... [详细]
  • 本文介绍了C#中数据集DataSet对象的使用及相关方法详解,包括DataSet对象的概述、与数据关系对象的互联、Rows集合和Columns集合的组成,以及DataSet对象常用的方法之一——Merge方法的使用。通过本文的阅读,读者可以了解到DataSet对象在C#中的重要性和使用方法。 ... [详细]
  • 本文介绍了OC学习笔记中的@property和@synthesize,包括属性的定义和合成的使用方法。通过示例代码详细讲解了@property和@synthesize的作用和用法。 ... [详细]
  • Mac OS 升级到11.2.2 Eclipse打不开了,报错Failed to create the Java Virtual Machine
    本文介绍了在Mac OS升级到11.2.2版本后,使用Eclipse打开时出现报错Failed to create the Java Virtual Machine的问题,并提供了解决方法。 ... [详细]
  • 在说Hibernate映射前,我们先来了解下对象关系映射ORM。ORM的实现思想就是将关系数据库中表的数据映射成对象,以对象的形式展现。这样开发人员就可以把对数据库的操作转化为对 ... [详细]
  • 本文介绍了在SpringBoot中集成thymeleaf前端模版的配置步骤,包括在application.properties配置文件中添加thymeleaf的配置信息,引入thymeleaf的jar包,以及创建PageController并添加index方法。 ... [详细]
  • 知识图谱——机器大脑中的知识库
    本文介绍了知识图谱在机器大脑中的应用,以及搜索引擎在知识图谱方面的发展。以谷歌知识图谱为例,说明了知识图谱的智能化特点。通过搜索引擎用户可以获取更加智能化的答案,如搜索关键词"Marie Curie",会得到居里夫人的详细信息以及与之相关的历史人物。知识图谱的出现引起了搜索引擎行业的变革,不仅美国的微软必应,中国的百度、搜狗等搜索引擎公司也纷纷推出了自己的知识图谱。 ... [详细]
  • 本文讲述了作者通过点火测试男友的性格和承受能力,以考验婚姻问题。作者故意不安慰男友并再次点火,观察他的反应。这个行为是善意的玩人,旨在了解男友的性格和避免婚姻问题。 ... [详细]
  • 本文详细介绍了Linux中进程控制块PCBtask_struct结构体的结构和作用,包括进程状态、进程号、待处理信号、进程地址空间、调度标志、锁深度、基本时间片、调度策略以及内存管理信息等方面的内容。阅读本文可以更加深入地了解Linux进程管理的原理和机制。 ... [详细]
  • 1,关于死锁的理解死锁,我们可以简单的理解为是两个线程同时使用同一资源,两个线程又得不到相应的资源而造成永无相互等待的情况。 2,模拟死锁背景介绍:我们创建一个朋友 ... [详细]
  • 后台获取视图对应的字符串
    1.帮助类后台获取视图对应的字符串publicclassViewHelper{将View输出为字符串(注:不会执行对应的ac ... [详细]
  • 《数据结构》学习笔记3——串匹配算法性能评估
    本文主要讨论串匹配算法的性能评估,包括模式匹配、字符种类数量、算法复杂度等内容。通过借助C++中的头文件和库,可以实现对串的匹配操作。其中蛮力算法的复杂度为O(m*n),通过随机取出长度为m的子串作为模式P,在文本T中进行匹配,统计平均复杂度。对于成功和失败的匹配分别进行测试,分析其平均复杂度。详情请参考相关学习资源。 ... [详细]
  • 本文介绍了通过ABAP开发往外网发邮件的需求,并提供了配置和代码整理的资料。其中包括了配置SAP邮件服务器的步骤和ABAP写发送邮件代码的过程。通过RZ10配置参数和icm/server_port_1的设定,可以实现向Sap User和外部邮件发送邮件的功能。希望对需要的开发人员有帮助。摘要长度:184字。 ... [详细]
  • 动态规划算法的基本步骤及最长递增子序列问题详解
    本文详细介绍了动态规划算法的基本步骤,包括划分阶段、选择状态、决策和状态转移方程,并以最长递增子序列问题为例进行了详细解析。动态规划算法的有效性依赖于问题本身所具有的最优子结构性质和子问题重叠性质。通过将子问题的解保存在一个表中,在以后尽可能多地利用这些子问题的解,从而提高算法的效率。 ... [详细]
author-avatar
鱼儿没有水会死_543
这个家伙很懒,什么也没留下!
PHP1.CN | 中国最专业的PHP中文社区 | DevBox开发工具箱 | json解析格式化 |PHP资讯 | PHP教程 | 数据库技术 | 服务器技术 | 前端开发技术 | PHP框架 | 开发工具 | 在线工具
Copyright © 1998 - 2020 PHP1.CN. All Rights Reserved | 京公网安备 11010802041100号 | 京ICP备19059560号-4 | PHP1.CN 第一PHP社区 版权所有