結果

問題 No.132 点と平面との距離
コンテスト
ユーザー kuuso1
提出日時 2015-01-21 00:12:31
言語 C#(csc)
(csc 3.9.0)
結果
AC  
実行時間 546 ms / 5,000 ms
コード長 1,992 bytes
コンパイル時間 2,228 ms
コンパイル使用メモリ 107,904 KB
実行使用メモリ 19,456 KB
最終ジャッジ日時 2024-06-22 23:27:43
合計ジャッジ時間 3,060 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 3
権限があれば一括ダウンロードができます
コンパイルメッセージ
Microsoft (R) Visual C# Compiler version 3.9.0-6.21124.20 (db94f4cc)
Copyright (C) Microsoft Corporation. All rights reserved.

ソースコード

diff # 

using System;
using System.Collections;
using System.Collections.Generic;
 
class TEST{
	static void Main(){
		Sol mySol =new Sol();
		mySol.Solve();
	}
}

class Sol{
	public void Solve(){
		
		double ans=0.0;
		for(int i=0;i<N;i++){
			for(int j=i+1;j<N;j++){
				for(int k=j+1;k<N;k++){
					double V=vol(Q[i],Q[j],Q[k],P);
					double Ar=area(Q[i],Q[j],Q[k]);
//Console.WriteLine(V);Console.WriteLine(Ar);
					ans+=V/Ar*3.0;
				}
			}
		}
		Console.WriteLine(ans);
		
	}
	
	double vol(double[] A,double[] B,double[] C,double[] O){
		double ret= (A[0]-O[0])*(B[1]-O[1])*(C[2]-O[2])
			   +(B[0]-O[0])*(C[1]-O[1])*(A[2]-O[2])
			   +(C[0]-O[0])*(A[1]-O[1])*(B[2]-O[2])
			   -(A[0]-O[0])*(C[1]-O[1])*(B[2]-O[2])
			   -(B[0]-O[0])*(A[1]-O[1])*(C[2]-O[2])
			   -(C[0]-O[0])*(B[1]-O[1])*(A[2]-O[2]);
		ret=Math.Abs(ret);
		return ret/6;
	}
	
	double area(double[] A,double[] B,double[] O){
		double ip=(A[0]-O[0])*(B[0]-O[0])+(A[1]-O[1])*(B[1]-O[1])+(A[2]-O[2])*(B[2]-O[2]);
		double ra=(A[0]-O[0])*(A[0]-O[0])+(A[1]-O[1])*(A[1]-O[1])+(A[2]-O[2])*(A[2]-O[2]);
		double rb=(B[0]-O[0])*(B[0]-O[0])+(B[1]-O[1])*(B[1]-O[1])+(B[2]-O[2])*(B[2]-O[2]);
		double theta=Math.Acos(Math.Abs(ip)/(Math.Sqrt(ra)*Math.Sqrt(rb)));
		return Math.Sqrt(ra)*Math.Sqrt(rb)*Math.Sin(theta)/2;
	}
	
	
	
	
	int N;
	double[] P;
	double[][] Q;
	public Sol(){
		N=ri();
		P=rda();
		Q=new double[N][];
		for(int i=0;i<N;i++)Q[i]=rda();
	}



	static String rs(){return Console.ReadLine();}
	static int ri(){return int.Parse(Console.ReadLine());}
	static long rl(){return long.Parse(Console.ReadLine());}
	static double rd(){return double.Parse(Console.ReadLine());}
	static String[] rsa(){return Console.ReadLine().Split(' ');}
	static int[] ria(){return Array.ConvertAll(Console.ReadLine().Split(' '),e=>int.Parse(e));}
	static long[] rla(){return Array.ConvertAll(Console.ReadLine().Split(' '),e=>long.Parse(e));}
	static double[] rda(){return Array.ConvertAll(Console.ReadLine().Split(' '),e=>double.Parse(e));}
}
0