Длинная арифметика: вычислить 100! - 2^100 - C++ - Ответ 1511853

#include <iostream>
#include <conio.h>
#include <locale.h>
using  namespace std;
void revers(char * s)
{
    int n = strlen(s);
    char c;
    for(int i = 0 ; i< int((n+1)/2); i++)
    {
        c = s[i];
        s[i] = s[n-1-i];
        s[n-i-1]= c;    
    }
}
 
void sum(char * str,const char * str2)
{
    char s2[1000] = "";
    char s1[1000] = "";
    strcat(s1,  str );
    strcat(s2,  str2);
    int i,n1 = strlen(s1);
    int sum,n2 = strlen(s2);
    revers(&s1[0]);
    revers(&s2[0]);
    bool f = false;
    while(true)
    {
        if(i == n1)
        {
            if(i == n2) break;
            strcat(s1, "0");
            n1++;
        }
        if(i == n2)
        {
            if(i == n1) break;
            strcat(s2, "0");
            n2++;
        }
        sum = s1[i] + s2[i] - 96;
        sum +=(f)?1:0;
        if (sum >= 10)       f = true;
        else                 f = false;
        if (f)              sum -=10;
        s1[i] = sum+48; 
        i++;
    }
    if(f) strcat(s1,"1");
    revers(&s1[0]);
    strcpy(str, s1);
}
void mul (char * str,const char * str2)
{
    char s2[1000] = "";
    char s1[1000] = "";
    char s4[1000] = "0";
    char s3[1000] = "0";
    char c [1000] = "";
    char c1 = '1';
    strcat(s1,  str );
    strcat(s2,  str2);
    revers(&s2[0]);
    int i = 0,
        n1 = strlen(s1);
    int n2 = strlen(s2);
    while(i<n2)
    {
        char ch;
        ch = s2[i];
        while(ch != '0')
        {
            sum(s3, s1);
            ch = ch - c1 + 48;          
        }
        strcat(s3,c);
        sum(s4, s3);
        strcpy(s3, "");     
        strcat(c,"0");
        i++;
    }
    strcpy(str, s4);
}
void fact(char * str, const char * str2)
{
    char s1[1000] = "1";
    char s2[4] = "1";
    while(strcmp(s2, str2))
    {
        mul(s1, s2);
        sum(s2, "1");
    }
    mul(s1, str2);
    strcpy(str, s1);
}
void power(char * str, int N)
{
    char t[1000] = "";
    char s1[1000] = "";
    char s2[1000] = "";
    strcat(t,   str );
    int i = 0,n;
    int M = N;
    while (M!=0)
    {
        if(M%2!=0)
            strcat(s2,"1");
        else
            strcat(s2,"0");
        M = (int) M/2;
    }
    strcpy(s1, "1");
    n = strlen(s2);
    while (i<n) 
    {    
        if (s2[i]=='1') mul(s1, t);
        mul(t, t);
        i++;
    }
    strcpy(str, s1);
}
void sub(char * str,const char * str2)
{
    char s2[1000] = "";
    char s1[1000] = "";
        strcat(s1,str);
        strcat(s2, str2) ;
        int i,n1 = strlen(s1);
        int subject,n2 = strlen(s2);
        revers(&s1[0]);
        revers(&s2[0]);
        bool f = false;
        while(true){
            if(i == n1){
                if(i == n2) break;
                strcat(s1, "0");
                n1++;
            }
            if(i == n2){
                if(i == n1) break;
                strcat(s2, "0");
                n2++;
            }
            subject = s1[i] - s2[i];
            subject -=(f)?1:0;
            if (subject < 0)       f = true;
            else                 f = false;
            if (f)              subject +=10;
            s1[i] = subject+48; 
            i++;
        }
        i = n1-1;
        while(s1[i] == '0' && i > 0){
            s1[i] = '\0';
            i--;
        }
        revers(&s1[0]);
        strcpy(str, s1);
}
void main()
{
    setlocale(0 ,"russian");
    cout<<"100! - 2^100\n\n";
    char str[1000] = "";
    char str2[1000] = "2";
    
    getch();
    fact(str, "100"); //   100!
    power(str2, 100); //   2^100
    sub(str, str2);   //   100! - 2^100
    cout<<str;
    getch();
    return;
}