题目链接:
http://acm.nbu.edu.cn/v1.0/Problems/Problem.php?pid=2296
题目大意:求1~m(1<=m<=10000)中有几个数可以表示成p*q的形式,且p,q均为质数(p可以等于q).
题目思路:
第一种方法:纯暴力(时间4000MS).
第二种方法:先把筛出质数,然后利用 dp[i]=dp[i-1]+{所有p*q>i-1的个数} ,O(10000)的复杂度预处理出所有答案.
题目大意:求1~m(1<=m<=10000)中有几个数可以表示成p*q的形式,且p,q均为质数(p可以等于q).
题目思路:
第一种方法:纯暴力(时间4000MS).
第二种方法:先把筛出质数,然后利用 dp[i]=dp[i-1]+{所有p*q>i-1的个数} ,O(10000)的复杂度预处理出所有答案.
代码:
#include <stdlib.h>
#include <string.h>
#include <stdio.h>
#include <ctype.h>
#include <math.h>
#include <time.h>
#include <stack>
#include <queue>
#include <map>
#include <set>
#include <vector>
#include <string>
#include <iostream>
#include <algorithm>
using namespace std;
#define ll __int64
//#define ll long long
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define middle l+r>>1
#define INF 0x3F3F3F3F
#define esp (1e-8)
#define MOD 100000007
#define type int
//const double pi=acos(-1.0);
const int M=10000 +5;
type min(type x,type y){return x<y? x:y;}
type max(type x,type y){return x>y? x:y;}
void swap(type& x,type& y){type t=x;x=y;y=t;}
int T,cas=0;
int m,f[M],ans[M];
void select(){
memset(f,0,sizeof(f));
int MM=sqrt(1.0*M),i,j,k;
for(i=3;i<=MM;i+=2) if(!f[i])
for(j=i;i*j<=M;j+=2)
f[i*j]=1;
for(i=4;i<=M;i+=2) f[i]=1;
ans[1]=ans[2]=ans[3]=0;
for(i=4;i<=M;i++){
ans[i]=ans[i-1];
for(j=2;j<=(i>>1);j++)if(!f[j]){
for(k=i/j;k>=j;k--)if(!f[k]){
if(j*k<=i-1) break;
ans[i]++;
}
}
}
};
void run(){
scanf("%d",&m);
printf("%d\n",ans[m]);
}
int main(){
//freopen("1.in","r",stdin);
//freopen("1.out","w",stdout);
//run();
select();
for(scanf("%d",&T),cas=1;cas<=T;cas++) run();
//while(~scanf("%d",&n) && n) run();
//system("pause");
return 0;
}