ACM-ICPC/CCPC/JSCPC板子总结（不断更新ing）

ACM-ICPC/CCPC/JSCPC板子

Team: Three alchemists（三个炼金师）
@zhoubo,18CS,Suzhou University of Science and Technology
@huangyangbang,20ICS,Suzhou University of Science and Technology
@shiweichun,20CS,Suzhou University of Science and Technology

一、数论

1.试除法判定质数

bool is_prime(int x)
{
    if (x < 2) return false;
    for (int i = 2; i <= x / i; i ++ )
        if (x % i == 0)
            return false;
    return true;
}

2.试除法分解质因数

void divide(int x)
{
    for (int i = 2; i <= x / i; i ++ )
        if (x % i == 0)
        {
            int s = 0;
            while (x % i == 0) x /= i, s ++ ;
            cout << i << ' ' << s << endl;
        }
    if (x > 1) cout << x << ' ' << 1 << endl;
    cout << endl;
}

3.朴素筛法求素数

int primes[N], cnt;     // primes[]存储所有素数
bool st[N];         // st[x]存储x是否被筛掉
void get_primes(int n)
{
    for (int i = 2; i <= n; i ++ )
    {
        if (st[i]) continue;
        primes[cnt ++ ] = i;
        for (int j = i + i; j <= n; j += i)
            st[j] = true;
    }
}

4.线性筛求素数

int primes[N], cnt;     // primes[]存储所有素数
bool st[N];         // st[x]存储x是否被筛掉
void get_primes(int n)
{
    for (int i = 2; i <= n; i ++ )
    {
        if (!st[i]) primes[cnt ++ ] = i;
        for (int j = 0; primes[j] <= n / i; j ++ )
        {
            st[primes[j] * i] = true;
            if (i % primes[j] == 0) break;
        }
    }
}

5.试除法求所有约数

vector<int> get_divisors(int x)
{
    vector<int> res;
    for (int i = 1; i <= x / i; i ++ )
        if (x % i == 0)
        {
            res.push_back(i);
            if (i != x / i) res.push_back(x / i);
        }
    sort(res.begin(), res.end());
    return res;
}

6.约数个数和约数之和

如果 N = p1^c1 * p2^c2 * ... *pk^ck
约数个数： (c1 + 1) * (c2 + 1) * ... * (ck + 1)
约数之和： (p1^0 + p1^1 + ... + p1^c1) * ... * (pk^0 + pk^1 + ... + pk^ck)

7.欧几里德算法

int gcd(int a, int b) {return b ? gcd(b, a % b) : a;}

8.求欧拉函数

$$欧拉函数公式：\varphi (x)=x\prod _{i=1}^n(1-\frac{1}{pi})$$

int phi(int x)
{
    int res = x;
    for (int i = 2; i <= x / i; i ++ )
        if (x % i == 0)
        {
            res = res / i * (i - 1);
            while (x % i == 0) x /= i;
        }
    if (x > 1) res = res / x * (x - 1);

    return res;
}

9.线性筛求欧拉函数

int primes[N], cnt;     // primes[]存储所有素数
int euler[N];           // 存储每个数的欧拉函数
bool st[N];         // st[x]存储x是否被筛掉
void get_eulers(int n)
{
    euler[1] = 1;
    for (int i = 2; i <= n; i ++ )
    {
        if (!st[i])
        {
            primes[cnt ++ ] = i;
            euler[i] = i - 1; 
        }
        for (int j = 0; primes[j] <= n / i; j ++ )
        {
            int t = primes[j] * i;
            st[t] = true;
            if (i % primes[j] == 0)
            {
                euler[t] = euler[i] * primes[j];
                break;
            }
            euler[t] = euler[i] * (primes[j] - 1);
        }
    }
}

10.快速幂

//求a^k mod p
typedef long long LL;
int qmi(int a, int k, int p)
{
    int res = 1 % p;
    while(k)
    {	
        if(k & 1) res = (LL)res * a % p;
        k >>= 1;
        a = (LL)a * a % p;
    }
    return res;
}c

11.费马小定理

$$如果p是一个质数，而整数a不是p的倍数，则有a^{p-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}m)$$

12.逆元

$$\begin{gathered} b\ast x\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p) \\ b存在乘法逆元的充要条件是b与模数p互质。当模数m为质数时，b^{p-2}即为b的乘法逆元 \\ b是p的倍数的时候，显然这个式子无解 \end{gathered}$$

13.裴蜀定理

$$对于任意正整数a,b,一定存在非零整数x,y,使得ax+by=gcd(a,b)$$

14.扩展欧几里德

//求x,y使得ax + by = gcd(a,b)
int exgcd(int a, int b, int &x, int &y)
{
    if(!b)
    {
        x = 1, y = 0;
        return a;
    }
    int d = exgcd(b, a % b, y, x);
    y -= (a / b) * x;
   return d;
}

$$可用来求线性同余方程ax\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m)$$

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
int exgcd(int a, int b, int &x, int &y) {
   if(!b) {
       x = 1, y = 0;
       return a;
   }
   int d = exgcd(b, a % b, y, x);
   y -= a / b * x;
   return d;
}
int main() {
    int n;
    scanf("%d", &n);
    while(n --) {
        int a, b, m;
        scanf("%d%d%d", &a, &b, &m);
        int x, y;
        int d = exgcd(a, m, x, y);
        if(b % d) puts("impossible");
        else printf("%d\n", (LL) x * (b / d) % m);
    }
    return 0;
}

16.高斯消元

// a[N][N]是增广矩阵
int gauss()
{
    int c, r;
    for (c = 0, r = 0; c < n; c ++ )
    {
        int t = r;
        for (int i = r; i < n; i ++ )   // 找到绝对值最大的行
            if (fabs(a[i][c]) > fabs(a[t][c]))
                t = i;
        if (fabs(a[t][c]) < eps) continue;
        for (int i = c; i <= n; i ++ ) swap(a[t][i], a[r][i]);      // 将绝对值最大的行换到最顶端
        for (int i = n; i >= c; i -- ) a[r][i] /= a[r][c];      // 将当前上的首位变成1
        for (int i = r + 1; i < n; i ++ )       // 用当前行将下面所有的列消成0
            if (fabs(a[i][c]) > eps)
                for (int j = n; j >= c; j -- )
                    a[i][j] -= a[r][j] * a[i][c];
        r ++ ;
    }
    if (r < n)
    {
        for (int i = r; i < n; i ++ )
            if (fabs(a[i][n]) > eps)
                return 2; // 无解
        return 1; // 有无穷多组解
    }
    for (int i = n - 1; i >= 0; i -- )
        for (int j = i + 1; j < n; j ++ )
            a[i][n] -= a[i][j] * a[j][n];
    return 0; // 有唯一解
}

二、组合数学

1.递归法求组合数

$$n组询问，每组询问包括一组a和b,1\le n\le 10000,1\le b\le a\le 2000O(n^2)$$

const int mod = 1e9 + 7;
// c[a][b] 表示从a个苹果中选b个的方案数
for (int i = 0; i < N; i ++ )
    for (int j = 0; j <= i; j ++ )
        if (!j) c[i][j] = 1;
        else c[i][j] = (c[i - 1][j] + c[i - 1][j - 1]) % mod;

2.通过预处理逆元的方式求组合数

$$n组询问，每组询问包括一组a和b,1\le n\le 10000,1\le b\le$$