// The Priest Mathematician (牧师数学家) // PC/UVa IDs: 110606/10254, Popularity: C, Success rate: high Level: 2 // Verdict: Accepted // Submission Date: 2011-06-07 // UVa Run Time: 0.048s // // 版权所有(C)2011,邱秋。metaphysis # yeah dot net // // 设 N 个盘子在 4 柱情况下的移动最小步数为 T(N),在 3 柱子情况下的最小移动步数为 H(N),则 // 通过下列方法可以得到最小移动步数 T(N)。假设第一次先将 1 个盘子通过 4 根柱子移动到非目标柱子, // 然后将 N - 1 个盘子通过 3 根柱子移动到目标柱子,最后把最先移动的盘子移动到目标柱子上,则总步 // 数为 H(N - 1) + 2 * T(1),若第一次移动两个盘子,则步数为 H(N - 2) + 2 * T(2),依 // 此类推,所求的就是当 N 一定时,2 * T(M) + H(N - M)(0 <= M <= N) 的最小值。由于最多 // 可有 10000 个盘子,移动方法数将是一个大整数,超出 long long 型的表示范围。实际上,也可以不 // 需要计算比较最小移动步数,可以观察到 T(N) 的后一项与前一项的差形成以下数列 (N >= 1): // // 1 2 2 4 4 4 8 8 8 8 16 16 16 16 16 32 32 32 32 32 32 …… // // 2 的 m 次幂出现 (m + 1) 次,m >= 0。 #include <iostream> #include <vector> #include <iomanip> #include <algorithm> using namespace std; class integer { friend ostream& operator<<(ostream&, const integer&); public: integer() { }; // 将无符号整数转换为大整数。 integer(unsigned int a) { if (a == 0) digits.push_back(0); else { while (a) { digits.push_back(a % base); a /= base; } } }; ~integer() { }; integer operator*(const integer&); integer& operator*=(unsigned int); integer& operator*=(const integer&); integer operator+(integer&); private: void zero_justify(void); vector < unsigned int > digits; // 数位。 static unsigned int const base = 10000; // 基数。 static unsigned int const width = 4; // 数位宽度。 }; // 重载输出符号 <<。 ostream& operator<<(ostream& os, const integer& number) { os << number.digits[number.digits.size() - 1]; for (int i = number.digits.size() - 2; i >= 0; i--) os << setw(number.width) << setfill('0') << number.digits[i]; return os; } // 移除大数运算产生的前导 0。 void integer::zero_justify(void) { for (int i = digits.size() - 1; i >= 1; i--) { if (digits[i] == 0) digits.erase(digits.begin() + i); else break; } } integer& integer::operator*=(const integer& b) { return *this = *this * b; } // 大数乘以无符号整数。 integer& integer::operator*=(unsigned int b) { if (b >= base) { *this *= integer(b); return *this; } int carry = 0; for (int i = 0; i < digits.size(); i++) { digits[i] = digits[i] * b + carry; carry = digits[i] / base; digits[i] %= base; } if (carry) digits.push_back(carry); return *this; } // 乘法。 integer integer::operator*(const integer& b) { integer c; // 预先分配足够空间。 c.digits.resize(digits.size() + b.digits.size()); fill(c.digits.begin(), c.digits.end(), 0); // 一行一行算,逐行相加。 for (int i = 0; i < b.digits.size(); i++) for (int j = 0; j < digits.size(); j++) { c.digits[i + j] += digits[j] * b.digits[i]; c.digits[i + j + 1] += c.digits[i + j] / base; c.digits[i + j] %= base; } // 去掉前导0。 c.zero_justify(); return c; } // 加法。 integer integer::operator+(integer& b) { integer c; c.digits.resize(max(digits.size(), b.digits.size()) + 1); fill(c.digits.begin(), c.digits.end(), 0); int carry = 0, x, y, t; for (int i = 0; i < c.digits.size(); i++) { x = y = 0; if (i < digits.size()) x = digits[i]; if (i < b.digits.size()) y = b.digits[i]; t = x + y + carry; carry = t / base; c.digits[i] = t % base; } // 清除前导 0。 c.zero_justify(); return c; } #define MAXDISKS 10001 // 利用最少步骤数前后项相差所形成的数列特点。2 的 m 次幂出现 (m + 1) 次。 void hanoi(vector < int > disks) { // 求出相应盘子数的最少移动步骤数。 vector < integer > steps(MAXDISKS); steps[0] = 0; integer exponent(1); for (int index = 1, sequence = 1; index < MAXDISKS; sequence++) { for (int start = 1; start <= sequence && index < MAXDISKS; start++, index++) steps[index] = steps[index - 1] + exponent; exponent *= 2; } // 输出结果。 for (int i = 0; i < disks.size(); i++) cout << steps[disks[i]] << endl; } int main(int ac, char *av[]) { vector < int > disks; int d; // 读入盘子数。 while (cin >> d) disks.push_back(d); // 求最少移动步骤数。 hanoi(disks); return 0; }