結果

問題 No.132 点と平面との距離
コンテスト
ユーザー akakimidori
提出日時 2016-12-12 20:14:02
言語 C90
(gcc 15.2.0)
コンパイル:
gcc-15 -O2 -std=c90 -DONLINE_JUDGE -o a.out _filename_ -lm
実行:
./a.out
結果
AC  
実行時間 38 ms / 5,000 ms
コード長 972 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 277 ms
コンパイル使用メモリ 40,948 KB
最終ジャッジ日時 2026-02-23 23:47:26
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 3
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<stdio.h>
#include<stdlib.h>
#include<math.h>

#define POS(i,j) ((i)*3+(j))

void run(void){
  int n;
  scanf("%d",&n);

  double p[3];
  scanf("%lf%lf%lf",p,p+1,p+2);

  double *q=(double *)malloc(sizeof(double)*3*n);
  int i;
  for(i=0;i<n;i++){
    int k;
    for(k=0;k<3;k++){
      scanf("%lf",q+i*3+k);
      q[POS(i,k)]-=p[k];
    }
  }

  double sum=0;
  int j,k,l;
  for(i=0;i<n;i++){
    double I[3];
    for(l=0;l<3;l++){
      I[l]=q[POS(i,l)];
    }
    for(j=i+1;j<n;j++){
      double J[3];
      for(l=0;l<3;l++){
	J[l]=q[POS(j,l)]-I[l];
      }
      for(k=j+1;k<n;k++){
	double K[3];
	for(l=0;l<3;l++){
	  K[l]=q[POS(k,l)]-I[l];
	}
	double t[3];
	for(l=0;l<3;l++){
	  t[l]=J[(l+1)%3]*K[(l+2)%3]-J[(l+2)%3]*K[(l+1)%3];
	}
	double nu=0;
	double de=0;
	for(l=0;l<3;l++){
	  nu+=t[l]*I[l];
	  de+=t[l]*t[l];
	}
	nu=nu>=0?nu:-nu;
	de=sqrt(de);
	sum+=nu/de;
      }
    }
  }
  printf("%.9lf\n",sum);
  return;
}

int main(void){
  run();
  return 0;
}
0