agc007C - Pushing Balls(期望 等差数列)

#include<bits/stdc++.h>
using namespace std;
long double N, d1, x, ans;
int main() {
    cin >> N >> d1 >> x;
    for(int i = N; i >= 1; i--) {
        long double s = d1 * 2 * N + N * (2 * N - 1) * x;
        ans += s / 2 / N;
        s = s - (4 * d1 + 4 * N * x - 2 * x)  / 2 / N;
        d1 = ((2 * N - 2) * d1 + d1 + 2 * x + 3 * d1 + 3 * x) / 2 / N;
        N--;
        x = (s - 2 * N * d1) / N / (2 * N - 1);
    //    if(i > 990) printf("%.10lf\n", (double)x);
    }
    printf("%.15lf", (double)ans);
    return 0;
}