for row in data:
    if row[2] == dong_code:
        population[row[1]] += row[3]

반복문이 끝나면 과정 3에서 population 리스트를 31로 나누어 시간대별 평균인구를 구합니다. 이때도 리스트 내포를 활용합니다. population의 각 요소를 p에 차례대로 담은 후 p를 31로 나눈 값을 리스트에 다시 저장하면 됩니다.

population = [p/31 for p in population]
print(population)

실행결과
[37098.63627419355, 35891.815432258074 … 43570.433745161296, 39358.195206451615]

population의 모든 요소를 31로 나누니 숫자가 확 줄었습니다. 0시 인구는 3만7천 명, 1시 인구는 3만6천 명 정도입니다. 오후에 사람이 많을 때는 대략 6만3천 명까지 있고요.