#include <iostream>
#include <cmath>
using namespace std;
// x * min + y * max + xi = b * sqrt(a)
// x + y = m - 1

int main() {
    int m, p, a, b, cantMax, cantMin;
    double max, min, xi, sumaCuad;

    while (cin >> m >> p >> a >> b) {
        max = sqrt(a);
        min = -1 / max;
        cantMax = (int) ((b * max - m * min) / (max - min));
        cantMin = m - 1 - cantMax;
        xi = b * max - cantMax * max - cantMin * min;
        sumaCuad = pow(xi, p) + cantMax * pow(max, p) + cantMin * pow(min, p);
        cout << (int) (sumaCuad + 0.5) << endl;
    }
    return 0;
}