-nan(ind)

#define _CRT_SECURE_NO_WARNINGS
#include <iostream>
#include <fstream>
#include <iomanip>
#include <cmath>
 
long double powDouble(long double a, long double b) {
    long double k = a;
    for(int i = 2; i <= b; ++i) {
        a *= k;
    }
    return a;
}
 
int main() {
    std::ofstream out("c:\\out.txt");
    long double lambda = 2., tetta = 400., PiConst = 0.04, beta = 0.6, So = 0.02, L = 750., alpha = 0.5;
    long double y, PiSave, S = 1000000.;
    int VSave, ESave;
    for(int V = 1; V <= 200; ++V) {
        y = lambda * tetta / (long double)V;
        for(int E = 1; E <= 200; ++E) {
            long double Pi;
            if(y != 1.L) {
                Pi = powDouble(y, (long double)E) * (1.L - y) / (1.L - powDouble(y, long double(E + 1)));
            } else {
                Pi = 1.L / (long double)(E + 1);
            }
            out  << std::setprecision(10) << E << "     " << V << "     " << Pi << "   " << PiConst << std::endl;
            if(Pi >= PiConst) {
                
                continue;
            }
            long double STemp = alpha * powDouble((long double)V, (long double)beta) + So * E * L; // S = Sн + Sп = alpha * V^beta + So*E*L
            if(STemp < S) {
                S = STemp;
                PiSave = Pi;
                VSave = V;
                ESave = E;
            }
            //out << E << " " << V << std::endl;
        }
    }
    //out << ESave << " " << VSave << " " << PiSave << " " << S << std::endl;
    return 0;
}