Egyptian Fractions (HARD version)（UVA 12558）

网址如下：

Egyptian Fractions (HARD version) - UVA 12558 - Virtual Judge (vjudge.net)

（第三方网站）

唐完啦！唐完啦！

这题《埃及数字》甚至在算法书里面有题解，虽然说是hard version，但也没有harder到哪里去

就加了个限制条件：有k个数字不能用

然后我就按照印象中的思路开始写，结果直接裂了，把循环结束条件按习惯写成剪枝条件，分析了一会，加了个全局变量改改，感觉还是不行，回过头去看源代码才知道自己有多唐

总的来说就是用迭代加深搜索，因为有用到剪枝的方法，实际上是IDA*

代码如下：

#include<cstdio>
#include<cstring>
#include<algorithm>
typedef long long LL;

const int maxk = 5, maxn = 1000;

LL num[maxk], ans[maxn], tmp[maxn];
int T, a, b, k, cnt, maxd;

int get_first(LL aa, LL bb){
    for(int i = 2; ; i++) if(bb <= aa * i) return i;
}
bool is_k(LL n){
    for(int i = 0; i < k; i++){
        if(num[i] == n) return true;
    }
    return false;
}
bool better(int d){
    for(int i = d; i >= 0; i--) if(ans[i] != tmp[i]) return ans[i] == 0 || tmp[i] < ans[i];
    return false;
}
bool dfs(int d, int from, LL aa, LL bb){
    if(d == maxd){
        if(bb % aa || is_k(bb / aa)) return false;
        tmp[d] = bb / aa;
        if(better(d)) memcpy(ans, tmp, sizeof(LL) * (maxd + 1));
        return true;
    }
    bool ok = false;
    from = std::max(from, get_first(aa, bb));
    for(int i = from; ; i++){
        if(bb * (maxd + 1 - d) <= i * aa) break;
        if(!is_k(i)){
            tmp[d] = i;
            LL a2 = aa * i - bb;
            LL b2 = bb * i;
            LL g = std::__gcd(a2, b2);
            if(dfs(d + 1, i + 1, a2 / g, b2 / g)) ok = true;
        }
    }
    return ok;
}

int main(void)
{
    scanf("%d", &T);
    for(int kase = 1; kase <= T; kase++){
        scanf("%d%d%d", &a, &b, &k);
        for(int i = 0; i < k; i++) scanf("%lld", &num[i]);
        for(maxd = 1; ; maxd++){
            memset(ans, 0, sizeof(ans));
            if(dfs(0, get_first(a, b), a, b)) break;
        }
        //输出
        printf("Case %d: %lld/%lld=", kase, a, b);
        for(int i = 0; i <= maxd; i++){
            if(i) putchar('+');
            printf("1/%lld", ans[i]);
        }
        putchar('\n');
    }

    return 0;
}

其中并没有用到浮点数，直接用两个整型来表示一个分数

循环结束条件：

if(bb * (maxd + 1 - d) <= i * aa) break;

即1/i * (maxd + 1 - d) <= aa / bb

因为后面枚举的分数都不会比1/i大了，且剩下的分数个数为maxd + 1 - d，在最好情况下（就是都等于1/i），剩余分数的总和如果不能比aa/bb大，就没必要继续了，直接剪枝，又因为实际上分数不能相等，所以应该加一个等于号

其他就是迭代加深搜索的正常操作了