题目
解答
#include <stdio.h>
#include <stdlib.h> //提供exit()原型
#define maxSize 1000 //最大点个数
typedef struct{
float x;
float y;
}Point;
int main(){
int i, j, k, m=0;
int count, max, maxi, maxj;
FILE* fp;
Point pset[maxSize]; //点集
/*读取数据*/
if((fp = fopen("2007-06.txt", "r")) == NULL)
exit(0);
while(fscanf(fp, "%d%d", &i, &j) != EOF){
pset[m].x = i;
pset[m].y = j;
m++; //m代表点的个数
}
max = 0; //先将最大值置为0
/*计算任意两点构成的直线通过的点数,并返回最大值*/
/*用i和j代表任意两点在pset中的位置*/
for(i = 0;i<m;i++){
for(j = 0; j < m; j++){
/*通过上面那两重循环,可以使任意两点相互交叉组合形成许多条线段*/
/*【假设有N个点,则行成的直线个数为N*N-N,这里把P1P2和P2P1看成了两条线段】*/
if(i!=j){ //排除掉自己跟自己组合的情况【同一个点不构成线段】
count = 0; //计数器置为0
/*第三层循环遍历所有的点,判断该点是否在线段PiPj上*/
for(k = 0; k < m; k++){
if(k != i && k != j){ //排除掉该点在直线PiPj两个端点的情况【Pi或Pj】
/*线段PkPj与线段PiPj无论什么情况下都不可能平行【因为有个共同的端点Pj】*/
/*若点Pk在线段PiPj的延长线上,则线段PkPj与线段PiPj共线*/
/*若向量Pi-Pj×向量Pk-Pj=0,则可以证明线段PiPj与线段PkPj平行或共线,但二者不可能平行→→推出二者共线*/
if((pset[j].x - pset[k].x)*(pset[j].y - pset[i].y) == (pset[j].y - pset[k].y)*(pset[j].x - pset[i].x))
count++; //若果点Pk在线段PiPj上,点的个数++
}
}
/*若当前线段上点数大,就更新最大点个数【不包括当前线段两端的端点】以及该线段两端端点对应的下标*/
if(count>max){
max = count;
maxi = i;
maxj = j;
}
}
}
}
/*循环结束后最终max记录的是最大点个数【不包括线段PiPj两端的端点】*/
/*maxi和maxj记录的是线段PiPj两端的端点在pset中的下标【该线段所在直线可以通过最大的点数】*/
/*max要加上2,因为还要加上线段两端的端点*/
printf("同一条直线所能通过的最多点数为:%d\n", max+2);
/*找出该直线上的点,并打印这些点*/
for(k = 0; k < m; ++k){
/*若向量Pmaxi-Pmaxj×向量Pk-Pmaxj=0,则可以证明线段PmaxiPmaxj与线段PkPmaxj平行或共线,但二者不可能平行→→推出二者共线*/
/*如果k==maxi或maxj,则直线PkPmaxi或直线PkPmaxj都是点,但同样符合下面的判断条件*/
if((pset[maxj].x - pset[k].x)*(pset[maxj].y - pset[maxi].y) == (pset[maxj].y - pset[k].y)*(pset[maxj].x - pset[maxi].x))
printf("第%d个点(%.2f,%.2f)\n", k+1,pset[k].x,pset[k].y);
}
return 0;
}