Hex Artifact Content
Not logged in

Artifact 27e70b189824397209fadfababb1007c82345954:


0000: 23 69 6e 63 6c 75 64 65 20 3c 69 6f 73 74 72 65  #include <iostre
0010: 61 6d 3e 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 73  am>..#include <s
0020: 73 74 72 65 61 6d 3e 0d 0a 23 69 6e 63 6c 75 64  stream>..#includ
0030: 65 20 3c 69 6f 6d 61 6e 69 70 3e 0d 0a 23 69 6e  e <iomanip>..#in
0040: 63 6c 75 64 65 20 3c 76 65 63 74 6f 72 3e 0d 0a  clude <vector>..
0050: 23 69 6e 63 6c 75 64 65 20 3c 73 74 72 69 6e 67  #include <string
0060: 3e 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 6d 61 70  >..#include <map
0070: 3e 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 73 65 74  >..#include <set
0080: 3e 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 61 6c 67  >..#include <alg
0090: 6f 72 69 74 68 6d 3e 0d 0a 23 69 6e 63 6c 75 64  orithm>..#includ
00a0: 65 20 3c 6e 75 6d 65 72 69 63 3e 0d 0a 23 69 6e  e <numeric>..#in
00b0: 63 6c 75 64 65 20 3c 69 74 65 72 61 74 6f 72 3e  clude <iterator>
00c0: 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 66 75 6e 63  ..#include <func
00d0: 74 69 6f 6e 61 6c 3e 0d 0a 23 69 6e 63 6c 75 64  tional>..#includ
00e0: 65 20 3c 63 6f 6d 70 6c 65 78 3e 0d 0a 23 69 6e  e <complex>..#in
00f0: 63 6c 75 64 65 20 3c 71 75 65 75 65 3e 0d 0a 23  clude <queue>..#
0100: 69 6e 63 6c 75 64 65 20 3c 73 74 61 63 6b 3e 0d  include <stack>.
0110: 0a 23 69 6e 63 6c 75 64 65 20 3c 63 6d 61 74 68  .#include <cmath
0120: 3e 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 63 61 73  >..#include <cas
0130: 73 65 72 74 3e 0d 0a 23 69 6e 63 6c 75 64 65 20  sert>..#include 
0140: 3c 63 73 74 72 69 6e 67 3e 0d 0a 75 73 69 6e 67  <cstring>..using
0150: 20 6e 61 6d 65 73 70 61 63 65 20 73 74 64 3b 0d   namespace std;.
0160: 0a 74 79 70 65 64 65 66 20 6c 6f 6e 67 20 6c 6f  .typedef long lo
0170: 6e 67 20 4c 4c 3b 0d 0a 74 79 70 65 64 65 66 20  ng LL;..typedef 
0180: 63 6f 6d 70 6c 65 78 3c 64 6f 75 62 6c 65 3e 20  complex<double> 
0190: 43 4d 50 3b 0d 0a 0d 0a 63 6c 61 73 73 20 44 65  CMP;....class De
01a0: 70 6f 73 69 74 20 7b 20 70 75 62 6c 69 63 3a 0d  posit { public:.
01b0: 0a 09 76 65 63 74 6f 72 3c 43 4d 50 3e 20 6d 61  ..vector<CMP> ma
01c0: 6b 65 5f 70 6f 6c 79 67 6f 6e 28 63 6f 6e 73 74  ke_polygon(const
01d0: 20 76 65 63 74 6f 72 3c 69 6e 74 3e 26 20 78 73   vector<int>& xs
01e0: 2c 20 63 6f 6e 73 74 20 76 65 63 74 6f 72 3c 69  , const vector<i
01f0: 6e 74 3e 26 20 79 73 29 0d 0a 09 7b 0d 0a 09 09  nt>& ys)...{....
0200: 76 65 63 74 6f 72 3c 43 4d 50 3e 20 76 73 3b 0d  vector<CMP> vs;.
0210: 0a 09 09 66 6f 72 28 69 6e 74 20 69 3d 30 3b 20  ...for(int i=0; 
0220: 69 3c 78 73 2e 73 69 7a 65 28 29 3b 20 2b 2b 69  i<xs.size(); ++i
0230: 29 0d 0a 09 09 09 76 73 2e 70 75 73 68 5f 62 61  ).....vs.push_ba
0240: 63 6b 28 20 43 4d 50 28 78 73 5b 69 5d 2c 20 79  ck( CMP(xs[i], y
0250: 73 5b 69 5d 29 20 29 3b 0d 0a 09 09 72 65 74 75  s[i]) );....retu
0260: 72 6e 20 76 73 3b 0d 0a 09 7d 0d 0a 0d 0a 09 64  rn vs;...}.....d
0270: 6f 75 62 6c 65 20 73 75 63 63 65 73 73 50 72 6f  ouble successPro
0280: 62 61 62 69 6c 69 74 79 28 76 65 63 74 6f 72 20  bability(vector 
0290: 3c 69 6e 74 3e 20 73 69 74 65 58 2c 20 76 65 63  <int> siteX, vec
02a0: 74 6f 72 20 3c 69 6e 74 3e 20 73 69 74 65 59 2c  tor <int> siteY,
02b0: 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 64 65   vector <int> de
02c0: 70 6f 73 69 74 58 2c 20 76 65 63 74 6f 72 20 3c  positX, vector <
02d0: 69 6e 74 3e 20 64 65 70 6f 73 69 74 59 29 20 0d  int> depositY) .
02e0: 0a 09 7b 0d 0a 09 09 63 6f 6e 73 74 20 76 65 63  ..{....const vec
02f0: 74 6f 72 3c 43 4d 50 3e 20 73 73 20 3d 20 6d 61  tor<CMP> ss = ma
0300: 6b 65 5f 70 6f 6c 79 67 6f 6e 28 73 69 74 65 58  ke_polygon(siteX
0310: 2c 20 73 69 74 65 59 29 3b 0d 0a 09 09 63 6f 6e  , siteY);....con
0320: 73 74 20 76 65 63 74 6f 72 3c 43 4d 50 3e 20 64  st vector<CMP> d
0330: 64 20 3d 20 6d 61 6b 65 5f 70 6f 6c 79 67 6f 6e  d = make_polygon
0340: 28 64 65 70 6f 73 69 74 58 2c 20 64 65 70 6f 73  (depositX, depos
0350: 69 74 59 29 3b 0d 0a 09 09 63 6f 6e 73 74 20 69  itY);....const i
0360: 6e 74 20 53 4e 20 3d 20 73 73 2e 73 69 7a 65 28  nt SN = ss.size(
0370: 29 3b 0d 0a 0d 0a 09 09 64 6f 75 62 6c 65 20 74  );......double t
0380: 6f 74 61 6c 6c 65 6e 3d 30 2c 20 68 69 74 6c 65  otallen=0, hitle
0390: 6e 3d 30 3b 0d 0a 09 09 66 6f 72 28 69 6e 74 20  n=0;....for(int 
03a0: 69 3d 30 3b 20 69 3c 53 4e 3b 20 2b 2b 69 29 0d  i=0; i<SN; ++i).
03b0: 0a 09 09 7b 0d 0a 09 09 09 63 6f 6e 73 74 20 43  ...{.....const C
03c0: 4d 50 20 70 30 20 3d 20 73 73 5b 69 5d 3b 0d 0a  MP p0 = ss[i];..
03d0: 09 09 09 63 6f 6e 73 74 20 43 4d 50 20 70 31 20  ...const CMP p1 
03e0: 3d 20 73 73 5b 28 69 2b 31 29 25 53 4e 5d 3b 0d  = ss[(i+1)%SN];.
03f0: 0a 09 09 09 74 6f 74 61 6c 6c 65 6e 20 2b 3d 20  ....totallen += 
0400: 61 62 73 28 70 31 2d 70 30 29 3b 0d 0a 0d 0a 0d  abs(p1-p0);.....
0410: 0a 09 09 09 43 4d 50 20 70 76 20 3d 20 70 31 20  ....CMP pv = p1 
0420: 2d 20 70 30 3b 20 70 76 20 2f 3d 20 61 62 73 28  - p0; pv /= abs(
0430: 70 76 29 3b 0d 0a 0d 0a 09 09 09 2f 2f 20 63 61  pv);.......// ca
0440: 6c 63 20 66 75 72 74 68 65 73 74 20 6f 6e 20 70  lc furthest on p
0450: 30 20 2d 20 70 31 0d 0a 09 09 09 64 6f 75 62 6c  0 - p1.....doubl
0460: 65 20 71 20 3d 20 31 2e 30 3b 0d 0a 09 09 09 66  e q = 1.0;.....f
0470: 6f 72 28 69 6e 74 20 6b 3d 69 2b 31 3b 20 6b 21  or(int k=i+1; k!
0480: 3d 69 3b 20 6b 3d 28 6b 2b 31 29 25 53 4e 29 0d  =i; k=(k+1)%SN).
0490: 0a 09 09 09 7b 0d 0a 09 09 09 09 63 6f 6e 73 74  ....{......const
04a0: 20 43 4d 50 20 71 30 20 3d 20 73 73 5b 6b 5d 3b   CMP q0 = ss[k];
04b0: 0d 0a 09 09 09 09 63 6f 6e 73 74 20 43 4d 50 20  ......const CMP 
04c0: 71 31 20 3d 20 73 73 5b 28 6b 2b 31 29 25 53 4e  q1 = ss[(k+1)%SN
04d0: 5d 3b 0d 0a 09 09 09 09 2f 2f 20 63 6f 6d 70 75  ];......// compu
04e0: 74 65 20 76 65 72 74 69 63 61 6c 2d 65 71 73 70  te vertical-eqsp
04f0: 6c 69 74 74 65 72 20 6f 66 20 6c 69 6e 65 73 65  litter of linese
0500: 67 20 71 30 20 2d 20 71 31 0d 0a 09 09 09 09 43  g q0 - q1......C
0510: 4d 50 20 6d 69 64 20 3d 20 28 71 30 2b 71 31 29  MP mid = (q0+q1)
0520: 2f 32 2e 30 3b 0d 0a 09 09 09 09 43 4d 50 20 76  /2.0;......CMP v
0530: 65 63 20 3d 20 28 71 30 2d 71 31 29 2a 43 4d 50  ec = (q0-q1)*CMP
0540: 28 30 2c 31 29 3b 3b 0d 0a 09 09 09 09 2f 2f 20  (0,1);;......// 
0550: 63 6f 6d 70 75 74 65 20 63 72 6f 73 73 69 6e 67  compute crossing
0560: 20 77 69 74 68 20 70 30 20 2d 20 70 31 0d 0a 09   with p0 - p1...
0570: 09 09 09 64 6f 75 62 6c 65 20 71 71 20 3d 20 73  ...double qq = s
0580: 6f 6c 76 65 28 70 30 2c 70 76 2c 6d 69 64 2c 76  olve(p0,pv,mid,v
0590: 65 63 29 3b 0d 0a 09 09 09 09 69 66 28 20 71 71  ec);......if( qq
05a0: 20 3c 20 71 20 29 0d 0a 09 09 09 09 7b 0d 0a 09   < q )......{...
05b0: 09 09 09 09 69 66 28 20 71 71 20 3c 3d 20 30 2e  ....if( qq <= 0.
05c0: 30 20 29 0d 0a 09 09 09 09 09 09 71 71 20 3d 20  0 )........qq = 
05d0: 30 2e 30 3b 0d 0a 09 09 09 09 09 2f 2f 20 69 6e  0.0;.......// in
05e0: 20 28 71 71 2c 71 29 2c 20 66 75 72 74 68 65 73   (qq,q), furthes
05f0: 74 20 69 73 20 71 31 0d 0a 09 09 09 09 09 2f 2f  t is q1.......//
0600: 20 5b 74 6f 64 6f 5d 20 63 6f 6d 70 75 74 65 20   [todo] compute 
0610: 74 68 65 20 66 72 61 63 74 69 6f 6e 20 6f 66 20  the fraction of 
0620: 74 68 65 20 72 65 67 69 6f 6e 20 63 72 6f 73 73  the region cross
0630: 69 6e 67 20 77 69 74 68 20 74 68 65 20 64 65 70  ing with the dep
0640: 6f 73 69 74 0d 0a 09 09 09 09 09 2f 2f 20 62 79  osit.......// by
0650: 20 74 65 72 6e 65 72 79 20 73 65 61 72 63 68 20   ternery search 
0660: 6f 72 20 73 6f 6d 65 74 68 69 6e 67 20 6c 69 6b  or something lik
0670: 65 20 74 68 61 74 0d 0a 09 09 09 09 09 64 6f 75  e that.......dou
0680: 62 6c 65 20 72 20 3d 20 63 6f 6d 70 75 74 65 46  ble r = computeF
0690: 72 61 63 74 69 6f 6e 4f 66 43 72 6f 73 73 28 70  ractionOfCross(p
06a0: 30 2b 71 71 2a 70 76 2c 20 70 30 2b 71 2a 70 76  0+qq*pv, p0+q*pv
06b0: 2c 20 71 31 2c 20 64 64 29 3b 0d 0a 0d 0a 09 09  , q1, dd);......
06c0: 09 09 09 68 69 74 6c 65 6e 20 2b 3d 20 61 62 73  ...hitlen += abs
06d0: 28 70 30 2d 70 31 29 2a 28 71 2d 71 71 29 2a 72  (p0-p1)*(q-qq)*r
06e0: 3b 0d 0a 09 09 09 09 09 71 20 3d 20 71 71 3b 0d  ;.......q = qq;.
06f0: 0a 09 09 09 09 7d 0d 0a 09 09 09 09 69 66 28 20  .....}......if( 
0700: 71 20 3d 3d 20 30 2e 30 20 29 0d 0a 09 09 09 09  q == 0.0 )......
0710: 09 62 72 65 61 6b 3b 0d 0a 09 09 09 7d 0d 0a 09  .break;.....}...
0720: 09 7d 0d 0a 09 09 72 65 74 75 72 6e 20 68 69 74  .}....return hit
0730: 6c 65 6e 20 2f 20 74 6f 74 61 6c 6c 65 6e 3b 0d  len / totallen;.
0740: 0a 09 7d 0d 0a 09 64 6f 75 62 6c 65 20 73 6f 6c  ..}...double sol
0750: 76 65 28 43 4d 50 20 61 2c 20 43 4d 50 20 62 2c  ve(CMP a, CMP b,
0760: 20 43 4d 50 20 63 2c 20 43 4d 50 20 64 29 0d 0a   CMP c, CMP d)..
0770: 09 7b 0d 0a 09 09 2f 2f 20 73 6f 6c 76 65 3a 20  .{....// solve: 
0780: 61 20 2b 20 3c 72 65 74 75 72 6e 20 3f 3e 2a 62  a + <return ?>*b
0790: 20 3d 3d 20 63 20 2b 20 3f 3f 2a 64 0d 0a 09 09   == c + ??*d....
07a0: 72 65 74 75 72 6e 20 2d 31 3b 0d 0a 09 7d 0d 0a  return -1;...}..
07b0: 09 64 6f 75 62 6c 65 20 63 6f 6d 70 75 74 65 46  .double computeF
07c0: 72 61 63 74 69 6f 6e 4f 66 43 72 6f 73 73 28 43  ractionOfCross(C
07d0: 4d 50 20 70 30 2c 20 43 4d 50 20 70 31 2c 20 43  MP p0, CMP p1, C
07e0: 4d 50 20 71 2c 20 63 6f 6e 73 74 20 76 65 63 74  MP q, const vect
07f0: 6f 72 3c 43 4d 50 3e 26 20 64 29 0d 0a 09 7b 0d  or<CMP>& d)...{.
0800: 0a 09 09 2f 2f 20 72 65 74 75 72 20 74 68 65 20  ...// retur the 
0810: 66 72 61 63 74 69 6f 6e 20 6f 66 20 6b 20 69 6e  fraction of k in
0820: 20 5b 30 2c 31 5d 20 20 73 2e 74 2e 20 20 6c 69   [0,1]  s.t.  li
0830: 6e 65 20 28 70 30 20 2b 20 6b 28 70 31 2d 70 30  ne (p0 + k(p1-p0
0840: 29 29 20 74 6f 20 71 20 63 72 6f 73 73 65 73 20  )) to q crosses 
0850: 64 2e 0d 0a 09 09 72 65 74 75 72 6e 20 2d 31 3b  d.....return -1;
0860: 0d 0a 09 7d 0d 0a 7d 3b 0d 0a 0d 0a 2f 2f 20 42  ...}..};....// B
0870: 45 47 49 4e 20 43 55 54 20 48 45 52 45 0d 0a 23  EGIN CUT HERE..#
0880: 69 6e 63 6c 75 64 65 20 3c 63 74 69 6d 65 3e 0d  include <ctime>.
0890: 0a 64 6f 75 62 6c 65 20 73 74 61 72 74 5f 74 69  .double start_ti
08a0: 6d 65 3b 20 73 74 72 69 6e 67 20 74 69 6d 65 72  me; string timer
08b0: 28 29 0d 0a 20 7b 20 6f 73 74 72 69 6e 67 73 74  ().. { ostringst
08c0: 72 65 61 6d 20 6f 73 3b 20 6f 73 20 3c 3c 20 22  ream os; os << "
08d0: 20 28 22 20 3c 3c 20 69 6e 74 28 28 63 6c 6f 63   (" << int((cloc
08e0: 6b 28 29 2d 73 74 61 72 74 5f 74 69 6d 65 29 2f  k()-start_time)/
08f0: 43 4c 4f 43 4b 53 5f 50 45 52 5f 53 45 43 2a 31  CLOCKS_PER_SEC*1
0900: 30 30 30 29 20 3c 3c 20 22 20 6d 73 65 63 29 22  000) << " msec)"
0910: 3b 20 72 65 74 75 72 6e 20 6f 73 2e 73 74 72 28  ; return os.str(
0920: 29 3b 20 7d 0d 0a 74 65 6d 70 6c 61 74 65 3c 74  ); }..template<t
0930: 79 70 65 6e 61 6d 65 20 54 3e 20 6f 73 74 72 65  ypename T> ostre
0940: 61 6d 26 20 6f 70 65 72 61 74 6f 72 3c 3c 28 6f  am& operator<<(o
0950: 73 74 72 65 61 6d 26 20 6f 73 2c 20 63 6f 6e 73  stream& os, cons
0960: 74 20 76 65 63 74 6f 72 3c 54 3e 26 20 76 29 0d  t vector<T>& v).
0970: 0a 20 7b 20 6f 73 20 3c 3c 20 22 7b 20 22 3b 0d  . { os << "{ ";.
0980: 0a 20 20 20 66 6f 72 28 74 79 70 65 6e 61 6d 65  .   for(typename
0990: 20 76 65 63 74 6f 72 3c 54 3e 3a 3a 63 6f 6e 73   vector<T>::cons
09a0: 74 5f 69 74 65 72 61 74 6f 72 20 69 74 3d 76 2e  t_iterator it=v.
09b0: 62 65 67 69 6e 28 29 3b 20 69 74 21 3d 76 2e 65  begin(); it!=v.e
09c0: 6e 64 28 29 3b 20 2b 2b 69 74 29 0d 0a 20 20 20  nd(); ++it)..   
09d0: 6f 73 20 3c 3c 20 27 5c 22 27 20 3c 3c 20 2a 69  os << '\"' << *i
09e0: 74 20 3c 3c 20 27 5c 22 27 20 3c 3c 20 28 69 74  t << '\"' << (it
09f0: 2b 31 3d 3d 76 2e 65 6e 64 28 29 20 3f 20 22 22  +1==v.end() ? ""
0a00: 20 3a 20 22 2c 20 22 29 3b 20 6f 73 20 3c 3c 20   : ", "); os << 
0a10: 22 20 7d 22 3b 20 72 65 74 75 72 6e 20 6f 73 3b  " }"; return os;
0a20: 20 7d 0d 0a 76 6f 69 64 20 76 65 72 69 66 79 5f   }..void verify_
0a30: 63 61 73 65 28 63 6f 6e 73 74 20 64 6f 75 62 6c  case(const doubl
0a40: 65 26 20 45 78 70 65 63 74 65 64 2c 20 63 6f 6e  e& Expected, con
0a50: 73 74 20 64 6f 75 62 6c 65 26 20 52 65 63 65 69  st double& Recei
0a60: 76 65 64 29 20 7b 0d 0a 20 62 6f 6f 6c 20 6f 6b  ved) {.. bool ok
0a70: 20 3d 20 28 61 62 73 28 45 78 70 65 63 74 65 64   = (abs(Expected
0a80: 20 2d 20 52 65 63 65 69 76 65 64 29 20 3c 20 31   - Received) < 1
0a90: 65 2d 39 29 3b 0d 0a 20 69 66 28 6f 6b 29 20 63  e-9);.. if(ok) c
0aa0: 65 72 72 20 3c 3c 20 22 50 41 53 53 45 44 22 20  err << "PASSED" 
0ab0: 3c 3c 20 74 69 6d 65 72 28 29 20 3c 3c 20 65 6e  << timer() << en
0ac0: 64 6c 3b 20 20 65 6c 73 65 20 7b 20 63 65 72 72  dl;  else { cerr
0ad0: 20 3c 3c 20 22 46 41 49 4c 45 44 22 20 3c 3c 20   << "FAILED" << 
0ae0: 74 69 6d 65 72 28 29 20 3c 3c 20 65 6e 64 6c 3b  timer() << endl;
0af0: 0d 0a 20 63 65 72 72 20 3c 3c 20 22 5c 74 6f 3a  .. cerr << "\to:
0b00: 20 5c 22 22 20 3c 3c 20 45 78 70 65 63 74 65 64   \"" << Expected
0b10: 20 3c 3c 20 27 5c 22 27 20 3c 3c 20 65 6e 64 6c   << '\"' << endl
0b20: 20 3c 3c 20 22 5c 74 78 3a 20 5c 22 22 20 3c 3c   << "\tx: \"" <<
0b30: 20 52 65 63 65 69 76 65 64 20 3c 3c 20 27 5c 22   Received << '\"
0b40: 27 20 3c 3c 20 65 6e 64 6c 3b 20 7d 20 7d 0d 0a  ' << endl; } }..
0b50: 23 64 65 66 69 6e 65 20 43 41 53 45 28 4e 29 20  #define CASE(N) 
0b60: 7b 63 65 72 72 20 3c 3c 20 22 54 65 73 74 20 43  {cerr << "Test C
0b70: 61 73 65 20 23 22 20 3c 3c 20 4e 20 3c 3c 20 22  ase #" << N << "
0b80: 2e 2e 2e 22 20 3c 3c 20 66 6c 75 73 68 3b 20 73  ..." << flush; s
0b90: 74 61 72 74 5f 74 69 6d 65 3d 63 6c 6f 63 6b 28  tart_time=clock(
0ba0: 29 3b 0d 0a 23 64 65 66 69 6e 65 20 45 4e 44 09  );..#define END.
0bb0: 20 76 65 72 69 66 79 5f 63 61 73 65 28 5f 2c 20   verify_case(_, 
0bc0: 44 65 70 6f 73 69 74 28 29 2e 73 75 63 63 65 73  Deposit().succes
0bd0: 73 50 72 6f 62 61 62 69 6c 69 74 79 28 73 69 74  sProbability(sit
0be0: 65 58 2c 20 73 69 74 65 59 2c 20 64 65 70 6f 73  eX, siteY, depos
0bf0: 69 74 58 2c 20 64 65 70 6f 73 69 74 59 29 29 3b  itX, depositY));
0c00: 7d 0d 0a 69 6e 74 20 6d 61 69 6e 28 29 7b 0d 0a  }..int main(){..
0c10: 0d 0a 43 41 53 45 28 30 29 0d 0a 09 69 6e 74 20  ..CASE(0)...int 
0c20: 73 69 74 65 58 5f 5b 5d 20 3d 20 7b 30 2c 34 2c  siteX_[] = {0,4,
0c30: 34 2c 30 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72  4,0};...  vector
0c40: 20 3c 69 6e 74 3e 20 73 69 74 65 58 28 73 69 74   <int> siteX(sit
0c50: 65 58 5f 2c 20 73 69 74 65 58 5f 2b 73 69 7a 65  eX_, siteX_+size
0c60: 6f 66 28 73 69 74 65 58 5f 29 2f 73 69 7a 65 6f  of(siteX_)/sizeo
0c70: 66 28 2a 73 69 74 65 58 5f 29 29 3b 20 0d 0a 09  f(*siteX_)); ...
0c80: 69 6e 74 20 73 69 74 65 59 5f 5b 5d 20 3d 20 7b  int siteY_[] = {
0c90: 30 2c 30 2c 34 2c 34 7d 3b 0d 0a 09 20 20 76 65  0,0,4,4};...  ve
0ca0: 63 74 6f 72 20 3c 69 6e 74 3e 20 73 69 74 65 59  ctor <int> siteY
0cb0: 28 73 69 74 65 59 5f 2c 20 73 69 74 65 59 5f 2b  (siteY_, siteY_+
0cc0: 73 69 7a 65 6f 66 28 73 69 74 65 59 5f 29 2f 73  sizeof(siteY_)/s
0cd0: 69 7a 65 6f 66 28 2a 73 69 74 65 59 5f 29 29 3b  izeof(*siteY_));
0ce0: 20 0d 0a 09 69 6e 74 20 64 65 70 6f 73 69 74 58   ...int depositX
0cf0: 5f 5b 5d 20 3d 20 7b 31 2c 32 2c 32 2c 31 7d 3b  _[] = {1,2,2,1};
0d00: 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74  ...  vector <int
0d10: 3e 20 64 65 70 6f 73 69 74 58 28 64 65 70 6f 73  > depositX(depos
0d20: 69 74 58 5f 2c 20 64 65 70 6f 73 69 74 58 5f 2b  itX_, depositX_+
0d30: 73 69 7a 65 6f 66 28 64 65 70 6f 73 69 74 58 5f  sizeof(depositX_
0d40: 29 2f 73 69 7a 65 6f 66 28 2a 64 65 70 6f 73 69  )/sizeof(*deposi
0d50: 74 58 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 64 65  tX_)); ...int de
0d60: 70 6f 73 69 74 59 5f 5b 5d 20 3d 20 7b 31 2c 31  positY_[] = {1,1
0d70: 2c 32 2c 32 7d 3b 0d 0a 09 20 20 76 65 63 74 6f  ,2,2};...  vecto
0d80: 72 20 3c 69 6e 74 3e 20 64 65 70 6f 73 69 74 59  r <int> depositY
0d90: 28 64 65 70 6f 73 69 74 59 5f 2c 20 64 65 70 6f  (depositY_, depo
0da0: 73 69 74 59 5f 2b 73 69 7a 65 6f 66 28 64 65 70  sitY_+sizeof(dep
0db0: 6f 73 69 74 59 5f 29 2f 73 69 7a 65 6f 66 28 2a  ositY_)/sizeof(*
0dc0: 64 65 70 6f 73 69 74 59 5f 29 29 3b 20 0d 0a 09  depositY_)); ...
0dd0: 64 6f 75 62 6c 65 20 5f 20 3d 20 30 2e 36 36 36  double _ = 0.666
0de0: 36 36 36 36 36 36 36 36 36 36 36 36 36 3b 20 0d  6666666666666; .
0df0: 0a 45 4e 44 0d 0a 43 41 53 45 28 31 29 0d 0a 09  .END..CASE(1)...
0e00: 69 6e 74 20 73 69 74 65 58 5f 5b 5d 20 3d 20 7b  int siteX_[] = {
0e10: 30 2c 34 2c 34 2c 30 7d 3b 0d 0a 09 20 20 76 65  0,4,4,0};...  ve
0e20: 63 74 6f 72 20 3c 69 6e 74 3e 20 73 69 74 65 58  ctor <int> siteX
0e30: 28 73 69 74 65 58 5f 2c 20 73 69 74 65 58 5f 2b  (siteX_, siteX_+
0e40: 73 69 7a 65 6f 66 28 73 69 74 65 58 5f 29 2f 73  sizeof(siteX_)/s
0e50: 69 7a 65 6f 66 28 2a 73 69 74 65 58 5f 29 29 3b  izeof(*siteX_));
0e60: 20 0d 0a 09 69 6e 74 20 73 69 74 65 59 5f 5b 5d   ...int siteY_[]
0e70: 20 3d 20 7b 30 2c 30 2c 34 2c 34 7d 3b 0d 0a 09   = {0,0,4,4};...
0e80: 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 73    vector <int> s
0e90: 69 74 65 59 28 73 69 74 65 59 5f 2c 20 73 69 74  iteY(siteY_, sit
0ea0: 65 59 5f 2b 73 69 7a 65 6f 66 28 73 69 74 65 59  eY_+sizeof(siteY
0eb0: 5f 29 2f 73 69 7a 65 6f 66 28 2a 73 69 74 65 59  _)/sizeof(*siteY
0ec0: 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 64 65 70 6f  _)); ...int depo
0ed0: 73 69 74 58 5f 5b 5d 20 3d 20 7b 31 2c 33 2c 33  sitX_[] = {1,3,3
0ee0: 2c 31 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20  ,1};...  vector 
0ef0: 3c 69 6e 74 3e 20 64 65 70 6f 73 69 74 58 28 64  <int> depositX(d
0f00: 65 70 6f 73 69 74 58 5f 2c 20 64 65 70 6f 73 69  epositX_, deposi
0f10: 74 58 5f 2b 73 69 7a 65 6f 66 28 64 65 70 6f 73  tX_+sizeof(depos
0f20: 69 74 58 5f 29 2f 73 69 7a 65 6f 66 28 2a 64 65  itX_)/sizeof(*de
0f30: 70 6f 73 69 74 58 5f 29 29 3b 20 0d 0a 09 69 6e  positX_)); ...in
0f40: 74 20 64 65 70 6f 73 69 74 59 5f 5b 5d 20 3d 20  t depositY_[] = 
0f50: 7b 31 2c 31 2c 33 2c 33 7d 3b 0d 0a 09 20 20 76  {1,1,3,3};...  v
0f60: 65 63 74 6f 72 20 3c 69 6e 74 3e 20 64 65 70 6f  ector <int> depo
0f70: 73 69 74 59 28 64 65 70 6f 73 69 74 59 5f 2c 20  sitY(depositY_, 
0f80: 64 65 70 6f 73 69 74 59 5f 2b 73 69 7a 65 6f 66  depositY_+sizeof
0f90: 28 64 65 70 6f 73 69 74 59 5f 29 2f 73 69 7a 65  (depositY_)/size
0fa0: 6f 66 28 2a 64 65 70 6f 73 69 74 59 5f 29 29 3b  of(*depositY_));
0fb0: 20 0d 0a 09 64 6f 75 62 6c 65 20 5f 20 3d 20 31   ...double _ = 1
0fc0: 2e 30 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28  .0; ..END..CASE(
0fd0: 32 29 0d 0a 09 69 6e 74 20 73 69 74 65 58 5f 5b  2)...int siteX_[
0fe0: 5d 20 3d 20 7b 35 2c 32 2c 2d 31 2c 30 2c 34 7d  ] = {5,2,-1,0,4}
0ff0: 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e  ;...  vector <in
1000: 74 3e 20 73 69 74 65 58 28 73 69 74 65 58 5f 2c  t> siteX(siteX_,
1010: 20 73 69 74 65 58 5f 2b 73 69 7a 65 6f 66 28 73   siteX_+sizeof(s
1020: 69 74 65 58 5f 29 2f 73 69 7a 65 6f 66 28 2a 73  iteX_)/sizeof(*s
1030: 69 74 65 58 5f 29 29 3b 20 0d 0a 09 69 6e 74 20  iteX_)); ...int 
1040: 73 69 74 65 59 5f 5b 5d 20 3d 20 7b 33 2c 34 2c  siteY_[] = {3,4,
1050: 32 2c 30 2c 30 7d 3b 0d 0a 09 20 20 76 65 63 74  2,0,0};...  vect
1060: 6f 72 20 3c 69 6e 74 3e 20 73 69 74 65 59 28 73  or <int> siteY(s
1070: 69 74 65 59 5f 2c 20 73 69 74 65 59 5f 2b 73 69  iteY_, siteY_+si
1080: 7a 65 6f 66 28 73 69 74 65 59 5f 29 2f 73 69 7a  zeof(siteY_)/siz
1090: 65 6f 66 28 2a 73 69 74 65 59 5f 29 29 3b 20 0d  eof(*siteY_)); .
10a0: 0a 09 69 6e 74 20 64 65 70 6f 73 69 74 58 5f 5b  ..int depositX_[
10b0: 5d 20 3d 20 7b 33 2c 33 2c 34 7d 3b 0d 0a 09 20  ] = {3,3,4};... 
10c0: 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 64 65   vector <int> de
10d0: 70 6f 73 69 74 58 28 64 65 70 6f 73 69 74 58 5f  positX(depositX_
10e0: 2c 20 64 65 70 6f 73 69 74 58 5f 2b 73 69 7a 65  , depositX_+size
10f0: 6f 66 28 64 65 70 6f 73 69 74 58 5f 29 2f 73 69  of(depositX_)/si
1100: 7a 65 6f 66 28 2a 64 65 70 6f 73 69 74 58 5f 29  zeof(*depositX_)
1110: 29 3b 20 0d 0a 09 69 6e 74 20 64 65 70 6f 73 69  ); ...int deposi
1120: 74 59 5f 5b 5d 20 3d 20 7b 33 2c 32 2c 31 7d 3b  tY_[] = {3,2,1};
1130: 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74  ...  vector <int
1140: 3e 20 64 65 70 6f 73 69 74 59 28 64 65 70 6f 73  > depositY(depos
1150: 69 74 59 5f 2c 20 64 65 70 6f 73 69 74 59 5f 2b  itY_, depositY_+
1160: 73 69 7a 65 6f 66 28 64 65 70 6f 73 69 74 59 5f  sizeof(depositY_
1170: 29 2f 73 69 7a 65 6f 66 28 2a 64 65 70 6f 73 69  )/sizeof(*deposi
1180: 74 59 5f 29 29 3b 20 0d 0a 09 64 6f 75 62 6c 65  tY_)); ...double
1190: 20 5f 20 3d 20 30 2e 36 31 31 32 37 30 30 32 30   _ = 0.611270020
11a0: 39 38 35 35 34 32 33 3b 20 0d 0a 45 4e 44 0d 0a  9855423; ..END..
11b0: 43 41 53 45 28 33 29 0d 0a 09 69 6e 74 20 73 69  CASE(3)...int si
11c0: 74 65 58 5f 5b 5d 20 3d 20 7b 32 30 30 2c 2d 39  teX_[] = {200,-9
11d0: 39 2c 2d 34 30 35 2c 2d 36 30 31 2c 2d 37 30 38  9,-405,-601,-708
11e0: 2c 2d 34 39 34 2c 2d 33 30 30 2c 2d 38 38 7d 3b  ,-494,-300,-88};
11f0: 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74  ...  vector <int
1200: 3e 20 73 69 74 65 58 28 73 69 74 65 58 5f 2c 20  > siteX(siteX_, 
1210: 73 69 74 65 58 5f 2b 73 69 7a 65 6f 66 28 73 69  siteX_+sizeof(si
1220: 74 65 58 5f 29 2f 73 69 7a 65 6f 66 28 2a 73 69  teX_)/sizeof(*si
1230: 74 65 58 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 73  teX_)); ...int s
1240: 69 74 65 59 5f 5b 5d 20 3d 20 7b 35 32 30 2c 35  iteY_[] = {520,5
1250: 31 36 2c 34 30 37 2c 33 32 31 2c 31 30 34 2c 2d  16,407,321,104,-
1260: 39 37 2c 2d 32 32 31 2c 2d 31 30 31 7d 3b 0d 0a  97,-221,-101};..
1270: 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20  .  vector <int> 
1280: 73 69 74 65 59 28 73 69 74 65 59 5f 2c 20 73 69  siteY(siteY_, si
1290: 74 65 59 5f 2b 73 69 7a 65 6f 66 28 73 69 74 65  teY_+sizeof(site
12a0: 59 5f 29 2f 73 69 7a 65 6f 66 28 2a 73 69 74 65  Y_)/sizeof(*site
12b0: 59 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 64 65 70  Y_)); ...int dep
12c0: 6f 73 69 74 58 5f 5b 5d 20 3d 20 7b 2d 31 30 31  ositX_[] = {-101
12d0: 2c 2d 32 30 31 2c 2d 32 39 36 2c 2d 34 30 30 2c  ,-201,-296,-400,
12e0: 2d 34 30 32 7d 3b 0d 0a 09 20 20 76 65 63 74 6f  -402};...  vecto
12f0: 72 20 3c 69 6e 74 3e 20 64 65 70 6f 73 69 74 58  r <int> depositX
1300: 28 64 65 70 6f 73 69 74 58 5f 2c 20 64 65 70 6f  (depositX_, depo
1310: 73 69 74 58 5f 2b 73 69 7a 65 6f 66 28 64 65 70  sitX_+sizeof(dep
1320: 6f 73 69 74 58 5f 29 2f 73 69 7a 65 6f 66 28 2a  ositX_)/sizeof(*
1330: 64 65 70 6f 73 69 74 58 5f 29 29 3b 20 0d 0a 09  depositX_)); ...
1340: 69 6e 74 20 64 65 70 6f 73 69 74 59 5f 5b 5d 20  int depositY_[] 
1350: 3d 20 7b 33 31 38 2c 33 39 36 2c 34 30 32 2c 33  = {318,396,402,3
1360: 30 35 2c 32 30 30 7d 3b 0d 0a 09 20 20 76 65 63  05,200};...  vec
1370: 74 6f 72 20 3c 69 6e 74 3e 20 64 65 70 6f 73 69  tor <int> deposi
1380: 74 59 28 64 65 70 6f 73 69 74 59 5f 2c 20 64 65  tY(depositY_, de
1390: 70 6f 73 69 74 59 5f 2b 73 69 7a 65 6f 66 28 64  positY_+sizeof(d
13a0: 65 70 6f 73 69 74 59 5f 29 2f 73 69 7a 65 6f 66  epositY_)/sizeof
13b0: 28 2a 64 65 70 6f 73 69 74 59 5f 29 29 3b 20 0d  (*depositY_)); .
13c0: 0a 09 64 6f 75 62 6c 65 20 5f 20 3d 20 30 2e 34  ..double _ = 0.4
13d0: 39 38 39 32 37 35 36 32 30 37 31 30 30 37 34 37  9892756207100747
13e0: 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 34 29  ; ..END..CASE(4)
13f0: 0d 0a 09 69 6e 74 20 73 69 74 65 58 5f 5b 5d 20  ...int siteX_[] 
1400: 3d 20 7b 2d 33 2c 36 2c 31 30 2c 31 30 2c 2d 32  = {-3,6,10,10,-2
1410: 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69  };...  vector <i
1420: 6e 74 3e 20 73 69 74 65 58 28 73 69 74 65 58 5f  nt> siteX(siteX_
1430: 2c 20 73 69 74 65 58 5f 2b 73 69 7a 65 6f 66 28  , siteX_+sizeof(
1440: 73 69 74 65 58 5f 29 2f 73 69 7a 65 6f 66 28 2a  siteX_)/sizeof(*
1450: 73 69 74 65 58 5f 29 29 3b 20 0d 0a 09 69 6e 74  siteX_)); ...int
1460: 20 73 69 74 65 59 5f 5b 5d 20 3d 20 7b 2d 31 2c   siteY_[] = {-1,
1470: 2d 31 30 2c 2d 39 2c 37 2c 33 7d 3b 0d 0a 09 20  -10,-9,7,3};... 
1480: 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 73 69   vector <int> si
1490: 74 65 59 28 73 69 74 65 59 5f 2c 20 73 69 74 65  teY(siteY_, site
14a0: 59 5f 2b 73 69 7a 65 6f 66 28 73 69 74 65 59 5f  Y_+sizeof(siteY_
14b0: 29 2f 73 69 7a 65 6f 66 28 2a 73 69 74 65 59 5f  )/sizeof(*siteY_
14c0: 29 29 3b 20 0d 0a 09 69 6e 74 20 64 65 70 6f 73  )); ...int depos
14d0: 69 74 58 5f 5b 5d 20 3d 20 7b 33 2c 35 2c 36 2c  itX_[] = {3,5,6,
14e0: 34 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c  4};...  vector <
14f0: 69 6e 74 3e 20 64 65 70 6f 73 69 74 58 28 64 65  int> depositX(de
1500: 70 6f 73 69 74 58 5f 2c 20 64 65 70 6f 73 69 74  positX_, deposit
1510: 58 5f 2b 73 69 7a 65 6f 66 28 64 65 70 6f 73 69  X_+sizeof(deposi
1520: 74 58 5f 29 2f 73 69 7a 65 6f 66 28 2a 64 65 70  tX_)/sizeof(*dep
1530: 6f 73 69 74 58 5f 29 29 3b 20 0d 0a 09 69 6e 74  ositX_)); ...int
1540: 20 64 65 70 6f 73 69 74 59 5f 5b 5d 20 3d 20 7b   depositY_[] = {
1550: 2d 35 2c 2d 37 2c 2d 36 2c 2d 34 7d 3b 0d 0a 09  -5,-7,-6,-4};...
1560: 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 64    vector <int> d
1570: 65 70 6f 73 69 74 59 28 64 65 70 6f 73 69 74 59  epositY(depositY
1580: 5f 2c 20 64 65 70 6f 73 69 74 59 5f 2b 73 69 7a  _, depositY_+siz
1590: 65 6f 66 28 64 65 70 6f 73 69 74 59 5f 29 2f 73  eof(depositY_)/s
15a0: 69 7a 65 6f 66 28 2a 64 65 70 6f 73 69 74 59 5f  izeof(*depositY_
15b0: 29 29 3b 20 0d 0a 09 64 6f 75 62 6c 65 20 5f 20  )); ...double _ 
15c0: 3d 20 30 2e 31 36 36 34 31 36 36 32 32 35 31 35  = 0.166416622515
15d0: 32 31 35 33 38 3b 20 0d 0a 45 4e 44 0d 0a 43 41  21538; ..END..CA
15e0: 53 45 28 35 29 0d 0a 09 69 6e 74 20 73 69 74 65  SE(5)...int site
15f0: 58 5f 5b 5d 20 3d 20 7b 2d 32 2c 34 2c 2d 32 7d  X_[] = {-2,4,-2}
1600: 0d 0a 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c  ..;...  vector <
1610: 69 6e 74 3e 20 73 69 74 65 58 28 73 69 74 65 58  int> siteX(siteX
1620: 5f 2c 20 73 69 74 65 58 5f 2b 73 69 7a 65 6f 66  _, siteX_+sizeof
1630: 28 73 69 74 65 58 5f 29 2f 73 69 7a 65 6f 66 28  (siteX_)/sizeof(
1640: 2a 73 69 74 65 58 5f 29 29 3b 20 0d 0a 09 69 6e  *siteX_)); ...in
1650: 74 20 73 69 74 65 59 5f 5b 5d 20 3d 20 7b 2d 32  t siteY_[] = {-2
1660: 2c 2d 32 2c 38 7d 0d 0a 3b 0d 0a 09 20 20 76 65  ,-2,8}..;...  ve
1670: 63 74 6f 72 20 3c 69 6e 74 3e 20 73 69 74 65 59  ctor <int> siteY
1680: 28 73 69 74 65 59 5f 2c 20 73 69 74 65 59 5f 2b  (siteY_, siteY_+
1690: 73 69 7a 65 6f 66 28 73 69 74 65 59 5f 29 2f 73  sizeof(siteY_)/s
16a0: 69 7a 65 6f 66 28 2a 73 69 74 65 59 5f 29 29 3b  izeof(*siteY_));
16b0: 20 0d 0a 09 69 6e 74 20 64 65 70 6f 73 69 74 58   ...int depositX
16c0: 5f 5b 5d 20 3d 20 7b 2d 31 2c 30 2c 30 7d 0d 0a  _[] = {-1,0,0}..
16d0: 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e  ;...  vector <in
16e0: 74 3e 20 64 65 70 6f 73 69 74 58 28 64 65 70 6f  t> depositX(depo
16f0: 73 69 74 58 5f 2c 20 64 65 70 6f 73 69 74 58 5f  sitX_, depositX_
1700: 2b 73 69 7a 65 6f 66 28 64 65 70 6f 73 69 74 58  +sizeof(depositX
1710: 5f 29 2f 73 69 7a 65 6f 66 28 2a 64 65 70 6f 73  _)/sizeof(*depos
1720: 69 74 58 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 64  itX_)); ...int d
1730: 65 70 6f 73 69 74 59 5f 5b 5d 20 3d 20 7b 30 2c  epositY_[] = {0,
1740: 2d 31 2c 30 7d 0d 0a 3b 0d 0a 09 20 20 76 65 63  -1,0}..;...  vec
1750: 74 6f 72 20 3c 69 6e 74 3e 20 64 65 70 6f 73 69  tor <int> deposi
1760: 74 59 28 64 65 70 6f 73 69 74 59 5f 2c 20 64 65  tY(depositY_, de
1770: 70 6f 73 69 74 59 5f 2b 73 69 7a 65 6f 66 28 64  positY_+sizeof(d
1780: 65 70 6f 73 69 74 59 5f 29 2f 73 69 7a 65 6f 66  epositY_)/sizeof
1790: 28 2a 64 65 70 6f 73 69 74 59 5f 29 29 3b 20 0d  (*depositY_)); .
17a0: 0a 09 64 6f 75 62 6c 65 20 5f 20 3d 20 30 2e 30  ..double _ = 0.0
17b0: 34 35 31 38 38 35 30 32 31 39 30 37 32 32 39 31  4518850219072291
17c0: 3b 20 0d 0a 45 4e 44 0d 0a 2f 2a 0d 0a 43 41 53  ; ..END../*..CAS
17d0: 45 28 36 29 0d 0a 09 69 6e 74 20 73 69 74 65 58  E(6)...int siteX
17e0: 5f 5b 5d 20 3d 20 3b 0d 0a 09 20 20 76 65 63 74  _[] = ;...  vect
17f0: 6f 72 20 3c 69 6e 74 3e 20 73 69 74 65 58 28 73  or <int> siteX(s
1800: 69 74 65 58 5f 2c 20 73 69 74 65 58 5f 2b 73 69  iteX_, siteX_+si
1810: 7a 65 6f 66 28 73 69 74 65 58 5f 29 2f 73 69 7a  zeof(siteX_)/siz
1820: 65 6f 66 28 2a 73 69 74 65 58 5f 29 29 3b 20 0d  eof(*siteX_)); .
1830: 0a 09 69 6e 74 20 73 69 74 65 59 5f 5b 5d 20 3d  ..int siteY_[] =
1840: 20 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69   ;...  vector <i
1850: 6e 74 3e 20 73 69 74 65 59 28 73 69 74 65 59 5f  nt> siteY(siteY_
1860: 2c 20 73 69 74 65 59 5f 2b 73 69 7a 65 6f 66 28  , siteY_+sizeof(
1870: 73 69 74 65 59 5f 29 2f 73 69 7a 65 6f 66 28 2a  siteY_)/sizeof(*
1880: 73 69 74 65 59 5f 29 29 3b 20 0d 0a 09 69 6e 74  siteY_)); ...int
1890: 20 64 65 70 6f 73 69 74 58 5f 5b 5d 20 3d 20 3b   depositX_[] = ;
18a0: 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74  ...  vector <int
18b0: 3e 20 64 65 70 6f 73 69 74 58 28 64 65 70 6f 73  > depositX(depos
18c0: 69 74 58 5f 2c 20 64 65 70 6f 73 69 74 58 5f 2b  itX_, depositX_+
18d0: 73 69 7a 65 6f 66 28 64 65 70 6f 73 69 74 58 5f  sizeof(depositX_
18e0: 29 2f 73 69 7a 65 6f 66 28 2a 64 65 70 6f 73 69  )/sizeof(*deposi
18f0: 74 58 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 64 65  tX_)); ...int de
1900: 70 6f 73 69 74 59 5f 5b 5d 20 3d 20 3b 0d 0a 09  positY_[] = ;...
1910: 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 64    vector <int> d
1920: 65 70 6f 73 69 74 59 28 64 65 70 6f 73 69 74 59  epositY(depositY
1930: 5f 2c 20 64 65 70 6f 73 69 74 59 5f 2b 73 69 7a  _, depositY_+siz
1940: 65 6f 66 28 64 65 70 6f 73 69 74 59 5f 29 2f 73  eof(depositY_)/s
1950: 69 7a 65 6f 66 28 2a 64 65 70 6f 73 69 74 59 5f  izeof(*depositY_
1960: 29 29 3b 20 0d 0a 09 64 6f 75 62 6c 65 20 5f 20  )); ...double _ 
1970: 3d 20 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28  = ; ..END..CASE(
1980: 37 29 0d 0a 09 69 6e 74 20 73 69 74 65 58 5f 5b  7)...int siteX_[
1990: 5d 20 3d 20 3b 0d 0a 09 20 20 76 65 63 74 6f 72  ] = ;...  vector
19a0: 20 3c 69 6e 74 3e 20 73 69 74 65 58 28 73 69 74   <int> siteX(sit
19b0: 65 58 5f 2c 20 73 69 74 65 58 5f 2b 73 69 7a 65  eX_, siteX_+size
19c0: 6f 66 28 73 69 74 65 58 5f 29 2f 73 69 7a 65 6f  of(siteX_)/sizeo
19d0: 66 28 2a 73 69 74 65 58 5f 29 29 3b 20 0d 0a 09  f(*siteX_)); ...
19e0: 69 6e 74 20 73 69 74 65 59 5f 5b 5d 20 3d 20 3b  int siteY_[] = ;
19f0: 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74  ...  vector <int
1a00: 3e 20 73 69 74 65 59 28 73 69 74 65 59 5f 2c 20  > siteY(siteY_, 
1a10: 73 69 74 65 59 5f 2b 73 69 7a 65 6f 66 28 73 69  siteY_+sizeof(si
1a20: 74 65 59 5f 29 2f 73 69 7a 65 6f 66 28 2a 73 69  teY_)/sizeof(*si
1a30: 74 65 59 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 64  teY_)); ...int d
1a40: 65 70 6f 73 69 74 58 5f 5b 5d 20 3d 20 3b 0d 0a  epositX_[] = ;..
1a50: 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20  .  vector <int> 
1a60: 64 65 70 6f 73 69 74 58 28 64 65 70 6f 73 69 74  depositX(deposit
1a70: 58 5f 2c 20 64 65 70 6f 73 69 74 58 5f 2b 73 69  X_, depositX_+si
1a80: 7a 65 6f 66 28 64 65 70 6f 73 69 74 58 5f 29 2f  zeof(depositX_)/
1a90: 73 69 7a 65 6f 66 28 2a 64 65 70 6f 73 69 74 58  sizeof(*depositX
1aa0: 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 64 65 70 6f  _)); ...int depo
1ab0: 73 69 74 59 5f 5b 5d 20 3d 20 3b 0d 0a 09 20 20  sitY_[] = ;...  
1ac0: 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 64 65 70  vector <int> dep
1ad0: 6f 73 69 74 59 28 64 65 70 6f 73 69 74 59 5f 2c  ositY(depositY_,
1ae0: 20 64 65 70 6f 73 69 74 59 5f 2b 73 69 7a 65 6f   depositY_+sizeo
1af0: 66 28 64 65 70 6f 73 69 74 59 5f 29 2f 73 69 7a  f(depositY_)/siz
1b00: 65 6f 66 28 2a 64 65 70 6f 73 69 74 59 5f 29 29  eof(*depositY_))
1b10: 3b 20 0d 0a 09 64 6f 75 62 6c 65 20 5f 20 3d 20  ; ...double _ = 
1b20: 3b 20 0d 0a 45 4e 44 0d 0a 2a 2f 0d 0a 7d 0d 0a  ; ..END..*/..}..
1b30: 2f 2f 20 45 4e 44 20 43 55 54 20 48 45 52 45 0d  // END CUT HERE.
1b40: 0a                                               .