RSA с длинными ключами

#include <iostream>
using namespace std;
 
string reduce(string const &a)
{
    string sign;
    string z;
    int p;
    
    if (a[0]=='-') 
    {
       sign=string("-");
       z=a.substr(1);
    }
    else
    {
       z=a;    
       sign=string("");
    }
 
    p=-1;
    for (int i=0; i<z.size(); i++)
        if (z[i] != '0')
        {
            p=i;
            break;
        }
    if (p==-1) return string("0");    
    return sign+z.substr(p);
}
 
 
string longSub(string const &a1, string const &a2)
{
     int sz1=a1.size();
     int sz2=a2.size();
     int sz;
     string aa1,aa2,sgn;
     int ii,q1,q2,b;
     
     if (sz1==sz2)
     {
         aa1=a1;
         aa2=a2;
         sz=sz1+1;
     }
     else if (sz1 > sz2)
     {
         aa1=a1;
         aa2=string(sz1-sz2,'0')+a2;
         sz=sz1+1;
     }
     else
     {
         aa1=string(sz2-sz1,'0')+a1;
         aa2=a2;
         sz=sz2+1;
     }
     
     if (aa1 > aa2)
        sgn="";
     else
     {
         swap(aa1,aa2);
         sgn="-";
     }
 
    char tmp[sz];   
 
    for(ii=0; ii<sz; ii++) tmp[ii]=' ';
    tmp[sz-1]=0;
    
    ii=sz-2;
    b=0;
    
    while(1)
    {
        q1=aa1[ii]-'0'-b;
        q2=aa2[ii]-'0';
        if (q1>=q2)
        {
           tmp[ii]=(q1-q2)+'0';
           b=0;
        }
        else
        {
           tmp[ii]=(q1+10-q2)+'0';
           b=1;
        }
        ii--;
        if (ii<0) break;
    }
    
    return sgn+string(tmp);
}
 
string longAdd(string const &a1, string const &a2)
{
    string res;
    int sz1=a1.size();
    int sz2=a2.size();
    int sz=max(sz1,sz2)+2;
    char tmp[sz];
    int i1,i2,ii,q1,q2,z,p;
    
    for(ii=0; ii<sz; ii++) tmp[ii]=' ';
    tmp[sz-1]=0;
 
    i1=sz1-1;
    i2=sz2-1;
    ii=sz-2;
    
    p=0;
    while(1)
    {
        if (i1>=0)
           q1=a1[i1]-'0';
        else
           q1=0;
 
        if (i2>=0)
           q2=a2[i2]-'0';
        else
           q2=0;
        
        z=p+q1+q2;
        tmp[ii]=z%10+'0';
        p=z/10;
        
        ii--;
        i1--;
        i2--;
        if (i1<0 && i2<0) break;
    }
    
    if (p>0) tmp[ii]=p+'0';
    
    for (ii=0; ii<sz; ii++)
       if (tmp[ii] != ' ')
       {
           res=string(tmp+ii);
           break;
       }
 
    return res;
    
}
 
string longDiv(string const &a1, string const &a2)
{
    if (a2=="0") throw "Длинное деление на нуль";
    
    string aa1=a1;
    string aa2=a2;
    
    int l1=aa1.size();
    int l2=aa2.size();
    
    if (l1>l2)
       aa2=string(l1-l2,'0')+aa2;
    
    if (l2>l1)
       aa1=string(l2-l1,'0')+aa1;
    
    string res("0");
 
    while(aa1>=aa2)
    {
        aa1=longSub(aa1,aa2);
        res=longAdd(res,"1");
    }
    
    return reduce(res);
}
 
string longRem(string const &a1, string const &a2)
{
    if (a2=="0") throw "Длинное деление на нуль";
    
    string aa1=a1;
    string aa2=a2;
    
    int l1=aa1.size();
    int l2=aa2.size();
    
    if (l1>l2)
       aa2=string(l1-l2,'0')+aa2;
    
    if (l2>l1)
       aa1=string(l2-l1,'0')+aa1;
    
    //string res("0");
 
    while(aa1>=aa2)
    {
        aa1=longSub(aa1,aa2);
        //res=longAdd(res,"1");
    }
    
    return reduce(aa1);
}
 
string longGcd(string const &a1, string const &a2)
{
 
    string aa1=a1;
    string aa2=a2;
    string tmp;
    string zero=string("0");
    
    while (aa2 != zero)
    {
        tmp=longRem(aa1,aa2);
        aa1=aa2;
        aa2=tmp;
    }
    
    return aa1;
}
 
string longMult(string const &a1, string const &a2)
{
    string add("");
    string res("0");
    string tmp;
    int k,c,p,u;
    
    int sz2=a2.size();
    int sz1=a1.size();
 
    for (int i=sz2-1; i>=0; i--)
    {
        tmp=string("");
 
        k=(a2[i]-'0');
 
        if (k > 0)
        {
            p=0;
            for (int j=sz1-1; j>=0; j--)
            {
                c=(a1[j]-'0')*k+p;
                u=c%10;
                p=c/10;
                tmp=(char)(u+'0')+tmp;
            }
            
            if (p>0) tmp=(char)(p+'0')+tmp;
            
            tmp=tmp+add;
            
            res=longAdd(res,tmp);
 
        }
        
        add=add+"0";
    }
    
    return res;
}
 
string factorial(int n)
{
    string f("1");
    string a("1");
    
    for (int i=2; i<=n; i++)
    {
        a=longAdd(a,string("1"));
        f=longMult(f,a);
    }
    
    return f;
    
}
 
 
int main()
{
 
    // 300 первых чисел Фибоначчи
    
    string c="0", p="1", t;
    
    for (int i=1; i<=300; i++)
    {
        t=longAdd(c,p);
        
        cout << i << " " << t << endl;
        
        p=c;
        c=t;
        
    }
 
    return 0;
}