【重口味。。】CTSC2011杀菌计划

在奋战了7个小时后才AC了这道题的人只能瞻仰在考场上AC的FHQ了。。。

代码很恐怖，写了350行8KB，最可恶的是一页草稿纸的伪代码和各种特殊情况。

 

涉及的几个关键字：三维叉积、半平面交、梯形剖分

 

三维叉积的一个核心操作就是：三维旋系，另外还需要用来写：面面交，面线交，判三维空间上下左右，二维投影。

 

具体做法是：

      一束光线首先在每一个面上投一个影，然后这个投影需要和这个面求一个交集（这里我用的是半平面交），之后再反射递归。

 

      计算答案时将每个面上曾经经过的所有投影求一个面积并（梯形剖分）。

 

当然细节很多，不过还好栋哥没有出细节太多的数据，本来我4.30小时过了样例之后直接去测就已经过了7个点，但令我哭笑不得的是第6个点（全部面积之和）和第7个点（0）错了。。。之后又发现疏漏了几种特殊情况。。

 

很郁闷的是半平面交又错了老地方：最开始的时候忘记判平行和判无解

 

扣的：

program ctscbox; uses math; const e=1e-7; type point=record x,y,z:real end; jid=record x,y,z:point end; side=record n:longint; o:array[0..20] of point; d:jid; end; fs=record g:point; p:side; end; cm=object n:longint; m:array[0..20] of side; function cal:real; procedure new(a:side); end; var dd:array[0..100] of point; q:array[0..100] of cm; a:array[0..100] of side; m,n,i,j,k,lim:longint; ans:real; f:fs; operator -(a,b:point)c:point;begin c.x:=a.x-b.x;c.y:=a.y-b.y;c.z:=a.z-b.z end; operator +(a,b:point)c:point;begin c.x:=a.x+b.x;c.y:=a.y+b.y;c.z:=a.z+b.z end; operator /(a,b:point)c:real;begin c:=a.x*b.x+a.y*b.y+a.z*b.z end; function len(a:point):real;begin with a do len:=sqrt(x*x+y*y+z*z) end; operator /(a:point;b:real)c:point;begin c.x:=a.x/b;c.y:=a.y/b;c.z:=a.z/b end; operator *(a:point;b:real)c:point;begin c.x:=a.x*b;c.y:=a.y*b;c.z:=a.z*b end; procedure cm.new(a:side);begin inc(n);m[n]:=a end; operator *(a,b:point)c:point; begin c.x:=a.y*b.z-a.z*b.y; c.y:=a.z*b.x-a.x*b.z; c.z:=a