for (long double N = 1; N < sqrt(2*S)+3; N++) {
}