编个程序求二元一次不定方程最小正整数解

C/C++ > C++ 语言
收藏 回复
[问题点数：100分，结帖人jandy123]
更多帖子
私信 空间 博客
jandy123
等级
本版专家分：25
结帖率 100%
编个程序求二元一次不定方程最小正整数解
已知 ：ax+by=c
求:x最小的正整数解。
输入 a b c
输出 x y （或impossible)

(a,b,c都是整数)
只看楼主 引用 举报 楼主 收起
展开阅读全文
更多帖子
私信 空间 博客
iatsbg
等级
本版专家分：1298
用辗转相除法。参考http://community.csdn.net/Expert/topic/5734/5734611.xml?temp=.7709619
只看TA 引用 举报 #1得分 40
更多帖子
私信 空间 博客
冷月清晖
等级
本版专家分：10345
初等数论 先学习下
只看TA 引用 举报 #2得分 10
更多帖子
私信 空间 博客
syy64
等级
本版专家分：154
勋章
Blank
黄花 2007年3月 VC/MFC大版内专家分月排行榜第二
2007年2月 VC/MFC大版内专家分月排行榜第二
2007年1月 VC/MFC大版内专家分月排行榜第二
x = (c-by)/a

当Y大时，X就小，且c〉by
只看TA 引用 举报 #3得分 10
更多帖子
私信 空间 博客
expter
等级
本版专家分：1899
对于LS在枚举？
只看TA 引用 举报 #4得分 10
更多帖子
私信 空间 博客
lin_style
等级
本版专家分：3972
昨天刚出了个取模的。。
只看TA 引用 举报 #5得分 10
更多帖子
私信 空间 博客
jandy123
等级
本版专家分：25
#include <stdio.h>
#include <math.h>

long gcd(long a,long b)
{
long c=0,d=0;
if(a<b){c=a;a=b;b=c;}
while(b!=0)
{
c=a%b;
a=b;
b=c;
}
return a;
}

bool findans(long a,long b,long c,long *x1,long *y1)
{
if(a%b==0){
if((a-c)%b!=0)return 0;
else 
{
x1[0]=b;
y1[0]=(c-a)/b;
return 1;
}
}
if(c%gcd(a,b)!=0)return 0;
_int64 s[3]={0},t[3]={0},q[3]={0},x[3]={0},y[3]={0};
int k=0;
s[0]=a;t[0]=s[1]=b;t[1]=a%b;
q[0]=a/b;q[1]=s[1]/t[1];
x[0]=0;y[0]=1;
x[1]=1;y[1]=q[0];
while(t[1]!=0)
{
t[2]=s[1]%t[1];
s[0]=s[1];s[1]=t[0]=t[1];t[1]=t[2];
q[0]=q[1];
if(t[1]==0)break;
q[1]=s[1]/t[1];
x[2]=x[1]*q[0]+x[0];
y[2]=y[1]*q[0]+y[0];
x[0]=x[1];
y[0]=y[1];
x[1]=x[2];
y[1]=y[2];
k++;
}
x1[0]=(long)(pow(-1,k)*x[1]*c/gcd(a,b));
y1[0]=(long)(pow(-1,k-1)*y[1]*c/gcd(a,b));
return 1;
}

int main()
{
long a,b,c,x,y;
bool l=0,q=0;
scanf("%ld%ld%ld",&a,&b,&c);
if(a<0){l=1;a=-a;}
if(b<0){q=1;b=-b;}
findans(a,b,c,&x,&y);
if(l){x=-x;a=-a;}
if(q){y=-y;b=-b;}
if(x>0){y+=x/b*a;x=x%b;}
else if(x<0){y-=(1-(x/b))*a;x=(x%b)+b;}
if(x==0){x+=b;y-=a;}
printf("%d %d\n",x,y);
return 0;
}

