#include <stdio.h>  
#include <time.h>  

#define KERNEL_SIZE 3  

float kernel[KERNEL_SIZE][KERNEL_SIZE] = {
    {1 / 16.0, 2 / 16.0, 1 / 16.0},
    {2 / 16.0, 4 / 16.0, 2 / 16.0},
    {1 / 16.0, 2 / 16.0, 1 / 16.0}
};

void applyGaussianBlur(float src[][5], float dst[][5], int h, int w, float kernel[KERNEL_SIZE][KERNEL_SIZE]) {
    for (int i = 1; i < h - 1; i++) {  
        for (int j = 1; j < w - 1; j++) { 
            float sum = 0.0f;
            for (int ki = 0; ki < KERNEL_SIZE; ki++) {
                for (int kj = 0; kj < KERNEL_SIZE; kj++) {
                    sum += src[i + ki - 1][j + kj - 1] * kernel[ki][kj];
                }
            }
            dst[i][j] = sum;
        }
    }
}

int main() {
    float src[5][5] = {
        {1, 2, 3, 4, 5},
        {6, 7, 8, 9, 10},
        {11, 12, 13, 14, 15},
        {16, 17, 18, 19, 20},
        {21, 22, 23, 24, 25}
    };

    float dst[5][5] = { 0 }; 

    clock_t start = clock();  
    applyGaussianBlur(src, dst, 5, 5, kernel);
    clock_t end = clock();  

    printf("模糊后的图像矩阵:\n");
    for (int i = 0; i < 5; i++) {
        for (int j = 0; j < 5; j++) {
            printf("%.2f ", dst[i][j]);
        }
        printf("\n");
    }

    double time_taken = (double)(end - start) / CLOCKS_PER_SEC;
    printf("运行时间: %e 秒\n", time_taken);

    return 0;
}