#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;
}
}