題目大意:
給定一個H-number,問有多少H-semi-primes:
定義一種數叫H-numbers,它是所有能除以四餘一的數。
在H-numbers中分三種數:
1、H-primes,這種數只能被1和它本身整除,不能被其他的H-number整除,例如9是一個H-number,能被1,3,9整除,但3不是H-number,所以他是H-primes。
2、H-semi-primes是由兩個H-primes相乘得出的。
3、剩下的是H-composite。
問給一個數,求1到這個數之間有多少個H-semi-primes。
解題思路:
像打素數表一樣打出這個semi表;
AC代碼;
#include <map>
#include <set>
#include <stack>
#include <cmath>
#include <queue>
#include <bitset>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <fstream>
#include <cstdlib>
#include <sstream>
#include <cstring>
#include <iostream>
#include <algorithm>
#pragma comment(linker, "/STACK:1024000000,1024000000")
using namespace std;
#define maxn 1010
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define ms(x,y) memset(x,y,sizeof(x))
#define rep(i,n) for(int i=0;i<(n);i++)
#define repf(i,a,b) for(int i=(a);i<=(b);i++)
#define PI pair<int,int>
//#define mp make_pair
#define FI first
#define SE second
#define IT iterator
#define PB push_back
#define Times 10
typedef long long ll;
typedef unsigned long long ull;
typedef long double ld;
typedef pair<int ,int > P;
//#define N 100
const double eps = 1e-10;
const double pi = acos(-1.0);
const ll mod = 1e9+7;
const int inf = 0x3f3f3f3f;
const ll INF = (ll)1e18+300;
const int maxd = 1000000 + 10;
int vis[maxd];
int pri[maxd];
void init() {
int n = maxd;
ms(vis, 0);
ms(pri, 0);
for (int i = 5; i < n; i += 4) {
for (int j = 5; j < n; j += 4) {
if(i * j > n) {
break;
}
if(!vis[j] && !vis[i]) {
vis[i*j] = 1;
}
else{
vis[i*j] = -1;
}
}
}
int num = 0;
for (int i = 0; i < n; i++) {
if(vis[i] == 1) {
num ++;
}
pri[i] = num;
}
}
int main() {
int n;
init();
while(cin >> n && n) {
cout << n << " " << pri[n] << endl;
}
}