FFmpeg coverage


Directory: ../../../ffmpeg/
File: src/libavcodec/aom_film_grain.c
Date: 2024-04-19 17:50:32
Exec Total Coverage
Lines: 3 205 1.5%
Functions: 1 5 20.0%
Branches: 1 129 0.8%

Line Branch Exec Source
1 /*
2 * AOM film grain synthesis
3 * Copyright (c) 2023 Niklas Haas <ffmpeg@haasn.xyz>
4 *
5 * This file is part of FFmpeg.
6 *
7 * FFmpeg is free software; you can redistribute it and/or
8 * modify it under the terms of the GNU Lesser General Public
9 * License as published by the Free Software Foundation; either
10 * version 2.1 of the License, or (at your option) any later version.
11 *
12 * FFmpeg is distributed in the hope that it will be useful,
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15 * Lesser General Public License for more details.
16 *
17 * You should have received a copy of the GNU Lesser General Public
18 * License along with FFmpeg; if not, write to the Free Software
19 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
20 */
21
22 /**
23 * @file
24 * AOM film grain synthesis.
25 * @author Niklas Haas <ffmpeg@haasn.xyz>
26 */
27
28 #include "libavutil/avassert.h"
29 #include "libavutil/imgutils.h"
30
31 #include "aom_film_grain.h"
32 #include "get_bits.h"
33
34 // Common/shared helpers (not dependent on BIT_DEPTH)
35 static inline int get_random_number(const int bits, unsigned *const state) {
36 const int r = *state;
37 unsigned bit = ((r >> 0) ^ (r >> 1) ^ (r >> 3) ^ (r >> 12)) & 1;
38 *state = (r >> 1) | (bit << 15);
39
40 return (*state >> (16 - bits)) & ((1 << bits) - 1);
41 }
42
43 static inline int round2(const int x, const uint64_t shift) {
44 return (x + ((1 << shift) >> 1)) >> shift;
45 }
46
47 enum {
48 GRAIN_WIDTH = 82,
49 GRAIN_HEIGHT = 73,
50 SUB_GRAIN_WIDTH = 44,
51 SUB_GRAIN_HEIGHT = 38,
52 FG_BLOCK_SIZE = 32,
53 };
54
55 static const int16_t gaussian_sequence[2048];
56
57 #define BIT_DEPTH 16
58 #include "aom_film_grain_template.c"
59 #undef BIT_DEPTH
60
61 #define BIT_DEPTH 8
62 #include "aom_film_grain_template.c"
63 #undef BIT_DEPTH
64
65
66 int ff_aom_apply_film_grain(AVFrame *out, const AVFrame *in,
67 const AVFilmGrainParams *params)
68 {
69 const AVFilmGrainAOMParams *const data = &params->codec.aom;
70 const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(out->format);
71 const int subx = desc->log2_chroma_w, suby = desc->log2_chroma_h;
72 const int pxstep = desc->comp[0].step;
73
74 av_assert0(out->format == in->format);
75 av_assert0(params->type == AV_FILM_GRAIN_PARAMS_AV1);
76
77 // Copy over the non-modified planes
78 if (!params->codec.aom.num_y_points) {
79 av_image_copy_plane(out->data[0], out->linesize[0],
80 in->data[0], in->linesize[0],
81 out->width * pxstep, out->height);
82 }
83 for (int uv = 0; uv < 2; uv++) {
84 if (!data->num_uv_points[uv]) {
85 av_image_copy_plane(out->data[1+uv], out->linesize[1+uv],
86 in->data[1+uv], in->linesize[1+uv],
87 AV_CEIL_RSHIFT(out->width, subx) * pxstep,
88 AV_CEIL_RSHIFT(out->height, suby));
89 }
90 }
91
92 switch (in->format) {
93 case AV_PIX_FMT_GRAY8:
94 case AV_PIX_FMT_YUV420P:
95 case AV_PIX_FMT_YUV422P:
96 case AV_PIX_FMT_YUV444P:
97 case AV_PIX_FMT_YUVJ420P:
98 case AV_PIX_FMT_YUVJ422P:
99 case AV_PIX_FMT_YUVJ444P:
100 return apply_film_grain_8(out, in, params);
101 case AV_PIX_FMT_GRAY9:
102 case AV_PIX_FMT_YUV420P9:
103 case AV_PIX_FMT_YUV422P9:
104 case AV_PIX_FMT_YUV444P9:
105 return apply_film_grain_16(out, in, params, 9);
106 case AV_PIX_FMT_GRAY10:
107 case AV_PIX_FMT_YUV420P10:
108 case AV_PIX_FMT_YUV422P10:
109 case AV_PIX_FMT_YUV444P10:
110 return apply_film_grain_16(out, in, params, 10);
111 case AV_PIX_FMT_GRAY12:
112 case AV_PIX_FMT_YUV420P12:
113 case AV_PIX_FMT_YUV422P12:
114 case AV_PIX_FMT_YUV444P12:
115 return apply_film_grain_16(out, in, params, 12);
116 }
117
118 /* The AV1 spec only defines film grain synthesis for these formats */
119 return AVERROR_INVALIDDATA;
120 }
121
122 int ff_aom_parse_film_grain_sets(AVFilmGrainAFGS1Params *s,
123 const uint8_t *payload, int payload_size)
124 {
125 GetBitContext gbc, *gb = &gbc;
126 AVFilmGrainAOMParams *aom;
127 AVFilmGrainParams *fgp, *ref = NULL;
128 int ret, num_sets, n, i, uv, num_y_coeffs, update_grain, luma_only;
129
130 ret = init_get_bits8(gb, payload, payload_size);
131 if (ret < 0)
132 return ret;
133
134 s->enable = get_bits1(gb);
135 if (!s->enable)
136 return 0;
137
138 skip_bits(gb, 4); // reserved
139 num_sets = get_bits(gb, 3) + 1;
140 for (n = 0; n < num_sets; n++) {
141 int payload_4byte, payload_size, set_idx, apply_units_log2, vsc_flag;
142 int predict_scaling, predict_y_scaling, predict_uv_scaling[2];
143 int payload_bits, start_position;
144
145 start_position = get_bits_count(gb);
146 payload_4byte = get_bits1(gb);
147 payload_size = get_bits(gb, payload_4byte ? 2 : 8);
148 set_idx = get_bits(gb, 3);
149 fgp = &s->sets[set_idx];
150 aom = &fgp->codec.aom;
151
152 fgp->type = get_bits1(gb) ? AV_FILM_GRAIN_PARAMS_AV1 : AV_FILM_GRAIN_PARAMS_NONE;
153 if (!fgp->type)
154 continue;
155
156 fgp->seed = get_bits(gb, 16);
157 update_grain = get_bits1(gb);
158 if (!update_grain)
159 continue;
160
161 apply_units_log2 = get_bits(gb, 4);
162 fgp->width = get_bits(gb, 12) << apply_units_log2;
163 fgp->height = get_bits(gb, 12) << apply_units_log2;
164 luma_only = get_bits1(gb);
165 if (luma_only) {
166 fgp->subsampling_x = fgp->subsampling_y = 0;
167 } else {
168 fgp->subsampling_x = get_bits1(gb);
169 fgp->subsampling_y = get_bits1(gb);
170 }
171
172 fgp->bit_depth_luma = fgp->bit_depth_chroma = 0;
173 fgp->color_primaries = AVCOL_PRI_UNSPECIFIED;
174 fgp->color_trc = AVCOL_TRC_UNSPECIFIED;
175 fgp->color_space = AVCOL_SPC_UNSPECIFIED;
176 fgp->color_range = AVCOL_RANGE_UNSPECIFIED;
177
178 vsc_flag = get_bits1(gb); // video_signal_characteristics_flag
179 if (vsc_flag) {
180 int cicp_flag;
181 fgp->bit_depth_luma = get_bits(gb, 3) + 8;
182 if (!luma_only)
183 fgp->bit_depth_chroma = fgp->bit_depth_luma;
184 cicp_flag = get_bits1(gb);
185 if (cicp_flag) {
186 fgp->color_primaries = get_bits(gb, 8);
187 fgp->color_trc = get_bits(gb, 8);
188 fgp->color_space = get_bits(gb, 8);
189 fgp->color_range = get_bits1(gb) ? AVCOL_RANGE_JPEG : AVCOL_RANGE_MPEG;
190 if (fgp->color_primaries > AVCOL_PRI_NB ||
191 fgp->color_primaries == AVCOL_PRI_RESERVED ||
192 fgp->color_primaries == AVCOL_PRI_RESERVED0 ||
193 fgp->color_trc > AVCOL_TRC_NB ||
194 fgp->color_trc == AVCOL_TRC_RESERVED ||
195 fgp->color_trc == AVCOL_TRC_RESERVED0 ||
196 fgp->color_space > AVCOL_SPC_NB ||
197 fgp->color_space == AVCOL_SPC_RESERVED)
198 goto error;
199 }
200 }
201
202 predict_scaling = get_bits1(gb);
203 if (predict_scaling && (!ref || ref == fgp))
204 goto error; // prediction must be from valid, different set
205
206 predict_y_scaling = predict_scaling ? get_bits1(gb) : 0;
207 if (predict_y_scaling) {
208 int y_scale, y_offset, bits_res;
209 y_scale = get_bits(gb, 9) - 256;
210 y_offset = get_bits(gb, 9) - 256;
211 bits_res = get_bits(gb, 3);
212 if (bits_res) {
213 int res[14], pred, granularity;
214 aom->num_y_points = ref->codec.aom.num_y_points;
215 for (i = 0; i < aom->num_y_points; i++)
216 res[i] = get_bits(gb, bits_res);
217 granularity = get_bits(gb, 3);
218 for (i = 0; i < aom->num_y_points; i++) {
219 pred = ref->codec.aom.y_points[i][1];
220 pred = ((pred * y_scale + 8) >> 4) + y_offset;
221 pred += (res[i] - (1 << (bits_res - 1))) * granularity;
222 aom->y_points[i][0] = ref->codec.aom.y_points[i][0];
223 aom->y_points[i][1] = av_clip_uint8(pred);
224 }
225 }
226 } else {
227 aom->num_y_points = get_bits(gb, 4);
228 if (aom->num_y_points > 14) {
229 goto error;
230 } else if (aom->num_y_points) {
231 int bits_inc, bits_scaling;
232 int y_value = 0;
233 bits_inc = get_bits(gb, 3) + 1;
234 bits_scaling = get_bits(gb, 2) + 5;
235 for (i = 0; i < aom->num_y_points; i++) {
236 y_value += get_bits(gb, bits_inc);
237 if (y_value > UINT8_MAX)
238 goto error;
239 aom->y_points[i][0] = y_value;
240 aom->y_points[i][1] = get_bits(gb, bits_scaling);
241 }
242 }
243 }
244
245 if (luma_only) {
246 aom->chroma_scaling_from_luma = 0;
247 aom->num_uv_points[0] = aom->num_uv_points[1] = 0;
248 } else {
249 aom->chroma_scaling_from_luma = get_bits1(gb);
250 if (aom->chroma_scaling_from_luma) {
251 aom->num_uv_points[0] = aom->num_uv_points[1] = 0;
252 } else {
253 for (uv = 0; uv < 2; uv++) {
254 predict_uv_scaling[uv] = predict_scaling ? get_bits1(gb) : 0;
255 if (predict_uv_scaling[uv]) {
256 int uv_scale, uv_offset, bits_res;
257 uv_scale = get_bits(gb, 9) - 256;
258 uv_offset = get_bits(gb, 9) - 256;
259 bits_res = get_bits(gb, 3);
260 aom->uv_mult[uv] = ref->codec.aom.uv_mult[uv];
261 aom->uv_mult_luma[uv] = ref->codec.aom.uv_mult_luma[uv];
262 aom->uv_offset[uv] = ref->codec.aom.uv_offset[uv];
263 if (bits_res) {
264 int res[10], pred, granularity;
265 aom->num_uv_points[uv] = ref->codec.aom.num_uv_points[uv];
266 for (i = 0; i < aom->num_uv_points[uv]; i++)
267 res[i] = get_bits(gb, bits_res);
268 granularity = get_bits(gb, 3);
269 for (i = 0; i < aom->num_uv_points[uv]; i++) {
270 pred = ref->codec.aom.uv_points[uv][i][1];
271 pred = ((pred * uv_scale + 8) >> 4) + uv_offset;
272 pred += (res[i] - (1 << (bits_res - 1))) * granularity;
273 aom->uv_points[uv][i][0] = ref->codec.aom.uv_points[uv][i][0];
274 aom->uv_points[uv][i][1] = av_clip_uint8(pred);
275 }
276 }
277 } else {
278 int bits_inc, bits_scaling, uv_offset;
279 int uv_value = 0;
280 aom->num_uv_points[uv] = get_bits(gb, 4);
281 if (aom->num_uv_points[uv] > 10)
282 goto error;
283 bits_inc = get_bits(gb, 3) + 1;
284 bits_scaling = get_bits(gb, 2) + 5;
285 uv_offset = get_bits(gb, 8);
286 for (i = 0; i < aom->num_uv_points[uv]; i++) {
287 uv_value += get_bits(gb, bits_inc);
288 if (uv_value > UINT8_MAX)
289 goto error;
290 aom->uv_points[uv][i][0] = uv_value;
291 aom->uv_points[uv][i][1] = get_bits(gb, bits_scaling) + uv_offset;
292 }
293 }
294 }
295 }
296 }
297
298 aom->scaling_shift = get_bits(gb, 2) + 8;
299 aom->ar_coeff_lag = get_bits(gb, 2);
300 num_y_coeffs = 2 * aom->ar_coeff_lag * (aom->ar_coeff_lag + 1);
301 if (aom->num_y_points) {
302 int ar_bits = get_bits(gb, 2) + 5;
303 for (i = 0; i < num_y_coeffs; i++)
304 aom->ar_coeffs_y[i] = get_bits(gb, ar_bits) - (1 << (ar_bits - 1));
305 }
306 for (uv = 0; uv < 2; uv++) {
307 if (aom->chroma_scaling_from_luma || aom->num_uv_points[uv]) {
308 int ar_bits = get_bits(gb, 2) + 5;
309 for (i = 0; i < num_y_coeffs + !!aom->num_y_points; i++)
310 aom->ar_coeffs_uv[uv][i] = get_bits(gb, ar_bits) - (1 << (ar_bits - 1));
311 }
312 }
313 aom->ar_coeff_shift = get_bits(gb, 2) + 6;
314 aom->grain_scale_shift = get_bits(gb, 2);
315 for (uv = 0; uv < 2; uv++) {
316 if (aom->num_uv_points[uv] && !predict_uv_scaling[uv]) {
317 aom->uv_mult[uv] = get_bits(gb, 8) - 128;
318 aom->uv_mult_luma[uv] = get_bits(gb, 8) - 128;
319 aom->uv_offset[uv] = get_bits(gb, 9) - 256;
320 }
321 }
322 aom->overlap_flag = get_bits1(gb);
323 aom->limit_output_range = get_bits1(gb);
324
325 // use first set as reference only if it was fully transmitted
326 if (n == 0)
327 ref = fgp;
328
329 payload_bits = get_bits_count(gb) - start_position;
330 if (payload_bits > payload_size * 8)
331 goto error;
332 skip_bits(gb, payload_size * 8 - payload_bits);
333 }
334 return 0;
335
336 error:
337 memset(s, 0, sizeof(*s));
338 return AVERROR_INVALIDDATA;
339 }
340
341 34457 int ff_aom_attach_film_grain_sets(const AVFilmGrainAFGS1Params *s, AVFrame *frame)
342 {
343 AVFilmGrainParams *fgp;
344
1/2
✓ Branch 0 taken 34457 times.
✗ Branch 1 not taken.
34457 if (!s->enable)
345 34457 return 0;
346
347 for (int i = 0; i < FF_ARRAY_ELEMS(s->sets); i++) {
348 if (s->sets[i].type != AV_FILM_GRAIN_PARAMS_AV1)
349 continue;
350 fgp = av_film_grain_params_create_side_data(frame);
351 if (!fgp)
352 return AVERROR(ENOMEM);
353 memcpy(fgp, &s->sets[i], sizeof(*fgp));
354 }
355
356 return 0;
357 }
358
359 // Taken from the AV1 spec. Range is [-2048, 2047], mean is 0 and stddev is 512
360 static const int16_t gaussian_sequence[2048] = {
361 56, 568, -180, 172, 124, -84, 172, -64, -900, 24, 820,
362 224, 1248, 996, 272, -8, -916, -388, -732, -104, -188, 800,
363 112, -652, -320, -376, 140, -252, 492, -168, 44, -788, 588,
364 -584, 500, -228, 12, 680, 272, -476, 972, -100, 652, 368,
365 432, -196, -720, -192, 1000, -332, 652, -136, -552, -604, -4,
366 192, -220, -136, 1000, -52, 372, -96, -624, 124, -24, 396,
367 540, -12, -104, 640, 464, 244, -208, -84, 368, -528, -740,
368 248, -968, -848, 608, 376, -60, -292, -40, -156, 252, -292,
369 248, 224, -280, 400, -244, 244, -60, 76, -80, 212, 532,
370 340, 128, -36, 824, -352, -60, -264, -96, -612, 416, -704,
371 220, -204, 640, -160, 1220, -408, 900, 336, 20, -336, -96,
372 -792, 304, 48, -28, -1232, -1172, -448, 104, -292, -520, 244,
373 60, -948, 0, -708, 268, 108, 356, -548, 488, -344, -136,
374 488, -196, -224, 656, -236, -1128, 60, 4, 140, 276, -676,
375 -376, 168, -108, 464, 8, 564, 64, 240, 308, -300, -400,
376 -456, -136, 56, 120, -408, -116, 436, 504, -232, 328, 844,
377 -164, -84, 784, -168, 232, -224, 348, -376, 128, 568, 96,
378 -1244, -288, 276, 848, 832, -360, 656, 464, -384, -332, -356,
379 728, -388, 160, -192, 468, 296, 224, 140, -776, -100, 280,
380 4, 196, 44, -36, -648, 932, 16, 1428, 28, 528, 808,
381 772, 20, 268, 88, -332, -284, 124, -384, -448, 208, -228,
382 -1044, -328, 660, 380, -148, -300, 588, 240, 540, 28, 136,
383 -88, -436, 256, 296, -1000, 1400, 0, -48, 1056, -136, 264,
384 -528, -1108, 632, -484, -592, -344, 796, 124, -668, -768, 388,
385 1296, -232, -188, -200, -288, -4, 308, 100, -168, 256, -500,
386 204, -508, 648, -136, 372, -272, -120, -1004, -552, -548, -384,
387 548, -296, 428, -108, -8, -912, -324, -224, -88, -112, -220,
388 -100, 996, -796, 548, 360, -216, 180, 428, -200, -212, 148,
389 96, 148, 284, 216, -412, -320, 120, -300, -384, -604, -572,
390 -332, -8, -180, -176, 696, 116, -88, 628, 76, 44, -516,
391 240, -208, -40, 100, -592, 344, -308, -452, -228, 20, 916,
392 -1752, -136, -340, -804, 140, 40, 512, 340, 248, 184, -492,
393 896, -156, 932, -628, 328, -688, -448, -616, -752, -100, 560,
394 -1020, 180, -800, -64, 76, 576, 1068, 396, 660, 552, -108,
395 -28, 320, -628, 312, -92, -92, -472, 268, 16, 560, 516,
396 -672, -52, 492, -100, 260, 384, 284, 292, 304, -148, 88,
397 -152, 1012, 1064, -228, 164, -376, -684, 592, -392, 156, 196,
398 -524, -64, -884, 160, -176, 636, 648, 404, -396, -436, 864,
399 424, -728, 988, -604, 904, -592, 296, -224, 536, -176, -920,
400 436, -48, 1176, -884, 416, -776, -824, -884, 524, -548, -564,
401 -68, -164, -96, 692, 364, -692, -1012, -68, 260, -480, 876,
402 -1116, 452, -332, -352, 892, -1088, 1220, -676, 12, -292, 244,
403 496, 372, -32, 280, 200, 112, -440, -96, 24, -644, -184,
404 56, -432, 224, -980, 272, -260, 144, -436, 420, 356, 364,
405 -528, 76, 172, -744, -368, 404, -752, -416, 684, -688, 72,
406 540, 416, 92, 444, 480, -72, -1416, 164, -1172, -68, 24,
407 424, 264, 1040, 128, -912, -524, -356, 64, 876, -12, 4,
408 -88, 532, 272, -524, 320, 276, -508, 940, 24, -400, -120,
409 756, 60, 236, -412, 100, 376, -484, 400, -100, -740, -108,
410 -260, 328, -268, 224, -200, -416, 184, -604, -564, -20, 296,
411 60, 892, -888, 60, 164, 68, -760, 216, -296, 904, -336,
412 -28, 404, -356, -568, -208, -1480, -512, 296, 328, -360, -164,
413 -1560, -776, 1156, -428, 164, -504, -112, 120, -216, -148, -264,
414 308, 32, 64, -72, 72, 116, 176, -64, -272, 460, -536,
415 -784, -280, 348, 108, -752, -132, 524, -540, -776, 116, -296,
416 -1196, -288, -560, 1040, -472, 116, -848, -1116, 116, 636, 696,
417 284, -176, 1016, 204, -864, -648, -248, 356, 972, -584, -204,
418 264, 880, 528, -24, -184, 116, 448, -144, 828, 524, 212,
419 -212, 52, 12, 200, 268, -488, -404, -880, 824, -672, -40,
420 908, -248, 500, 716, -576, 492, -576, 16, 720, -108, 384,
421 124, 344, 280, 576, -500, 252, 104, -308, 196, -188, -8,
422 1268, 296, 1032, -1196, 436, 316, 372, -432, -200, -660, 704,
423 -224, 596, -132, 268, 32, -452, 884, 104, -1008, 424, -1348,
424 -280, 4, -1168, 368, 476, 696, 300, -8, 24, 180, -592,
425 -196, 388, 304, 500, 724, -160, 244, -84, 272, -256, -420,
426 320, 208, -144, -156, 156, 364, 452, 28, 540, 316, 220,
427 -644, -248, 464, 72, 360, 32, -388, 496, -680, -48, 208,
428 -116, -408, 60, -604, -392, 548, -840, 784, -460, 656, -544,
429 -388, -264, 908, -800, -628, -612, -568, 572, -220, 164, 288,
430 -16, -308, 308, -112, -636, -760, 280, -668, 432, 364, 240,
431 -196, 604, 340, 384, 196, 592, -44, -500, 432, -580, -132,
432 636, -76, 392, 4, -412, 540, 508, 328, -356, -36, 16,
433 -220, -64, -248, -60, 24, -192, 368, 1040, 92, -24, -1044,
434 -32, 40, 104, 148, 192, -136, -520, 56, -816, -224, 732,
435 392, 356, 212, -80, -424, -1008, -324, 588, -1496, 576, 460,
436 -816, -848, 56, -580, -92, -1372, -112, -496, 200, 364, 52,
437 -140, 48, -48, -60, 84, 72, 40, 132, -356, -268, -104,
438 -284, -404, 732, -520, 164, -304, -540, 120, 328, -76, -460,
439 756, 388, 588, 236, -436, -72, -176, -404, -316, -148, 716,
440 -604, 404, -72, -88, -888, -68, 944, 88, -220, -344, 960,
441 472, 460, -232, 704, 120, 832, -228, 692, -508, 132, -476,
442 844, -748, -364, -44, 1116, -1104, -1056, 76, 428, 552, -692,
443 60, 356, 96, -384, -188, -612, -576, 736, 508, 892, 352,
444 -1132, 504, -24, -352, 324, 332, -600, -312, 292, 508, -144,
445 -8, 484, 48, 284, -260, -240, 256, -100, -292, -204, -44,
446 472, -204, 908, -188, -1000, -256, 92, 1164, -392, 564, 356,
447 652, -28, -884, 256, 484, -192, 760, -176, 376, -524, -452,
448 -436, 860, -736, 212, 124, 504, -476, 468, 76, -472, 552,
449 -692, -944, -620, 740, -240, 400, 132, 20, 192, -196, 264,
450 -668, -1012, -60, 296, -316, -828, 76, -156, 284, -768, -448,
451 -832, 148, 248, 652, 616, 1236, 288, -328, -400, -124, 588,
452 220, 520, -696, 1032, 768, -740, -92, -272, 296, 448, -464,
453 412, -200, 392, 440, -200, 264, -152, -260, 320, 1032, 216,
454 320, -8, -64, 156, -1016, 1084, 1172, 536, 484, -432, 132,
455 372, -52, -256, 84, 116, -352, 48, 116, 304, -384, 412,
456 924, -300, 528, 628, 180, 648, 44, -980, -220, 1320, 48,
457 332, 748, 524, -268, -720, 540, -276, 564, -344, -208, -196,
458 436, 896, 88, -392, 132, 80, -964, -288, 568, 56, -48,
459 -456, 888, 8, 552, -156, -292, 948, 288, 128, -716, -292,
460 1192, -152, 876, 352, -600, -260, -812, -468, -28, -120, -32,
461 -44, 1284, 496, 192, 464, 312, -76, -516, -380, -456, -1012,
462 -48, 308, -156, 36, 492, -156, -808, 188, 1652, 68, -120,
463 -116, 316, 160, -140, 352, 808, -416, 592, 316, -480, 56,
464 528, -204, -568, 372, -232, 752, -344, 744, -4, 324, -416,
465 -600, 768, 268, -248, -88, -132, -420, -432, 80, -288, 404,
466 -316, -1216, -588, 520, -108, 92, -320, 368, -480, -216, -92,
467 1688, -300, 180, 1020, -176, 820, -68, -228, -260, 436, -904,
468 20, 40, -508, 440, -736, 312, 332, 204, 760, -372, 728,
469 96, -20, -632, -520, -560, 336, 1076, -64, -532, 776, 584,
470 192, 396, -728, -520, 276, -188, 80, -52, -612, -252, -48,
471 648, 212, -688, 228, -52, -260, 428, -412, -272, -404, 180,
472 816, -796, 48, 152, 484, -88, -216, 988, 696, 188, -528,
473 648, -116, -180, 316, 476, 12, -564, 96, 476, -252, -364,
474 -376, -392, 556, -256, -576, 260, -352, 120, -16, -136, -260,
475 -492, 72, 556, 660, 580, 616, 772, 436, 424, -32, -324,
476 -1268, 416, -324, -80, 920, 160, 228, 724, 32, -516, 64,
477 384, 68, -128, 136, 240, 248, -204, -68, 252, -932, -120,
478 -480, -628, -84, 192, 852, -404, -288, -132, 204, 100, 168,
479 -68, -196, -868, 460, 1080, 380, -80, 244, 0, 484, -888,
480 64, 184, 352, 600, 460, 164, 604, -196, 320, -64, 588,
481 -184, 228, 12, 372, 48, -848, -344, 224, 208, -200, 484,
482 128, -20, 272, -468, -840, 384, 256, -720, -520, -464, -580,
483 112, -120, 644, -356, -208, -608, -528, 704, 560, -424, 392,
484 828, 40, 84, 200, -152, 0, -144, 584, 280, -120, 80,
485 -556, -972, -196, -472, 724, 80, 168, -32, 88, 160, -688,
486 0, 160, 356, 372, -776, 740, -128, 676, -248, -480, 4,
487 -364, 96, 544, 232, -1032, 956, 236, 356, 20, -40, 300,
488 24, -676, -596, 132, 1120, -104, 532, -1096, 568, 648, 444,
489 508, 380, 188, -376, -604, 1488, 424, 24, 756, -220, -192,
490 716, 120, 920, 688, 168, 44, -460, 568, 284, 1144, 1160,
491 600, 424, 888, 656, -356, -320, 220, 316, -176, -724, -188,
492 -816, -628, -348, -228, -380, 1012, -452, -660, 736, 928, 404,
493 -696, -72, -268, -892, 128, 184, -344, -780, 360, 336, 400,
494 344, 428, 548, -112, 136, -228, -216, -820, -516, 340, 92,
495 -136, 116, -300, 376, -244, 100, -316, -520, -284, -12, 824,
496 164, -548, -180, -128, 116, -924, -828, 268, -368, -580, 620,
497 192, 160, 0, -1676, 1068, 424, -56, -360, 468, -156, 720,
498 288, -528, 556, -364, 548, -148, 504, 316, 152, -648, -620,
499 -684, -24, -376, -384, -108, -920, -1032, 768, 180, -264, -508,
500 -1268, -260, -60, 300, -240, 988, 724, -376, -576, -212, -736,
501 556, 192, 1092, -620, -880, 376, -56, -4, -216, -32, 836,
502 268, 396, 1332, 864, -600, 100, 56, -412, -92, 356, 180,
503 884, -468, -436, 292, -388, -804, -704, -840, 368, -348, 140,
504 -724, 1536, 940, 372, 112, -372, 436, -480, 1136, 296, -32,
505 -228, 132, -48, -220, 868, -1016, -60, -1044, -464, 328, 916,
506 244, 12, -736, -296, 360, 468, -376, -108, -92, 788, 368,
507 -56, 544, 400, -672, -420, 728, 16, 320, 44, -284, -380,
508 -796, 488, 132, 204, -596, -372, 88, -152, -908, -636, -572,
509 -624, -116, -692, -200, -56, 276, -88, 484, -324, 948, 864,
510 1000, -456, -184, -276, 292, -296, 156, 676, 320, 160, 908,
511 -84, -1236, -288, -116, 260, -372, -644, 732, -756, -96, 84,
512 344, -520, 348, -688, 240, -84, 216, -1044, -136, -676, -396,
513 -1500, 960, -40, 176, 168, 1516, 420, -504, -344, -364, -360,
514 1216, -940, -380, -212, 252, -660, -708, 484, -444, -152, 928,
515 -120, 1112, 476, -260, 560, -148, -344, 108, -196, 228, -288,
516 504, 560, -328, -88, 288, -1008, 460, -228, 468, -836, -196,
517 76, 388, 232, 412, -1168, -716, -644, 756, -172, -356, -504,
518 116, 432, 528, 48, 476, -168, -608, 448, 160, -532, -272,
519 28, -676, -12, 828, 980, 456, 520, 104, -104, 256, -344,
520 -4, -28, -368, -52, -524, -572, -556, -200, 768, 1124, -208,
521 -512, 176, 232, 248, -148, -888, 604, -600, -304, 804, -156,
522 -212, 488, -192, -804, -256, 368, -360, -916, -328, 228, -240,
523 -448, -472, 856, -556, -364, 572, -12, -156, -368, -340, 432,
524 252, -752, -152, 288, 268, -580, -848, -592, 108, -76, 244,
525 312, -716, 592, -80, 436, 360, 4, -248, 160, 516, 584,
526 732, 44, -468, -280, -292, -156, -588, 28, 308, 912, 24,
527 124, 156, 180, -252, 944, -924, -772, -520, -428, -624, 300,
528 -212, -1144, 32, -724, 800, -1128, -212, -1288, -848, 180, -416,
529 440, 192, -576, -792, -76, -1080, 80, -532, -352, -132, 380,
530 -820, 148, 1112, 128, 164, 456, 700, -924, 144, -668, -384,
531 648, -832, 508, 552, -52, -100, -656, 208, -568, 748, -88,
532 680, 232, 300, 192, -408, -1012, -152, -252, -268, 272, -876,
533 -664, -648, -332, -136, 16, 12, 1152, -28, 332, -536, 320,
534 -672, -460, -316, 532, -260, 228, -40, 1052, -816, 180, 88,
535 -496, -556, -672, -368, 428, 92, 356, 404, -408, 252, 196,
536 -176, -556, 792, 268, 32, 372, 40, 96, -332, 328, 120,
537 372, -900, -40, 472, -264, -592, 952, 128, 656, 112, 664,
538 -232, 420, 4, -344, -464, 556, 244, -416, -32, 252, 0,
539 -412, 188, -696, 508, -476, 324, -1096, 656, -312, 560, 264,
540 -136, 304, 160, -64, -580, 248, 336, -720, 560, -348, -288,
541 -276, -196, -500, 852, -544, -236, -1128, -992, -776, 116, 56,
542 52, 860, 884, 212, -12, 168, 1020, 512, -552, 924, -148,
543 716, 188, 164, -340, -520, -184, 880, -152, -680, -208, -1156,
544 -300, -528, -472, 364, 100, -744, -1056, -32, 540, 280, 144,
545 -676, -32, -232, -280, -224, 96, 568, -76, 172, 148, 148,
546 104, 32, -296, -32, 788, -80, 32, -16, 280, 288, 944,
547 428, -484
548 };
549