题目链接:
http://acm.nyist.net/JudgeOnline/problem.php?pid=26
描述
写一个程序,找出给出素数范围内的所有孪生素数的组数。一般来说,孪生素数就是指两个素数距离为2,近的不能再近的相邻素数。有些童鞋一看到题就开始写程序,不仔细看题,咱们为了遏制一下读题不认真仔细的童鞋,规定,两个素数相邻为1的也成为孪生素数。
输入
第一行给出N(0<N<100)表示测试数据组数。
接下来组测试数据给出m,表示找出m之前的所有孪生素数。
(0<m<1000000)
输出
每组测试数据输出占一行,该行为m范围内所有孪生素数组数。
样例输入
1 14 |
样例输出
4 |
算法思想:
首先,生成所有的素数,并将其标识。然后从m开始循环遍历,判断a[i]与a[i - 2]是否均为素数,如果是,则ans++,当m>=3时,ans的初始值为1,因为3和2也是一组孪生素数。
源代码
#include <iostream>
#include <cmath>
using namespace std;
#define MAX_P 1000000
int a[MAX_P];
//生成素数算法
void getPrimes()
{
int i = 0, j = 0, n = sqrt(MAX_P) + 1;
for (i = 2; i <= n; i++) //筛选循环
for (j = 2; j <= MAX_P / i; j++)
a[j*i] = 1;
}
int main()
{
int N, m, ans;
getPrimes();
cin >> N;
while (N--)
{
ans = 0;
cin >> m;
if (m >= 3)
{
ans = 1;
}
for (int i = m; i >=4; i--)
{
if (!a[i - 2] && !a[i])
ans++;
}
cout << ans << endl;
}
cout << endl;
return 0;
}
最优源代码
#include <iostream>
#include<math.h>
using namespace std;
bool vis[1000010];
int main()
{
int n=1000010;
int m =sqrt(n+0.5);
int c=0;
for(int i =2;i<=m;i++)
if(!vis[i])
{
for (int j = i*i;j<=n;j+=i)
vis[j]=1;
}
cin>>n;
while(n--)
{
int count=0,m;
cin>>m;
for(int i=3;i<m-1;i++)
{
if(!vis[i] && !vis[i+2]) count++;
}
if(m>3)
cout<<count+1<<endl;
else if(m==3) cout<<"1"<<endl;
else cout<<"0"<<endl;
}
return 0;
}
算法复杂度:
由源代码可知,算法时间复杂度在于求出所有的素数,两层循环,故算法时间复杂度为O(sqrt(n) * sqrt(n)) = O(n)。