C语言分治法求最近对问题 运行一直报错 求高手修改

Python018

C语言分治法求最近对问题 运行一直报错 求高手修改,第1张

#include<stdio.h>

#include<math.h>

//using namespace std

//using namespace std

#define MAX 100000000

#define N 100000

typedef struct Point

{

double x,y

}Point

Point S[N*2],S1[N],S2[N],P1[N],P2[N]

double Dis(Point a,Point b)

{

return sqrt((a.x-b.x)*(a.x-b.y)+(a.y-b.y)*(a.y-b.y))

}

int CMP1(Point a,Point b)

{

return a.x<b.x

}

int CMP2(Point a,Point b)

{

return a.y <b.y

}

double min(double a,double b)

{

double c

if(a>b)

{

c=b

}

return c

}

void sort(Point *p1,Point *p2,int (* cmp)(Point a,Point b))

{

//这个函数你自己添加吧

}

double ClosePoint(Point S[],int n)

{

if(n<2)

return MAX

else

{

double d0=MAX

int i,j,m

int j1=0,j2=0

double d1,d2,d

int k1=0

int k2=0

sort(S,S+n,CMP1)

Point p=S[n/2]

m=p.x

for(i=0i<(n+1)/2i++)

{

S1[j1].x=S[i].x

S1[j1].y=S[i].y

j1++

}//构造S1中点坐标小于m

for(i=(n+1)/2i<ni++)

{

S2[j2].x=S[i].x

S2[j2].x=S[i].y

}//构造S2中点大于m

d1=ClosePoint(S1,j1)

d2=ClosePoint(S2,j2)

d=min(d1,d2)

for(i=0i<j1i++)

{

if(m-S1[i].x<d)

{

P1[k1].x=S1[i].x

P2[k2].y=S1[i].y

k1++

}

}

for(i=0i<j2i++)

if(S2[i].x-m<d)

{

P2[k2].x=S2[i].x

P2[k2].y=S2[i].y

k2++

}

sort(P1,P1+k1,CMP2)

sort(P2,P2+k2,CMP2)

for(i=0i<k1i++)

{

for(j=0j<k2j++)

{

double A=Dis(P1[i],P2[j])

d0=min(d0,A)

}

}

return min(d0,d)

}

}

int main()

{

int n

int i

double d

printf("请输入端点的个数:\n")

scanf("%d",&n)

printf("请输入端点的坐标:\n")

for(i=1i<=ni++)

{

scanf("%lf",&S[i].x)

scanf("%lf",&S[i].y)

}

d=ClosePoint(&S[i],n)

printf("最短距离为:%lf\n",d)

}

#include<stdio.h>

#include<stdlib.h>

#include<math.h>

#include<time.h>

typedef struct

{

    double x

    double y

}Building

// 创建所有建筑的坐标

Building* creatCoordinate(int iBuilNum)

{

    int i, j

    double x, y

    srand((int)time(NULL))

    Building* bBuil = (Building*)malloc(sizeof(Building)*iBuilNum)

    for(i=0 i<iBuilNum i++)

    {

        x = (rand()%10000)*0.01

        y = (rand()%10000)*0.01

        for(j=0 j<i j++)

        {

            if(x == bBuil[j].x && y == bBuil[j].y)

            {

                i --

                continue

            }

        }

        bBuil[i].x = x

        bBuil[i].y = y

    }

    return bBuil

}

// 求两点间的距离

double getDis(Building bA, Building bB)

{

    double dDifferenceX = bA.x - bB.x

    double dDifferenceY = bA.y - bB.y

    return sqrt(dDifferenceX*dDifferenceX + dDifferenceY*dDifferenceY)

}

// 求所有点间的距离

double* getAllDis(Building* abCoordinate, int iBuilNum, int iDisNum, int* aiMinIndex)

{

    int i, j, k

    double dMin = 1000.0

    double* adDis = (double*)malloc(sizeof(double)*iDisNum)

    for(i=iBuilNum-1, k=0 i>0 i--)

    {

        for(j=i j>0 j--)

        {

            //printf("%d %d\n", iBuilNum-i-1, iBuilNum-j)

            adDis[k] = getDis(abCoordinate[iBuilNum-i-1], abCoordinate[iBuilNum-j])

            // 截取最小距离

            if(adDis[k] < dMin)

            {

                dMin = adDis[k]

                // 截取最小距离相应坐标数组的下标

                aiMinIndex[0] = iBuilNum-i-1

                aiMinIndex[1] = iBuilNum-j

            }

            k++

        }

    }

    return adDis

}

int main(void)

