//
// main.cpp
// Richard
//
// Created by 邵金杰 on 16/8/27.
// Copyright © 2016年 邵金杰. All rights reserved.
//
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long LL;
const int maxn=10000;
LL v[maxn],ans[maxn];
int maxd;
int get_first(LL a,LL b)
{
return (int)(b/a)+1;
}
int gcd(LL a,LL b)
{
return b==0?(int)a:gcd(b,a%b);
}
bool is_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,int from,LL aa,LL bb)
{
if(d==maxd)
{
if(bb%aa) return false;
v[d]=(bb/aa);
if(is_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++)
{
if(bb*(maxd-d+1)<=aa*i) break;
v[d]=i;
LL a2=aa*i-bb;
LL b2=bb*i;
LL g=gcd(a2,b2);
if(dfs(d+1,i+1,a2/g,b2/g)) ok=true;
}
return ok;
}
int main()
{
int a,b;
while(scanf("%d%d",&a,&b)!=EOF)
{
int ok=false;
for(maxd=1;;maxd++)
{
memset(ans,-1,sizeof(ans));
if(dfs(0,get_first(a,b),a,b)) {ok=true;break;}
}
if(ok)
{
printf("%d/%d=",a,b);
for(int i=0;i<=maxd;i++)
{
if(i==0) printf("1/%lld",ans[i]);
else printf("+1/%lld",ans[i]);
}
printf("\n");
}
}
return 0;
}
埃及分数(IDA*)
最新推荐文章于 2022-03-30 20:08:41 发布