#include <iostream>
#include <iomanip>
#include <cmath>
#include <vector>
#include <string>
#include <chrono>

using namespace std;

class muTimer
{
    using Clock = std::chrono::high_resolution_clock;
    bool active = false;
    Clock::duration   duration_;
    Clock::time_point start_ = Clock::now(), stop_ = Clock::now();

    muTimer(const muTimer&)             = delete;
    muTimer& operator=(const muTimer&)  = delete;
public:
    using ns       = std::chrono::nanoseconds;
    using mks      = std::chrono::microseconds;
    using ms       = std::chrono::milliseconds;
    muTimer() { reset(); start(); }
    ~muTimer() = default;
    muTimer& reset()
    {
        duration_ = std::chrono::nanoseconds(0);
        active    = false;
        return *this;
    }
    muTimer& start()
    {
        if (!active)
        {
            start_ = Clock::now();
            active = true;
        }
        return *this;
    }
    muTimer& stop()
    {
        if (active)
        {
            stop_      = Clock::now();
            duration_ += stop_ - start_;
            active     = false;
        }
        return *this;
    }
    template<typename T = mks>
        unsigned long long duration()
    {
        return static_cast<unsigned long long>
            (std::chrono::duration_cast<T>(stop_-start_).count());
    }
};


int main(int argc, const char * argv[])
{

    const int N = 10000000;
    {
        muTimer mt;
        double sum = 0;
        for(int i = 0; i < N; ++i)
        {
            double x = i/(1.0*N), y = N/double(N+i);
            sum += pow(x,y);
        }
        mt.stop();
        cout << mt.duration<>() << " mks\n";
        cout << sum << endl;
    }
    {
        muTimer mt;
        double sum = 0;
        for(int i = 0; i < N; ++i)
        {
            double x = i/(1.0*N), y = N/double(N+i);
            sum += exp(log(x)*y);
        }
        mt.stop();
        cout << mt.duration<>() << " mks\n";
        cout << sum << endl;
    }


}
