题目意思是找到2 ~n-1中最大公因数是d的个数
这道题可以转化一下。
转化的思路跟欧拉的素数筛法的原理有异曲同工之妙,so先看一下欧拉筛法的原理(先上代码)
for(long i = 2;i<=MAXN;i++){
if(!is_prime[i])//如果是质数,值为0
{
prime[n++] = i;//存放素数
}
for(int j = 0;j<n&&i*prime[j]<MAXN;j++){
is_prime[i*prime[j]] = 1;
if(!i%prime[j])
break;
}
欧拉筛法区别与埃式筛法的原理是它避免了重复计算,使得每个数只被计算了一次,时间复杂度为O(n)。
那么它是如何实现的呢?
首先如果一个合数能够写成X=ab;其中若a为合数 ,b为素数,那将a拆分成a1,a2,那么X就可以表示成X=a1(a2b)。这个筛法的关键在于
if(!i%prime[j])
break;
它是从小到大搜索,把一个数的素因子升序排列(a1a2a3…an),在扫描过程中,如果prime等于它的最小素数,继续扫描将会重复计算。
关于题意的转化,假如X= d * a,按照题目的意思,d为最大公因数,那么,a必为素数;反过来a为素数d不一定就是X的最大公因数,如果d是合数(a1,a2),a1<a时,就不成立。所以计算有多少个数是,应该统计的是有多少个
小于d的最小公因数。
if(!i%prime[j])
break;
上完整代码
#include<iostream>
#include<cstdio>
using namespace std;
const int MAXN = 1e+5;
int is_prime[MAXN+5];
int prime[MAXN+5];
long num = 0;
void _prime(){
for(long i = 2;i<MAXN;i++){
if(!is_prime[i])
{
prime[num++] = i;
}
for(int j = 0;j<num&&i*prime[j]<MAXN;j++){
is_prime[i*prime[j]] = 1;
if(!i%prime[j])
break;
}
}
}
int main(){
_prime();
long long N,n,d;
scanf("%d",&N);
for(int i = 0;i<N;i++){
scanf("%d%d",&n,&d);
int count = 0;
for(int j = 0;j<num;j++){
if(prime[j]>d ||d*prime[j]>=n)
break;
else if(d%prime[j]==0){
count++;
break;
}
count++;
}
printf("%d\n",count);
}
return 0;
}
第一次提交TLE,后来把cin,cout改成scanf,printf,就过了。‘