牛客寒假3 f k合因字数（埃氏筛的扩展应用）

这道题第一遍写的写的简单暴力，不给过。后面了解到要用埃氏筛算法。
在这我先介绍一下这种算法的思想，其实埃拉托斯特尼筛法是一种很古老的算法，它的原理很简单，先标记所有数是质数，从二开始，把所有二的倍数重新标记成合数，知道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