poj 3134 && LA 3621 Power Calculus (迭代加深深度优先搜索)

https://icpcarchive.ecs.baylor.edu/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=1622

3134 -- Power Calculus

　　又搞了几天才出一题。

　　这题我是用迭代加深深度优先搜索来做的。网上的题解多种多样，大部分都是BFS的，不过好多都是没有证明的方法。表示个人认为迭代加深是一个十分好的正解。

　　至于这题为何用迭代加深，个人认为原因是：

（1）BFS剪枝不容易加上去，即便是加上一两个剪枝也无法剪枝到系统可以接受的时间和空间。（网上有题解说搜索到多少个后停止搜索的方法，可惜没有详细的证明，故认为那是水过的）

（2）IDDFS与BFS相比，因为是在搜索的过程中直接记录路径，所以相对BFS节省了大量空间。

　　除此以外，表示打了好几个版本的代码，IDDFS是相当轻便的一种方法。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <vector>
#include <cmath>

using namespace std;

const int N = 1 << 11;
int ep2[15], lg2[N];
#define REP(i, n) for (int i = 0; i < (n); i++)
#define REP_1(i, n) for (int i = 1; i <= (n); i++)
#define DEC(i, a, b) for (int i = (a); i >= (b); i--)
#define PB push_back
#define SZ(x) ((int) (x).size())

void PRE() {
    REP(i, 15) {
        ep2[i] =  1 << i;
    }
    int t = 0;
    REP(i, N) {
        if (ep2[t + 1] <= i) t++;
        lg2[i] = t;
    }
//    REP(i, 10) cout << lg2[i] << endl;
}

int st[20], top;

bool dfs(int n, int depth) {
    if (st[top] == n) return true;
    if (top >= depth) return false;
    int last = st[top];
    if ((last << (depth - top)) < n) return false;
    if (st[top] <= 1200) {
        DEC(i, top, 0) {
            st[++top] = last + st[i];
            if (dfs(n, depth)) return true;
            top--;
        }
    }
    REP(i, top + 1) {
        st[++top] = last - st[i];
        if (st[top] > 1) {
            if (dfs(n, depth)) return true;
        }
        top--;
    }
    return false;
}

int main() {
//    freopen("out", "w", stdout);
    int n;
    PRE();
//    REP_1(n, 1000) {
//        int step = lg2[n];
//        while (step <= 12) {
//            st[top = 0] = 1;
//            if (dfs(n, step)) break;
//            step++;
//        }
//        cout << n << ' ' << step << endl;
//    }
//    return 0;
    while (cin >> n && n) {
        int step = lg2[n];
        while (step <= 12) {
            st[top = 0] = 1;
            if (dfs(n, step)) break;
            step++;
        }
        cout << step << endl;
    }
    return 0;
}

——written by Lyon