Динамический массив, много циклов и простые числа. Как ускорить работу программы ?

#include <iostream>
using namespace std;
 
int main()
{
    int n, k; 
    double q=0, h=0; //q-счетчик не нулевых элементов, h-счетчик не нулевых элементов с шагом k
    int counter=0; //побочный счетчик для нахождения простых чисел
    cout<<"Enter n=";
    cin>>n;
    cout<<"Enter k=";
    cin>>k;
    int *arr=new int[n]; //создаем динамический массив с размером n
    for(int i=0; i<n; i++) //цикл присвоения элементам массива чисел от 1 до n
        arr[i]=i+1;
    for(int i=0; i<n; i++) //цикл вывода массива на экран
cout<<arr[i]<<" ";
    cout<<endl;
for(int i=0; i<n; i++) //цикл перебора всех элементов массива
{
    for(int z=1; z!=arr[i]; z++) //цикл нахождения всех простых чисел от 2 до arr[i] (элемента массива с которым мы работаем)
    {
        counter=0;
        for(int j=1; j!=z+1; j++) //цикл проверки числа на простоту
        {
            if(z%j==0)
                counter++;
        }
        if(counter==2) //если число делится на единицу и на само себя (делилось всего 2 раза), значит оно простое
        {
            if(arr[i]%(z*z)==0) //если элемент массива делится на квадрат простого числа - выйти из цикла и проверять следующий элемент
            {
                arr[i]=0;
                break;
            }
        }
    }
}
for(int i=0; i<n; i++) //цикл подсчета не нулевых элементов
if(arr[i]!=0)
q++;
for(int i=0; i<n; i+=k) //цикл подсчета не нулевых элементов с шагом k
if(arr[i]!=0)
h++;
cout<<"--------"<<endl;
for(int i=0; i<n; i++) //цикл вывода массива на экран
cout<<arr[i]<<" ";
cout<<endl;
cout<<"h="<<h<<endl;
cout<<"q="<<q<<endl;
cout<<"h/q="<<h/q<<endl;
system("pause");
return 0;
}

#include <iostream>
#include <math.h>
using namespace std;
 
int main()
{
    int n, k; 
    double q=0, h=0; //q-счетчик не нулевых элементов, h-счетчик не нулевых элементов с шагом k
    int counter=0; //побочный счетчик для нахождения простых чисел
    cout<<"Enter n=";
    cin>>n;
    cout<<"Enter k=";
    cin>>k;
    int *arr=new int[n]; //создаем динамический массив с размером n
    for(int i=0; i<n; i++) //цикл присвоения элементам массива чисел от 1 до n
    {
        arr[i]=i+1;
        cout<<i+1<<" ";
    }
    cout<<endl;
for(int i=0; i<n; i++) //цикл перебора всех элементов массива
{
    for(int z=1; z!=arr[i]; z++) //цикл нахождения всех простых чисел от 2 до arr[i] (элемента массива с которым мы работаем)
    {
        counter=0;
        for(int j=1; j<sqrt((float)z); j++) //цикл проверки числа на простоту
        {
            if(z%j==0)
                counter++;
        }
        if(counter==1) //если число делится на единицу и на само себя (делилось всего 2 раза), значит оно простое
        {
            if(arr[i]%(z*z)==0) //если элемент массива делится на квадрат простого числа - выйти из цикла и проверять следующий элемент
            {
                arr[i]=0;
                break;
            }
        }
    }
}
for(int i=0; i<n; i++) //цикл подсчета не нулевых элементов
if(arr[i]!=0)
q++;
for(int i=0; i<n; i+=k) //цикл подсчета не нулевых элементов с шагом k
if(arr[i]!=0)
h++;
cout<<"--------"<<endl;
for(int i=0; i<n; i++) //цикл вывода массива на экран
cout<<arr[i]<<" ";
cout<<endl;
cout<<"h="<<h<<endl;
cout<<"q="<<q<<endl;
cout<<"h/q="<<h/q<<endl;
delete[] arr;
system("pause");
return 0;
}
}

