埃及分数问题

// by BNU_LZM
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;

int ans[1010], v[1010], maxd;

int gcd(int a, int b)
{
	return b == 0 ? a : gcd(b, a%b);
}

int get_from(int a, int b)
{
	int i;
	for(i = 2; ; i++)
	{
		if(b <= a*i) break;
	}
	return i;
}

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;
}

int dfs(int cur, int from, int a, int b)
{
	if(cur+1 == maxd)
	{
		if(b%a == 0)
		{
			b /= a;
			if(v[cur] == b) return 0;
			v[maxd] = b;
			if(better(maxd)) memcpy(ans, v, sizeof(v));
			return 1;
		}
		else return 0;
	}
	int ok = 0;
	from = max(from, get_from(a, b));
	for(int i = from; ; i++)
	{
		if(b*(maxd-cur) < a*i) break;
		v[cur+1] = i;
		int bb = b*i;
		int aa = a*i-b;
		int g = gcd(aa, bb);
		//printf("%d\n", cur);
		if(dfs(cur+1, i+1, aa/g, bb/g)) ok = 1;
	}
	return ok;
}

int main()
{
	//freopen("in.txt", "r", stdin);
	int a, b;
	cin >> a >> b;
	printf("%d/%d=", a, b);
	int ok = 0;
	for(maxd = 1; ; maxd++)
	{
		memset(ans, -1, sizeof(ans));
		if(dfs(0, get_from(a, b), a, b)) 
		{
			ok = 1;
			break;
		}
	}
	if(ok)
	{
		for(int i = 1; i <= maxd; i++) 
		{
			if(i == 1)
				printf("1/%d", ans[i]);
			else printf("+1/%d", ans[i]);
		}
	}
		return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值