P1082 [NOIP2012 提高组] 同余方程

合集 - 刷题(27)

1.uva253骰子涂色2023-11-302.洛谷uva2202023-11-303.uva12096集合栈计算机 The SetStack Computer2023-12-054.uva101The Blocks Problem2023-12-055.uva400 Unix Is命令 Unix ls2023-12-106.uva10391 复合词 Compound Words2023-12-107.小球下落 Dropping Balls （uva679）2023-12-198.树的层次遍历 Trees on the level uva1222023-12-199.树 Tree uva5482023-12-1910.天平 Not so Mobile uva8392023-12-1911.P1102 A-B 数对2024-02-1012.成最多水的容器2024-02-1213.P3367 【模板】并查集2024-02-1314.947. 移除最多的同行或同列石头2024-02-14

15.P1082 [NOIP2012 提高组] 同余方程2024-02-22

16.803 打砖块2024-02-2317.B3644 【模板】拓扑排序 / 家谱树2024-02-2318.U107394 拓扑排序模板2024-02-2319.P1137 旅行计划2024-02-2420.936 戳印序列2024-02-2721.P4017 最大食物链计数2024-02-2722.2050 并行课程 III2024-02-2723.P3366 【模板】最小生成树2024-03-0124.P2330 [SCOI2005] 繁忙的都市2024-03-0425.Abbott的复仇 Abbott's Revenge2024-03-0626.P1960 郁闷的记者2024-03-2227.P4391 [BOI2009] Radio Transmission 无线传输2024-05-03

收起

原题链接

扩展欧几里得算法的应用，关于原理性的讲解这里就略去了，这边给出学习链接即模板。

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

文章链接      视频链接

首先原题方程可以化简为 ax=by+1 即 ax-by=1（正负无所谓）；

那么根据裴蜀定理我们可以得到得到该方程的一个特解，该解不一定是题目要求的最小正整数解。

此时我们对x批量的进行b的增减，原因如下：

ax+by=1;

ax+by+k*ab-k*ab=1;

a(x-kb)+b(y+ka)=1;

显然此时x-kb，y+ka仍为整数。

因此为了保证x0为最小整数，我们进行如下操作：

x=（x%b+b）%b

code

 

#include<bits/stdc++.h>
using namespace std;
int exgcd(int a,int b,int &x,int &y){
    if (b==0){
        x=1;
        y=0;
        return x;
    }
    int d=exgcd(b,a%b,x,y);
    x=y;
    y=d-a/b*y;
    return x;
}
int main(){
    int x,y,a,b;
    cin>>a>>b;
    int m=exgcd(a,b,x,y);
    cout<<(m%b+b)%b;
    return 0;
}