博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
HDU1576 A/B【扩展欧几里得算法】
阅读量:6647 次
发布时间:2019-06-25

本文共 946 字,大约阅读时间需要 3 分钟。

问题链接

问题简述:参见上述链接。

问题分析这个问题可以用解整数的不定方程来解决,即使用扩展欧几里德算法。

根据题意,输入的n=A%9973(没有输入A),A%B=0(A必能被B整除),B与9973互素(GCD(B,9973)=1)。

解题过程首先是建立方程,然后才能编写程序。

设x=(A/B)%9973(x是最终想计算的值),则9973k+x=A/B(k为整数),得A=9973Bk+xB。

因为n=A%9973与A=9973Bk+xB,所以xB%9973=n,得xB=n+9973y。

故:(x/n)B+(-y/n)9973=1=GCD(B,9973),该方程有解。

要求x和y,先求X=x/n和Y=-y/n,即先解方程BX+9973Y=1。

最后,x=X*n。

需要注意的是,求得的x有可能是负值,需要进行调整。

不过,这个计算方法好像比较花时间。

程序说明:(略)。

AC的程序如下:

#include 
// 递推法实现扩展欧几里德算法long exgcd(long a, long b, long *x, long *y){ long x0=1, y0=0, x1=0, y1=1; long r, q; *x=0; *y=1; r = a % b; q = (a - r) / b; while(r) { *x = x0 - q * x1; *y = y0 - q * y1; x0 = x1; y0 = y1; x1 = *x; y1 = *y; a = b; b = r; r = a % b; q = (a - r) / b; } return b;}int main(void){ int t, i; long n, b, a=9973, x, y; scanf("%d", &t); for(i=0; i

转载于:https://www.cnblogs.com/tigerisland/p/7564930.html

你可能感兴趣的文章
IDEA 极速导包功能
查看>>
推荐子龙山人的emacs视频教程
查看>>
细说shiro之二:组件架构
查看>>
Linux---解压缩
查看>>
高可用Keepalived
查看>>
CentOS(RHEL)内核版本与系统版本对应关系
查看>>
Linux之本地yum源
查看>>
安装配置SVN服务-1
查看>>
vs2010开发安卓系统
查看>>
Splunk Forward简单部署_Win
查看>>
Oracle 双机热备:Oracle dataguard 和Oracle rac的区别和联系
查看>>
如何只显示不同字段值的行?
查看>>
挂载ISO文件
查看>>
DataGridView 经典用法总结(一)
查看>>
Java并发显式锁和显式条件队列
查看>>
云OS可国产替代
查看>>
try-catch 能否监听多线程中的错误?
查看>>
Android使用Token 实现单点登录
查看>>
模拟器可以,但是真机却不行
查看>>
CSS: hover选择器的使用
查看>>