for (j = 2; j < 2 * n; ++j) {
        for (i = 2 * n - 1; i >= j; --i) {
            if (i & 1)
                d[i >> 1] = (d[i >> 1] - y[i >> 1]) / (x[i >> 1] - x[(i - j) >> 1]);
            else
                y[i >> 1] = (y[i >> 1] - d[(i >> 1) - 1]) / (x[i >> 1] - x[(i - j) >> 1]);
        }
    }