oj.uz

Submission #962301

# Submission time Handle Problem Language Result Execution time Memory
962301 2024-04-13T10:17:00 Z danikoynov Art Class (IOI13_artclass) C++14
100 / 100
 64 ms 9704 KB 
#include "artclass.h"
#include<bits/stdc++.h>
using namespace std;

struct pixel
{
    double r, g, b;
    pixel(double _r = 0.0, double _g = 0.0, double _b = 0.0)
    {
         r = _r;
         g = _g;
         b = _b;
    }
};

pixel avg_st[5];
pixel image[510][510];
int h, w;

pixel avg_dev[5];

void init()
{
     avg_st[1] = pixel(189.852, 166.758, 147.354);
     avg_st[2] = pixel(102.854, 98.1113, 64.2302);
    // avg_st[3] = pixel(140.586, 128.459, 107.031);
    // avg_st[4] = pixel(158.108, 89.7725, 55.3149);

    avg_dev[1] = pixel(420.707, 429.94, 475.923);
    avg_dev[2] = pixel(255.938, 241.82, 189.427);
    avg_dev[3] = pixel(1563.29, 1532.6, 1435.09);
    avg_dev[4] = pixel(25.4556, 23.6549, 35.4133);
}

double penalty[5];

pixel average_pixel(int stx, int sty, int enx, int eny)
{
     double sr = 0.0;
     double sg = 0.0;
     double sb = 0.0;

     for (int i = stx; i <= enx; i ++)
          for (int j = sty; j <= eny; j ++)
     {
          sr += image[i][j].r;
          sg += image[i][j].g;
          sb += image[i][j].b;
     }

     double cnt = (enx - stx + 1) * (eny - sty + 1);
     sr /= cnt;
     sg /= cnt;
     sb /= cnt;
     return pixel(sr, sg, sb);
}

pixel region_deviation(int stx, int sty, int enx, int eny)
{
     pixel avg = average_pixel(stx, sty, enx, eny);

     double dr = 0.0;
     double dg = 0.0;
     double db = 0.0;

     for (int i = stx; i <= enx; i ++)
     {
          for (int j = sty; j <= eny; j ++)
          {
               double delta_red = abs(image[i][j].r - avg.r);
               double delta_green = abs(image[i][j].g - avg.g);
               double delta_blue = abs(image[i][j].b - avg.b);

               delta_red = delta_red * delta_red;
               delta_green = delta_green * delta_green;
               delta_blue = delta_blue * delta_blue;

               dr += delta_red;
               dg += delta_green;
               db += delta_blue;
          }
     }

     double cnt = (enx - stx + 1) * (eny - sty + 1);

     dr /= cnt;
     dg /= cnt;
     db /= cnt;

     return pixel(dr, dg, db);
}

const int block = 4;
pixel get_average_block_deviation()
{
     pixel avg;
     double cnt = 0.0;
         for (int i = 0; i < h - block + 1; i ++)
     for (int j = 0; j < w - block + 1; j ++)
     {
          int stx = i, sty = j;
          int enx = min(h, i + block) - 1;
          int eny = min(w, j + block) - 1;


          pixel cur = region_deviation(stx, sty, enx, eny);
          avg.r += cur.r;
          avg.g += cur.g;
          avg.b += cur.b;
          cnt ++;
     }

     avg.r /= cnt;
     avg.g /= cnt;
     avg.b /= cnt;
     return avg;
     //cout << "average deviation" << endl;
     //cout << avg.r << " " << avg.g << " " << avg.b << endl;
}

pixel get_average_pixel()
{
     pixel s;
     for (int i = 0; i < h; i ++)
          for (int j = 0; j < w; j ++)
     {
          s.r += image[i][j].r;
          s.g += image[i][j].g;
          s.b += image[i][j].b;
     }

     double cnt = (double)(h * w);
     s.r /= cnt;
     s.g /= cnt;
     s.b /= cnt;
     return s;
}


int further_check()
{
     pixel avg = get_average_pixel();

     penalty[1] = penalty[2] = 0;
     for (int i = 1; i <= 2; i ++)
     {
          double dr = abs(avg.r - avg_st[i].r);
          double dg = abs(avg.g - avg_st[i].g);
          double db = abs(avg.b - avg_st[i].b);

          dr = dr * dr;
          dg = dg * dg;
          db = db * db;
          penalty[i] = dr + dg + db;
     }

     if (penalty[1] < penalty[2])
          return 1;
     return 2;
}
int style(int H, int W, int R[500][500], int G[500][500], int B[500][500]) {
    init();
    h = H;
    w = W;
    for (int i = 0; i < h; i ++)
     for (int j = 0; j < w; j ++)
    {
         image[i][j] = pixel(R[i][j], G[i][j], B[i][j]);
    }

    pixel dev = get_average_block_deviation();
    for (int s = 1; s <= 4; s ++)
    {
         double dr = abs(dev.r - avg_dev[s].r);
         double dg = abs(dev.g - avg_dev[s].g);
         double db = abs(dev.b - avg_dev[s].b);

         dr = dr * dr;
         dg = dg * dg;
         db = db * db;
         penalty[s] += dr + dg + db;
    }
    int mx = 1;
    for (int i = 1; i <= 4; i ++)
    {
         if (penalty[i] < penalty[mx])
          mx = i;
    }

    if (mx == 1 || mx == 2)
     return further_check();
    return mx;
}

Subtask #1
100.0 / 100.0

