嗯...

题目链接:https://www.luogu.org/problem/P1082

这道题很明显涉及到了同余和exgcd的问题,下面推导一下:

首先证明有解情况:

    ax + by = m有解的必要条件是 m mod gcd(a, b) = 0

    a为gcd(a, b)的倍数,b为gcd(a, b)的倍数,x、y为整数,

    所以ax + by是gcd(a, b)的倍数,所以m是gcd(a, b)的倍数

然后证明a、b互质(下面会用到):

    本题中1 mod gcd(a, b) = 0,所以gcd(a, b)  = 1,所以a、b互质

同余:

    a≡b(mod n) --> 含义:a和b关于模n同余,即 a mod n = b mod n。

    所以不难推出,a≡b的充要条件是a-b是n的整数倍(a > b)。

    因为a-b是n的整数倍,所以a-b = ny(y为倍数)

所以,根据同余,我们可以把本题中的同余式转化为(注:这里的a.b与上文不同):

    ax≡1(modb)  --> ax % b = 1 % b --> ax - 1 = by --> ax - by = 1

下一步,便进行exgcd(关于exgcd的证明见:https://www.luogu.org/problemnew/solution/P1082),分别求出ax - by = 1中x和y的值。

最后进行答案处理:

    因为答案要求是x的最小正整数,所以我们进行一个答案处理:x = (x % b + b) % b

证明其正确性:

    设新得到的x为xn,x = kb + q(q < b)则:

    x % b = q ,把x % b = q带入 xn = (x % b + b) % b,得

    xn = (x % b + b) % b = (q + b) % b = (q % b + b % b) % b = q % b = q

     把xn = q带入x = kb + q,得,x = kb + xn, 所以xn = x - kb,然后再根据下面的推导得知它是正确的...

证明:    

    x批量地减去或加上 b,能保证 ax + by = ax + by = 1:

    ax + by = 1

    ax + by + k*ba - k*ba = 1

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

    1.显然这并不会把方程中任何整数变成非整数。

    2.“加上或减去 b”不会使 x 错过任何解。可以这么理解:

    已经求出一组整数 x,y 使得 ax+by =1 ,也就是(1 - ax) / b = y。y 是整数,可见目前 1−ax 是 b 的倍数。现在想改变 x并使得方程仍然成立。

    已知 a,b 互质,假若x的变化量Δx不是b的倍数,则1−ax 的变化量−a*Δx也不是 bb 的倍数,这会使得1-ax不再是b的倍数,则y不是整数了。

    仅当x的变化量是b的倍数时,1−ax能保持自己是b的倍数,此时就出现新的解了。

AC代码:

 #include<cstdio>
#include<iostream> using namespace std;
//ax % b == 1 % b --> ax - 1 = y * b --> ax - yb == 1 long long d, x, y; inline void exgcd(long long a, long long b, long long &d, long long &x, long long &y){
if(!b) {d = a; x = ; y = ;}
else{ exgcd(b, a % b, d, y, x); y -= x * (a / b);}
}
int main(){
long long a,b;
scanf("%lld%lld", &a, &b);
exgcd(a, b, d, x, y);
x = (x % b + b) % b;
printf("%lld", x);
return ;
}

AC代码

