Hex Artifact Content
Not logged in

Artifact dc68ca069cbe1c919585cb445002c044ab380f3b:


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 41 69  CMP;....class Ai
01a0: 72 63 72 61 66 74 20 7b 20 70 75 62 6c 69 63 3a  rcraft { public:
01b0: 0d 0a 09 64 6f 75 62 6c 65 20 66 28 63 6f 6e 73  ...double f(cons
01c0: 74 20 64 6f 75 62 6c 65 20 70 5b 5d 2c 20 63 6f  t double p[], co
01d0: 6e 73 74 20 64 6f 75 62 6c 65 20 76 5b 5d 2c 20  nst double v[], 
01e0: 64 6f 75 62 6c 65 20 74 29 0d 0a 09 7b 0d 0a 09  double t)...{...
01f0: 09 72 65 74 75 72 6e 20 73 71 72 74 28 20 28 70  .return sqrt( (p
0200: 5b 30 5d 2b 74 2a 76 5b 30 5d 29 2a 28 70 5b 30  [0]+t*v[0])*(p[0
0210: 5d 2b 74 2a 76 5b 30 5d 29 20 2b 20 28 70 5b 31  ]+t*v[0]) + (p[1
0220: 5d 2b 74 2a 76 5b 31 5d 29 2a 28 70 5b 31 5d 2b  ]+t*v[1])*(p[1]+
0230: 74 2a 76 5b 31 5d 29 20 2b 20 28 70 5b 32 5d 2b  t*v[1]) + (p[2]+
0240: 74 2a 76 5b 32 5d 29 2a 28 70 5b 32 5d 2b 74 2a  t*v[2])*(p[2]+t*
0250: 76 5b 32 5d 29 20 29 3b 0d 0a 09 7d 0d 0a 09 73  v[2]) );...}...s
0260: 74 72 69 6e 67 20 6e 65 61 72 4d 69 73 73 28 76  tring nearMiss(v
0270: 65 63 74 6f 72 20 3c 69 6e 74 3e 20 70 31 2c 20  ector <int> p1, 
0280: 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 76 31 2c  vector <int> v1,
0290: 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 70 32   vector <int> p2
02a0: 2c 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 76  , vector <int> v
02b0: 32 2c 20 69 6e 74 20 52 29 20 0d 0a 09 7b 0d 0a  2, int R) ...{..
02c0: 09 09 63 6f 6e 73 74 20 64 6f 75 62 6c 65 20 70  ..const double p
02d0: 5b 5d 20 3d 20 7b 70 31 5b 30 5d 2d 70 32 5b 30  [] = {p1[0]-p2[0
02e0: 5d 2c 20 70 31 5b 31 5d 2d 70 32 5b 31 5d 2c 20  ], p1[1]-p2[1], 
02f0: 70 31 5b 32 5d 2d 70 32 5b 32 5d 7d 3b 0d 0a 09  p1[2]-p2[2]};...
0300: 09 63 6f 6e 73 74 20 64 6f 75 62 6c 65 20 76 5b  .const double v[
0310: 5d 20 3d 20 7b 76 31 5b 30 5d 2d 76 32 5b 30 5d  ] = {v1[0]-v2[0]
0320: 2c 20 76 31 5b 31 5d 2d 76 32 5b 31 5d 2c 20 76  , v1[1]-v2[1], v
0330: 31 5b 32 5d 2d 76 32 5b 32 5d 7d 3b 0d 0a 0d 0a  1[2]-v2[2]};....
0340: 09 09 64 6f 75 62 6c 65 20 78 20 3d 20 30 2e 30  ..double x = 0.0
0350: 3b 0d 0a 09 09 64 6f 75 62 6c 65 20 77 20 3d 20  ;....double w = 
0360: 31 65 2b 31 30 3b 0d 0a 09 09 66 6f 72 28 69 6e  1e+10;....for(in
0370: 74 20 6b 3d 30 3b 20 6b 3c 31 30 30 30 30 30 30  t k=0; k<1000000
0380: 3b 20 2b 2b 6b 29 0d 0a 09 09 7b 0d 0a 09 09 09  ; ++k)....{.....
0390: 64 6f 75 62 6c 65 20 79 20 3d 20 32 2a 78 2f 33  double y = 2*x/3
03a0: 20 2b 20 77 2f 33 3b 0d 0a 09 09 09 64 6f 75 62   + w/3;.....doub
03b0: 6c 65 20 7a 20 3d 20 78 2f 33 20 2b 20 32 2a 77  le z = x/3 + 2*w
03c0: 2f 33 3b 0d 0a 09 09 09 64 6f 75 62 6c 65 20 66  /3;.....double f
03d0: 78 20 3d 20 66 28 70 2c 76 2c 78 29 3b 0d 0a 09  x = f(p,v,x);...
03e0: 09 09 64 6f 75 62 6c 65 20 66 79 20 3d 20 66 28  ..double fy = f(
03f0: 70 2c 76 2c 79 29 3b 0d 0a 09 09 09 64 6f 75 62  p,v,y);.....doub
0400: 6c 65 20 66 7a 20 3d 20 66 28 70 2c 76 2c 7a 29  le fz = f(p,v,z)
0410: 3b 0d 0a 09 09 09 64 6f 75 62 6c 65 20 66 77 20  ;.....double fw 
0420: 3d 20 66 28 70 2c 76 2c 77 29 3b 0d 0a 09 09 09  = f(p,v,w);.....
0430: 69 66 28 20 66 78 20 3c 20 66 79 20 29 0d 0a 09  if( fx < fy )...
0440: 09 09 09 77 20 3d 20 79 3b 0d 0a 09 09 09 65 6c  ...w = y;.....el
0450: 73 65 20 69 66 28 20 66 7a 20 3e 20 66 77 20 29  se if( fz > fw )
0460: 0d 0a 09 09 09 09 78 20 3d 20 7a 3b 0d 0a 09 09  ......x = z;....
0470: 09 65 6c 73 65 20 69 66 28 20 66 79 20 3c 20 66  .else if( fy < f
0480: 7a 20 29 0d 0a 09 09 09 09 77 20 3d 20 7a 3b 0d  z )......w = z;.
0490: 0a 09 09 09 65 6c 73 65 0d 0a 09 09 09 09 78 20  ....else......x 
04a0: 3d 20 79 3b 0d 0a 09 09 7d 0d 0a 0d 0a 09 09 72  = y;....}......r
04b0: 65 74 75 72 6e 20 66 28 70 2c 76 2c 78 29 2d 30  eturn f(p,v,x)-0
04c0: 2e 30 30 30 30 30 31 3c 3d 52 20 3f 20 22 59 45  .000001<=R ? "YE
04d0: 53 22 20 3a 20 22 4e 4f 22 3b 0d 0a 09 7d 0d 0a  S" : "NO";...}..
04e0: 7d 3b 0d 0a 0d 0a 2f 2f 20 42 45 47 49 4e 20 43  };....// BEGIN C
04f0: 55 54 20 48 45 52 45 0d 0a 23 69 6e 63 6c 75 64  UT HERE..#includ
0500: 65 20 3c 63 74 69 6d 65 3e 0d 0a 64 6f 75 62 6c  e <ctime>..doubl
0510: 65 20 73 74 61 72 74 5f 74 69 6d 65 3b 20 73 74  e start_time; st
0520: 72 69 6e 67 20 74 69 6d 65 72 28 29 0d 0a 20 7b  ring timer().. {
0530: 20 6f 73 74 72 69 6e 67 73 74 72 65 61 6d 20 6f   ostringstream o
0540: 73 3b 20 6f 73 20 3c 3c 20 22 20 28 22 20 3c 3c  s; os << " (" <<
0550: 20 69 6e 74 28 28 63 6c 6f 63 6b 28 29 2d 73 74   int((clock()-st
0560: 61 72 74 5f 74 69 6d 65 29 2f 43 4c 4f 43 4b 53  art_time)/CLOCKS
0570: 5f 50 45 52 5f 53 45 43 2a 31 30 30 30 29 20 3c  _PER_SEC*1000) <
0580: 3c 20 22 20 6d 73 65 63 29 22 3b 20 72 65 74 75  < " msec)"; retu
0590: 72 6e 20 6f 73 2e 73 74 72 28 29 3b 20 7d 0d 0a  rn os.str(); }..
05a0: 74 65 6d 70 6c 61 74 65 3c 74 79 70 65 6e 61 6d  template<typenam
05b0: 65 20 54 3e 20 6f 73 74 72 65 61 6d 26 20 6f 70  e T> ostream& op
05c0: 65 72 61 74 6f 72 3c 3c 28 6f 73 74 72 65 61 6d  erator<<(ostream
05d0: 26 20 6f 73 2c 20 63 6f 6e 73 74 20 76 65 63 74  & os, const vect
05e0: 6f 72 3c 54 3e 26 20 76 29 0d 0a 20 7b 20 6f 73  or<T>& v).. { os
05f0: 20 3c 3c 20 22 7b 20 22 3b 0d 0a 20 20 20 66 6f   << "{ ";..   fo
0600: 72 28 74 79 70 65 6e 61 6d 65 20 76 65 63 74 6f  r(typename vecto
0610: 72 3c 54 3e 3a 3a 63 6f 6e 73 74 5f 69 74 65 72  r<T>::const_iter
0620: 61 74 6f 72 20 69 74 3d 76 2e 62 65 67 69 6e 28  ator it=v.begin(
0630: 29 3b 20 69 74 21 3d 76 2e 65 6e 64 28 29 3b 20  ); it!=v.end(); 
0640: 2b 2b 69 74 29 0d 0a 20 20 20 6f 73 20 3c 3c 20  ++it)..   os << 
0650: 27 5c 22 27 20 3c 3c 20 2a 69 74 20 3c 3c 20 27  '\"' << *it << '
0660: 5c 22 27 20 3c 3c 20 28 69 74 2b 31 3d 3d 76 2e  \"' << (it+1==v.
0670: 65 6e 64 28 29 20 3f 20 22 22 20 3a 20 22 2c 20  end() ? "" : ", 
0680: 22 29 3b 20 6f 73 20 3c 3c 20 22 20 7d 22 3b 20  "); os << " }"; 
0690: 72 65 74 75 72 6e 20 6f 73 3b 20 7d 0d 0a 76 6f  return os; }..vo
06a0: 69 64 20 76 65 72 69 66 79 5f 63 61 73 65 28 63  id verify_case(c
06b0: 6f 6e 73 74 20 73 74 72 69 6e 67 26 20 45 78 70  onst string& Exp
06c0: 65 63 74 65 64 2c 20 63 6f 6e 73 74 20 73 74 72  ected, const str
06d0: 69 6e 67 26 20 52 65 63 65 69 76 65 64 29 20 7b  ing& Received) {
06e0: 0d 0a 20 62 6f 6f 6c 20 6f 6b 20 3d 20 28 45 78  .. bool ok = (Ex
06f0: 70 65 63 74 65 64 20 3d 3d 20 52 65 63 65 69 76  pected == Receiv
0700: 65 64 29 3b 0d 0a 20 69 66 28 6f 6b 29 20 63 65  ed);.. if(ok) ce
0710: 72 72 20 3c 3c 20 22 50 41 53 53 45 44 22 20 3c  rr << "PASSED" <
0720: 3c 20 74 69 6d 65 72 28 29 20 3c 3c 20 65 6e 64  < timer() << end
0730: 6c 3b 20 20 65 6c 73 65 20 7b 20 63 65 72 72 20  l;  else { cerr 
0740: 3c 3c 20 22 46 41 49 4c 45 44 22 20 3c 3c 20 74  << "FAILED" << t
0750: 69 6d 65 72 28 29 20 3c 3c 20 65 6e 64 6c 3b 0d  imer() << endl;.
0760: 0a 20 63 65 72 72 20 3c 3c 20 22 5c 74 6f 3a 20  . cerr << "\to: 
0770: 5c 22 22 20 3c 3c 20 45 78 70 65 63 74 65 64 20  \"" << Expected 
0780: 3c 3c 20 27 5c 22 27 20 3c 3c 20 65 6e 64 6c 20  << '\"' << endl 
0790: 3c 3c 20 22 5c 74 78 3a 20 5c 22 22 20 3c 3c 20  << "\tx: \"" << 
07a0: 52 65 63 65 69 76 65 64 20 3c 3c 20 27 5c 22 27  Received << '\"'
07b0: 20 3c 3c 20 65 6e 64 6c 3b 20 7d 20 7d 0d 0a 23   << endl; } }..#
07c0: 64 65 66 69 6e 65 20 43 41 53 45 28 4e 29 20 7b  define CASE(N) {
07d0: 63 65 72 72 20 3c 3c 20 22 54 65 73 74 20 43 61  cerr << "Test Ca
07e0: 73 65 20 23 22 20 3c 3c 20 4e 20 3c 3c 20 22 2e  se #" << N << ".
07f0: 2e 2e 22 20 3c 3c 20 66 6c 75 73 68 3b 20 73 74  .." << flush; st
0800: 61 72 74 5f 74 69 6d 65 3d 63 6c 6f 63 6b 28 29  art_time=clock()
0810: 3b 0d 0a 23 64 65 66 69 6e 65 20 45 4e 44 09 20  ;..#define END. 
0820: 76 65 72 69 66 79 5f 63 61 73 65 28 5f 2c 20 41  verify_case(_, A
0830: 69 72 63 72 61 66 74 28 29 2e 6e 65 61 72 4d 69  ircraft().nearMi
0840: 73 73 28 70 31 2c 20 76 31 2c 20 70 32 2c 20 76  ss(p1, v1, p2, v
0850: 32 2c 20 52 29 29 3b 7d 0d 0a 69 6e 74 20 6d 61  2, R));}..int ma
0860: 69 6e 28 29 7b 0d 0a 0d 0a 43 41 53 45 28 30 29  in(){....CASE(0)
0870: 0d 0a 09 69 6e 74 20 70 31 5f 5b 5d 20 3d 20 7b  ...int p1_[] = {
0880: 31 35 2c 35 30 2c 35 7d 3b 0d 0a 09 20 20 76 65  15,50,5};...  ve
0890: 63 74 6f 72 20 3c 69 6e 74 3e 20 70 31 28 70 31  ctor <int> p1(p1
08a0: 5f 2c 20 70 31 5f 2b 73 69 7a 65 6f 66 28 70 31  _, p1_+sizeof(p1
08b0: 5f 29 2f 73 69 7a 65 6f 66 28 2a 70 31 5f 29 29  _)/sizeof(*p1_))
08c0: 3b 20 0d 0a 09 69 6e 74 20 76 31 5f 5b 5d 20 3d  ; ...int v1_[] =
08d0: 20 7b 32 35 2c 31 2c 30 7d 3b 0d 0a 09 20 20 76   {25,1,0};...  v
08e0: 65 63 74 6f 72 20 3c 69 6e 74 3e 20 76 31 28 76  ector <int> v1(v
08f0: 31 5f 2c 20 76 31 5f 2b 73 69 7a 65 6f 66 28 76  1_, v1_+sizeof(v
0900: 31 5f 29 2f 73 69 7a 65 6f 66 28 2a 76 31 5f 29  1_)/sizeof(*v1_)
0910: 29 3b 20 0d 0a 09 69 6e 74 20 70 32 5f 5b 5d 20  ); ...int p2_[] 
0920: 3d 20 7b 31 36 31 2c 31 30 32 2c 39 7d 3b 0d 0a  = {161,102,9};..
0930: 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20  .  vector <int> 
0940: 70 32 28 70 32 5f 2c 20 70 32 5f 2b 73 69 7a 65  p2(p2_, p2_+size
0950: 6f 66 28 70 32 5f 29 2f 73 69 7a 65 6f 66 28 2a  of(p2_)/sizeof(*
0960: 70 32 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 76 32  p2_)); ...int v2
0970: 5f 5b 5d 20 3d 20 7b 2d 31 30 2c 2d 31 30 2c 2d  _[] = {-10,-10,-
0980: 31 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c  1};...  vector <
0990: 69 6e 74 3e 20 76 32 28 76 32 5f 2c 20 76 32 5f  int> v2(v2_, v2_
09a0: 2b 73 69 7a 65 6f 66 28 76 32 5f 29 2f 73 69 7a  +sizeof(v2_)/siz
09b0: 65 6f 66 28 2a 76 32 5f 29 29 3b 20 0d 0a 09 69  eof(*v2_)); ...i
09c0: 6e 74 20 52 20 3d 20 31 30 3b 20 0d 0a 09 73 74  nt R = 10; ...st
09d0: 72 69 6e 67 20 5f 20 3d 20 22 59 45 53 22 3b 20  ring _ = "YES"; 
09e0: 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 31 29 0d 0a  ..END..CASE(1)..
09f0: 09 69 6e 74 20 70 31 5f 5b 5d 20 3d 20 7b 30 2c  .int p1_[] = {0,
0a00: 30 2c 30 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72  0,0};...  vector
0a10: 20 3c 69 6e 74 3e 20 70 31 28 70 31 5f 2c 20 70   <int> p1(p1_, p
0a20: 31 5f 2b 73 69 7a 65 6f 66 28 70 31 5f 29 2f 73  1_+sizeof(p1_)/s
0a30: 69 7a 65 6f 66 28 2a 70 31 5f 29 29 3b 20 0d 0a  izeof(*p1_)); ..
0a40: 09 69 6e 74 20 76 31 5f 5b 5d 20 3d 20 7b 32 2c  .int v1_[] = {2,
0a50: 32 2c 30 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72  2,0};...  vector
0a60: 20 3c 69 6e 74 3e 20 76 31 28 76 31 5f 2c 20 76   <int> v1(v1_, v
0a70: 31 5f 2b 73 69 7a 65 6f 66 28 76 31 5f 29 2f 73  1_+sizeof(v1_)/s
0a80: 69 7a 65 6f 66 28 2a 76 31 5f 29 29 3b 20 0d 0a  izeof(*v1_)); ..
0a90: 09 69 6e 74 20 70 32 5f 5b 5d 20 3d 20 7b 39 2c  .int p2_[] = {9,
0aa0: 30 2c 35 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72  0,5};...  vector
0ab0: 20 3c 69 6e 74 3e 20 70 32 28 70 32 5f 2c 20 70   <int> p2(p2_, p
0ac0: 32 5f 2b 73 69 7a 65 6f 66 28 70 32 5f 29 2f 73  2_+sizeof(p2_)/s
0ad0: 69 7a 65 6f 66 28 2a 70 32 5f 29 29 3b 20 0d 0a  izeof(*p2_)); ..
0ae0: 09 69 6e 74 20 76 32 5f 5b 5d 20 3d 20 7b 2d 32  .int v2_[] = {-2
0af0: 2c 32 2c 30 7d 3b 0d 0a 09 20 20 76 65 63 74 6f  ,2,0};...  vecto
0b00: 72 20 3c 69 6e 74 3e 20 76 32 28 76 32 5f 2c 20  r <int> v2(v2_, 
0b10: 76 32 5f 2b 73 69 7a 65 6f 66 28 76 32 5f 29 2f  v2_+sizeof(v2_)/
0b20: 73 69 7a 65 6f 66 28 2a 76 32 5f 29 29 3b 20 0d  sizeof(*v2_)); .
0b30: 0a 09 69 6e 74 20 52 20 3d 20 35 3b 20 0d 0a 09  ..int R = 5; ...
0b40: 73 74 72 69 6e 67 20 5f 20 3d 20 22 59 45 53 22  string _ = "YES"
0b50: 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 32 29  ; ..END..CASE(2)
0b60: 0d 0a 09 69 6e 74 20 70 31 5f 5b 5d 20 3d 20 7b  ...int p1_[] = {
0b70: 30 2c 30 2c 30 7d 3b 0d 0a 09 20 20 76 65 63 74  0,0,0};...  vect
0b80: 6f 72 20 3c 69 6e 74 3e 20 70 31 28 70 31 5f 2c  or <int> p1(p1_,
0b90: 20 70 31 5f 2b 73 69 7a 65 6f 66 28 70 31 5f 29   p1_+sizeof(p1_)
0ba0: 2f 73 69 7a 65 6f 66 28 2a 70 31 5f 29 29 3b 20  /sizeof(*p1_)); 
0bb0: 0d 0a 09 69 6e 74 20 76 31 5f 5b 5d 20 3d 20 7b  ...int v1_[] = {
0bc0: 2d 32 2c 32 2c 30 7d 3b 0d 0a 09 20 20 76 65 63  -2,2,0};...  vec
0bd0: 74 6f 72 20 3c 69 6e 74 3e 20 76 31 28 76 31 5f  tor <int> v1(v1_
0be0: 2c 20 76 31 5f 2b 73 69 7a 65 6f 66 28 76 31 5f  , v1_+sizeof(v1_
0bf0: 29 2f 73 69 7a 65 6f 66 28 2a 76 31 5f 29 29 3b  )/sizeof(*v1_));
0c00: 20 0d 0a 09 69 6e 74 20 70 32 5f 5b 5d 20 3d 20   ...int p2_[] = 
0c10: 7b 39 2c 30 2c 35 7d 3b 0d 0a 09 20 20 76 65 63  {9,0,5};...  vec
0c20: 74 6f 72 20 3c 69 6e 74 3e 20 70 32 28 70 32 5f  tor <int> p2(p2_
0c30: 2c 20 70 32 5f 2b 73 69 7a 65 6f 66 28 70 32 5f  , p2_+sizeof(p2_
0c40: 29 2f 73 69 7a 65 6f 66 28 2a 70 32 5f 29 29 3b  )/sizeof(*p2_));
0c50: 20 0d 0a 09 69 6e 74 20 76 32 5f 5b 5d 20 3d 20   ...int v2_[] = 
0c60: 7b 32 2c 32 2c 30 7d 3b 0d 0a 09 20 20 76 65 63  {2,2,0};...  vec
0c70: 74 6f 72 20 3c 69 6e 74 3e 20 76 32 28 76 32 5f  tor <int> v2(v2_
0c80: 2c 20 76 32 5f 2b 73 69 7a 65 6f 66 28 76 32 5f  , v2_+sizeof(v2_
0c90: 29 2f 73 69 7a 65 6f 66 28 2a 76 32 5f 29 29 3b  )/sizeof(*v2_));
0ca0: 20 0d 0a 09 69 6e 74 20 52 20 3d 20 35 3b 20 0d   ...int R = 5; .
0cb0: 0a 09 73 74 72 69 6e 67 20 5f 20 3d 20 22 4e 4f  ..string _ = "NO
0cc0: 22 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 33  "; ..END..CASE(3
0cd0: 29 0d 0a 09 69 6e 74 20 70 31 5f 5b 5d 20 3d 20  )...int p1_[] = 
0ce0: 7b 2d 32 38 33 38 2c 2d 37 39 34 30 2c 2d 32 39  {-2838,-7940,-29
0cf0: 33 36 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20  36};...  vector 
0d00: 3c 69 6e 74 3e 20 70 31 28 70 31 5f 2c 20 70 31  <int> p1(p1_, p1
0d10: 5f 2b 73 69 7a 65 6f 66 28 70 31 5f 29 2f 73 69  _+sizeof(p1_)/si
0d20: 7a 65 6f 66 28 2a 70 31 5f 29 29 3b 20 0d 0a 09  zeof(*p1_)); ...
0d30: 69 6e 74 20 76 31 5f 5b 5d 20 3d 20 7b 31 2c 31  int v1_[] = {1,1
0d40: 2c 2d 32 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72  ,-2};...  vector
0d50: 20 3c 69 6e 74 3e 20 76 31 28 76 31 5f 2c 20 76   <int> v1(v1_, v
0d60: 31 5f 2b 73 69 7a 65 6f 66 28 76 31 5f 29 2f 73  1_+sizeof(v1_)/s
0d70: 69 7a 65 6f 66 28 2a 76 31 5f 29 29 3b 20 0d 0a  izeof(*v1_)); ..
0d80: 09 69 6e 74 20 70 32 5f 5b 5d 20 3d 20 7b 35 33  .int p2_[] = {53
0d90: 32 2c 33 38 35 30 2c 39 35 39 30 7d 3b 0d 0a 09  2,3850,9590};...
0da0: 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 70    vector <int> p
0db0: 32 28 70 32 5f 2c 20 70 32 5f 2b 73 69 7a 65 6f  2(p2_, p2_+sizeo
0dc0: 66 28 70 32 5f 29 2f 73 69 7a 65 6f 66 28 2a 70  f(p2_)/sizeof(*p
0dd0: 32 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 76 32 5f  2_)); ...int v2_
0de0: 5b 5d 20 3d 20 7b 31 2c 30 2c 2d 33 7d 3b 0d 0a  [] = {1,0,-3};..
0df0: 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20  .  vector <int> 
0e00: 76 32 28 76 32 5f 2c 20 76 32 5f 2b 73 69 7a 65  v2(v2_, v2_+size
0e10: 6f 66 28 76 32 5f 29 2f 73 69 7a 65 6f 66 28 2a  of(v2_)/sizeof(*
0e20: 76 32 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 52 20  v2_)); ...int R 
0e30: 3d 20 33 34 31 30 3b 20 0d 0a 09 73 74 72 69 6e  = 3410; ...strin
0e40: 67 20 5f 20 3d 20 22 59 45 53 22 3b 20 0d 0a 45  g _ = "YES"; ..E
0e50: 4e 44 0d 0a 43 41 53 45 28 34 29 0d 0a 09 69 6e  ND..CASE(4)...in
0e60: 74 20 70 31 5f 5b 5d 20 3d 20 7b 2d 38 35 30 39  t p1_[] = {-8509
0e70: 2c 39 35 36 30 2c 33 34 35 7d 3b 0d 0a 09 20 20  ,9560,345};...  
0e80: 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 70 31 28  vector <int> p1(
0e90: 70 31 5f 2c 20 70 31 5f 2b 73 69 7a 65 6f 66 28  p1_, p1_+sizeof(
0ea0: 70 31 5f 29 2f 73 69 7a 65 6f 66 28 2a 70 31 5f  p1_)/sizeof(*p1_
0eb0: 29 29 3b 20 0d 0a 09 69 6e 74 20 76 31 5f 5b 5d  )); ...int v1_[]
0ec0: 20 3d 20 7b 2d 38 39 2c 2d 33 33 2c 36 32 7d 3b   = {-89,-33,62};
0ed0: 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74  ...  vector <int
0ee0: 3e 20 76 31 28 76 31 5f 2c 20 76 31 5f 2b 73 69  > v1(v1_, v1_+si
0ef0: 7a 65 6f 66 28 76 31 5f 29 2f 73 69 7a 65 6f 66  zeof(v1_)/sizeof
0f00: 28 2a 76 31 5f 29 29 3b 20 0d 0a 09 69 6e 74 20  (*v1_)); ...int 
0f10: 70 32 5f 5b 5d 20 3d 20 7b 2d 35 31 38 35 2c 2d  p2_[] = {-5185,-
0f20: 31 34 31 37 2c 32 38 34 36 7d 3b 0d 0a 09 20 20  1417,2846};...  
0f30: 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 70 32 28  vector <int> p2(
0f40: 70 32 5f 2c 20 70 32 5f 2b 73 69 7a 65 6f 66 28  p2_, p2_+sizeof(
0f50: 70 32 5f 29 2f 73 69 7a 65 6f 66 28 2a 70 32 5f  p2_)/sizeof(*p2_
0f60: 29 29 3b 20 0d 0a 09 69 6e 74 20 76 32 5f 5b 5d  )); ...int v2_[]
0f70: 20 3d 20 7b 2d 35 38 2c 32 34 2c 32 36 7d 3b 0d   = {-58,24,26};.
0f80: 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e  ..  vector <int>
0f90: 20 76 32 28 76 32 5f 2c 20 76 32 5f 2b 73 69 7a   v2(v2_, v2_+siz
0fa0: 65 6f 66 28 76 32 5f 29 2f 73 69 7a 65 6f 66 28  eof(v2_)/sizeof(
0fb0: 2a 76 32 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 52  *v2_)); ...int R
0fc0: 20 3d 20 38 33 34 34 3b 20 0d 0a 09 73 74 72 69   = 8344; ...stri
0fd0: 6e 67 20 5f 20 3d 20 22 59 45 53 22 3b 20 0d 0a  ng _ = "YES"; ..
0fe0: 45 4e 44 0d 0a 43 41 53 45 28 35 29 0d 0a 09 69  END..CASE(5)...i
0ff0: 6e 74 20 70 31 5f 5b 5d 20 3d 20 7b 2d 37 31 36  nt p1_[] = {-716
1000: 33 2c 2d 33 37 31 2c 2d 32 34 35 39 7d 3b 0d 0a  3,-371,-2459};..
1010: 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20  .  vector <int> 
1020: 70 31 28 70 31 5f 2c 20 70 31 5f 2b 73 69 7a 65  p1(p1_, p1_+size
1030: 6f 66 28 70 31 5f 29 2f 73 69 7a 65 6f 66 28 2a  of(p1_)/sizeof(*
1040: 70 31 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 76 31  p1_)); ...int v1
1050: 5f 5b 5d 20 3d 20 7b 2d 35 39 2c 2d 34 31 2c 2d  _[] = {-59,-41,-
1060: 31 34 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20  14};...  vector 
1070: 3c 69 6e 74 3e 20 76 31 28 76 31 5f 2c 20 76 31  <int> v1(v1_, v1
1080: 5f 2b 73 69 7a 65 6f 66 28 76 31 5f 29 2f 73 69  _+sizeof(v1_)/si
1090: 7a 65 6f 66 28 2a 76 31 5f 29 29 3b 20 0d 0a 09  zeof(*v1_)); ...
10a0: 69 6e 74 20 70 32 5f 5b 5d 20 3d 20 7b 2d 32 33  int p2_[] = {-23
10b0: 39 38 2c 2d 34 32 36 2c 2d 35 34 38 37 7d 3b 0d  98,-426,-5487};.
10c0: 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e  ..  vector <int>
10d0: 20 70 32 28 70 32 5f 2c 20 70 32 5f 2b 73 69 7a   p2(p2_, p2_+siz
10e0: 65 6f 66 28 70 32 5f 29 2f 73 69 7a 65 6f 66 28  eof(p2_)/sizeof(
10f0: 2a 70 32 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 76  *p2_)); ...int v
1100: 32 5f 5b 5d 20 3d 20 7b 2d 34 33 2c 32 37 2c 36  2_[] = {-43,27,6
1110: 37 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c  7};...  vector <
1120: 69 6e 74 3e 20 76 32 28 76 32 5f 2c 20 76 32 5f  int> v2(v2_, v2_
1130: 2b 73 69 7a 65 6f 66 28 76 32 5f 29 2f 73 69 7a  +sizeof(v2_)/siz
1140: 65 6f 66 28 2a 76 32 5f 29 29 3b 20 0d 0a 09 69  eof(*v2_)); ...i
1150: 6e 74 20 52 20 3d 20 35 34 31 30 3b 20 0d 0a 09  nt R = 5410; ...
1160: 73 74 72 69 6e 67 20 5f 20 3d 20 22 4e 4f 22 3b  string _ = "NO";
1170: 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 36 29 0d   ..END..CASE(6).
1180: 0a 09 69 6e 74 20 70 31 5f 5b 5d 20 3d 20 7b 31  ..int p1_[] = {1
1190: 37 37 34 2c 2d 34 34 39 31 2c 37 38 31 30 7d 3b  774,-4491,7810};
11a0: 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74  ...  vector <int
11b0: 3e 20 70 31 28 70 31 5f 2c 20 70 31 5f 2b 73 69  > p1(p1_, p1_+si
11c0: 7a 65 6f 66 28 70 31 5f 29 2f 73 69 7a 65 6f 66  zeof(p1_)/sizeof
11d0: 28 2a 70 31 5f 29 29 3b 20 0d 0a 09 69 6e 74 20  (*p1_)); ...int 
11e0: 76 31 5f 5b 5d 20 3d 20 7b 2d 31 32 2c 31 39 2c  v1_[] = {-12,19,
11f0: 2d 32 34 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72  -24};...  vector
1200: 20 3c 69 6e 74 3e 20 76 31 28 76 31 5f 2c 20 76   <int> v1(v1_, v
1210: 31 5f 2b 73 69 7a 65 6f 66 28 76 31 5f 29 2f 73  1_+sizeof(v1_)/s
1220: 69 7a 65 6f 66 28 2a 76 31 5f 29 29 3b 20 0d 0a  izeof(*v1_)); ..
1230: 09 69 6e 74 20 70 32 5f 5b 5d 20 3d 20 7b 32 33  .int p2_[] = {23
1240: 32 32 2c 33 37 39 33 2c 39 38 39 37 7d 3b 0d 0a  22,3793,9897};..
1250: 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20  .  vector <int> 
1260: 70 32 28 70 32 5f 2c 20 70 32 5f 2b 73 69 7a 65  p2(p2_, p2_+size
1270: 6f 66 28 70 32 5f 29 2f 73 69 7a 65 6f 66 28 2a  of(p2_)/sizeof(*
1280: 70 32 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 76 32  p2_)); ...int v2
1290: 5f 5b 5d 20 3d 20 7b 2d 31 32 2c 31 39 2c 2d 32  _[] = {-12,19,-2
12a0: 34 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c  4};...  vector <
12b0: 69 6e 74 3e 20 76 32 28 76 32 5f 2c 20 76 32 5f  int> v2(v2_, v2_
12c0: 2b 73 69 7a 65 6f 66 28 76 32 5f 29 2f 73 69 7a  +sizeof(v2_)/siz
12d0: 65 6f 66 28 2a 76 32 5f 29 29 3b 20 0d 0a 09 69  eof(*v2_)); ...i
12e0: 6e 74 20 52 20 3d 20 31 30 30 30 30 3b 20 0d 0a  nt R = 10000; ..
12f0: 09 73 74 72 69 6e 67 20 5f 20 3d 20 22 59 45 53  .string _ = "YES
1300: 22 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 37  "; ..END..CASE(7
1310: 29 0d 0a 09 69 6e 74 20 70 31 5f 5b 5d 20 3d 20  )...int p1_[] = 
1320: 7b 33 37 33 31 2c 38 35 33 37 2c 35 36 36 31 7d  {3731,8537,5661}
1330: 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e  ;...  vector <in
1340: 74 3e 20 70 31 28 70 31 5f 2c 20 70 31 5f 2b 73  t> p1(p1_, p1_+s
1350: 69 7a 65 6f 66 28 70 31 5f 29 2f 73 69 7a 65 6f  izeof(p1_)/sizeo
1360: 66 28 2a 70 31 5f 29 29 3b 20 0d 0a 09 69 6e 74  f(*p1_)); ...int
1370: 20 76 31 5f 5b 5d 20 3d 20 7b 2d 37 30 2c 37 31   v1_[] = {-70,71
1380: 2c 33 32 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72  ,32};...  vector
1390: 20 3c 69 6e 74 3e 20 76 31 28 76 31 5f 2c 20 76   <int> v1(v1_, v
13a0: 31 5f 2b 73 69 7a 65 6f 66 28 76 31 5f 29 2f 73  1_+sizeof(v1_)/s
13b0: 69 7a 65 6f 66 28 2a 76 31 5f 29 29 3b 20 0d 0a  izeof(*v1_)); ..
13c0: 09 69 6e 74 20 70 32 5f 5b 5d 20 3d 20 7b 38 37  .int p2_[] = {87
13d0: 30 31 2c 2d 31 38 38 36 2c 2d 35 31 31 35 7d 3b  01,-1886,-5115};
13e0: 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74  ...  vector <int
13f0: 3e 20 70 32 28 70 32 5f 2c 20 70 32 5f 2b 73 69  > p2(p2_, p2_+si
1400: 7a 65 6f 66 28 70 32 5f 29 2f 73 69 7a 65 6f 66  zeof(p2_)/sizeof
1410: 28 2a 70 32 5f 29 29 3b 20 0d 0a 09 69 6e 74 20  (*p2_)); ...int 
1420: 76 32 5f 5b 5d 20 3d 20 7b 32 38 2c 2d 31 33 2c  v2_[] = {28,-13,
1430: 37 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c  7};...  vector <
1440: 69 6e 74 3e 20 76 32 28 76 32 5f 2c 20 76 32 5f  int> v2(v2_, v2_
1450: 2b 73 69 7a 65 6f 66 28 76 32 5f 29 2f 73 69 7a  +sizeof(v2_)/siz
1460: 65 6f 66 28 2a 76 32 5f 29 29 3b 20 0d 0a 09 69  eof(*v2_)); ...i
1470: 6e 74 20 52 20 3d 20 39 37 36 36 3b 20 0d 0a 09  nt R = 9766; ...
1480: 73 74 72 69 6e 67 20 5f 20 3d 20 22 4e 4f 22 3b  string _ = "NO";
1490: 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 38 29 0d   ..END..CASE(8).
14a0: 0a 09 69 6e 74 20 70 31 5f 5b 5d 20 3d 20 7b 30  ..int p1_[] = {0
14b0: 2c 30 2c 30 7d 3b 0d 0a 09 20 20 76 65 63 74 6f  ,0,0};...  vecto
14c0: 72 20 3c 69 6e 74 3e 20 70 31 28 70 31 5f 2c 20  r <int> p1(p1_, 
14d0: 70 31 5f 2b 73 69 7a 65 6f 66 28 70 31 5f 29 2f  p1_+sizeof(p1_)/
14e0: 73 69 7a 65 6f 66 28 2a 70 31 5f 29 29 3b 20 0d  sizeof(*p1_)); .
14f0: 0a 09 69 6e 74 20 76 31 5f 5b 5d 20 3d 20 7b 30  ..int v1_[] = {0
1500: 2c 30 2c 30 7d 3b 0d 0a 09 20 20 76 65 63 74 6f  ,0,0};...  vecto
1510: 72 20 3c 69 6e 74 3e 20 76 31 28 76 31 5f 2c 20  r <int> v1(v1_, 
1520: 76 31 5f 2b 73 69 7a 65 6f 66 28 76 31 5f 29 2f  v1_+sizeof(v1_)/
1530: 73 69 7a 65 6f 66 28 2a 76 31 5f 29 29 3b 20 0d  sizeof(*v1_)); .
1540: 0a 09 69 6e 74 20 70 32 5f 5b 5d 20 3d 20 7b 31  ..int p2_[] = {1
1550: 2c 30 2c 30 7d 3b 0d 0a 09 20 20 76 65 63 74 6f  ,0,0};...  vecto
1560: 72 20 3c 69 6e 74 3e 20 70 32 28 70 32 5f 2c 20  r <int> p2(p2_, 
1570: 70 32 5f 2b 73 69 7a 65 6f 66 28 70 32 5f 29 2f  p2_+sizeof(p2_)/
1580: 73 69 7a 65 6f 66 28 2a 70 32 5f 29 29 3b 20 0d  sizeof(*p2_)); .
1590: 0a 09 69 6e 74 20 76 32 5f 5b 5d 20 3d 20 7b 30  ..int v2_[] = {0
15a0: 2c 30 2c 30 7d 3b 0d 0a 09 20 20 76 65 63 74 6f  ,0,0};...  vecto
15b0: 72 20 3c 69 6e 74 3e 20 76 32 28 76 32 5f 2c 20  r <int> v2(v2_, 
15c0: 76 32 5f 2b 73 69 7a 65 6f 66 28 76 32 5f 29 2f  v2_+sizeof(v2_)/
15d0: 73 69 7a 65 6f 66 28 2a 76 32 5f 29 29 3b 20 0d  sizeof(*v2_)); .
15e0: 0a 09 69 6e 74 20 52 20 3d 20 31 3b 20 0d 0a 09  ..int R = 1; ...
15f0: 73 74 72 69 6e 67 20 5f 20 3d 20 22 59 45 53 22  string _ = "YES"
1600: 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 39 29  ; ..END..CASE(9)
1610: 0d 0a 09 69 6e 74 20 70 31 5f 5b 5d 20 3d 20 7b  ...int p1_[] = {
1620: 30 2c 30 2c 30 7d 3b 0d 0a 09 20 20 76 65 63 74  0,0,0};...  vect
1630: 6f 72 20 3c 69 6e 74 3e 20 70 31 28 70 31 5f 2c  or <int> p1(p1_,
1640: 20 70 31 5f 2b 73 69 7a 65 6f 66 28 70 31 5f 29   p1_+sizeof(p1_)
1650: 2f 73 69 7a 65 6f 66 28 2a 70 31 5f 29 29 3b 20  /sizeof(*p1_)); 
1660: 0d 0a 09 69 6e 74 20 76 31 5f 5b 5d 20 3d 20 7b  ...int v1_[] = {
1670: 30 2c 30 2c 30 7d 3b 0d 0a 09 20 20 76 65 63 74  0,0,0};...  vect
1680: 6f 72 20 3c 69 6e 74 3e 20 76 31 28 76 31 5f 2c  or <int> v1(v1_,
1690: 20 76 31 5f 2b 73 69 7a 65 6f 66 28 76 31 5f 29   v1_+sizeof(v1_)
16a0: 2f 73 69 7a 65 6f 66 28 2a 76 31 5f 29 29 3b 20  /sizeof(*v1_)); 
16b0: 0d 0a 09 69 6e 74 20 70 32 5f 5b 5d 20 3d 20 7b  ...int p2_[] = {
16c0: 32 2c 30 2c 30 7d 3b 0d 0a 09 20 20 76 65 63 74  2,0,0};...  vect
16d0: 6f 72 20 3c 69 6e 74 3e 20 70 32 28 70 32 5f 2c  or <int> p2(p2_,
16e0: 20 70 32 5f 2b 73 69 7a 65 6f 66 28 70 32 5f 29   p2_+sizeof(p2_)
16f0: 2f 73 69 7a 65 6f 66 28 2a 70 32 5f 29 29 3b 20  /sizeof(*p2_)); 
1700: 0d 0a 09 69 6e 74 20 76 32 5f 5b 5d 20 3d 20 7b  ...int v2_[] = {
1710: 30 2c 30 2c 30 7d 3b 0d 0a 09 20 20 76 65 63 74  0,0,0};...  vect
1720: 6f 72 20 3c 69 6e 74 3e 20 76 32 28 76 32 5f 2c  or <int> v2(v2_,
1730: 20 76 32 5f 2b 73 69 7a 65 6f 66 28 76 32 5f 29   v2_+sizeof(v2_)
1740: 2f 73 69 7a 65 6f 66 28 2a 76 32 5f 29 29 3b 20  /sizeof(*v2_)); 
1750: 0d 0a 09 69 6e 74 20 52 20 3d 20 31 3b 20 0d 0a  ...int R = 1; ..
1760: 09 73 74 72 69 6e 67 20 5f 20 3d 20 22 4e 4f 22  .string _ = "NO"
1770: 3b 20 0d 0a 45 4e 44 0d 0a 0d 0a 7d 0d 0a 2f 2f  ; ..END....}..//
1780: 20 45 4e 44 20 43 55 54 20 48 45 52 45 0d 0a      END CUT HERE..