В чем причина глюка вычисления квадратного корня вычитанием нечётных чисел?

#include <stdio.h>
int h,w,k,z,e,g,j,i,y;                              //  c - остаток ; d - вычитаемое ;
int main() {
    char q,a[3101],c[1510],d[1510];
    while ((g=getchar())<49 || g>57); a[k++]=g-48;          // k - размер числа a (к-во цифр)
    while ((g=getchar())>47 && g<58) a[k++]=g-48;
    if (k&1) g=a[z++]; else g=a[z++]*10+a[z++];     // z - позиция следующей цифры (a[z])
    for (j=1; j<=g; j+=2) { g-=j; w++; }
    printf("%d",w); j--;
    if (z==k) return 0;
    if (j<9) { d[0]=j; d[1]=1; y=2; }                       // y - размер числа d (вычитаемое)
    else { d[0]=j/10; d[1]=j%10; d[2]=1; y=3; }
    c[0]=g/10; c[1]=g%10;                           // е - размер числа С с ведущими нулями
    c[2]=a[z++]; c[3]=a[z++]; e=g=4;                // g - чистый размер числа С (остаток)
    for (j=0; j<e; j++)                                // h - начало числа С без ведущих нулей
        if (c[j]) j=e; else { h++; g--; }
    while (1) {
        w=0;                                          // w - к-во вычитаний (след. цифра)
        while (1) {
            if (g<y) break;                          // если остаток < вычитаемого, то выход
            q=0;
            if (g==y) {
                i=0;
                for (j=h; j<e;) {
                    if (c[j]<d[i]) { q=1; j=e; }
                    if (c[j++]>d[i++]) j=e;
                }
            }
            if (q) break;                           // --------------//-------------
            i=y; w++;
            for (j=e; j-->h;) {                     // вычитаем : c = c - d ; 
                if (i) c[j]-=d[--i];
                if (c[j]<0) { c[j]+=10; c[j-1]--; }
            }
            for (j=h; j<e; j++) if (c[j]) j=e; else { h++; g--; }
            i=y-1;                              // убираем ведущие нули
            if (d[i]==9) { d[i]=1; d[i-1]++; }
            else d[i]+=2;                           // увеличиваем вычитаемое d на 2
        }
        printf("%d",w);
        if (z==k) break;
        c[e++]=a[z++]; c[e++]=a[z++]; g+=2;
        for (j=h; j<e; j++) if (c[j]) j=e; else { h++; g--; }
        d[y-1]--; d[y++]=1;
    }
    return 0;
}