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、乘法逆元、快速幂、快速乘、筛素数、快速求逆元、组合数的更多相关文章

  1. hdu-5667 Sequence(矩阵快速幂+费马小定理+快速幂)

    题目链接: Sequence Time Limit: 2000/1000 MS (Java/Others)     Memory Limit: 65536/65536 K (Java/Others) ...

  2. 【BZOJ 2186】 2186: [Sdoi2008]沙拉公主的困惑 (欧拉筛,线性求逆元)

    2186: [Sdoi2008]沙拉公主的困惑 Description 大富翁国因为通货膨胀,以及假钞泛滥,政府决定推出一项新的政策:现有钞票编号范围为1到N的阶乘,但是,政府只发行编号与M!互质的钞 ...

  3. 51Nod 3的幂的和(扩展欧几里德求逆元)

    求:3^0 + 3^1 +...+ 3^(N) mod 1000000007 Input 输入一个数N(0 <= N <= 10^9) Output 输出:计算结果 Input示例 3 O ...

  4. 整数快速乘法/快速幂+矩阵快速幂+Strassen算法

    快速幂算法可以说是ACM一类竞赛中必不可少,并且也是非常基础的一类算法,鉴于我一直学的比较零散,所以今天用这个帖子总结一下 快速乘法通常有两类应用:一.整数的运算,计算(a*b) mod c  二.矩 ...

  5. 矩阵快速幂(入门) 学习笔记hdu1005, hdu1575, hdu1757

    矩阵快速幂是基于普通的快速幂的一种扩展,如果不知道的快速幂的请参见http://www.cnblogs.com/Howe-Young/p/4097277.html.二进制这个东西太神奇了,好多优秀的算 ...

  6. P3390 【模板】矩阵快速幂

    题目背景 矩阵快速幂 题目描述 给定n*n的矩阵A,求A^k 输入输出格式 输入格式: 第一行,n,k 第2至n+1行,每行n个数,第i+1行第j个数表示矩阵第i行第j列的元素 输出格式: 输出A^k ...

  7. poj3613 Cow Relays【好题】【最短路】【快速幂】

    Cow Relays Time Limit: 1000MS   Memory Limit: 65536K Total Submissions:9207   Accepted: 3604 Descrip ...

  8. [C/C++] 快速幂讲解

    转自:http://www.cnblogs.com/CXCXCXC/p/4641812.html 快速幂这个东西比较好理解,但实现起来到不老好办,记了几次老是忘,今天把它系统的总结一下防止忘记. 首先 ...

  9. [学习笔记]快速幂&&快速乘

    本质:二进制拆分(你说倍增我也没脾气).然后是一个配凑. 合起来就是边二进制拆分,边配凑. 快速乘(其实龟速):把乘数二进制拆分.利用乘法分配率. 用途:防止爆long long 代码: ll qk( ...

  10. 矩阵快速幂--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, ...

随机推荐

  1. 转: 使用Eclipse的Working Set,界面清爽多了

    from: http://iyuanbo.iteye.com/blog/1158136 使用Eclipse的Working Set,界面清爽多了   想必大家的Eclipse里也会有这么多得工程... ...

  2. box-shadow 的一些使用

    1.只有左侧有阴影 box-shadow: -10px 0px 3px 1px #aaaaaa;

  3. 见鬼了,swiper

    1.今天不知怎么swiper的onInit函数不起作用,怎么弄都不行: 把以前能行的案例的包都导进去还是不行,但是onSlideChangeEnd可以触发,晕死了.... 不,它触发了一次onInit ...

  4. PAT 1013. 数素数 (20)

    令Pi表示第i个素数.现任给两个正整数M <= N <= 104,请输出PM到PN的所有素数. 输入格式: 输入在一行中给出M和N,其间以空格分隔. 输出格式: 输出从PM到PN的所有素数 ...

  5. React问题总结与归纳

    欢迎大家指导与讨论 : ) [持续更新]本文主要记录笔者在学习中遇到的问题,并作出相应总结.有错误的地方希望各位能够指出. 一.在es6中getInitialState( 摘要:  construct ...

  6. (原创)解决Excel 互操作错误"检索COML类工厂中 CLSID为 {00024500-0000-0000-C000-000000000046}的组件时失败,原因是出现以下错误: 80070005"

    最近在.net中处理Excel文件数据导入时报出以下错误: 检索COML类工厂中 CLSID为 {00024500-0000-0000-C000-000000000046}的组件时失败,原因是出现以下 ...

  7. Qt Creator 常用快捷键

    多行注释模式                                                                                            Ct ...

  8. 是什么时候开始学习gulp了

    转自:http://www.ydcss.com/archives/18 简介: gulp是前端开发过程中对代码进行构建的工具,是自动化项目的构建利器:她不仅能对网站资源进行优化,而且在开发过程中很多重 ...

  9. .NET MVC AjaxHelper

    我们首先必须开启 非入侵式 Ajax:导入Jquery和unobtrusiveAjax文件 已经默认开启客户端验证 和 非侵入式js <add key="ClientValidatio ...

  10. js队列

    用指针和数组模拟基本队列 http://blog.csdn.net/zhuwq585/article/details/53177192 js下的事件队列,或者异步队列 http://www.jb51. ...