素数个数

今天这都什么破题啊……

首先是统计素数个数的破题

这是zrt讲的线性筛法

//O(n)
// mindiv 一个数字的最小素因子
// prime  
for(int i=2;i<=n;i++){
    if(!mindiv[i]){    //没有就存为自己 
        mindiv[i]=i;
        prime[tot++]=i;    //判断为素数 
        //printf("%d ",i);
    }
    for(int j=0;prime[j]*i<=n;j++){    //去掉不是素数的
        mindiv[prime[j]*i]=prime[j]; // prime[j]<=mindiv[i]
        //printf("# %d=%d*%d
",i*prime[j],i,prime[j]);
        if(prime[j]==mindiv[i]) break;    //后面的就不用了，免得重复 
    }
}
// N = p1^a1 p2^a2 p3^a3 ....
// 1) a1 == 1 i= p2^a2 p3^a3 .... prime[j] = p1
// 2) a1 >1  i=p1^(a1-1) p2^a2 p3^a3 .... prime[j] = p1
// O(n)

听有人说线性筛法好像不行（

反正我不知道（

这是正解给出的代码

#include <cstdio> 
#include <cstring> 
#include <iostream>
using namespace std;  
  
bool isPrime[60000001];    //是不是素数 
int primeCount, primes[5761455];    //存储素数 
   
int main() {
   freopen("prime.in","r",stdin);
   freopen("prime.out","w",stdout); 
    int n; 
    scanf("%d",&n);
    memset(isPrime, true, sizeof(isPrime));    //其实不这样也能做，不过有点别扭（ 
    primeCount = 0; 
    for (int i = 2; i <= n; ++ i) { 
        if (isPrime[i]) { 
            primes[primeCount ++] = i; 
        } 
        for (int j = 0; j < primeCount && i * primes[j] <= n; ++ j) { 
            isPrime[i * primes[j]] = false; 
            if (i % primes[j] == 0) {    //同线性筛法，防止重复 
                break; 
            } 
        } 
    } 
    printf("%d",primeCount); 
    return 0; 
}

下面这个是……这是我写的暴力……

//我为什么要把我的破代码贴上来？

顺便说一句这个是老师让我们画的奇怪的图，我不想画，所以

mmp的我改题都懒得改，原题都懒得放