結果

問題 No.132 点と平面との距離
ユーザー kuuso1kuuso1
提出日時 2015-01-21 00:12:31
言語 C#(csc)
(csc 3.9.0)
結果
AC  
実行時間 658 ms / 5,000 ms
コード長 1,992 bytes
コンパイル時間 5,982 ms
コンパイル使用メモリ 105,108 KB
実行使用メモリ 23,268 KB
最終ジャッジ日時 2023-09-05 02:38:48
合計ジャッジ時間 6,911 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge14
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 91 ms
23,268 KB
testcase_01 AC 243 ms
23,132 KB
testcase_02 AC 658 ms
23,032 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
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