# Verdict Execution time Memory Grader output
1 Incorrect 55 ms 9456 KB Output isn't correct
2 Correct 49 ms 9308 KB Output is correct
3 Correct 15 ms 9308 KB Output is correct
4 Correct 41 ms 9304 KB Output is correct
5 Correct 46 ms 9304 KB Output is correct
6 Correct 24 ms 9476 KB Output is correct
7 Correct 47 ms 9460 KB Output is correct
8 Correct 44 ms 9304 KB Output is correct
9 Correct 12 ms 9308 KB Output is correct
10 Correct 42 ms 9320 KB Output is correct
11 Correct 41 ms 9456 KB Output is correct
12 Correct 38 ms 9304 KB Output is correct
13 Incorrect 59 ms 9308 KB Output isn't correct
14 Correct 45 ms 9308 KB Output is correct
15 Correct 38 ms 9308 KB Output is correct
16 Correct 45 ms 9308 KB Output is correct
17 Correct 61 ms 9456 KB Output is correct
18 Correct 29 ms 9304 KB Output is correct
19 Correct 42 ms 9704 KB Output is correct
20 Correct 43 ms 9308 KB Output is correct
21 Correct 43 ms 9452 KB Output is correct
22 Correct 56 ms 9308 KB Output is correct
23 Correct 47 ms 9308 KB Output is correct
24 Correct 44 ms 9304 KB Output is correct
25 Correct 53 ms 9324 KB Output is correct
26 Correct 39 ms 9460 KB Output is correct
27 Correct 55 ms 9452 KB Output is correct
28 Correct 32 ms 9480 KB Output is correct
29 Correct 49 ms 9484 KB Output is correct
30 Correct 34 ms 9372 KB Output is correct
31 Correct 33 ms 9308 KB Output is correct
32 Correct 48 ms 9304 KB Output is correct
33 Correct 51 ms 9420 KB Output is correct
34 Correct 49 ms 9308 KB Output is correct
35 Correct 28 ms 9304 KB Output is correct
36 Correct 52 ms 9304 KB Output is correct
37 Correct 40 ms 9304 KB Output is correct
38 Correct 45 ms 9308 KB Output is correct
39 Correct 39 ms 9308 KB Output is correct
40 Correct 50 ms 9308 KB Output is correct
41 Correct 41 ms 9304 KB Output is correct
42 Correct 43 ms 9308 KB Output is correct
43 Correct 42 ms 9560 KB Output is correct
44 Correct 40 ms 9456 KB Output is correct
45 Correct 44 ms 9320 KB Output is correct
46 Correct 45 ms 9308 KB Output is correct
47 Correct 23 ms 9308 KB Output is correct
48 Correct 44 ms 9460 KB Output is correct
49 Correct 43 ms 9304 KB Output is correct
50 Correct 52 ms 9304 KB Output is correct
51 Correct 42 ms 9456 KB Output is correct
52 Correct 33 ms 9304 KB Output is correct
53 Correct 47 ms 9308 KB Output is correct
54 Correct 52 ms 9308 KB Output is correct
55 Correct 52 ms 9452 KB Output is correct
56 Correct 41 ms 9304 KB Output is correct
57 Correct 53 ms 9308 KB Output is correct
58 Correct 41 ms 9304 KB Output is correct
59 Incorrect 42 ms 9560 KB Output isn't correct
60 Correct 50 ms 9308 KB Output is correct
61 Correct 47 ms 9308 KB Output is correct
62 Correct 45 ms 9308 KB Output is correct
63 Correct 45 ms 9308 KB Output is correct
64 Correct 46 ms 9452 KB Output is correct
65 Correct 64 ms 9308 KB Output is correct
66 Correct 34 ms 9308 KB Output is correct
67 Correct 39 ms 9456 KB Output is correct
68 Correct 41 ms 9456 KB Output is correct
69 Correct 33 ms 9308 KB Output is correct
70 Correct 28 ms 9484 KB Output is correct
71 Correct 45 ms 9472 KB Output is correct
72 Correct 53 ms 9308 KB Output is correct
73 Correct 46 ms 9308 KB Output is correct
74 Correct 42 ms 9648 KB Output is correct
75 Correct 45 ms 9308 KB Output is correct
76 Correct 43 ms 9452 KB Output is correct
77 Correct 57 ms 9460 KB Output is correct
78 Correct 42 ms 9424 KB Output is correct
79 Incorrect 30 ms 9308 KB Output isn't correct
80 Correct 28 ms 9308 KB Output is correct
81 Correct 53 ms 9320 KB Output is correct
82 Correct 34 ms 9644 KB Output is correct
83 Correct 33 ms 9304 KB Output is correct
84 Correct 35 ms 9304 KB Output is correct
85 Incorrect 51 ms 9464 KB Output isn't correct
86 Correct 43 ms 9308 KB Output is correct
87 Correct 45 ms 9472 KB Output is correct
88 Correct 28 ms 9328 KB Output is correct
89 Correct 38 ms 9476 KB Output is correct
90 Correct 48 ms 9308 KB Output is correct
91 Incorrect 52 ms 9308 KB Output isn't correct
92 Correct 32 ms 9304 KB Output is correct
93 Correct 22 ms 9308 KB Output is correct
94 Correct 52 ms 9308 KB Output is correct
95 Correct 41 ms 9452 KB Output is correct
96 Correct 25 ms 9308 KB Output is correct
97 Correct 35 ms 9476 KB Output is correct
98 Incorrect 42 ms 9456 KB Output isn't correct
99 Correct 45 ms 9308 KB Output is correct
100 Correct 41 ms 9456 KB Output is correct
101 Correct 54 ms 9308 KB Output is correct
102 Correct 42 ms 9304 KB Output is correct