UVa 12558 - Egyptian Fractions (HARD version)（IDA*）

最新推荐文章于 2021-07-27 11:47:59 发布

wcr1996

最新推荐文章于 2021-07-27 11:47:59 发布

阅读量2.3k

点赞数

分类专栏：算法竞赛入门经典（第二版）第七章暴力求解法文章标签： uva

本文链接：https://blog.csdn.net/wcr1996/article/details/43632113

版权

算法竞赛入门经典（第二版）第七章暴力求解法专栏收录该内容

17 篇文章 0 订阅

订阅专栏

给出一个真分数，把它分解成最少的埃及分数的和。同时给出了k个数，不能作为分母出现，要求解的最小的分数的分母尽量大。

使用IDA*算法解决，在例题求埃及分数的基础上，加上禁用限制就可以了，一开始禁用限制没加好，一直WA。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<set>
using namespace std;
typedef long long LL;
const int maxn=10010;
int maxd,t,tt;
set<LL> sk;
LL ans[maxn],v[maxn];
LL gcd(LL a,LL b){
    return b?gcd(b,a%b):a;
}
LL get_first(LL a,LL b){
    return b/a+1;
}
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;
}
bool dfs(int d,LL from,LL a,LL b){
    if(d==maxd){
        if(b%a) return false;
        v[d]=b/a;
        if(sk.count(b/a)) return false;
        if(better(d)) memcpy(ans,v,sizeof(LL)*(d+1));
        return true;
    }
    bool ok=false;
    for(LL i=max(from,get_first(a,b));;++i){
        if(b*(maxd+1-d)<=i*a) break;
        if(sk.count(i)) continue;
        v[d]=i;
        LL b2=b*i;
        LL a2=a*i-b;
        LL g=gcd(a2,b2);
        if(dfs(d+1,i+1,a2/g,b2/g)) ok=true;
    }
    return ok;
}
int main(){
    scanf("%d",&t);
    while(t--){
        LL a,b,k,sk0;
        sk.clear();
        scanf("%lld%lld%lld",&a,&b,&k);
        while(k--) scanf("%lld",&sk0),sk.insert(sk0);
        for(maxd=0;;++maxd){
            memset(ans,-1,sizeof(ans));
            if(dfs(0,get_first(a,b),a,b)) break;
        }
        printf("Case %d: %lld/%lld=",++tt,a,b);
        for(int i=0;i<=maxd;++i){
            if(i) printf("+");
            printf("1/%lld",ans[i]);
        }
        printf("\n");
    }
    return 0;
}