XTU oj 1279 Dual Prime

Dual Prime

题目描述

如果一个合数x=p⋅q,p,q是素数且p≠q，我们称x是双素数。 现给你一个区间[a,b]，求区间内的的双素数个数。

输入

第一行是一个整数T(1≤T≤30000)，为样例的数目。 以后每行一个样例，为两个整数a,b(1≤a≤b≤106)

输出

依次每行输出一个样例的结果。

样例输入

3
1 10
1 100
1 1000000

样例输出

2
30
209867

        思路：和之前的各种素数是一样的，只是这题关键在于统计区间内Dual Prime个数。这题时间限制很松（老师说的），所以可以先判断，再枚举计数。

        数组要是开到1000000而且记录Dual Prime的数组不先清零的话，最后一组样例会多一个，而且oj报错了，估计p*q没算到1000000。开到1000500也行，而且不用初始化这个数组了。

#include <stdio.h>
#include <stdbool.h>//包含布尔变量的头文件
#include <memory.h>//包含memset的头文件，也可以用stdlib.h
#define MAX 1000000
bool isprime[MAX];
int prime[MAX]={0};
int num[MAX]={0};//记录Dual Prime
int ans[MAX]={0};//记下Dual Prime的个数
void eular();
void count();
int main (){
	int T;
	int a,b;
	eular ();
	count ();
	scanf ("%d",&T);
	while (T--){
		scanf ("%d %d",&a,&b);
		printf ("%d\n",ans[b]-ans[a-1]);//区间内Dual Prime	
	}
	return 0;
} 

void eular(){//欧拉筛
	int i,j;
	memset (isprime,true,sizeof(isprime));
	isprime[1]=false;
	for (i=2;i<MAX;i++){
		if (isprime[i]) prime[++prime[0]]=i;
		for (j=1;j<=prime[0]&&i*prime[j]<=MAX;j++){
			isprime[i*prime[j]]=false;
			if (i%prime[j]==0) break;
		}
	}
}

void count(){
	int p,q,cnt=0,i;
	for (p=2;p<=MAX;p++){//判断Dual Prime
		if (isprime[p]){//别把if写到for里面，否则碰到第一个不满足条件的就会跳出，导致漏数
			for (q=2;p*q<=MAX;q++){
				if (isprime[q]&&p!=q) num[p*q]=1;
			}
		}
	}
	for (i=6;i<=MAX;i++){//计数
		if (num[i]) ans[i]=ans[i-1]+1;
		else ans[i]=ans[i-1];
	}
	
}