也是復(fù)習(xí)題了?傻逼衣叉還我一等

兩個野人無法相遇當(dāng)
Ci+Pi*x=Cj+Pj*x (mod m)
(Pi-pj)*x=Cj-Ci (mod m)
無解 或 最小正整數(shù)解>min(Li,Lj)

轉(zhuǎn)換一下
(Pi-pj)*x+m*y=Cj-Ci

要約一個最大公因數(shù)（不約居然會T囧）

?

#include<cstdio> #include<iostream> #include<cstring> #include<cstdlib> #include<algorithm> #include<cmath> using namespace std;int exgcd(int a,int b,int &x,int &y) {if(a==0){x=0;y=1;return b;}else{int tx,ty;int d=exgcd(b%a,a,tx,ty);x=ty-(b/a)*tx;y=tx;return d;} } int gcd(int a,int b) {if(a==0)return b;else return gcd(b%a,a); }int c[20],p[20],l[20]; int main() {int n,m=0;scanf("%d",&n);for(int i=1;i<=n;i++){scanf("%d%d%d",&c[i],&p[i],&l[i]);m=max(m,c[i]);}while(1){bool bk=true;for(int i=1;i<=n;i++){for(int j=i+1;j<=n;j++){int A=p[i]-p[j],B=m,K=c[j]-c[i];int t=gcd(A,B);if(K%t==0){int x,y;A/=t,B/=t,K/=t;int d=exgcd(A,B,x,y);B=abs(B);x=((x*K/d)%(B/d)+(B/d))%(B/d);if(x<=min(l[i],l[j])){bk=false;break;}}}if(bk==false)break;}if(bk==true){printf("%d\n",m);break;}m++;}return 0; }

?

轉(zhuǎn)載于:https://www.cnblogs.com/AKCqhzdy/p/8783172.html

總結(jié)

以上是生活随笔為你收集整理的bzoj1407: [Noi2002]Savage的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。