日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

POJ1061 青蛙的约会(线性同余方程)

發布時間:2023/12/20 编程问答 44 豆豆
生活随笔 收集整理的這篇文章主要介紹了 POJ1061 青蛙的约会(线性同余方程) 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

線性同余方程$ ax \equiv b \pmod n$可以用擴展歐幾里得算法求解。

這一題假設青蛙們跳t次后相遇,則可列方程:

$$ Mt+X \equiv Nt+Y \pmod L$$ $$ (M-N)t \equiv Y-X \pmod L$$

于是就構造出一個線性同余方程,即可對t求解,解出最小非負整數解。

1 #include<cstdio> 2 #include<cstring> 3 using namespace std; 4 #define mod(x,y) (((x)%(y)+(y))%(y)) 5 #define lld long long 6 //a*x+b*y=gcd(a,b) 7 lld exgcd(lld a,lld b,lld &x,lld &y){ 8 if(b==0){ 9 x=1; y=0; 10 return a; 11 } 12 lld d=exgcd(b,mod(a,b),x,y); 13 lld t=y; 14 y=x-a/b*y; 15 x=t; 16 return d; 17 } 18 //ax??b (mod n) 19 lld MLES(lld a,lld b,lld n){ 20 lld x,y; 21 lld d=exgcd(a,n,x,y); 22 if(b%d) return -1; 23 return mod(x*(b/d),n/d); 24 } 25 26 int main(){ 27 lld x,y,m,n,l; 28 scanf("%lld%lld%lld%lld%lld",&x,&y,&m,&n,&l); 29 lld res=MLES(m-n,y-x,l); 30 if(res==-1) puts("Impossible"); 31 else printf("%lld",res); 32 return 0; 33 }

?

轉載于:https://www.cnblogs.com/WABoss/p/5179995.html

總結

以上是生活随笔為你收集整理的POJ1061 青蛙的约会(线性同余方程)的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。