Длинная арифметика. Умножение двух длинных чисел.

#include <iostream>
#include <vector>
#include <algorithm>
void readlong(std::vector<int> &a); //чтение
void writelong(std::vector<int> a){  //вывод
    for (size_t i=0; i< a.size();i++) std::cout << a[i];}
short compare(std::vector<int> a,std::vector<int> b); //сравнение
void nol(std::vector<int> &a,std::vector<int> &b); //добавление нулей в начало(для суммы)
std::vector<int> sum(std::vector<int> a, std::vector <int> b); //сумма
std::vector<int> smult(std::vector<int> a,int b); //умножение длинного на короткое
void ten(std::vector<int> &a,int b);//добавление b нулей в конец
long long vectoint(std::vector<int> a);//перевод вектора в int
std::vector<int> mult (std::vector<int> a, std::vector<int> b); //умножение длинного на длинное
std::vector<int> factorial(std::vector<int> b); //факториал
int main(){
 
    std::vector<int> a,b,c,rez;
    readlong(a);
    writelong(factorial(a));
    return 0;
}
void readlong(std::vector<int> &a){
    char c;
    for(;;){
        std::cin.get(c);
        if ((c==' ')||(c=='\n'))break;
        a.push_back(c-48);
    }
}
short compare(std::vector<int> a,std::vector<int> b){
    if (a.size()>b.size()) return 1;
    if (a.size()<b.size()) return 2;
    for (size_t i=0; i<a.size();i++){
    if (a[i]>b[i]) return 1;
    if (a[i]<b[i]) return 2;
    }
    return 0;
 
}
void nol(std::vector<int> &a,std::vector<int> &b){
    if (a.size()==b.size()) return;
    if (a.size()>b.size()) std::swap(a,b);
    for (size_t i=0; i< a.size()-b.size(); i++){
        a.insert(a.begin(),0);
    }
}
std::vector<int> sum(std::vector<int> a, std::vector<int> b){
    std::vector<int>c;
    short p=0;
    nol(a,b);
    c.assign(a.size(),0);
    for (short i = a.size()-1; i >= 0; i--) {
        c[i]=a[i]+b[i]+p;
        p=0;
        if (c[i]>9) {c[i]-=10; p=1;}
    }
    if (p) c.insert(c.begin(),1);
    return c;
}
std::vector<int> smult(std::vector<int> a,int b){
    std::vector<int> c;
    if (!b||!a[0]) {c.push_back(0); return c; }
    for (size_t i = 0; i <a.size(); i++) {
        if (i==c.size()) c.push_back(0);
        c[i]=a[i]*b;
    }
    for (size_t i=c.size()-1; i > 0; i--) {
        if (c[i]>9) {
            c[i-1]+=c[i]/10;
            c[i]%=10;
            }
        }
    if (c[0]<10) return c;
    short p=c[0];
    c[0]%=10;
    for(;;){
        if (!(p/10)) break;
        c.insert(c.begin(),p/10);
        p/=10;
    }
    return c;
}
void ten(std::vector<int>& a,int b){
    for (; b; b--)
        a.push_back(0);
}
std::vector<int> factorial(std::vector<int>(b)){
    std::vector<int> a,o,c;
    a.push_back(1);
    o.push_back(1);
    c.push_back(1);
    for(long int i=0;i<vectoint(b)-1;i++){
        c=sum(c,o);
        a=mult(c,a);
    }
    return a;
}
long long vectoint(std::vector<int> a){
    long long c=0,q=1;
    for (short i = a.size()-1; i >=0; i--) {
        c+=a[i]*q;
        q*=10;
    }
    return c;
}
std::vector<int> mult (std::vector<int> a, std::vector<int> b)
{
    reverse(a.begin(),a.end());
    reverse(b.begin(),b.end());
    std::vector<int> c(a.size()+b.size());
    for (int i = 0; i < a.size(); ++i)
        for (int j = 0, carry = 0; j < b.size() || carry; ++j)
        {
            long long cur = c[i+j] + a[i] * (j < b.size() ? b[j] : 0) + carry;
            c[i+j] = cur % 10;
            carry = cur / 10;
        }
 
    while (c.size() > 1 && c.back() == 0)
        c.pop_back();
        reverse(c.begin(),c.end());
        return c;
}

