codevs 1288 埃及分数 (迭代加深搜索)

题目大意：给你一个分数$a/b$，把它拆解成$sum_{i=1}^{n}1/ai$的形式，必须保证$ai$互不相同的情况下，尽量保证n最小，其次保证分母最大的分数的分母最小

什么鬼玄学题！！！

因为要保证$n$最小，所以从小到大遍历最大层数$n$，令$a'/b'$表示当前剩余的需要被拆解的数

如果当前分数的数量等于$n$，更新最优解并返回，继续搜总层数为$n$的所有情况，然后输出答案，因为这一层的最优解一定是全局最优解

每次选择一个数$i$，表示当前层的数选择了$i$，把剩余的数改成$a'/b'-1/i$进入下一层

显然暴力枚举$i$会$T$，那就加一些玄学的剪枝

1.i要大于等于$left lceil b/a ight ceil$，否则就剩余负数了

2.为了减少枚举次数，可以保证前几层选择的数升序排列，所以$i$要大于$now[dep-1]$

3.$(n-dep)/i+1/i>=a'/b'$，因为后几层选择的数都小于等于$1/i$，所以后几层和的最大值要大于等于$a'/b'$

有了这几个玄学的剪枝，我们就能以$O($能过$)$的复杂度通过本题！

#include <map>
#include <queue>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define NN 5010
#define MM 2000
#define ll long long 
#define uint unsigned int
#define ull unsigned long long 
using namespace std;

ll A,B;
int K=500;
ll now[NN],ans[NN];
ll gcd(ll x,ll y){if(!y)return x;return gcd(y,x%y);}
const ll maxn=666666666666ll;
void dfs(int dep,ll a,ll b,int ma)
{
    if(dep==ma+1){
        if(a>0) return;
        if(now[ma]<ans[ma])
            for(int i=1;i<=ma;i++)
            ans[i]=now[i];
        return;
    }
    ll na,nb,g;
    for(ll i=max(now[dep-1]+1,(b/a)+(b%a==0?0:1));a*i<=(ma-dep+1)*b;i++){
        na=a*i-b,nb=b*i;
        g=gcd(na,nb);
        now[dep]=i;
        dfs(dep+1,na/g,nb/g,ma);
        now[dep]=0;
    }
}

int main()
{
    //freopen("t2.in","r",stdin);
    scanf("%lld%lld",&A,&B);
    for(int k=1;k<=K;k++)
        ans[k]=maxn;
    for(int k=1;k<=K;k++)
    {
        dfs(1,A,B,k);
        if(ans[k]<maxn){
            for(int i=1;i<=k;i++)
                printf("%lld ",ans[i]);
            puts("");
            break;
        }
    }
    return 0;
}