【板子】gcd、exgcd、乘法逆元、快速幂、快速乘、筛素数、快速求逆元、组合数
1.gcd
int gcd(int a,int b){
return b?gcd(b,a%b):a;
}
2.扩展gcd )extend great common divisor
ll exgcd(ll l,ll r,ll &x,ll &y)
{
if(r==){x=;y=;return l;}
else
{
ll d=exgcd(r,l%r,y,x);
y-=l/r*x;
return d;
}
}
3.求a关于m的乘法逆元
ll mod_inverse(ll a,ll m){
ll x,y;
if(exgcd(a,m,x,y)==)//ax+my=1
return (x%m+m)%m;
return -;//不存在
}
补充:求逆元还可以用$$ans = \frac{a}{b} \bmod m = (a \bmod (m\cdot b)) /b $$
4.快速幂quick power
ll qpow(ll a,ll b,ll m){
ll ans=;
ll k=a;
while(b){
if(b&)ans=ans*k%m;
k=k*k%m;
b>>=;
}
return ans;
}
5.快速乘,直接乘会爆ll时需要它,也叫二分乘法。
ll qmul(ll a,ll b,ll m){
ll ans=;
ll k=a;
ll f=;//f是用来存负号的
if(k<){f=-;k=-k;}
if(b<){f*=-;b=-b;}
while(b){
if(b&)
ans=(ans+k)%m;
k=(k+k)%m;
b>>=;
}
return ans*f;
}
6.中国剩余定理CRT (x=ai mod mi)
ll china(ll n, ll *a,ll *m) {
ll M=,y,x=,d;
for(ll i = ; i <= n; i++) M *= m[i];
for(ll i = ; i <= n; i++) {
ll w = M /m[i];
exgcd(m[i], w, d, y);//m[i]*d+w*y=1
x = (x + y*w*a[i]) % M;
}
return (x+M)%M;
}
7.筛素数,全局:int cnt,prime[N],p[N];
void isprime()
{
cnt = ;
memset(prime,true,sizeof(prime));
for(int i=; i<N; i++)
{
if(prime[i])
{
p[cnt++] = i;
for(int j=i+i; j<N; j+=i)
prime[j] = false;
}
}
}
8.快速计算逆元
补充:>>关于快速算逆元的递推式的证明<<
void inverse(){
inv[] = ;
for(int i=;i<N;i++)
{
if(i >= M) break;
inv[i] = (M-M/i)*inv[M%i]%M;
}
}
9.组合数取模
n和m 10^5时,预处理出逆元和阶乘
ll fac[N]={,},inv[N]={,},f[N]={,};
ll C(ll a,ll b){
if(b>a)return ;
return fac[a]*inv[b]%M*inv[a-b]%M;
}
void init(){//快速计算阶乘的逆元
for(int i=;i<N;i++){
fac[i]=fac[i-]*i%M;
f[i]=(M-M/i)*f[M%i]%M;
inv[i]=inv[i-]*f[i]%M;
}
}
n较大10^9,但是m较小10^5时,
ll C(ll n,ll m){
if(m>n)return ;
ll ans=;
for(int i=;i<=m;i++)
ans=ans*(n-i+)%M*qpow(i,M-,M)%M;
return ans;
}
n和m特别大10^18时但是p较小10^5时用lucas
10.Lucas大组合取模
#define N 100005
#define M 100007
ll n,m,fac[N]={};
ll C(ll n,ll m){
if(m>n)return ;
return fac[n]*qpow(fac[m],M-,M)%M*qpow(fac[n-m],M-,M)%M;//费马小定理求逆元
}
ll lucas(ll n,ll m){
if(!m)return ;
return(C(n%M,m%M)*lucas(n/M,m/M))%M;
}
void init(){
for(int i=;i<=M;i++)
fac[i]=fac[i-]*i%M;
}
to be continued...
【板子】gcd、exgcd、乘法逆元、快速幂、快速乘、筛素数、快速求逆元、组合数的更多相关文章
- hdu-5667 Sequence(矩阵快速幂+费马小定理+快速幂)
题目链接: Sequence Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others) ...
- 【BZOJ 2186】 2186: [Sdoi2008]沙拉公主的困惑 (欧拉筛,线性求逆元)
2186: [Sdoi2008]沙拉公主的困惑 Description 大富翁国因为通货膨胀,以及假钞泛滥,政府决定推出一项新的政策:现有钞票编号范围为1到N的阶乘,但是,政府只发行编号与M!互质的钞 ...
- 51Nod 3的幂的和(扩展欧几里德求逆元)
求:3^0 + 3^1 +...+ 3^(N) mod 1000000007 Input 输入一个数N(0 <= N <= 10^9) Output 输出:计算结果 Input示例 3 O ...
- 整数快速乘法/快速幂+矩阵快速幂+Strassen算法
快速幂算法可以说是ACM一类竞赛中必不可少,并且也是非常基础的一类算法,鉴于我一直学的比较零散,所以今天用这个帖子总结一下 快速乘法通常有两类应用:一.整数的运算,计算(a*b) mod c 二.矩 ...
- 矩阵快速幂(入门) 学习笔记hdu1005, hdu1575, hdu1757
矩阵快速幂是基于普通的快速幂的一种扩展,如果不知道的快速幂的请参见http://www.cnblogs.com/Howe-Young/p/4097277.html.二进制这个东西太神奇了,好多优秀的算 ...
- P3390 【模板】矩阵快速幂
题目背景 矩阵快速幂 题目描述 给定n*n的矩阵A,求A^k 输入输出格式 输入格式: 第一行,n,k 第2至n+1行,每行n个数,第i+1行第j个数表示矩阵第i行第j列的元素 输出格式: 输出A^k ...
- poj3613 Cow Relays【好题】【最短路】【快速幂】
Cow Relays Time Limit: 1000MS Memory Limit: 65536K Total Submissions:9207 Accepted: 3604 Descrip ...
- [C/C++] 快速幂讲解
转自:http://www.cnblogs.com/CXCXCXC/p/4641812.html 快速幂这个东西比较好理解,但实现起来到不老好办,记了几次老是忘,今天把它系统的总结一下防止忘记. 首先 ...
- [学习笔记]快速幂&&快速乘
本质:二进制拆分(你说倍增我也没脾气).然后是一个配凑. 合起来就是边二进制拆分,边配凑. 快速乘(其实龟速):把乘数二进制拆分.利用乘法分配率. 用途:防止爆long long 代码: ll qk( ...
- 矩阵快速幂--51nod-1242斐波那契数列的第N项
斐波那契额数列的第N项 斐波那契数列的定义如下: F(0) = 0 F(1) = 1 F(n) = F(n - 1) + F(n - 2) (n >= 2) (1, 1, 2, 3, 5, 8, ...
随机推荐
- jsp的九大内置对象和四大作用域
定义:可以不加声明就在JSP页面脚本(Java程序片和Java表达式)中使用的成员变量? JSP共有以下9种基本内置组件(可与ASP的6种内部组件相对应):? 1.request对象(作用域)? 客户 ...
- CSS3之文本阴影text-shadow
- PAT 1003. 我要通过!(20)
"答案正确"是自动判题系统给出的最令人欢喜的回复.本题属于PAT的"答案正确"大派送 -- 只要读入的字符串满足下列条件,系统就输出"答案正确&quo ...
- 表单验证——JavaScript和Jquery版
1.轻量级的JavaScript表单验证 在应用中引用 validator.min.js 文件 <script type="text/javascript" src=&quo ...
- usb驱动开发22之驱动生命线
我们总是很喜欢高潮,不是吗?那就好好对待她哦.我们来看一下linux中的高潮部分设备是怎么从Address进入Configured的. usb_set_configuration函数的代码就不贴了,可 ...
- 解决Ehcache缓存警告问题
警告: Creating a new instance of CacheManager using the diskStorePath "D:\Apache Tomcat 6.0.18\te ...
- Gruntjs: grunt-contrib-jst
预编译Underscore模板到JST文件(Underscore:JS工具库) generate JavaScript template functions Gruntfile的配置实例: modul ...
- TinyFrame升级之三:逻辑访问部分
在上一篇,我们打造了自己的数据访问部分,这篇,我们准备讲解如何打造逻辑访问部分. 在上一篇中,我们利用Repository模式构建了基于泛型的操作合集.由于这些操作的合集都是原子性的操作,也就是针对单 ...
- JavaScript函数劫持
一.为什么我会写这篇文章 这篇文章其实是在一个偶然的机会下发现了居然有JavaScript劫持这种东西,虽然这种东西在平时用的比较少,而且一般实用价值不高,但是在一些特殊的情况下还是要使用到的,所以在 ...
- Android手机截屏
刚开始打算做一个简单的截屏程序时,以为很轻松就能搞定. 在Activity上放一个按钮,点击完成截屏操作,并将数据以图片形式保存在手机中. 动手之前,自然是看书和网上各种查资料.结果发现了解的知识越多 ...