Возведение больших чисел в степень

#include <iostream>
#include <math.h>
#include <time.h>
using namespace std;
#define DIM 50
 
long int            modPow(int a, int b, int n); //àëãîðèòì áûñòðîãî âîçâåäåíèÿ â ñòåïåíü
int         sD(int a, int b, int i);
int         SmpDig(int num); // ïðîâåðÿåì ÷èñëî íà ïðîñòîòó
long int    SD(void);
void        setD(long int *d); // çàäàåì ðàíäîìíûå çíà÷åíèÿ â p è q
void        printD(long int *q, long int *p); // ñîîòâåòñòâåííî âûâîäè èõ íà ýêðàí
 
void printD(long int *q, long int *p)
{
    for(int i=0;i<9;i++)
        printf("*");
    printf("\n");
 
    printf("q = %ld ",*q);
    printf("p = %ld ",*p);
 
    printf("\n");
 
    for(int i=0;i<9;i++)
        printf("*");
    printf("\n");
 
}
 
void setD(long int *d)
{
    *d = SD();
}
 
long int SD(void) 
{
 
    long int d = rand()%DIM;
 
    while(SmpDig(d)!=1)
    {
        d = rand()%DIM;
    }
 
    return d;
 
}
 
int SmpDig(int num){
 
    int pr = 0;
 
    for(int i=1; i<=num;i++)
    {
        if((num % i) == 0)pr++;
    }
 
    if(pr ==2 )
        return 1;
    else
        return 0;
}
 
int gcd(int a, int b)
{
    int r;
    while(b!=0)
    {
        r = a%b;
        a = b;
        b = r;
    }
    return a;
}
 
int algev(int a, int b) 
{
    int i;
 
    for (i=2; i<b; i++) 
    {
        if (gcd(a,i) == 1) return i;
    }
    return 0;
}
 
int oalgev(int a, int b) 
{
    int u[]={a, 0};
    int v[]={b, 1};
    int t[]={0, 0};
    int q;
    while (v[0] != 0) {
        q = u[0]/v[0];
        t[0] = u[0]%v[0];t[1]=u[1]-q*v[1];
        u[0] = v[0]; u[1] = v[1];
        v[0] = t[0]; v[1] = t[1];
    }
    if (u[1] < 0) u[1] += a;
    return u[1];
}
 
long int modpow(long int x, long int e, int n) 
{
    int r = 1;
    while( e>0 ) 
    {
        if( (e%2)==1 ) { r = (r*x) % n; }
        e = e/2;
        x = (x * x) % n;
    }
    return r;
}
 
long sign (long msg, int d, int n)  // ñîçäàíèå öèôðîâîé ïîäïèñè
{   
    return modpow(msg,d,n);
}
 
long verify ( long msg,long sign, int e, int n ) // è å¸ âåðèôèêàöèÿ
{   
    if (msg==modpow(sign,e,n)) 
    {
        printf("VERIFICATION OK\n");
        return 1;
    }
    else printf("VERIFICATION FAILED, GOD DAMMIT\n");
    return 0;
 
}
    
       int main(int argc, const char * argv[])
{
 
    srand(time(NULL));
 
    long int p,q,n,m,d,e,msg,code, dcode, msg2,  s, ver ;
 
    setD(&p);
    setD(&q);
 
    n = p*q;
    m = ((p-1)*(q-1));
 
    e = algev(m,m);
    d = oalgev(m, e);
    printf ("Enter ur message for encodiing");
    cin>>msg;                
    code = modpow(msg, e, n);
    dcode = modpow(code, d, n);
    
    printf("Open keys: %ld %ld\n", n, e);
    printf("Close keys: %ld %ld\n", n, d);
    printf("Msgs: %ld %ld\n", code, dcode);
    printD(&q,&p);
        
    printf("Make ur digital signature. Enter ur msg\n");
    cin>>msg2;
    s=sign(msg2,d,n);
    printf ("Ur digital signature =%d\n",s);
    
    printf("So, you can proof my signature, enter your message and signature\n");
    cin>>msg2>>s;
    ver=verify(msg2,s,e,n);
        
    
    
    system("PAUSE");
    return 0;
}