这个程序不知行不？
只看TA 引用 举报 #6得分 0
更多帖子
私信 空间 博客
spofmy
等级
本版专家分：825
顶了
只看TA 引用 举报 #7得分 10
更多帖子
私信 空间 博客
jandy123
等级
本版专家分：25
#include <stdio.h>
#include <math.h>
#include <stdlib.h>

long gcd(long a,long b)
{
long c;
if(a<b){a+=b;b=a-b;a-=b;}
while(b!=0)
{
c=b;
b=a%b;
a=c;
}
return a;
}

long findans(long a,long b,long c)
{
long k=0;
if(c%gcd(a,b)!=0)
return -1;
if(a%b==0)
return 1;
bool l=0;
if(a<0){a=-a;l=1;}
if(b<0)b=-b;
_int64 s[3],t[3],r[3],x[3];
x[0]=0;
s[1]=b;t[1]=a%b;
r[1]=s[1]/t[1];
t[2]=s[1]%t[1];
x[1]=1;
while(t[2]!=0)
{
x[2]=x[1]*r[1]+x[0];
x[0]=x[1];
x[1]=x[2];
s[1]=t[1];
t[1]=t[2];
r[1]=s[1]/t[1];
t[2]=s[1]%t[1];
k++;
}
x[1]*=c/gcd(a,b);
if(l ==( k%2==0)){
x[1]=-x[1];
x[1]=x[1]%labs(b)+labs(b);
}
else x[1]%=labs(b);
if(x[1]==0)x[1]=labs(b);
return (long)x[1];
}

int main()
{
long a,b,c,x,y;
scanf("%ld%ld%ld",&a,&b,&c);
x=findans(a,b,c);
if(x==-1){printf("Impossible\n");return 0;}
y=-(a*x-c)/b;
printf("%ld %ld\n",x,y);
return 0;
}
这个算法改了一下。
只看TA 引用 举报 #8得分 0
更多帖子
私信 空间 博客
xiaciping
等级
本版专家分：20
学习
只看TA 引用 举报 #9得分 0
更多帖子
私信 空间 博客
chai2010
等级
本版专家分：1131
勋章
Blank
技术圈认证
Blank
GitHub
扩展的gcd算法啊
只看TA 引用 举报 #10得分 10
收起
收藏 回复
jandy123

等级：

关注 私信 TA的帖子
POJ 1061（不定方程的最小正整数解）

青蛙的约会 Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 108469   Accepted: 21745 ...它们很高兴地发现它们住在同一条纬度线上，于是它们约定各自朝西跳，直到碰

poj 2115 C Looooops（解二元一次不定方程）

题意： 。。。 思路： 和 青蛙的约会 差不多。。 d=B−A,m=2kd = B - A, m = 2^k 方程Cx−my=dCx - my = d 先利用拓展gcd求出 Cx0−my0=g=gcd(C,m)Cx_0-my_0 = g = gcd(C, m) ...最后求出最小正整数的解

二元一次不定方程的快速解法

二元一次不定方程（形如 a * x + b * y = c的方程，又叫丢番图方程，下简称不定方程），是初等数论经典的研究对象。二元一次不定方程应用广泛，如经典的找换问题和装箱问题（下面我将通过一系列的文章来解释我如何...

二元一次不定方程的整数解（扩展欧几里得算法）

二元一次不定方程的整数解（扩展欧几里得算法） （不得不说这是一堂数学*信竞课） 整数解解法 c(mod b)或ax+by=c有整数解当且仅当（a,b）|c 一般意义下的解法： 欧拉函数 扩展欧几里得算法 代码实现 ...

扩展欧几里德算法解二元一次不定方程

