请教前辈一个关于算法的问题、。。！

/*

 *说明方格的构造顶点坐标(x, y) 宽度(width) 高度(height)

 *程序中顶点坐标是使用一个间接的方式得出 即利用map当中的

 *二位数组的下标分别计算 因为10*10个方格的形状相同.

 */
#include <stdio.h>
#include <math.h>

#define WIDTH    1//宽度
#define HEIGHT    1//高度
#define DIM        10

double p1_x, p1_y, p2_x, p2_y;//两点的坐标
/*线性方程中的参数*/
double k;//斜率
double b;//常量

int map[DIM][DIM];//虚拟位置

/*

 *从终端输入两个点 得出一条直线

 *形如y = k*x + b;

 *1)若返回为1 则表示斜率存在

 *2)若返回为0 则表示斜率不存在

 */
int cacu_k_and_b()
{
    printf("输入第一个点的横纵坐标(例如: x, y): ");
    scanf("%lf, %lf", &p1_x, &p1_y);
    printf("输入第二个点的横纵坐标(例如: x, y): ");
    scanf("%lf, %lf", &p2_x, &p2_y);

    printf("输入的两个点分别为(%.3f, %.3f) 和 (%.3f, %.3f)\n",
        p1_x, p1_y, p2_x, p2_y );

    //斜率
    if (fabs(p1_x - p2_x) < 0.0001)
    {//表示斜率不存在 x = b;
        return 0;
    }
    else
    {
        k = (p1_y-p2_y)/(p1_x-p2_x);
        b = p1_y - k*p1_x;
        //printf("y = %.3fx + %.3f\n", k, b);
        return 1;
    }
}

int isx_in_scope(int x)
{
    if ( p1_x > p2_x )
    {
        if ( p1_x > x*WIDTH && p2_x < (1+x)*WIDTH )
        {
            return 1;
        }
        else
        {
            return 0;
        }
    }
    else
    {
        if ( p1_x < (x+1)*WIDTH && p2_x > x*WIDTH )
        {
            return 1;
        }
        else
        {
            return 0;
        }
    }
}
int isy_in_scope(int y)
{
    if ( p1_y > p2_y )
    {
        if ( p1_y > y*HEIGHT && p2_y < (1+y)*HEIGHT )
        {
            return 1;
        }
        else
        {
            return 0;
        }
    }
    else
    {
        if ( p1_y < (1+y)*HEIGHT && p2_y > y*HEIGHT )
        {
            return 1;
        }
        else
        {
            return 0;
        }
    }
}
/*

 *判断map[y][x]所标记的方框 是否在线段的影响范围内

 *1) 若在则返回 1

 *2) 若不在则返回 0

 */
int is_in_scope(int y, int x)
{
    return isx_in_scope(x) && isy_in_scope(y);
}
/*

 *判断直线与map[y][x]是否存在交点

 *1)若存在 返回1

 *2)若不存在 返回0

 */
int has_common_point( int y, int x)
{
    double p_x, p_y;

    //代入横坐标 x*WIDTH
    //1
    p_y = k * x * WIDTH + b;
    if ( p_y > y*HEIGHT && p_y < (y+1)*HEIGHT )
    {
        return 1;
    }
    //2
    p_y = k * (x+1) * WIDTH + b;
    if ( p_y > y*HEIGHT && p_y < (y+1)*HEIGHT )
    {
        return 1;
    }

    //代入纵坐标 y*WIDHT
    //1
    p_x = (y*HEIGHT-b)/k;
    if ( p_x > x*WIDTH && p_x < (x+1)*WIDTH )
    {
        return 1;
    }
    //2
    p_x = ((y+1)*HEIGHT-b)/k;
    if ( p_x > x*WIDTH && p_x < (x+1)*WIDTH )
    {
        return 1;
    }

    return 0;
}
/*

 *判断直线y = k*x + b 是否经过 虚拟点map[y][x]标记的实际方块

 *1)若是进过 则返回1

 *2)若是不经过 则返回0

 */
int is_cross_block(int y, int x, int flag)
{
    //1.贯穿整个方框则和边界(不同的边)有两个交点

    //2.一个点为线段的终点或者起点 驻留在方块里面 则两端和方块只有一个
    // 交点

    //3. 若斜率为 0  或者不存在的时候 存在这样的情况 直线落在边上
    if (flag)
    {
        if (0 != k)
        {
            if ( !is_in_scope(y, x) )
            {
                return 0;
            }

            if ( has_common_point(y, x) )
            {//如果有公共点则返回1
                return 1;
            }
        }
        else
        {
            if ( !isx_in_scope(x) )
            {
                return 0;
            }

            if ( (b-y*HEIGHT) > (10e-6) &&
                (((y+1)*HEIGHT)-b) > (10e-6) )
            {
                return 1;
            }      
        }
    }
    else
    {//斜率不存在
        if ( !isy_in_scope(y) )
        {
            return 0;
        }

        if ( (b-x*WIDTH) > (10e-6) &&
            (((x+1)*WIDTH)-b) > (10e-6) )
        {
            return 1;
        }
    }

    return 0;
}

void deal()
{
    int flag;//标记斜率是否存在标志
    int i, j;

    flag = cacu_k_and_b();

    for (i=0; i<DIM; ++i)
    {
        for (j=0; j<DIM; ++j)
        {
            if ( is_cross_block(i, j, flag) )
            {//若果经过则把值置为1
                map[i][j] = 1;
            }
        }
    }

    for (i=0; i<DIM; ++i)
    {
        for (j=0; j<DIM; ++j)
        {
            printf("%d ", map[i][j]);
        }
        printf("\n");
    }
    printf("\n");
}

int main(void)
{
    deal();

    return 0;
}