結果

問題 No.622 点と三角柱の内外判定
ユーザー ikdikd
提出日時 2017-12-25 16:48:00
言語 D
(dmd 2.105.2)
結果
AC  
実行時間 2 ms / 1,500 ms
コード長 2,164 bytes
コンパイル時間 1,983 ms
コンパイル使用メモリ 96,368 KB
実行使用メモリ 4,384 KB
最終ジャッジ日時 2023-09-03 17:48:05
合計ジャッジ時間 2,479 ms
ジャッジサーバーID
(参考情報)
judge13 / judge15
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,380 KB
testcase_01 AC 1 ms
4,380 KB
testcase_02 AC 2 ms
4,376 KB
testcase_03 AC 2 ms
4,376 KB
testcase_04 AC 1 ms
4,380 KB
testcase_05 AC 2 ms
4,376 KB
testcase_06 AC 2 ms
4,380 KB
testcase_07 AC 1 ms
4,376 KB
testcase_08 AC 1 ms
4,380 KB
testcase_09 AC 1 ms
4,376 KB
testcase_10 AC 1 ms
4,376 KB
testcase_11 AC 1 ms
4,380 KB
testcase_12 AC 2 ms
4,376 KB
testcase_13 AC 1 ms
4,376 KB
testcase_14 AC 1 ms
4,376 KB
testcase_15 AC 1 ms
4,380 KB
testcase_16 AC 2 ms
4,380 KB
testcase_17 AC 1 ms
4,380 KB
testcase_18 AC 2 ms
4,384 KB
testcase_19 AC 1 ms
4,380 KB
testcase_20 AC 2 ms
4,376 KB
testcase_21 AC 1 ms
4,376 KB
testcase_22 AC 2 ms
4,376 KB
testcase_23 AC 1 ms
4,380 KB
testcase_24 AC 1 ms
4,376 KB
testcase_25 AC 1 ms
4,380 KB
testcase_26 AC 1 ms
4,376 KB
testcase_27 AC 1 ms
4,380 KB
testcase_28 AC 1 ms
4,380 KB
testcase_29 AC 1 ms
4,380 KB
testcase_30 AC 1 ms
4,376 KB
testcase_31 AC 1 ms
4,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

void main(){
  import std.stdio, std.string, std.conv, std.algorithm;
  import std.typecons, std.math;

  double xx, yy, zz; rd(xx, yy, zz);
  alias P=Tuple!(double, "x", double, "y", double, "z");
  P va, vb, vc;
  rd(va.x, va.y, va.z); va.x-=xx; va.y-=yy; va.z-=zz;
  rd(vb.x, vb.y, vb.z); vb.x-=xx; vb.y-=yy; vb.z-=zz;
  rd(vc.x, vc.y, vc.z); vc.x-=xx; vc.y-=yy; vc.z-=zz;

  auto vn=P(va.y*vb.z-va.z*vb.y, va.z*vb.x-va.x*vb.z, va.x*vb.y-va.y*vb.x);
  if(sqrt(vn.x*vn.x+vn.y*vn.y+vn.z*vn.z)<1e-9){writeln("NO"); return;}
  auto T=new double[][](3, 3);
  T[0][0]=va.x; T[0][1]=vb.x; T[0][2]=vn.x;
  T[1][0]=va.y; T[1][1]=vb.y; T[1][2]=vn.y;
  T[2][0]=va.z; T[2][1]=vb.z; T[2][2]=vn.z;
  auto invT=invMatrix(T);
  auto x=invT[0][0]*vc.x+invT[0][1]*vc.y+invT[0][2]*vc.z,
       y=invT[1][0]*vc.x+invT[1][1]*vc.y+invT[1][2]*vc.z;
  
  if(x>=0 && y>=0 && x+y<=1) writeln("YES");
  else writeln("NO");
}

double[][] invMatrix(const double[][] A){
  auto invA=new double[][](3, 3);
  auto det=A[0][0]*A[1][1]*A[2][2]+A[0][1]*A[1][2]*A[2][0]+A[1][0]*A[2][1]*A[0][2]
          -A[0][0]*A[1][2]*A[2][1]-A[0][1]*A[1][0]*A[2][2]-A[0][2]*A[1][1]*A[2][0];
  invA[0][0]=(A[1][1]*A[2][2]-A[1][2]*A[2][1])/det;
  invA[0][1]=(A[0][2]*A[2][1]-A[0][1]*A[2][2])/det;
  invA[0][2]=(A[0][1]*A[1][2]-A[0][2]*A[1][1])/det;
  invA[1][0]=(A[1][2]*A[2][0]-A[1][0]*A[2][2])/det;
  invA[1][1]=(A[0][0]*A[2][2]-A[0][2]*A[2][0])/det;
  invA[1][2]=(A[0][2]*A[1][0]-A[0][0]*A[1][2])/det;
  invA[2][0]=(A[1][0]*A[2][1]-A[1][1]*A[2][0])/det;
  invA[2][1]=(A[0][1]*A[2][0]-A[0][0]*A[2][1])/det;
  invA[2][2]=(A[0][0]*A[1][1]-A[0][1]*A[1][0])/det;
  return invA;
}

unittest{
  import std.stdio;
  auto A=new double[][](3, 3);
  A[0][0]=1; A[0][1]=2; A[0][2]=1;
  A[1][0]=2; A[1][1]=1; A[1][2]=0;
  A[2][0]=1; A[2][1]=1; A[2][2]=2;
  auto invA=invMatrix(A);
  writeln(invA);
  // -0.4  0.6  0.2
  //  0.8 -0.2 -0.4
  // -0.2 -0.2  0.6
}

void rd(T...)(ref T x){
  import std.stdio, std.string, std.conv;
  auto l=readln.split;
  assert(l.length==x.length);
  foreach(i, ref e; x){
    e=l[i].to!(typeof(e));
  }
}
void wr(T...)(T x){
  import std.stdio;
  foreach(e; x) write(e, " ");
  writeln();
}
0