sicily 1099 Packing Passengers

/**
**  题意:  求x,y 满足 x*pa+y*pb=n 同时使得 p = x*ca+y*cb的值最小,若有多种可能,则选择最大x值的组合
**  分析:  x*pa+y*pb=n 可以用 线性同余方程 求得各组解: X=x+(pb/q)*t, Y=y-(pa/q)*t 
**   t为整数,显然 X,Y>=0,所以 -x/(pb/q) <= t <= y/(pa/q)
**   欲使 x*ca+y*cb 最小, 注意到 其中 [(pb/q)*t]*ca + [-(pa/q)*t]*cb = (pb*ca-pa*cb)/q*t
**   当pb*ca-pa*cb<0,t应取最大值;
**   当pb*ca-pa*cb>0,t应取最小值;
**   当pb*ca-pa*cb=0,应选择较大的X值,显然 pb/q总是大于0的,所以t 选择最大值,X取得更大值
**   http://zh.wikipedia.org/wiki/%E6%89%A9%E5%B1%95%E6%AC%A7%E5%87%A0%E9%87%8C%E5%BE%97%E7%AE%97%E6%B3%95
***/

 #include<iostream>
 #include<cmath>
 using namespace std;

 /**
 **    扩展欧几里德算法(辗转相除法), x*a+y*b=gcd(a,b)=q
 **/
 long long x,y,q;
 void extend_eulid(long long a,long long b){
     if(b==){
         x=;
         y=;
         q=a;
     }
     else{
         extend_eulid(b,a%b);
         //收集辗转相除法中产生的式子，倒回去，可以得到ax+by=gcd(a,b)的整数解
         long long temp=x;
         x = y;
         y = temp-a/b*y;
     }
 }
 int main()
 {
     int n,ca,pa,cb,pb,cnt=;
     long long upper,lower,t;
     while(cin>>n&&n){
         cin>>ca>>pa>>cb>>pb;
         extend_eulid(pa,pb);
         if(n%q!=)
             cout<<"Data set "<<cnt++<<": cannot be flown"<<endl;
         else{
             //求特解
             x = x*n/q;
             y = y*n/q;
             upper = (long long)floor((double)y/(pa/q));     //返回小于或者等于指定表达式的最大整数头
             lower = (long long)ceil((double)-x/(pb/q));     //返回大于或者等于指定表达式的最小整数头
             t = pb*ca-pa*cb<= ? upper : lower;
             cout<< "Data set "<<cnt++<<": "<<x+(pb/q)*t<<" aircraft A, "<<y-(pa/q)*t<<" aircraft B"<<endl;
         }
     }
     return ;
 }