数论 MillerRabin hdu 3792

本文介绍了一个使用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;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值