202009-4星际旅行 ccf java 90分

202009-4星际旅行 ccf java 90分

先定义类（坐标点）：

class Coordinate{
	double co;//点和o距离
	int [] dis;//各个维度坐标
	Coordinate(int dis[],int o[],int n){
		this.dis = dis;
		double temp = 0;
		for(int i=0;i<n;i++){
			temp += (dis[i]-o[i])*(dis[i]-o[i])*1.0;
		}
		this.co = Math.sqrt(temp*1.0);
	}	
}

处理两点间最短距离的方法：

public static double getDis(Coordinate a,Coordinate b){
	//如果o到ab距离<=r 则直接就是弧长
	double ab = 0;
	for(int i=0;i<n;i++){
		ab += (a.dis[i]-b.dis[i])*(a.dis[i]-b.dis[i])*1.0;
	}
	ab = Math.sqrt(ab*1.0);
	
	
	double cosa = (a.co*a.co+ab*ab-b.co*b.co)/(2*a.co*ab);
	double cosb = (b.co*b.co+ab*ab-a.co*a.co)/(2*b.co*ab);
	
	if(cosa<=0 || cosb<=0){//a或b为钝角或直角 则ab在黑洞外
		return ab;
	}else{
		double p = (ab+a.co+b.co)*1.0/2*1.0;
		double s = Math.sqrt(p*(p-a.co)*(p-b.co)*(p-ab));
		double h = 2*s/ab;
		
		if(h<r){//在内部 返回 ae+ef+fb 其中ef为切点
			double ae = Math.sqrt(a.co*a.co-r*r*1.0);
			double fb = Math.sqrt(b.co*b.co-r*r*1.0);
			double coso = (a.co*a.co+b.co*b.co-ab*ab)/(2*a.co*b.co*1.0);
			double cosaoe = r*1.0/a.co;
			double cosfob = r*1.0/b.co;
			double arfa = Math.acos(coso)-Math.acos(cosaoe)-Math.acos(cosfob); 
			double ef = r*arfa;
			return ae+fb+ef;
		}
		return ab;
	}

剩下的看注释就好啦
感觉计算距离那个地方完全就是数学推导

import java.util.*;

public class _4 {
public static int n,m,r;
static int[] o;//储存o点坐标
static Coordinate ms[];//储存所有旅行点坐标
static double[] diss;//计算第i个点和其他点之间的最短距离


public static void main (String args[]){
	Scanner sc = new Scanner(System.in);
	n = sc.nextInt();m  = sc.nextInt();r = sc.nextInt();
	o = new int[n];ms = new Coordinate[m];diss = new double[m];
	for(int i=0;i<n;i++){//o点坐标		
		o[i] =sc.nextInt();
	}
	for(int i=0;i<m;i++){
		int[] temp = new int[m];
		for(int j=0;j<n;j++){			
			temp[j] = sc.nextInt();
		}	
		ms[i] = new Coordinate(temp,o,n);
	}//输入处理完成
	//需要计算每两点之间
	for(int i=0;i<m;i++){
		for(int j=i+1;j<m;j++){
			double d = getDis(ms[i],ms[j]);
			diss[i] += d;
			diss[j] += d;
		}		
	}
	for(int i=0;i<m;i++){
		System.out.println(String.format("%.12f", diss[i]));
	}
}

}
}

不知道为啥超时了 还没想好不超时的算法
90分已经满足了（不是）