Kevin_Zhang — Thu, 19 Aug 2010 17:14:00 GMT

#include"iostream"
#include"stdio.h"
#include"stdlib.h"
#include"time.h" //�pȝ��旉��函数的头文�g
using namespace std;
int T;//the number of tests
int j,d,N;

void computejd(int n)//求j和d,�q�个是正��的�Q�d为odd
{int m=n-1;
j=0;
while(1)
{
   if(m%2==0){j++;m=m/2;}
   else {d=m;break;}

}
return ;
}//right

int GetRand()//求随机数�Q�种子�ؓ�pȝ��旉��
{
srand(time(NULL)) ;
int n = rand()%(N-1)+1;

return n ;
}//right

int main()
{
int k;
printf("输入��试�ơ数T�Q?/span>");
scanf("%d",&T);//right
for(int x=1;x<=T;x++)
  {scanf("%d",&N);
   computejd(N);
   for(k=1000;k>0;k--)
    {int b=GetRand();
     int r0, r,r1,i;
     r0=N-1;i=0;
     r=1; r1=b%N;
          for(int q=0; q<=31;q++)
          { if((d&1)==1) r=(r*r1)%N;
             r1=(r1*r1)%N; d>>=1;
          }

           while(r!=1&&i<j)
          {
              r0=r;r=(r*r)%N;i++;

          }
         if(r==1||i>=j)
          {
             if(r==1&&r0==N-1)continue;
             else break;
          }


       }

   if(k==0)printf("prime\n");
   else printf("not prime\n");

  }//outside cycle
return 0;
}

Kevin_Zhang 2010-08-20 01:14 发表评论

欧美午夜精品理论片a级大开眼界 欧美午夜精品理论片a级按摩 ,免费久久精品视频,欧美gay视频激情

欧美午夜精品理论片a级大开眼界欧美午夜精品理论片a级按摩 ,免费久久精品视频,欧美gay视频激情