#226063 (D) No.622 点と三角柱の内外判定

提出ソース

結果

問題	No.622 点と三角柱の内外判定
ユーザー	ikd
提出日時	2017-12-25 16:49:59
言語	D (dmd 2.106.1)
結果	AC
実行時間	1 ms / 1,500 ms
コード長	2,161 bytes
コンパイル時間	763 ms
コンパイル使用メモリ	107,696 KB
実行使用メモリ	6,944 KB
最終ジャッジ日時	2024-06-12 23:16:18
合計ジャッジ時間	1,523 ms
ジャッジサーバーID （参考情報）	judge2 / judge4

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	1 ms 6,812 KB
testcase_01	AC	1 ms 6,940 KB
testcase_02	AC	1 ms 6,940 KB
testcase_03	AC	1 ms 6,944 KB
testcase_04	AC	1 ms 6,944 KB
testcase_05	AC	1 ms 6,940 KB
testcase_06	AC	1 ms 6,944 KB
testcase_07	AC	1 ms 6,940 KB
testcase_08	AC	1 ms 6,944 KB
testcase_09	AC	1 ms 6,944 KB
testcase_10	AC	1 ms 6,940 KB
testcase_11	AC	1 ms 6,944 KB
testcase_12	AC	1 ms 6,940 KB
testcase_13	AC	1 ms 6,944 KB
testcase_14	AC	1 ms 6,940 KB
testcase_15	AC	1 ms 6,944 KB
testcase_16	AC	1 ms 6,940 KB
testcase_17	AC	1 ms 6,940 KB
testcase_18	AC	1 ms 6,940 KB
testcase_19	AC	1 ms 6,940 KB
testcase_20	AC	1 ms 6,944 KB
testcase_21	AC	1 ms 6,940 KB
testcase_22	AC	1 ms 6,940 KB
testcase_23	AC	1 ms 6,944 KB
testcase_24	AC	1 ms 6,944 KB
testcase_25	AC	1 ms 6,944 KB
testcase_26	AC	1 ms 6,940 KB
testcase_27	AC	1 ms 6,940 KB
testcase_28	AC	1 ms 6,940 KB
testcase_29	AC	1 ms 6,940 KB
testcase_30	AC	1 ms 6,944 KB
testcase_31	AC	1 ms 6,940 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

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();
}

yukicoder

結果

テストケース

ソースコード