Возведение числа в квадрат (длинная арифметика)

#include <iostream>
#include <string>
using namespace std;
 
const int NMAX = 10000;
 
void zero(int *A)
{
    for (int i = 0; i < NMAX; i++)
        A[i] = 0;
}
 
void TakeChislo(string x, int &nA, int *A)
{
    int len;
    len = x.length();
    for (int i = len - 1; i >= 0; i--)
        A[nA++] = x[i] - '0';
}
 
void Sum(int nA, int *A, int nB, int *B, int &nC, int *C)
{
    int max, i;
    max = nA>nB ? nA : nB;
    for (i = 0; i < max; i++)
        C[i] = A[i] + B[i];
    for (i = 0; i < max; i++)
    {
        C[i + 1] += C[i] / 10;
        C[i] %= 10;
    }
    if (C[max] >0) nC = max + 1;
    else nC = max;
}
 
void sdvig(int &nA, int *A, int k)
{
    for (int i = nA - 1; i >= 0; i--)
        A[i + k] = A[i];
    for (int i = 0; i < k; i++)
        A[i] = 0;
    nA += k;
}
 
void Ymnozh(int nA, int *A, int nB, int *B, int &nC, int *C)
{
    int max, i;
    int *X1, *X2, D[NMAX];
    int nX1, nX2, nD;
    zero(C); zero(D);
    max = nA>nB ? nA : nB;
    if (max == nA)
    {
        X1 = A; nX1 = nA;
        X2 = B; nX2 = nB;
    }
    else
    {
        X1 = B; nX1 = nB;
        X2 = A; nX2 = nA;
    }
    for (int j = 0; j < nX2; j++)
    {
        for (i = 0; i < nX1; i++)
            D[i] = X1[i] * X2[j];
 
        for (i = 0; i < nX1; i++)
        {
            D[i + 1] += D[i] / 10;
            D[i] %= 10;
        }
        if (D[nX1] != 0) nD = nX1 + 1;
        else nD = nX1;
        sdvig(nD, D, j);
        Sum(nC, C, nD, D, nC, C);
    }
}
 
void Print(int nA, int *A)
{
    for (int i = nA - 1; i >= 0; i--)
        cout << A[i];
    cout << endl;
}
 
void main()
{
    int A[NMAX], B[NMAX], C[NMAX];
    zero(A); zero(B); zero(C);
    int N,nA = 0, nB = 0, nC = 0;
    string s;
    cin >> s;
    TakeChislo(s, nA, A);
    for (int i = 0; i < nA; i++)
    {
        B[i] = A[i];
        nB = nA;
    }
    Ymnozh(nA, A, nB, B, nC, C);
    Print(nC, C);
    system("pause");
}