char *My_Factorial(int f)
{
    char *name;
    char C[4],A[10000],B[10000],D[10000];
    char E[]="123456789";
    int i,j,k,l,m,n,o,p,q,r,s,t;
    for(i=0;i<9999;i++)
    {
        A[i]='0';
        B[i]='0';
        D[i]='0';
    }
    A[9998]='1';
    A[9999]='\0';
    B[9999]='\0';
    D[9999]='\0';
    //f -> scanf
    for(i=1;i<=f;i++)
    {
        itoa(i,C,10);
        k=strlen(C);
        for(j=k;j>0;j--)
        {
            n=(C[j-1]-48);
            s=(abs(j-k));
            for(l=9998;l>=3;l--)
            {
                m=A[l]-48;
                o=B[l-s]-48;
                p=n*m+o;
                r=p%10;
                q=(int)(p-r)/10;
                if(p<10)
                {
                    r=p;
                    q=0;
                }
                B[l-s]=r+48;
                B[l-s-1]=q+48;
            }
            for(l=9999;l>0;l--)
            {
                m=B[l-1]-48;
                o=D[l-1]-48;
                p=m+o;
                r=p%10;
                q=(int)(p-r)/10;
                if(p<10){r=p;q=0;}
                D[l-1]=r+48;
                D[l-2]+=q;
            }
            for(l=0;l<9999;l++)B[l]='0';
        }
        strcpy(A,D);
        for(l=0;l<9999;l++)D[l]='0';
    }
    int xz=strlen(A);
    name=new char[xz];
    name[0]=NULL;
    name=strpbrk(A,E);
    name[0]=NULL;
    return name;
}

#include <iostream>
using std::cin;
using std::cout;
using std::endl;
 
#include <string>
using std::string;
 
#include <algorithm>
 
#include <windows.h>
 
string& multOneDigit(const string& value1, const char& value2)
{
    static string result;
    result = "";
    unsigned digit(value2 - '0');
    switch (digit)
    {
    case 0:
        return result = "0";
        break;
    case 1:
        return result = value1;
        break;
    }
    unsigned temp(0);
    for (unsigned i(0); value1[i]; ++i)
    {
        temp += (value1[i] - '0') * digit;
        result.push_back(temp % 10 + '0');
        temp /= 10;
    }
    if (temp)
        result.push_back(temp % 10 + '0');
    return result;
}
 
string& multDigitWithZeros(const string& value1, const string& value2)
{
    static string result;
    result = multOneDigit(value1, value2[0]);
    if (value2.length() == 1)
        return result;
    return result = string(value2.begin() + 1, value2.end()) + result;
}
 
string& toAdd(string& value1, string& value2)
{
    static string result;
    result = "";
    unsigned i(value1.length()), j(value2.length());
    if (i < j)
        for (; i < j; ++i)
            value1.push_back('0');
    else if (j < i)
        for (; j < i; ++j)
            value2.push_back('0');
    unsigned temp(0);
    for (i = 0; i < j; ++i)
    {
        temp += value1[i] + value2[i] - 2 * '0';
        result.push_back(temp % 10 + '0');
        temp /= 10;
    }
    if (temp)
        result.push_back(temp % 10 + '0');
    return result;
}
 
const string& toMult(string& value1, string& value2)
{
    static string result;
    result = "0";
    if (value1.length() < value2.length())
        swap(value1, value2);
    reverse(value1.begin(), value1.end()),
        reverse(value2.begin(), value2.end());
    for (unsigned i(0), count(value2.length()); i < count; ++i)
    {
        string temp(&value2[i], 1);
        for (unsigned j(0); j < i; ++j)
            temp.push_back('0');
        result = toAdd(multDigitWithZeros(value1, temp), result);
    }
    reverse(result.begin(), result.end());
    return result;
}
 
int main(void)
{
    SetConsoleCP(1251);
    SetConsoleOutputCP(1251);
    string value1("470827"), value2("60013");
    cout << toMult(value1, value2) << endl;
    system("pause");
    return 0;
}