#include <iostream>
#include <cmath> // for std::sqrt

int main(int argc, char** argv)
{
    if (argc < 2)
    {
	std::cerr << "Usage: " << argv[0] << " FILE" << std::endl;
	return 1;
    }

    // filename is stored in argv[1]

    // implement the steps here
}