博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
POJ 1061 青蛙的约会 扩展欧几里得
阅读量:7052 次
发布时间:2019-06-28

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

扩展欧几里得模板套一下就A了,不过要注意刚好整除的时候,代码中有注释

#include 
#include
#include
#include
#include
using namespace std;typedef long long ll;ll exgcd(ll a, ll b, ll&x, ll&y) { if (b == 0) { x = 1; y = 0; return a; } ll r = exgcd(b, a%b, y, x); ll t = x; y = y - a/b*t; return r;}bool modular_linear_equation(ll a, ll b, ll n) { ll x, y, x0, i; ll d = exgcd(a, n, x, y); if (b%d) { printf("Impossible\n"); return false; } x0 = x*(b/d)%n; //x0为方程的一个特解,可以为正也可以为负。题目要求的是最小的非负数 ll ans; if(x0<0) { ans=x0; for(i = 0;ans<0; i++) ans=(x0 + i*(n/d))%n; } else if(x0>0) { ans=x0; ll temp; for(i=0;ans>=0;i++) { temp=ans; ans=(x0 - i*(n/d))%n; } ans=temp; } else ans=n; //此时x0=0,但是结果一定会大于0,所以要加一个n printf("%I64d\n",ans); return true;}int main(){ //freopen("in.txt","r",stdin); //freopen("out.txt","w",stdout); ll x,y,m,n,L; while(~scanf("%I64d%I64d%I64d%I64d%I64d",&x,&y,&m,&n,&L)) { ll temp1=m-n,temp2=y-x; if(temp1<0) { temp1=-1*temp1; temp2=-1*temp2; } modular_linear_equation(temp1,temp2,L); } return 0;}

 

转载于:https://www.cnblogs.com/pach/p/5924896.html

你可能感兴趣的文章
171. Excel Sheet Column Number
查看>>
简单深搜
查看>>
java-http请求
查看>>
WMI VS. ExplorerOM
查看>>
问题的定义和管理复杂度《Code Complete 2》
查看>>
Android init.rc解析【转】
查看>>
算法(Algorithms)第4版 练习 2.2.11(2)
查看>>
云解放了计算机这台机器,让计算的能力彻底从一个箱子里释放出来,回归了计算的本质。...
查看>>
Java ---- baidu评价抽取关键词-商品评论
查看>>
排序算法之选择排序
查看>>
三门问题
查看>>
ES6之主要知识点(十)Proxy
查看>>
UITableView
查看>>
list详解
查看>>
oracle学习篇六:子查询
查看>>
关于获取客户端Mac地址
查看>>
紫书 例题 10-9 UVa 1636 (概率计算)
查看>>
51nod 01背包
查看>>
outlook anywhere 配置
查看>>
关于android帮助文档打开慢
查看>>