bzoj1407: [Noi2002]Savage
生活随笔
收集整理的這篇文章主要介紹了
bzoj1407: [Noi2002]Savage
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
也是復習題了?傻逼衣叉還我一等
兩個野人無法相遇當
Ci+Pi*x=Cj+Pj*x (mod m)
(Pi-pj)*x=Cj-Ci (mod m)
無解 或 最小正整數解>min(Li,Lj)
轉換一下
(Pi-pj)*x+m*y=Cj-Ci
要約一個最大公因數(不約居然會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; }?
轉載于:https://www.cnblogs.com/AKCqhzdy/p/8783172.html
總結
以上是生活随笔為你收集整理的bzoj1407: [Noi2002]Savage的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 前端开发的正确姿势——各种文件的目录结构
- 下一篇: 6大设计原则总结