欧几里得算法

就是求最大公约数的辗转相除法。

数学公式

$$\gcd(a, b)= \begin{cases} \gcd(b,a\bmod b) &,b\neq 0\\ a &,b=0 \end{cases}$$

模板

1
2
3
4
5

int gcd(int a,int b)
{
    if(b==0) return a;
    return gcd(b,a%b);
}

拓展欧几里得(exgcd)

适用问题

给定正整数 $a,b$，求出 $ax+by=\gcd(a,b)$ 的一组正整数解。

算法思路

由欧几里得算法得出 $ax+by=\gcd(a,b)=\gcd(b,a\bmod b)$，代回原式，即为：

$$bx^{\prime}+a\bmod b\times y^{\prime}=\gcd(b,a\bmod b)$$

假设我们已知上式的一组解 $x^{\prime},y^{\prime}$，可以得到：

$$\begin{aligned} ax+by=\gcd(a, b) &=b x^{\prime}+a \bmod b\times y^{\prime} \\ &=b x^{\prime}+(a-\lfloor \frac{a}{b} \rfloor \times b) y^{\prime}\\ &=a y^{\prime}+b(x^{\prime}-\lfloor \frac{a}{b}\rfloor \times y^{\prime}) \end{aligned}$$

所以

$$\begin{aligned} ax+by&=ay^{\prime}+b(x^{\prime}-\lfloor \frac{a}{b}\rfloor \times y^{\prime})\\ x&=y^{\prime}\\ y&=x^{\prime}-\frac{a}{b}\times y^{\prime} \end{aligned}$$

我们发现，当 $b=0$ 时，$a\times 1+b\times 0=a$，有一组解

$$\begin{cases} x=1\\ y=0 \end{cases}$$

模板

直接递归求解即可

1
2
3
4
5
6
7
8
9
10
11
12

int exgcd(int a,int b,int &x,int &y)
{
    if(b==0)
    {
        x=1;
        y=0;
        return a;
    }
    int g=exgcd(b,a%b,y,x);
    y=y-a/b*x;
    return g;
}

这时，我们已经得到了一组解，下一步是求解 $ax+by=c$。

求解 $ax+by=c$

接下来我们需要分类讨论。

设 $g=\gcd(a,b)$。

当 $c\bmod b\neq0$ 时，无解。

否则，令 $k=\frac{c}{g}$，

则

$$ax^{\prime}k+by^{\prime}k=g\times k=c$$

得

$$\begin{cases} x=x^{\prime}k\\ y=y^{\prime}k \end{cases}$$

设 $x$ 的周期 $T=\frac{b}{g}$，则最小正整数解为 $(x\bmod T+T)\bmod T$。

乘法逆元

定义

已知 $a,p$，求 $ax\equiv 1\pmod{p}$ 中的 $x$，称 $x$ 为 $a$ 关于 $p$ 的逆元。

求解

直接求解 $ax+py=1$ 即可。

最终答案为 $(x\bmod p+p)\bmod p$。

费马小定理

当 $p$ 为素数时，$a$ 关于 $p$ 的逆元是 $a^{p-2}$。

数论分块

balabala

$a\bmod b=a-b\times\lfloor\frac{a}{b}\rfloor$

一些题

P1516 青蛙的约会

Problem Link

这是一道比较裸的拓欧题。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29

#define int ll
int x,y,m,n,l;
int p,q;
int exgcd(int a,int b,int &x,int &y) {
    if(b==0) {
        x=1;
        y=0;
        return a;
    }
    int g=exgcd(b,a%b,y,x);
    y=y-a/b*x;
    return g;
}
signed main() {
    x=read();y=read();m=read();n=read();l=read();
    int a=n-m,b=x-y;
    if(a<0) {
        a=-a;
        b=-b;
    }
    int ans=exgcd(a,l);
    if(b%ans!=0) {
        cout<<"Impossible"<<endl;
    }
    else {
        cout<<((p*(b/ans))%(l/ans)+l/ans)%(l/ans)<<endl;
    }
    return 0;
}

CF1114C Trailing Loves

