C语言求解同一条直线所能通过的最多点数

题目

解答

	#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;
	}