#include <iostream>
#include <math.h>
#include <ctime>
using namespace std;
int main()
{
     int limit = 10001;
     int sqr_lim;
     bool is_prime[10002];
     int x2, y2;
     int i, j;
     int n;
     sqr_lim = (int)sqrt((long double)limit);
     for (i = 0; i <= limit; i++) 
         is_prime[i] = false; // Инициализация решета
     is_prime[2] = true;
     is_prime[3] = true;
// Предположительно простые - это целые с нечетным числом
// представлений в данных квадратных формах.
// x2 и y2 - это квадраты i и j (оптимизация).
     x2 = 0;
     for (i = 1; i <= sqr_lim; i++) 
     {
         x2 += 2 * i - 1;
         y2 = 0;
         for (j = 1; j <= sqr_lim; j++) 
         {
             y2 += 2 * j - 1;
             n = 4 * x2 + y2;
             if ((n <= limit) && (n % 12 == 1 || n % 12 == 5))
                 is_prime[n] = !is_prime[n];
             // n = 3 * x2 + y2; 
             n -= x2; // Оптимизация
             if ((n <= limit) && (n % 12 == 7))
                 is_prime[n] = !is_prime[n];
             // n = 3 * x2 - y2;
             n -= 2 * y2; // Оптимизация
             if ((i > j) && (n <= limit) && (n % 12 == 11))
                 is_prime[n] = !is_prime[n];
         }
     }
// Отсеиваем кратные квадратам простых чисел в интервале [5, sqrt(limit)].
// (основной этап не может их отсеять)
     for (i = 5; i <= sqr_lim; i++) 
     {
         if (is_prime[i]) 
         {
             n = i * i;
             for (j = n; j <= limit; j += n) 
             {
                 is_prime[j] = false;
             }
         }
     }
// Вывод списка простых чисел в консоль.
     printf("2 3 5"); 
     for (i = 6; i <= limit; i++) // проверка делимости на 3 и 5
     { 
         if ((is_prime[i]) && (i % 3 != 0) && (i % 5 != 0))
         {
             printf(" %d", i);
         }
     } 
     cin>>n;
}

#include <iostream>
#include <cmath>
using namespace std;
 
int main()
{
    int n, k; 
    double q=0, h=0; //q-счетчик не нулевых элементов, h-счетчик не нулевых элементов с шагом k
    int counter=0; //побочный счетчик для нахождения простых чисел
    cout<<"Enter n=";
    cin>>n;
    cout<<"Enter k=";
    cin>>k;
    int *arr=new int[n]; //создаем динамический массив с размером n
    for(int i=0; i<n; i++) //цикл присвоения элементам массива чисел от 1 до n
        arr[i]=i+1;
    for(int i=0; i<n; i++) //цикл вывода массива на экран
        cout<<arr[i]<<" ";
    cout<<endl;
for(int i=0; i<n; i++) //цикл перебора всех элементов массива
{
    for(int z=1; z!=arr[i]; z++) //цикл нахождения всех простых чисел от 2 до arr[i] (элемента массива с которым мы работаем)
    {
        counter=0;
        int kor=int(sqrt(float(z+1)));
        for(int j=1; j!=kor; j++) //цикл проверки числа на простоту
        {
            if(z%j==0)
            {
                counter++;
                if(counter > 2)
                break;
            }
        }
        if(counter==2) //если число делится на единицу и на само себя (делилось всего 2 раза), значит оно простое
        {
            if(arr[i]%(z*z)==0) //если элемент массива делится на квадрат простого числа - выйти из цикла и проверять следующий элемент
            {
                arr[i]=0;
                break;
            }
        }
    }
}
cout<<"--------"<<endl;
for(int i=0,j=0; i<n; i++,j+=k) //цикл подсчета не нулевых элементов
{
    cout<<arr[i]<<" ";
    if(arr[i]!=0)
        q++;
    if(j<n && arr[i]!=0)
        h++;
}
cout<<endl;
cout<<"h="<<h<<endl;
cout<<"q="<<q<<endl;
cout<<"h/q="<<h/q<<endl;
system("pause");
return 0;
}