TP-02/isqrt/isqrt.c

#include <stdio.h>
#include "isqrt.h"

unsigned long isqrt(unsigned long n)
{
    unsigned long r = n;
    while (r*r > n)
    {
        r = r + n/r;
        r = r/2;
    }
    return r;
}