首先从 $10$ 进制开始分析。求 $p$ 末尾 $0$ 的数量即为求 $p$ 的因数中含有 $2$ 和 $5$ 的数量的最小值。

例如设 $x_1=2^3\times 5^2$，则 $x_1$ 的末尾 $0$ 的数量即为 $\min(3,2)=2$。

接下来考虑 $(n!)_b$ 的末尾 $0$ 的数量。我们先枚举 $b$ 的因数 $b=p_1^{a_1}\times p_2^{a_2}\times \ldots \times p_k^{a_k}$，若 $p_i^{m}$ 是 $n!$ 的一个因数，则其对答案的贡献为 $\lfloor \dfrac{m}{a_i}\rfloor$。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26

const int N=2e6+10;
int n,b;
int pri[N],cnt=0;
int t[N];
int ans=INF;
signed main() {
    n=read();b=read();
    for(int i=2;i*i<=b;i++) {
        if(b%i==0) {
            pri[++cnt]=i;
            while(b%i==0) b/=i,t[cnt]++;
        }
    }
    if(b>1) pri[++cnt]=b,t[cnt]++;
    for(int i=1;i<=cnt;i++) {
        int nt=n;
        int sum=0;
        while(nt) {
            sum+=nt/pri[i];
            nt/=pri[i];
        }
        ans=min(sum/t[i],ans);
    }
    printf("%lld\n",ans);
    return 0;
}

CF1493D GCD of an Array

由于一个集合的最大公因数的一个质因数的指数是这个集合中该质数指数的最小值，可以想到每次动态维护集合 gcd 时，考虑将 $y$ 分解质因数 $y=p_1^{e_1}\times p_2^{e_2}\times \ldots \times p_k^{e_k}$，如果集合中除 $a_x$ 以外的其他元素均有可以对答案产生贡献的因数 $p_i$，那么答案就可以加入 $p_i$ 产生的贡献。

形式化的，如果我们定义 $v_p(a_i)$ 为 $a_i$ 中质因数 $p$ 的指数，那么若满足

$$\min(v_p(a_1), \ldots, v_p(a_{x-1}), v_p(y), v_p(a_{x+1}), \ldots, v_p(a_n)) > \min(v_p(a_1), v_p(a_2), \ldots, v_p(a_n))$$

设 $A$ 为上述两式的差值，则对答案的贡献为 $p^A$。

接下来考虑如何求。可以用线性筛预处理质数，每次操作用预处理的最小质因子分解。

我们用 multiset 记录质数 $p$ 出现在的元素下标，$t$ 数组记录 $p$ 在多少个不同的元素中出现。当执行一次操作时，判断是否能够给答案产生贡献。如果能，那么插入当前的下标，随后暴力将 multiset 中的 $1\sim n$ 下标删除一次。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50

const int N=2e5+10;
int n,m;
int a[N];
multiset<int> st[N];
int b[N];
int prime[N];
int cnt=0;
int t[N];
int ans=1;
void pri(int n) {
    for(int i=2;i<=n;i++) {   
        if(!b[i]) b[i]=i,prime[++cnt]=i;
        for(int j=1;j<=cnt;j++) {
            if(prime[j]>b[i] || prime[j]>n/i) break;
            b[i*prime[j]]=prime[j];
        }
    }
}
void modify(int p,int x) {
    while(x!=1) {
        int minn=b[x];
        if(st[minn].find(p)==st[minn].end()) t[minn]++;
        st[minn].insert(p);
        if(t[minn]==n) {
            for(int i=1;i<=n;i++) {
                st[minn].erase(st[minn].find(i));
                if(st[minn].find(i)==st[minn].end()) t[minn]--;
            }
            ans=(ans*minn)%mod;
        }
        x/=minn;
    }
}
signed main()
{
    pri(N-1);
    n=read();m=read();
    for(int i=1;i<=n;i++) {
        a[i]=read();
    }
    for(int i=1;i<=n;i++) {
        modify(i,a[i]);
    }
    while(m--) {
        int p=read(),x=read();
        modify(p,x);
        printf("%lld\n",ans);
    }
    return 0;
}