marbles-ACM

算法：设A为弹珠的总数，x,y分别为两种盒子的数目。Z为买这些盒子的花销。那么就有下面的等式：n1*x+n2*y=A,c1*x+c2*y=Z. 如果先救y的话，就要先求出n1与c1的最小公倍数。求最小公倍数的公式有lcm(x,y)=x*y/gcd(x,y),因为x*y可能会超出范围，还有一个公式是（x/gcd(x,y)）*y. 利用第二个公式可以求出n1与c1的最小公倍数。设lcm(c1,n1)=a*n1=b*c1. 两式相减得，y=(b*Z-a*A)/(b*c2-a*n2). 因为y>=0, 所以b*Z-a*A>=0, b*c2-a*n2>0. 如果第二个等式不成立的话，便可判定没有解决方案。如果成立的话，Z>=a*A/b. 这样就可以求出Z，然后代回等式，就可救出y,如果y不是整数的话，也说明没有解决方案。再求出x, 如果x不是整数，也可说明没有解决方案。

 

// Marbles.cpp : Defines the entry point for the console application. // #include "stdafx.h" using namespace std; int _tmain(int argc, _TCHAR* argv[]) { unsigned int MOD(unsigned int x,unsigned int y); unsigned int LCM(unsigned int c,unsigned int n); unsigned int marble=0; while(cin>>marble) { if(marble==0)break; else { unsigned int c1,c2,n1,n2,b1,b2; unsigned int lcm=0; unsigned int pay=0; bool flag=true; cin>>c1>>n1; cin>>c2>>n2; lcm=LCM(c1,n1); unsigned int a,b; a=lcm/n1; b=lcm/c1; if((b*c2-a*n2)<=0)flag=false; else { if((a*marble)%b!=0)pay=(a*marble)/b+1; else pay=(a*marble)/b; int tmp1=(b*pay-a*marble)%(b*c2-a*n2); if(tmp1==0) { b2=(b*pay-a*marble)/(b*c2-a*n2); int tmp2=(pay-c2*b2)%c1; if(tmp2==0)b1=(pay-c2*b2)/c1; else flag=false; } else flag=false; } if(flag==true)cout<<b1<<" "<<b2<<"/n"; else cout<<"failed"<<"/n"; } } system("pause"); return 0; } unsigned int MOD(unsigned int x,unsigned int y) { unsigned int i=x; unsigned int j=y; unsigned int r=x%y; while(j!=0) { r=i; i=j; j=r%j; } return i; } unsigned int LCM(unsigned int c,unsigned int n) { if(c>=n) return c/MOD(c,n)*n; else return n/MOD(n,c)*c; }