vijod1308(IDA*)

题意：在古埃及，人们使用单位分数的和(形如1/a的, a是自然数)表示一切有理数。如：2/3=1/2+1/6,但不允许2/3=1/3+1/3,因为加数中有相同的。对于一个分数a/b,表示方法有很多种，但是哪种最好呢？首先，加数少的比加数多的好，其次，加数个数相同的，最小的分数越大越好。
如：19/45=1/3 + 1/12 + 1/180
19/45=1/3 + 1/15 + 1/45
19/45=1/3 + 1/18 + 1/30,
19/45=1/4 + 1/6 + 1/180
19/45=1/5 + 1/6 + 1/18. 
最好的是最后一种，因为1/18比1/180,1/45,1/30,1/180都大。
给出a,b(0<a<b<1000),编程计算最好的表达方式。

链接：点击打开链接

代码：

#include <map>
#include <set>
#include <queue>
#include <string>
#include <math.h>
#include <vector>
#include <stdio.h>
#include <stdlib.h>
#include <iostream>
#include <string.h>
#include <algorithm>
using namespace std;
int dep,sig;
long long ans[1005],tmp[1005];
long long gcd(long long x,long long y){
    if(y==0)
    return x;
    return gcd(y,x%y);
}
long long get(long long x,long long y){
    if((y%x)==0)                                //找到小于等于x/y的最大单位分数
    return y/x;                                 //不要暴力去找
    return y/x+1;
}
int judge(int d){
    int i;
    if(ans[d]==-1||ans[d]>tmp[d]){              //因为是要保证最后一个分数最大，所以
        for(i=1;i<=d;i++)                       //每一个深度的所有答案要再检查一遍
        ans[i]=tmp[i];
        return 1;
    }
    return 0;
}
void dfs(int d,long long u,long long aa,long long bb){
    long long i,x,y,z;
    if(d==dep){
        if(aa!=1)
        return;
        tmp[d]=bb;
        if(judge(d)){
            sig=1;
            return;
        }
    }
    u=max(u,get(aa,bb));                
    for(i=u;;i++){                              //如果全选当前都不够直接break
        if((aa*i)>=(bb*(dep-d+1)))
        break;
        tmp[d]=i;
        x=aa*i-bb;
        y=bb*i;
        z=gcd(x,y);
        dfs(d+1,i+1,x/z,y/z);
    }
}
int main(){
    int i;
    long long a,b;
    while(scanf("%lld%lld",&a,&b)!=EOF){
        sig=0;
        memset(ans,-1,sizeof(ans));
        for(dep=1;;dep++){
            dfs(1,get(a,b),a,b);
            if(sig==1){
                for(i=1;i<=dep;i++){
                    if(i!=1)
                    printf(" ");
                    printf("%lld",ans[i]);
                }
                printf("\n");
                break;
            }
        }
    }
    return 0;
}