迭代加深搜索 紫书7-6

迭代加深搜索：  迭代+DFS（循环+DFS）

以埃及分数问题为例：

      埃及分数：分子是1的分数，叫单位分数。古代埃及人在进行分数运算时。只使用分子是1的分数。因此这种分数也叫做埃及分数，或者叫单分子分数。

紫书例题：

题目描述 Description

在古埃及，人们使用单位分数的和(形如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),编程计算最好的表达方式。

输入描述 Input Description

a b

输出描述 Output Description

若干个数，自小到大排列，依次是单位分数的分母。

样例输入 Sample Input

19 45

样例输出 Sample Output

5 6 18

分析：从小到大枚举深度上限maxd，每次执行只考虑不超过maxd的结点。只要深度有限，一定可以在有限时间内枚举到。

代码如下：

#include <bits/stdc++.h>


using namespace std;


typedef long long ll;
const int maxn = 1000;
ll ans[maxn],v[maxn];
int maxd;
//求最大公约数
ll gcd(ll a,ll b){
    if(b==0) return a;
    else return gcd(b,a%b);
}
//求满足1/c <= a/b 的最小c
int get_first(ll a,ll b){
    for(int i=1;;i++)
        if(b<=a*i) return i;
}
//如果当前解优与最优解，更新ans
bool better(int d){
    for(int i=d;i>=0;i--){
        if(v[i]!=ans[i]) return ans[i]==-1||v[i]<ans[i];
    }
    return false;
}
//当前深度为d，from为搜索起点，分数和为aa/bb
bool dfs(int d,int from,ll aa,ll bb){
    if(d==maxd){
        if(bb%aa) return false;
        v[d]=bb/aa;
        if(better(d)) memcpy(ans,v,sizeof(ll)*(d+1));//更新
        return true;
    }
    bool ok= false;
    from=max(from,get_first(aa,bb));
    for(int i=from;;i++){
        //剪枝:如果剩下的maxd+1-d个分数全部都是1/i,加起来仍然小于aa/bb,无解,剪
        if(bb*(maxd+1-d)<=i*aa) break;
        v[d]=i;
        //计算aa/bb-1/i
        ll b2=bb*i;
        ll a2=aa*i-bb;
        ll g=gcd(a2,b2);
        if(dfs(d+1,i+1,a2/g,b2/g)) ok=true;
    }
    return ok;
}
int main(){
    int a,b,kase=1;
    while(cin>>a>>b){
        bool ok = false;
        for(maxd=1;;maxd++){
            memset(ans,-1,sizeof(ans));
            if(dfs(0,get_first(a,b),a,b)) {ok=true;break;}
        }
        cout<<"Case "<<kase++<<" : ";
        if(ok){
            for(int i=0;i<=maxd;i++){
                if(i>0) cout<<"+";
                printf("1/%lld",ans[i]);
            }
            cout<<endl;
        }else cout<<-1<<endl;
    }
}