B、生日蛋糕

减枝优化
1、如果直接暴力枚举全部可能的长度与半径，递归深度太大了，走到很远才发现不能到终点，
2、考虑减枝优化，最容易想到的就是当前如果大于了我们假设的直接即可。
3、再进一步减枝，首先预处理出每一层最少的表面积与体积，比如最顶层一定是，那么体积和面积公式算到最小值。如果再过程中，现在的 或者 加上最小的体积已经大于等于可以直接
4、）做到上面减枝居然还要T……最最最厉害的一个减枝，，可以发现变两倍，缩小4倍，而我们可以从上面式子直接算到后序最小的表面积，即当取到下面一层的的时候算到的表面积最小，原因是已经缩小最多倍，只是放大根号倍，如果这个时候加上也要大于的话，也可以直接

减枝……好累阿

#include <bits/stdc++.h>
#pragma GCC optimize(2)
#pragma GCC optimize(3)
using namespace std;
#define js ios::sync_with_stdio(false);cin.tie(0); cout.tie(0)
typedef long long ll;
const ll MOD = 1e9 + 7;
inline ll read() { ll s = 0, w = 1; char ch = getchar(); while (ch < 48 || ch > 57) { if (ch == '-') w = -1; ch = getchar(); }    while (ch >= 48 && ch <= 57) s = (s << 1) + (s << 3) + (ch ^ 48), ch = getchar();    return s * w; }
inline void write(ll x) { if (!x) { putchar('0'); return; } char F[200]; ll tmp = x > 0 ? x : -x; if (x < 0)putchar('-');    int cnt = 0;    while (tmp > 0) { F[cnt++] = tmp % 10 + '0';        tmp /= 10; }    while (cnt > 0)putchar(F[--cnt]); }
inline ll gcd(ll x, ll y) { return y ? gcd(y, x % y) : x; }
ll qpow(ll a, ll b) { ll ans = 1;    while (b) { if (b & 1)    ans *= a;        b >>= 1;        a *= a; }    return ans; }    ll qpow(ll a, ll b, ll mod) { ll ans = 1; while (b) { if (b & 1)(ans *= a) %= mod; b >>= 1; (a *= a) %= mod; }return ans % mod; }
inline int lowbit(int x) { return x & (-x); }

const int INF = 0x3f3f3f3f;
const int N = 25;
int n, m, ans = INF;
int mins[N], minv[N];

void dfs(int t, int s, int v, int lr, int lh) {
    if (!t) {
        if (v == n)    ans = min(ans, s);
        return;
    }
    if (v + minv[t] > n)    return;
    if (s + mins[t] >= ans)    return;
    if (s + 2 * (n - v) / lr >= ans)    return; // s=pi*2*r*h <-  v=pi*r*r*h
    for (int i = lr - 1; i >= t; --i) {
        if (t == m)    s = i * i; //计数最底下底面积
        int maxh = min(lh - 1, (n - v - minv[t - 1]) / (i * i));
        for (int j = maxh; j >= t; --j)
            dfs(t - 1, s + 2 * i * j, v + i * i * j, i, j);
    }
}

int main() {
    n = read(), m = read();
    for (int i = 1; i <= m; ++i) {
        mins[i] = 2 * i * i;
        mins[i] += mins[i - 1];
        minv[i] = i * i * i;
        minv[i] += minv[i - 1];
    }
    dfs(m, 0, 0, sqrt(n), n);
    if (ans == INF)    return puts("0"), 0;
    printf("%d\n", ans);
    return 0;
}