开始之前先介绍两种求素数的方法,埃氏筛和线性筛
一 . 埃式筛法主要用于筛选素数,合数不用筛,思想就是如果一个数是质数,那么它的倍数肯定非质。先把所有整数列出来,然后把2的倍数全部剔除,然后是三的,以此类推,遍历所有素数,把倍数全部划去。
二 . 线性筛法,和埃氏筛法的区别是对于每一个要筛除的数,欧拉筛法只筛除一次,而埃氏筛法会重复筛除,比如8和16同时被2和4筛去,推荐使用欧拉筛法,当数据小于10^6时两者都可以
题目:
DD 在数学课上学习了素数,在信息课上学习了回文数,她现在很想知道,对于区间 [L,R],其中有多少个数既是素数又是回文数
输入格式
第一行两个整数分别表示 L,R
输出格式
一个整数表示有多少个数满足条件
数据范围
对于 20%的数据,L=R
对于 50% 的数据,L≤R≤10^3
对于 100% 的数据, L≤R≤10^6
打表(法一):
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
const int N = 1e6 + 10;
int l, r, n, m;
int cnt, res;
int primes[N];
bool vis[N];
int arr[N], pos;
void get_primes() //打表,用线性筛法把1--10^6的素数全部算出来储存在表里
{
for (int i = 2; i < N; i++) //i<N,N这时已知,打表
{
if (vis[i])
continue;
primes[cnt++] = i;
for (int j = i + i; j < N; j += i)
vis[j] = true;
}
}
bool is_palindrome(int x)
//判断回文数,现在写的是int类型下如何判断回文,下面会介绍string类型下回文的判断
{
int t = x;
int y = 0;
while (t)
{
y = y * 10 + t % 10;
t /= 10;
}
if (y == x)
return true;
else
return false;
}
int main()
{
get_primes(); //引入素数的打表
for (int i = 0; i < cnt; i++)
if (is_palindrome(primes[i]))
arr[pos++] = primes[i];//将所有的素数回文数都找出来存入数组
scanf("%d%d", &n, &m);
int l = lower_bound(arr, arr + pos, n) - arr;
//lower_bound()的参数是lower_bound(begin,end,value),作用是从begin位置到end-1位置二分查找第一个大于或等于value的数字,找到返回该数字的地址,不存在则返回end。
for (int i = l; i < pos; i++)
{
if (arr[i] > m)
break;
res++; //寻找个数
}
cout << res << endl;
return 0;
}
打表(法二):
#include<bits/stdc++.h>
using namespace std;
const int N=1e+6;
bool isprime[N];
int primes[N],arr[N];
int cnt,res,pos;
void is_prime() //打表,埃氏筛法
{
for(int i=0; i<=N; i++)
isprime[i]=true;
isprime[0]=isprime[1]=false;
for(int i=2; i<=N; i++)
{
if(isprime[i])
{
for(int j=i*2; j<=N; j+=i)
{
isprime[j]=false;
}
}
if(isprime[i]) primes[cnt++]=i;
}
}
bool isHuiwen(int x) //字符串式判断
{
int i,length,flag=1;
string s;
s=to_string(x); //将数字转为字符串
int k=s.length();
for(i=0;i<=k/2;i++) {
if(s[i]!=s[k-i-1]) {
flag=0;
break;
}
}
if(flag==1) return 1;
else return 0;
}
int main()
{
is_prime();
for (int i = 0; i < cnt; i++)
if (isHuiwen(primes[i]))
arr[pos++] = primes[i];//将所有的素数回文数都找出来存入数组
scanf("%d%d", &n, &m);
int l = lower_bound(arr, arr + pos, n) - arr;
for (int i = l; i < pos; i++)
{
if (arr[i] > m)
break;
res++; //寻找个数
}
cout << res << endl;
return 0;
}
还有一种非打表方法,不过主要内容基本一样,不再注释
#include<bits/stdc++.h>
using namespace std;
bool isPrime(int n)
{
if(n==1) return 0;
if(n==2||n==3) return 1;
if(n%6!=1&&n%6!=5) return 0;
for(int i=5; i*i<=n; i+=6)
{
if(n%i==0||n%(i+2)==0) return 0;
}
return 1;
}
bool isHuiwen(int x)
{
int t=x;
int y=0;
while(t)
{
y=y*10+t%10;
t=t/10;
}
if(y==x) return 1;
else return 0;
}
int main()
{
int n,m,res=0;
cin>>n>>m;
for(int i=n; i<=m; i++)
{
if(isHuiwen(i))
{
if(isPrime(i))
res++;
}
}
cout<<res;
return 0;
}