{

    int iBuilNum, iDisNum = 0, aiMinIndex[2], i

    printf("%s\n", "请输入建筑物数量:")

    scanf("%d", &iBuilNum)

    for(i=1 i<iBuilNum i++)

        iDisNum += i

    Building* abCoordinate = creatCoordinate(iBuilNum)

    double* adDis = getAllDis(abCoordinate, iBuilNum, iDisNum, aiMinIndex)

    // 输出所有坐标

    printf("%s\n", "所有建筑物的坐标为:")

    for(i=0 i<iBuilNum i++)

        printf("x: %0.2lf    y: %0.2lf\n", abCoordinate[i].x, abCoordinate[i].y)

    // 输出距离最近两个点的坐标

    printf("%s\n", "距离最近两个建筑物的坐标为:")

    printf("x: %0.2lf    y: %0.2lf\n", abCoordinate[aiMinIndex[0]].x, abCoordinate[aiMinIndex[0]].y)

    printf("x: %0.2lf    y: %0.2lf\n", abCoordinate[aiMinIndex[1]].x, abCoordinate[aiMinIndex[1]].y)

    // 输出距离最近两个点的距离

    printf("%s\n", "距离最近两个建筑物的距离为:")

    printf("%0.2lf\n", getDis(abCoordinate[aiMinIndex[0]], abCoordinate[aiMinIndex[1]]))

    // 输出所有距离

    //printf("%s\n", "所有建筑物的距离为:")

    //for(i=0 i<iDisNum i++)

    //    printf("%0.2lf\n", adDis[i])

    free(abCoordinate)

    free(adDis)

    return 0

}

1.书写标识符时,忽略了大小写字母的区别。

main()

{

int

a=5

printf("%d",A)

}

编译程序把a和A认为是两个不同的变量名,而显示出错信息。C认为大写字母和小写字母是两个不同的字符。习惯上,符号常量名用大写,变量名用小写表示,以增加可读性。

2.忽略了变量的类型,进行了不合法的运算。

main()

{

float

a,b

printf("%d",a%b)

}

%是求余运算,得到a/b的整余数。整型变量a和b可以进行求余运算,而实型变量则不允许进行“求余”运算。

3.将字符常量与字符串常量混淆。

char

c

c="a"

在这里就混淆了字符常量与字符串常量,字符常量是由一对单引号括起来的单个字符,字符串常量是一对双引号括起来的字符序列。C规定以“\”作字符串结束标志,它是由系统自动加上的,所以字符串“a”实际上包含两个字符:‘a'和‘\',而把它赋给一个字符变量是不行的。

4.忽略了“=”与“==”的区别。

在许多高级语言中,用“=”符号作为关系运算符“等于”。如在BASIC程序中可以写

if

(a=3)

then

但C语言中,“=”是赋值运算符,“==”是关系运算符。如:

if

(a==3)

a=b

前者是进行比较,a是否和3相等,后者表示如果a和3相等,把b值赋给a。由于习惯问题,初学者往往会犯这样的错误。

5.忘记加分号

分号是C语句中不可缺少的一部分,语句末尾必须有分号。

a=1

b=2

编译时,编译程序在“a=1”后面没发现分号,就把下一行“b=2”也作为上一行语句的一部分,这就会出现语法错误。改错时,有时在被指出有错的一行中未发现错误,就需要看一下上一行是否漏掉了分号。

{

z=x+y

t=z/100

printf("%f",t)

}

对于复合语句来说,最后一个语句中最后的分号不能忽略不写(这是和PASCAL不同的)。

6.多加分号。

对于一个复合语句,如:

{

z=x+y

t=z/100

printf("%f",t)

}

复合语句的花括号后不应再加分号,否则将会画蛇添足。

又如:

if

(a%3==0)

I++

本是如果3整除a,则I加1。但由于if

(a%3==0)后多加了分号,则if语句到此结束,程序将执行I++语句,不论3是否整除a,I都将自动加1。

再如:

for

(I=0I<5I++)

{scanf("%d",&x)printf("%d",x)}

本意是先后输入5个数,每输入一个数后再将它输出。由于for()后多加了一个分号,使循环体变为空语句,此时只能输入一个数并输出它。

7.输入变量时忘记加地址运算符“&”。

int

a,b

scanf("%d%d",a,b)

这是不合法的。Scanf函数的作用是:按照a、b在内存的地址将a、b的值存进去。“&a”指a在内存中的地址。

8.输入数据的方式与要求不符。①scanf("%d%d",&a,&b)

输入时,不能用逗号作两个数据间的分隔符,如下面输入不合法:

3,4

输入数据时,在两个数据之间以一个或多个空格间隔,也可用回车键,跳格键tab。

②scanf("%d,%d",&a,&b)

C规定:如果在“格式控制”字符串中除了格式说明以外还有其它字符,则在输入数据时应输入与这些字符相同的字符。下面输入是合法的:

3,4

此时不用逗号而用空格或其它字符是不对的。

3

4

3:4

又如:

scanf("a=%d,b=%d",&a,&b)

输入应如以下形式:

a=3,b=4

9.输入字符的格式与要求不一致。

在用“%c”格式输入字符时,“空格字符”和“转义字符”都作为有效字符输入。

scanf("%c%c%c",&c1,&c2,&c3)

如输入a

b

c

字符“a”送给c1,字符“

”送给c2,字符“b”送给c3,因为%c只要求读入一个字符,后面不需要用空格作为两个字符的间隔。

10.输入输出的数据类型与所用格式说明符不一致。

例如,a已定义为整型,b定义为实型

a=3b=4.5

printf("%f%d\n",a,b)

编译时不给出出错信息,但运行结果将与原意不符。这种错误尤其需要注意。