Submission #969251

#	Submission time	Handle	Problem	Language	Result	Execution time	Memory
969251	2024-04-24T19:45:15 Z	anango	Art Class (IOI13_artclass)	C++17	100 / 100	62 ms	9816 KB

artclass

#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	Correct	44 ms	9304 KB	Output is correct
2	Correct	23 ms	9308 KB	Output is correct
3	Correct	40 ms	9304 KB	Output is correct
4	Correct	52 ms	9552 KB	Output is correct
5	Correct	52 ms	9308 KB	Output is correct
6	Correct	29 ms	9308 KB	Output is correct
7	Correct	35 ms	9304 KB	Output is correct
8	Correct	33 ms	9304 KB	Output is correct
9	Correct	52 ms	9308 KB	Output is correct
10	Incorrect	52 ms	9456 KB	Output isn't correct
11	Incorrect	41 ms	9560 KB	Output isn't correct
12	Correct	58 ms	9308 KB	Output is correct
13	Correct	46 ms	9456 KB	Output is correct
14	Correct	49 ms	9456 KB	Output is correct
15	Correct	44 ms	9304 KB	Output is correct
16	Correct	42 ms	9308 KB	Output is correct
17	Correct	42 ms	9304 KB	Output is correct
18	Correct	45 ms	9304 KB	Output is correct
19	Correct	34 ms	9304 KB	Output is correct
20	Correct	54 ms	9308 KB	Output is correct
21	Correct	42 ms	9452 KB	Output is correct
22	Incorrect	31 ms	9308 KB	Output isn't correct
23	Correct	28 ms	9476 KB	Output is correct
24	Correct	48 ms	9308 KB	Output is correct
25	Correct	40 ms	9308 KB	Output is correct
26	Correct	41 ms	9304 KB	Output is correct
27	Correct	14 ms	9308 KB	Output is correct
28	Correct	47 ms	9312 KB	Output is correct
29	Correct	28 ms	9304 KB	Output is correct
30	Correct	45 ms	9308 KB	Output is correct
31	Correct	46 ms	9308 KB	Output is correct
32	Correct	40 ms	9308 KB	Output is correct
33	Correct	52 ms	9308 KB	Output is correct
34	Correct	41 ms	9308 KB	Output is correct
35	Correct	42 ms	9308 KB	Output is correct
36	Correct	12 ms	9412 KB	Output is correct
37	Correct	50 ms	9456 KB	Output is correct
38	Correct	37 ms	9308 KB	Output is correct
39	Correct	53 ms	9308 KB	Output is correct
40	Correct	44 ms	9308 KB	Output is correct
41	Correct	34 ms	9452 KB	Output is correct
42	Correct	31 ms	9480 KB	Output is correct
43	Correct	49 ms	9816 KB	Output is correct
44	Incorrect	44 ms	9308 KB	Output isn't correct
45	Correct	50 ms	9308 KB	Output is correct
46	Correct	47 ms	9476 KB	Output is correct
47	Correct	40 ms	9304 KB	Output is correct
48	Correct	43 ms	9304 KB	Output is correct
49	Correct	37 ms	9556 KB	Output is correct
50	Correct	45 ms	9308 KB	Output is correct
51	Correct	26 ms	9308 KB	Output is correct
52	Correct	44 ms	9384 KB	Output is correct
53	Correct	42 ms	9304 KB	Output is correct
54	Correct	39 ms	9308 KB	Output is correct
55	Correct	42 ms	9304 KB	Output is correct
56	Correct	57 ms	9452 KB	Output is correct
57	Correct	40 ms	9304 KB	Output is correct
58	Correct	33 ms	9308 KB	Output is correct
59	Correct	22 ms	9304 KB	Output is correct
60	Correct	48 ms	9452 KB	Output is correct
61	Correct	43 ms	9308 KB	Output is correct
62	Correct	46 ms	9308 KB	Output is correct
63	Correct	33 ms	9424 KB	Output is correct
64	Correct	42 ms	9308 KB	Output is correct
65	Correct	41 ms	9304 KB	Output is correct
66	Correct	48 ms	9308 KB	Output is correct
67	Correct	24 ms	9308 KB	Output is correct
68	Correct	42 ms	9308 KB	Output is correct
69	Correct	40 ms	9308 KB	Output is correct
70	Correct	44 ms	9308 KB	Output is correct
71	Incorrect	49 ms	9412 KB	Output isn't correct
72	Correct	52 ms	9308 KB	Output is correct
73	Correct	45 ms	9308 KB	Output is correct
74	Correct	41 ms	9308 KB	Output is correct
75	Correct	51 ms	9308 KB	Output is correct
76	Correct	41 ms	9308 KB	Output is correct
77	Correct	39 ms	9304 KB	Output is correct
78	Correct	38 ms	9308 KB	Output is correct
79	Correct	49 ms	9452 KB	Output is correct
80	Correct	49 ms	9308 KB	Output is correct
81	Correct	52 ms	9308 KB	Output is correct
82	Correct	45 ms	9456 KB	Output is correct
83	Correct	48 ms	9308 KB	Output is correct
84	Incorrect	62 ms	9460 KB	Output isn't correct
85	Correct	44 ms	9304 KB	Output is correct
86	Correct	35 ms	9308 KB	Output is correct
87	Correct	43 ms	9456 KB	Output is correct
88	Correct	45 ms	9304 KB	Output is correct
89	Correct	32 ms	9304 KB	Output is correct
90	Correct	54 ms	9308 KB	Output is correct
91	Correct	28 ms	9308 KB	Output is correct
92	Correct	40 ms	9408 KB	Output is correct
93	Correct	40 ms	9308 KB	Output is correct
94	Incorrect	50 ms	9456 KB	Output isn't correct
95	Correct	33 ms	9308 KB	Output is correct
96	Correct	47 ms	9452 KB	Output is correct
97	Correct	44 ms	9304 KB	Output is correct
98	Correct	57 ms	9456 KB	Output is correct
99	Correct	36 ms	9304 KB	Output is correct
100	Correct	36 ms	9304 KB	Output is correct
101	Correct	40 ms	9308 KB	Output is correct
102	Correct	46 ms	9452 KB	Output is correct