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

题目大意：

（来源：洛谷）

解题思路：

因为是第一次写IDA的博客，所以用oiwiki来简单介绍下IDA

• 即相当于dfs的改进，既有迭代加深又有估价函数

对于这道题我们先来想一个简单版，如果没有要求哪些分母不能用怎么办呢？

• 首先是不太可能用bfs来做，因为正整数有无数多个，所以方向是无数多个，所以这道题我们利用dfs来做
• 而dfs不限制深度是不行的，所以我们将其迭代加深来搜索
• 又因为不能有相同的分母，所以我们将分母从小到大来搜索
• 那么对于$\frac{a}{b}$的最小分母就是满足$\frac{1}{c}\le \frac{a}{b}$，所以$c\ge \frac{b}{a}$，所以我们往往向上取整$\frac{b}{a}$，就是直接 $\frac{b}{a}+1$（之后在代码还会有进一步解释，为什么不担心b整除a的情况）
• 最后再保持分母递增的情况下，到达指定层数时，再判断答案是否可行，以及是否最优即可
• 那么估价函数（不一定要写成函数形式，有这个体现即可）是啥呢，如果我当前要选$\frac{1}{i}$，如果余下$(maxd+1-d)\ast \frac{1}{i}<\frac{a}{b}$则说明继续下去没意义了

那么对于这题，如果要求了哪些分母不能用怎么办？

• 加个map保存一下不能用的以及用过的即可
• 还有要注意开long long

AC代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll maxn = 1e3 + 10;
map <ll, bool> mp;
ll v[maxn], ans[maxn];
//不担心b%a==0，是因为迭代加深的关系b%a==0一定会在最后一个才出现
ll get_first(ll a, ll b) {return b / a + 1;} 
ll gcd(ll a, ll b) {return b > 0 ? gcd(b, a % b) : a;}
bool better(ll d) {
  	for (ll i = d; i >= 0; i--)
    	if (v[i] != ans[i]) {
      		return ans[i] == -1 || v[i] < ans[i];
    	}
  	return false;
}
bool dfs(ll d, ll maxd, ll from, ll aa, ll bb) {
	//cout << from << endl;
	if (d == maxd) {
		if (bb % aa) return false; //如果最后bb%aa不是0的话，说明在该种选法下还要拓展层数，所以返回
		if (mp[bb / aa]) return false;
		v[d] = bb / aa;
		if (better(d)) memcpy(ans, v, sizeof(v));
		return true;
	}
	from = max(from, get_first(aa, bb));//取较大的那个
	bool ok = false;
	for (ll i = from, a2, b2; ; i++) {
		if (bb * (maxd + 1 - d) < aa * i) break;
		if (i <= 1000 && mp[i]) continue;
		a2 = i * aa - bb, b2 = bb * i;
		ll _d = gcd(a2, b2);
		a2 /= _d, b2 /= _d;
		mp[i] = true;
		v[d] = i;
		if (dfs(d + 1, maxd, i + 1, a2, b2)) ok = true;
		mp[i] = false;
	}
	return ok;
}
ll T, cas, a, b, k;
int main() {
	cin >> T;
	while (T--) {
		mp.clear();
		memset(v, 0, sizeof(v));
		memset(ans, 0, sizeof(ans));
		cin >> a >> b >> k;
		ll u;
		for (ll i = 1; i <= k; i++) {
			cin >> u;
			mp[u] = true;
		}
		bool ok = false;
		ll maxd;
		for (maxd = 2; ; maxd++) { //答案保证了有解
			memset(ans, -1, sizeof(ans));
			if (dfs(1, maxd, get_first(a, b), a, b)) {
				ok = true;
				break;
			}
		}
    	cout << "Case " << ++cas << ": ";
    	// if (ok) {
      	cout << a << "/" << b << "=";
      	for (ll i = 1; i <= maxd; i++) cout << "1/" << ans[i] << "+\n"[i == maxd];
    	// } 
		// else cout << "No solution.\n";
  	}
}