这道题第一遍写的写的简单暴力,不给过。后面了解到要用埃氏筛算法。
在这我先介绍一下这种算法的思想,其实埃拉托斯特尼筛法是一种很古老的算法,它的原理很简单,先标记所有数是质数,从二开始,把所有二的倍数重新标记成合数,知道n为止,每次都能保证我正在选择的质数没有比它本身和1以外其他的因子,也可以保证所有非质数都被排除。如果我们要求n以内的范围只需遍历到根号n
它的复杂度是O(n+log(log n))是一种很好用的算法
具体实现如下
int f[1000000]
cin>>n;
for(int i=2;i*i<=n;i++)
if(f[i]==0)for(int j=2*i;j<=n;j+=i)
f[j]=1;
到了这道题我们需要灵活改变一下思路。对于每一个数都用埃氏筛处理,如果是质数就正常选出所有质数,如果是合数就把它的所有倍数+1表示他有一个合数因子,最后在桶排序统计每一种和因子数的多少即可代码可以这样表示
#include<bits/c++std.h>
using namespace std;
typedef long long ll;
const ll mod=1e9+7;
bool prime[100005];
int f[100005],cnt[100005];
void judge(int n)
{
for(int i=2;i<=n;i++){
prime[i]=true;
}
for(int i=2;i<=n;i++){
if(prime[i]){
for(int j=i+i;j<=n;j+=i){
prime[j]=false;
}
}
else{
for(int j=i;j<=n;j+=i){
cnt[j]++;
}
}f[cnt[i]]++;
}
}
int main()
{
int n,m;
cin>>n>>m;
judge(n);
while(m--){
int k;
cin>>k;
cout<<f[k]<<endl;
}
return 0;
}
很简单的想法,所有合数的倍数都有这样一个合因子
后面无聊看其他人的代码,我发现了一种更简短的代码,它的原理和这个是一样的但实现时用了一种很巧妙的技巧
代码如下
#include <iostream>
using namespace std;
int a[100010],b[100010];
int main(){
int n,m,k;
cin>>n>>m;
for(int i=2;i<=n;i++)
for(int j=i+i;j<=n;j+=i)
if(a[j]==0 || a[i]!=0) a[j]++;
for(int i=0;i<=n;i++) b[a[i]]++;
while(m--){
cin>>k;
cout<<b[k]<<endl;
}
return 0;
}
其实这个代码主要有两个地方和上面不同,首先数组a[]的意义是当前这个数所有的合数个数,如果这个数是1那就代表他本身,这段代码的判断语句其实就是,如果这个数的倍数还是0,就将其标记为1,因为它已经是i的倍数了,它本身就是他自己的合因子,相当于去除质数的作用如果它本身已近被加过了就不用重复标记了,后半句意思就是如果这个因子本身就是合数,那么他的所有倍数的合因子数+1