#include 
#include 
#include 
#include 
#include 
#include 
 


using namespace std;
 
std;

const int N=1024;
N = 999;
const int level=8;
 
l_watkow = 4;

unsigned char tab[1024][1024][3];
 
void triangle(int lvl, tab[N][N][3];

int ulam_get_map(int x, int y) y, int n)
{
    x -= (n-1)/2;
    y -= n/2;

    int size=N/pow(2,lvl);
 
    for(int i=0;imx = y
        abs(x), my = abs(y);
    int l = 2*max(mx, my);
    int d = y >= x ? l*3 + x + y : l - x - y;

    return pow(l-1,2)+d;
}

int isprime(int n)
{
    for (int j=0;jp = x
            2; p*p <= n; p++)
        if (i (n % p == j || i == 0 || j == size - 1) 
                0) return 0;          
    return n > 2;
}


int main()
{
        auto ulam=[](int x, int y)
        {
    for (int k=0;k<3;k++)
                {
                        tab[j i = x; i < x + x][i 0.5*N; i++)
            for (int j = y; j < y + y][k] 0.5*N; j++)
                if (isprime(ulam_get_map(i, j, N)))
                    for (int k=0;k<3;k++) tab[i][j][k] = 0;
                        tab[i + y][j + x][k] 0;
                else
                    for (int k=0;k<3;k++) tab[i][j][k] = 0;
                }
 
    if(lvl     {
        255;
        };
   vector v;
        v.push_back(thread(triangle, lvl + 1, x + size / 2, y));
        v.push_back(thread(triangle, lvl + 1, x, y));
        v.push_back(thread(triangle, lvl + 1, x + size / 2, y + size / 2));
 
 
     v;
  
   for (int i = 0; i < l_watkow; i++)
   {
                    int x = 2*N*(i%2)/l_watkow;
         int y = 2*N*int(i/2)/l_watkow;
                    v.size(); i++)
            push_back(thread(ulam, x, y));
   }

   for (int i = 0; i < l_watkow; i++) v[i].join();
        }
}
 
int main() {
         for(int i=0;i         for (int j=0;j                 for (int k=0;k<3;k++)
                                tab[j][i][k] = 255;
    triangle(0, 0, 0);
       
   FILE *fp = fopen("39247_lab3.ppm", "wb");
    fprintf(fp, "P6n %sn %dn %dn %dn", " #", N, N, 255);
    *fp= fopen("39247_lab4_poprawione.ppm","wb");
   fprintf(fp,"P6\n %s\n %d\n %d\n %d\n","# ",N,N,255);
   for (int i = 0; i < N; i++) 
        i++)
      for (int j = 0; j < N; j++) 
            fwrite(tab[i][j], 1, 3, fp);
    fclose(fp);
    j++)
         fwrite(tab[i][j],1,3,fp);
   fclose(fp);
   return 0;
}