洛谷 P1082 同余方程(同余&&exgcd)的更多相关文章

  1. 洛谷——P1082 同余方程

    P1082 同余方程 题目描述 求关于 x 的同余方程 ax ≡ 1 (mod b)的最小正整数解. 输入输出格式 输入格式: 输入只有一行,包含两个正整数 a, b,用一个空格隔开. 输出格式: 输 ...

  2. 洛谷P1082 同余方程 [2012NOIP提高组D2T1] [2017年6月计划 数论06]

    P1082 同余方程 题目描述 求关于 x 的同余方程 ax ≡ 1 (mod b)的最小正整数解. 输入输出格式 输入格式: 输入只有一行,包含两个正整数 a, b,用一个空格隔开. 输出格式: 输 ...

  3. 洛谷 P1082 同余方程 —— exgcd

    题目:https://www.luogu.org/problemnew/show/P1082 用 exgcd 即可. 代码如下: #include<iostream> #include&l ...

  4. 洛谷 P1082 同余方程(exgcd)

    题目传送门 解题思路: 因为推导过程过于复杂,懒得写,所以题解传送门 AC代码: #include<iostream> #include<cstdio> using names ...

  5. [NOIP2012] 提高组 洛谷P1082 同余方程

    题目描述 求关于 x 的同余方程 ax ≡ 1 (mod b)的最小正整数解. 输入输出格式 输入格式: 输入只有一行,包含两个正整数 a, b,用一个空格隔开. 输出格式: 输出只有一行,包含一个正 ...

  6. 洛谷P1082 同余方程

    题目描述 求关于 x 的同余方程 ax ≡ 1 (mod b)的最小正整数解. 输入输出格式 输入格式: 输入只有一行,包含两个正整数 a, b,用一个空格隔开. 输出格式: 输出只有一行,包含一个正 ...

  7. 洛谷 P1082 同余方程

    题目描述 求关于 x 的同余方程 ax ≡ 1 (mod b)的最小正整数解. 输入输出格式 输入格式: 输入只有一行,包含两个正整数 a, b,用一个空格隔开. 输出格式: 输出只有一行,包含一个正 ...

  8. 洛谷P1082 同余方程 题解

    题目链接:https://www.luogu.com.cn/problem/P1082 题目大意: 求关于 \(x\) 的同余方程 ax≡1(mod b) 的最小正整数解. 告诉你 \(a,b\) 求 ...

  9. 洛谷 P1082 同余方程 题解

    每日一题 day31 打卡 Analysis 题目问的是满足 ax mod b = 1 的最小正整数 x.(a,b是正整数) 但是不能暴力枚举 x,会超时. 把问题转化一下.观察 ax mod b = ...

随机推荐

  1. php preg正则表达式的组成部分

    定界符号 : 多种都可以,常用为// 原子 : 最小的一个匹配单位 (放在定界符中),在一个正则表达式中,至少要有一个原子 1,打印字符(a-z A-Z 0-9 ~!@#$%^&*()_+.. ...

  2. python 第三方库安装

    1.首先安装pip 2.在cmd中找到pip的安装路径,(一般在python的scripts文件中) 3.pip install 第三方库名称

  3. 交换机的MAC地址?

    该示例中记录了Cat2960的MAC地址情况: Cat2960#sho mac address-table Mac Address Table----------------------------- ...

  4. Codeforces Round #624 (Div. 3) A. Add Odd or Subtract Even(水题)

    You are given two positive integers aa and bb . In one move, you can change aa in the following way: ...

  5. 实战【docker 镜像制作与使用】

    一.制作docker 镜像 使用spring boot 构建一个简单的web 项目,返回 “Hello,World ”字符串,使用 Maven 打成 jar 包,使用的Linux 环境是 Centos ...

  6. 多租户SaaS平台的数据库方案

    1.1 多租户是什么 多租户技术(Multi-TenancyTechnology)又称多重租赁技术:是一种软件架构技术,是实现如何在多用户环境下 (此处的多用户一般是面向企业用户)共用相同的系统或程序 ...

  7. pwnable.kr-cmd1-Writeup

    MarkdownPad Document *:first-child { margin-top: 0 !important; } body>*:last-child { margin-botto ...

  8. django model 高级进阶

    十.model高级用法: 10.1 ORM映射: Object Relational Mapping: orm映射的任务: 将list ,object 的操作转化为sql语句,根据对象生成数据表,将查 ...

  9. MindManager安装、激活

    MindManager安装 MindManager激活

  10. acm数论之旅--唯一分解定理

    题目: 给出n,问n = b^p中p符合该等式的最大值 分析: 先求出所有n的质因子,然后对这m个质因子分类统计,比如 n = 36时,可以分成 2个2,2个3,然后求出所有这些基数的 最大公因数gc ...