求 v=(a^e)mod p

#include<stdio.h>
unsigned long Func(unsigned long a,unsigned long e,unsigned long p);
int main()
{
 unsigned long p,a,e,v;
    //求 v=(a^e)mod p
    printf("The program is to compute v=(a^e)mod p\n\n");

 printf("Please input the arguments(form: p,a,e):\n");
 scanf("%lu,%lu,%lu",&p,&a,&e);    
 v=Func(a,e,p);
 printf("(%lu^%lu)mod %lu = %lu\n",a,e,p,v);
}
unsigned long Func(unsigned long a,unsigned long e,unsigned long p) //求 v=a(e次幂)mod p
{
 unsigned long r=1;//返回 r
 unsigned long i=0;
 while(i<e)
 {
  r=r*a;
  r=r % p; 
  i++;
 }
 return r%p;
}