已知两个不完全为 0 的非负整数 a，b，必然存在整数对 x，y ，使它们满足贝祖等式： 解一定存在，根据数论中的相关定理。下面给出代码： int extgcd(int a, int b, int&amp; x, int&amp; y) { int gcd...

二元一次不定方程的解法

二元一次不定方程的解法 时间:2008-12-17 14:47 点击: 147次 　我们知道，如果未知数的个数多于方程的个数，那么，一般来说，它的解往往是不确定的，例如方程 x-2y=3，　方程组　等，它们的解是不确定的．像这...

P5656 【模板】二元一次不定方程(exgcd)

ExgcdQwQ

poj 1061 青蛙的约会 --扩展欧几里得 解二元一次不定方程例题

它们很高兴地发现它们住在同一条纬度线上，于是它们约定各自朝西跳，直到碰面为止。可是它们出发之前忘记了一件很重要的事情，既没有问清楚对方的特征，也没有约定见面的具体位置。不过青蛙们都是很乐观的，它们觉得...

通过扩欧得出线性同余方程的通解以及x的最小正整数解

若 ax与b模m的余数相同（其中x为未知数，即所需要求的数），即ax%m=b%m，则这个式子可以记作成a≡b (mod m)。 设ax对m取模后的余数为r1，则有：ax=y1*m+r1。      ① 同理，设b对m...

java多元一次方程组求解_java 怎么解多元一次不定方程组

-同余筛数法uniqueleion6152018-10-24Python100例——第五章----不定方程的解wdt338516542013-07-19高斯消元求解多元一次方程组nikelong021632016-03-26二元一次不定方程的快速解法uniqueleion37972018-07-30JAVA解N...

python编写程序实现二元一次方程的求解_教你用python解析二元一次方程

二元一次函数的实现import cmath import math import sys 这里导入cmath包是在后面用来处理复数的情况导入math使用来处理 平方 根号等的运算而导入sys的意义是为了比较0 ，在python中float的精度值不够，所以在计算...

洛谷 P5656 【模板】二元一次不定方程(exgcd)

1、这里需要解方程 ax + by = c 的所有正整数解。 不过，有个前提，a, b, c 都是正整数， 后面方便处理很多。 2、记 d = gcd(a, b), 如果 c % d != 0, 说明方程无整数解。 3、 用 扩展欧几里得 exgcd 算出 ax + by =...

多元一次不定方程的强力算法---同余筛数法

1、求解多元一次不定方程 n元一次不定方程就是形如∑aixi = C的...举例: 3x+4y+6z = 7，为了将3元变2元，这里我们做一个假设，设4y+6z=w，由不定方程的性质可知 2 | w,即w是2的倍数，由此我们不妨假设4y+6z=2w，...

【模板】二元一次不定方程 (exgcd)

若方程有整数解，但没有正整数解，你需要输出所有整数解中 x 的最小正整数值， y 的最小正整数值。 正整数解即为 x, y 均为正整数的解，0 不是正整数。 整数解即为 x, y 均为整数的解。 x 的最小正整数值即所有 x 为...

扩展欧几里德解二元一次不定方程

http://blog.csdn.net/zhjchengfeng5/article/details/7786595对于二元一次不定方程: ax + by = c;有解的充要条件是 ： c % gcd(a,b) ==0; 可用扩展欧几里得算法求得 某一个解 x, 但不一定是最终解； 最终解应为X ...

求佩尔方程的解

求关于x y的二次不定方程的解 x2-ny2=1 Input 多组输入数据，先输入组数T 然后输入正整数n(n&lt;=100) Output 对于每组数据输出一行，求y&lt;=10000的最小正整数解 ，输出y的值，如果在此范围内没有解则输出...

裴蜀定理（针对于多元一次不定方程的求解）

裴蜀定理的基本内容是】 若a,b是整数,且gcd(a,b)...假设c是ax+by的最小正整数解，也即是ax+by=c，d=(a,b） 根据最大公约数的性质及同余方程的性质，可以知道：d|a,d|b,d|(ax+by) 假设r是a%c的余数，t是a/c的整除...

浅谈扩欧及exgcd对二元不定方程求解问题

今天模拟赛第一题是一道对欧几里得和扩欧算法的简单应用，可惜两种方法都不会的我只能用求导和最小矩阵来存不定方程在坐标系上的整数解，满打满算了七十多行代码，其实一个扩欧就能解决的问题，被我想的很复杂。...

题解 洛谷P5656 【【模板】二元一次不定方程(exgcd)】

否则，可以将a,b,ca,b,ca,b,c都除以gcd⁡(a,b)\gcd(a,b)gcd(a,b)，只要这个等式成立，原等式一定成立，这样也就使得gcd⁡(a,b)=1\gcd(a,b)=1gcd(a,b)=1了。 然后用exgcd⁡\operatorname{exgcd}...

乘法逆元、扩展欧几里得算法、二元一次方程、a的n次方取余

知识点：乘法逆元，逆元的求法，二元一次方程求通解，a的n次方求余数一，乘法逆元乘法逆元的概念类似于倒数（ax=1,a−1=xax=1,a^{-1}=x），不过是在取余数的情况下的倒数。 如果(a×x)%p=1(a\times x)\%p = 1，则称...

ZOJ 3759 几个二元二次不定方程的解法

如果存在，求出最小解 k边形数的定义 那么题目意思实际是求解关于x，y的不定方程 代入具体的k并变形可以得到 可以看出这几个方程具有相同的形式，我们考虑更一般的情况 配方可以得到 换元可以得到...

洛谷 P5656 【模板】二元一次不定方程 (exgcd+分类讨论)

传送门 刚刚写题以为是exgcdexgcdexgcd板子出问题了，就去洛谷上找到这个板子题搞一波 本题肯定是板子题，但是需要知道exgcdexgcdexgcd求出的最小整数解会有哪几种情况，因为a,b,ca,b,...先给出扩欧求二元一次不定方程

解不定方程

解不定方程及其应用

【初等数论】指数、原根与不定方程

【初等数论】指数、原根与不定方程 1、指数 现在我们就开始为剩余系建立“坐标”，完全剩余系是连续的，剩余类本身就是很好的坐标，所以这里我们只需讨论既约剩余系。前面已经知道(a,m)=1(a,m)=1(a,m)=1时，总存 d ...

HDU 2669 扩展欧几里德求二元不定方程，超水..

可是书上的第二种解二元一次不定方程方法，目前还在纠结中... 运用的定理中的推论得出一定存在整数 x , y 使得 ax + by =gcd( a,b); 成立. 同样有定理得到 如果 ax + by = n 有解，则 gcd( a,b ) |

codevs4290 二元一次不定方程~（三星）

(^q^/∠) 描述：已知a,b,c, 求满足ax+by=c的整数对(x,y)中x为正整数且最小；若此时0,y，则输出x y，若此时y则输出“sometimes naive”（没有双引号）；若没有整数对(x,y)满足条件则输出“too simple” 输入描述 ...

我的算法不可能这么简单—扩展欧几里得

文章目录

计算欧式距离的matlab程序

这是一个计算欧式距离的matlab程序，本人使用过，很好用。

Hibernate4视频教程_全面来袭

本Java视频教程基于 Hibernate4.x 录制。内容涵盖安装 Hibernatetools 插件、Session 核心方法、持久化对象生命周期、对象关系映射（1-n、1-1、n-n、继承映射）、检索策略、检索方式（对象导航图、OID 检索、HQL、QBC、本地SQL）、Hibernate 一\二级缓存、管理 Session、批量处理等 Hibernate 企业级开发的核心技术。 本Java视频教程将帮助掌握学习者编写出具有合理的软件架构，以及好的运行性能和并发性能的实用 Hibernate 应用。Java视频教程内容注重理论与实践相结合，列举大量具典型性和实用价值的 Hibernate应用实例，并提供详细的开发和部署步骤。

计算机设计大赛作品开发文档

参加的是2020年的计算机设计大赛，软件应用与开发赛道。我们的开发文档仅供参考。（20页）

相关热词 c# sdk 微信 c#后台管理模板 c# 或取目录下文件 c#操作word书签 c# 比较两个数 c#语言打开.srec c# mvc 特性验证 c#类的对象初始化 c# 字符串完全匹配 c# 串口