题目大意:
给定4n+1数(1、5、9、13、……)。将这些数分为unit(即为1)和prime(不是真正的素数),composite。规定一个semi-prime数为可为两个prime数乘积。题目给定一个4n+1数,要判断1~4n+1数之间(包含1和该4n+1数)的所有semi-prime个数。
(类似于筛选素数)
附AC代码:
#include <iostream> #include <algorithm> #include <cstdio> #include <cstring> #include <vector> #define maxn 1000010 typedef long long int ll; using namespace std; vector<int> v;//存放4n+1的数 bool isHP[maxn]; int ans[maxn]; void init() { memset(isHP,true,sizeof(isHP)); for(int i=5; i<maxn; i+=4) { if(isHP[i]) { v.push_back(i); for(int j=i+i;j<maxn;j+=i) { if(j%4==1) isHP[j]=false; } } } for(int i=0;i<v.size();i++) { if(v[i]>1010)//只需v[i]小于maxn的开方 break; for(int j=i;j<v.size();j++) { int temp=v[i]*v[j];//temp必定是4n+1形式(易证),无需判断 if(temp>=maxn) break; ans[temp]=1; } } for(int i=1;i<maxn;i++) { ans[i]+=ans[i-1]; } } int main() { init(); int n; while(scanf("%d",&n)&&n) { printf("%d %d\n",n,ans[n]); } return 0; }