给定 22 个整数 a,ba,b 求出它们之间(包括 a,ba,b)的所有质数的和。
输入格式
两个整数 a,b\ (1\le a,b\le 10^5)a,b (1≤a,b≤105)。
输出格式
一个整数,表示范围内的质数和。
输出时每行末尾的多余空格,不影响答案正确性
样例输入
999 10
样例输出
76110
思路一:暴力判断+统计即可;---数据范围小,更适合;<50ms
#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long LL;
bool prime(LL n){
if(n==2) return true;
if(n==1) return false;
for(int i=2;i<=sqrt(n);i++)
if(n%i==0) return false;
return true;
}
int main(){
LL a,b,ans=0;
scanf("%lld%lld",&a,&b);
if(a>=b) swap(a,b);
for(LL i=a;i<=b;i++)
if(prime(i))
ans+=i;
printf("%lld\n",ans);
return 0;
}
思路二:Miller_Rabin算法;数据范围更大适合,此题时间<500ms
#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long LL;
LL mul_mod(LL a,LL b,LL mod){
LL ret=0;
while(b){
if(b&1) ret=ret+a;
if(ret>=mod) ret-=mod;
a=a+a;
if(a>=mod) a-=mod;
b>>=1;
}
return ret;
}
LL pow_mod(LL a,LL b,LL mod){
LL ret=1;
while(b){
if(b&1)ret=mul_mod(ret,a,mod);
a=mul_mod(a,a,mod);
b>>=1;
}
return ret;
}
bool Miller_Rabin(LL n){//判断素数
LL u=n-1,pre,x;
int i,j,k=0;
if(n==2||n==3||n==5||n==7||n==11) return true;
if(n==1||(!(n%2))||(!(n%3))||(!(n%5))||(!(n%7))||(!(n%11))) return
false;//初始化判断素数
for(;!(u&1);k++,u>>=1);//按照要求转换形式
for(i=0;i<10;i++){
x=rand()%(n-2)+2;//生成随机数
x=pow_mod(x,u,n);
pre=x;
for(j=0;j<k;j++){
x=mul_mod(x,x,n);
if(x==1&&pre!=1&&pre!=(n-1))//二次探测判断
return false;
pre=x;
}
if(x!=1) return false;//用费马小定理判断
}
return true;
}
int main(){
LL a,b,ans=0;
scanf("%lld%lld",&a,&b);
if(a>=b) swap(a,b);
for(LL i=a;i<=b;i++)
if(Miller_Rabin(i))
ans+=i;
printf("%lld\n",ans);
return 0;
}