[시뮬레이션] BOJ 17144 미세먼지 안녕!

import sys
import copy
input = sys.stdin.readline
 
DELTAS = ((1, 0), (-1, 0), (0, 1), (0, -1))
R, C, T = map(int, input().split())
arr = [list(map(int, input().split())) for _ in range(R)]
 
# 시뮬레이션 문제
 
def diffuse(arr):
    new_arr = copy.deepcopy(arr)
    for y in range(R):
        for x in range(C):
            if arr[y][x] < 5: continue
            val = arr[y][x] // 5
            cnt = 0
            for dy, dx in DELTAS:
                ny = y + dy 
                nx = x + dx
                if ny < 0 or ny >= R or nx < 0 or nx >= C: continue 
                if arr[ny][nx] == -1: continue
                new_arr[ny][nx] += val
                cnt += 1
            new_arr[y][x] -= (val * cnt) 
    return new_arr
 
def clean_up(arr, y):
    prev = arr[y][1]
    arr[y][1] = 0
    for j in range(2, C-1):
        tmp =  arr[y][j]
        arr[y][j] = prev 
        prev = tmp
    for i in range(y, 0, -1):
        tmp = arr[i][C-1]
        arr[i][C-1] = prev 
        prev = tmp 
    for j in range(C-1, 0, -1):
        tmp = arr[0][j]
        arr[0][j] = prev 
        prev = tmp         
    for i in range(0, y):
        tmp = arr[i][0]
        arr[i][0] = prev 
        prev = tmp 
 
def clean_down(arr, y):
    prev = arr[y][1]
    arr[y][1] = 0
    for j in range(2, C-1):
        tmp =  arr[y][j]
        arr[y][j] = prev 
        prev = tmp
    for i in range(y, R-1):
        tmp = arr[i][C-1]
        arr[i][C-1] = prev 
        prev = tmp 
    for j in range(C-1, 0, -1):
        tmp = arr[R-1][j]
        arr[R-1][j] = prev 
        prev = tmp         
    for i in range(R-1, y, -1):
        tmp = arr[i][0]
        arr[i][0] = prev 
        prev = tmp 
 
            
clean_y = R
for i in range(R):
    if arr[i][0] == -1:
        clean_y = i 
        break
    
for t in range(T):
    arr = diffuse(arr)
    clean_up(arr, clean_y)
    clean_down(arr, clean_y+1)
 
print(sum([sum(a) for a in arr]) + 2)