UVA12558 Egyptian Fractions (HARD version)（埃及分数）

传送门

题目大意

给出一个真分数 a/b，要求出几个互不相同的埃及分数（从大到小），使得它们之和为 a/b

（埃及分数意思是分子为1的分数，详见百度百科）

如果有多组解，则分数数量少的优先

如果分数数量一样则分母最大的要尽量小，如果最大的分母同样大，则第二大的分母尽量小，以此类推

为了加大难度，会给出k个不能作为分母的数

（2<=a,b<=876，k<=5 并且 a,b 互质）

首先想的是数论，但是呢

推不出来...

然后发现a,b好像不大

貌似可以搜索

但是呢

不知道上界...

那就迭代加深搜索呗

然后想想怎么剪枝

如果知道 a/b，要怎么求出最小 k 使 1/k < a/b 呢（注意符号）

易知 a/b 为真分数 ，a,b,k均为整数
所以 a<b
所以必定有整数 x,y 使得 x*a+y=b 且y<a
所以可得 int(b/a)=int( (x*a+y) /a )=int(x*a/a)+int(y/a) = x
因为 int(b/a)<=b/a，int(b/a)=x
所以 x+ > b/a , 所以 int(b/a)+1>b/a，

所以1/(x+) < a/b 又 /int(b/a)=/x >=a/b 所以 x+1就是最小的 k 使得 /k < a/b 所以 k=int(b/a) +  证明完毕

知道了 k 的下界，那还要求 k 的上界

显然每次求的分母肯定要比上一次大（题目要求）

所以如果后面所有的分母都为 k ，加起来还没有当前要求的数大，那就不要再搜下去了，再下去也不会有结果的

然后就可以搜了，至于判断是否可用我用的是set

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<set>
using namespace std;
const int N=1e5+;
set <long long> p;
long long t,A,B,n,mxdep;
long long ans[N],w[N];
bool flag;
inline long long gcd(long long a,long long b)
{
    return b ? gcd(b,a%b) : a;
}//gcd是为了约分
inline bool pd()
{
    for(int i=mxdep;i;i--)
        if(ans[i]!=w[i])
            return ans[i]&&ans[i]<w[i];
    return ;
}//判断答案是否更优
inline void dfs(long long dep,long long a,long long b,long long mi)
{
    if(dep==mxdep)
    {
        if(b%a||(b<(mi-)*a+)||p.count(b/a)) return;//判断合法性
        //如果a/b可化为 1/k 的形式并且 1/k 可以选
        w[dep]=b/a;
        if(pd()) return;//判断答案是否更优
        memcpy(ans,w,sizeof(w)); flag=;//更新
        return;
    }
    mi=max(mi,b/a+);//求出下界
    for(long long i=mi;i;i++)
    {
        if( (mxdep-dep+)*b<=i*a ) return;//上界
        if(p.count(i)) continue;//判断合法
        w[dep]=i;//记录路径
        long long xa=a*i-b,xb=i*b;//通分
        long long z=gcd(xa,xb);
        dfs(dep+,xa/z,xb/z,i+);//约分，向下一层
    }
}
int main()
{
    cin>>t;
    for(long long i=;i<=t;i++)
    {
        flag=; mxdep=;
        cin>>A>>B>>n;
        long long c;
        p.clear();//细节
        while(n--) scanf("%lld",&c),p.insert(c);//存储不合法的数
        while(mxdep++)
        {
            memset(ans,,sizeof(ans));
            dfs(,A,B,B/A+);
            if(flag) break;
        }//迭代加深
        printf("Case %lld: %lld/%lld=1/%lld",i,A,B,ans[]);
        for(int i=;i<=mxdep;i++)
            printf("+1/%lld",ans[i]);
        printf("\n");
    }
    return ;
}