数论 MillerRabin hdu 3792

原创于 2011-03-14 08:53:00 发布 · 1k 阅读

0 ·

CC 4.0 BY-SA版权

Algorithm 专栏收录该内容

20 篇文章

订阅专栏

本文介绍了一个使用Miller-Rabin素性测试算法来生成素数表，并基于此表找出所有小于特定数值的孪生素数对的方法。通过__int64类型支持大整数运算，利用ModExp函数实现快速幂取模操作，提高素性测试效率。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include<iostream>
#include<cmath>
using namespace std;
int p[10002];
int ans[2002];
int len;
__int64 ModExp( __int64 a , __int64 b , __int64 n )
{
   __int64 m = 1;
   while( b )
   {
    if( b & 1 )
    { m = ( m % n ) * ( a % n ) % n; }
    b >>= 1;
    a = ( a % n ) * ( a % n ) % n;
   }
   return m;
}
bool miller_rabbin(__int64 n)
{
__int64 i , a;
for( i = 0 ; i < 50 ; i++ )
{
  a = rand() % ( n - 2 ) + 2;
  if( ModExp( a , n - 1 , n ) != 1 )
  { return false; }
}
return true;
}
void Prime()
{
__int64 i;
p[0] = 2;
for( i = 3 ; i <= 100000 ; i = i + 1 )
{
  if( miller_rabbin( i ) )
  { p[len++] = i; }
}
}
int main()
{
int i , j;
int n , k;
int sum;
len = 0;
k = 0;
Prime();
memset( ans , 0 , sizeof( ans ) );
for( j = 0 ; j < len - 1 ; ++j )
{
  if( p[j + 1] - p[j] == 2 )
  { ans[k++] = p[j + 1]; }
}
while( scanf( "%d" , &n ) )
{
  if( n < 0 )
  { break; }
  sum = 0;
  for( i = 0 ; i < k ; ++i )
  { if( ans[i] > n ) break; }

  printf( "%d/n" , i );
}
return 0;
}