埃及分数


#include<stdio.h>
#include<string.h>
const int N = 1000;
int maxd;
int ans[N];
int v[N];
int get_first(int a, int b)
{
	int k = b/a;
	return k*a >= b ? k : k+1;
} 
int gcd(int a, int b)
{
	return b==0? a : gcd(b, a%b);
}
bool better(int d)
{
	for(int i = d; i >= 0; --i)
	{
		if(ans[i] != v[i])
		{
			return ans[i] == -1 || v[i] < ans[i];
		}
	}
	return false;
}  
bool dfs(int cur, int from, int aa, int bb)
{
	if(cur == maxd)
	{
		if(bb%aa == 0)
		{
			v[cur] = bb/aa;
			if(better(cur))
			{
				memcpy(ans, v, sizeof(ans));
				return true;
			}
		}
		return false;
	}
	bool res = false;
	for(int i = from; true; ++i)
	{
		if((maxd+1-cur)*bb <= aa*i)
			break;
		v[cur] = i;
		int a = aa*i-bb;
		int b = bb*i;
		int g = gcd(a,b);
		a /= g;
		b /= g;
		int f = get_first(a,b) > i+1 ? get_first(a,b) : i+1;
		if(dfs(cur+1, f, a, b))
		{
			res = true;
		}
	}
	return res;

}
int main()
{
	//
	freopen("input.txt", "r", stdin);
	int a, b;
	while(scanf("%d%d", &a, &b) == 2)
	{
		bool res = false;
		if(a+b == 0)
			break;

		for(maxd = 1; !res; ++maxd)
		{
			memset(ans, -1, sizeof(ans));
			res = dfs(0, get_first(a, b), a, b);
		}
		maxd--;
		for(int i=0; i <= maxd; ++i)
		{
			printf("%d ", ans[i]);
		}
		putchar('\n');
	}
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值