Возведение в степень по модулю для большого числа

#include <vcl.h>
#pragma hdrstop
#include <iostream>
#include <math.h>
#include <conio.h>
using namespace std;
//---------------------------------------------------------------------------
#pragma argsused
long powmod(long a, long k9, long p)  // a - основание k9 степень p модуль
{
  long b=1;
  while (k9) {
    if (k9%2==0) {
      k9 /= 2;
      a = (a*a)%p;
      }
    else {
      k9--;
      b = (b*a)%p;
      }
  }
  return b;
}
int main()
{  long a,p;
 
  cout << "chislo dla testa: "; cin >> p ;
   a = rand() % (p-1);
  cout << "chislo a = " << a << endl ;
 //k= a ^((p-1) div 2) mod p;  - вычисляем по поэтапно
  double k0;
  k0=(p-1)/2;      // показатель степени
  int k01;
  k01 = (int) k0;       //  целая часть показателя степени (p-1) div 2
  cout<<"k01 = "<< k01 <<endl;
   //+++++
  long k9, k;
  k9= k01;
   k = powmod(a, k9, p);
   cout<<"!!!!! k = "<< k <<endl;
 
 
  int p01= p-1 ; // для условия
  if (k == 1) {
        cout<<"  Prostoe chislo "<< p  <<endl;
        }
        else if (k == p01 )
        {
        cout<<"  Prostoe chislo "<< p  <<endl;
        }
    else {
        cout<<" Ne prostoe chislo "<< p  <<endl;
    }
    getch(); getch();
    return 0;
}