4fd800b3a8 2011-02-23 kinaba: #include <iostream> 4fd800b3a8 2011-02-23 kinaba: #include <sstream> 4fd800b3a8 2011-02-23 kinaba: #include <iomanip> 4fd800b3a8 2011-02-23 kinaba: #include <vector> 4fd800b3a8 2011-02-23 kinaba: #include <string> 4fd800b3a8 2011-02-23 kinaba: #include <map> 4fd800b3a8 2011-02-23 kinaba: #include <set> 4fd800b3a8 2011-02-23 kinaba: #include <algorithm> 4fd800b3a8 2011-02-23 kinaba: #include <numeric> 4fd800b3a8 2011-02-23 kinaba: #include <iterator> 4fd800b3a8 2011-02-23 kinaba: #include <complex> 4fd800b3a8 2011-02-23 kinaba: #include <queue> 4fd800b3a8 2011-02-23 kinaba: #include <stack> 4fd800b3a8 2011-02-23 kinaba: #include <cmath> 4fd800b3a8 2011-02-23 kinaba: #include <cassert> 4fd800b3a8 2011-02-23 kinaba: using namespace std; 4fd800b3a8 2011-02-23 kinaba: typedef long long LL; 4fd800b3a8 2011-02-23 kinaba: 4fd800b3a8 2011-02-23 kinaba: class LaserShooting 4fd800b3a8 2011-02-23 kinaba: { 4fd800b3a8 2011-02-23 kinaba: public: 4fd800b3a8 2011-02-23 kinaba: double numberOfHits(vector <int> x, vector <int> y1, vector <int> y2) 4fd800b3a8 2011-02-23 kinaba: { 4fd800b3a8 2011-02-23 kinaba: double ans = 0.0; 4fd800b3a8 2011-02-23 kinaba: for(int i=0; i<x.size(); ++i) 4fd800b3a8 2011-02-23 kinaba: { 4fd800b3a8 2011-02-23 kinaba: double a1 = atan2(min(y1[i],y2[i]), x[i]); 4fd800b3a8 2011-02-23 kinaba: double a2 = atan2(max(y1[i],y2[i]), x[i]); 4fd800b3a8 2011-02-23 kinaba: ans += (a2-a1) / M_PI; 4fd800b3a8 2011-02-23 kinaba: } 4fd800b3a8 2011-02-23 kinaba: return ans; 4fd800b3a8 2011-02-23 kinaba: } 4fd800b3a8 2011-02-23 kinaba: }; 4fd800b3a8 2011-02-23 kinaba: 4fd800b3a8 2011-02-23 kinaba: // BEGIN CUT HERE 4fd800b3a8 2011-02-23 kinaba: #include <ctime> 4fd800b3a8 2011-02-23 kinaba: double start_time;string timer() { ostringstream os; os << " (" << int((clock()-start_time)/CLOCKS_PER_SEC*1000) << " msec)"; return os.str(); } 4fd800b3a8 2011-02-23 kinaba: 4fd800b3a8 2011-02-23 kinaba: template <typename T> string print_array(const vector<T> &V) { ostringstream os; os << "{ "; for (typename vector<T>::const_iterator iter = V.begin(); iter != V.end(); ++iter) os << '\"' << *iter << "\","; os << " }"; return os.str(); } 4fd800b3a8 2011-02-23 kinaba: int verify_case(const double &Expected, const double &Received) { double diff = Expected - Received; if (diff < 0) diff = -diff; if (diff < 1e-9) cerr << "PASSED" << timer() << endl; else { cerr << "FAILED" << timer() << endl; cerr << "\tExpected: \"" << Expected << '\"' << endl; cerr << "\tReceived: \"" << Received << '\"' << endl; } return 0;} 4fd800b3a8 2011-02-23 kinaba: 4fd800b3a8 2011-02-23 kinaba: template<int N> struct Case_ { Case_(){start_time=clock();} }; 4fd800b3a8 2011-02-23 kinaba: char Test_(...); 4fd800b3a8 2011-02-23 kinaba: int Test_(Case_<0>) { 4fd800b3a8 2011-02-23 kinaba: int x_[] = {1}; 4fd800b3a8 2011-02-23 kinaba: vector <int> x(x_, x_+sizeof(x_)/sizeof(*x_)); 4fd800b3a8 2011-02-23 kinaba: int y1_[] = {-1}; 4fd800b3a8 2011-02-23 kinaba: vector <int> y1(y1_, y1_+sizeof(y1_)/sizeof(*y1_)); 4fd800b3a8 2011-02-23 kinaba: int y2_[] = {1}; 4fd800b3a8 2011-02-23 kinaba: vector <int> y2(y2_, y2_+sizeof(y2_)/sizeof(*y2_)); 4fd800b3a8 2011-02-23 kinaba: double RetVal = 0.5; 4fd800b3a8 2011-02-23 kinaba: return verify_case(RetVal, LaserShooting().numberOfHits(x, y1, y2)); } 4fd800b3a8 2011-02-23 kinaba: int Test_(Case_<1>) { 4fd800b3a8 2011-02-23 kinaba: int x_[] = {1,2}; 4fd800b3a8 2011-02-23 kinaba: vector <int> x(x_, x_+sizeof(x_)/sizeof(*x_)); 4fd800b3a8 2011-02-23 kinaba: int y1_[] = {-1,-2}; 4fd800b3a8 2011-02-23 kinaba: vector <int> y1(y1_, y1_+sizeof(y1_)/sizeof(*y1_)); 4fd800b3a8 2011-02-23 kinaba: int y2_[] = {1,2}; 4fd800b3a8 2011-02-23 kinaba: vector <int> y2(y2_, y2_+sizeof(y2_)/sizeof(*y2_)); 4fd800b3a8 2011-02-23 kinaba: double RetVal = 1.0; 4fd800b3a8 2011-02-23 kinaba: return verify_case(RetVal, LaserShooting().numberOfHits(x, y1, y2)); } 4fd800b3a8 2011-02-23 kinaba: int Test_(Case_<2>) { 4fd800b3a8 2011-02-23 kinaba: int x_[] = {3,4,7,1}; 4fd800b3a8 2011-02-23 kinaba: vector <int> x(x_, x_+sizeof(x_)/sizeof(*x_)); 4fd800b3a8 2011-02-23 kinaba: int y1_[] = {1,2,3,4}; 4fd800b3a8 2011-02-23 kinaba: vector <int> y1(y1_, y1_+sizeof(y1_)/sizeof(*y1_)); 4fd800b3a8 2011-02-23 kinaba: int y2_[] = {4,3,2,1}; 4fd800b3a8 2011-02-23 kinaba: vector <int> y2(y2_, y2_+sizeof(y2_)/sizeof(*y2_)); 4fd800b3a8 2011-02-23 kinaba: double RetVal = 0.4623163952488826; 4fd800b3a8 2011-02-23 kinaba: return verify_case(RetVal, LaserShooting().numberOfHits(x, y1, y2)); } 4fd800b3a8 2011-02-23 kinaba: int Test_(Case_<3>) { 4fd800b3a8 2011-02-23 kinaba: int x_[] = {134,298,151,942}; 4fd800b3a8 2011-02-23 kinaba: vector <int> x(x_, x_+sizeof(x_)/sizeof(*x_)); 4fd800b3a8 2011-02-23 kinaba: int y1_[] = {-753,-76,19,568}; 4fd800b3a8 2011-02-23 kinaba: vector <int> y1(y1_, y1_+sizeof(y1_)/sizeof(*y1_)); 4fd800b3a8 2011-02-23 kinaba: int y2_[] = {440,689,-39,672}; 4fd800b3a8 2011-02-23 kinaba: vector <int> y2(y2_, y2_+sizeof(y2_)/sizeof(*y2_)); 4fd800b3a8 2011-02-23 kinaba: double RetVal = 1.444210260641501; 4fd800b3a8 2011-02-23 kinaba: return verify_case(RetVal, LaserShooting().numberOfHits(x, y1, y2)); } 4fd800b3a8 2011-02-23 kinaba: 4fd800b3a8 2011-02-23 kinaba: template<int N> void Run_() { cerr << "Test Case #" << N << "..." << flush; Test_(Case_<N>()); Run_<sizeof(Test_(Case_<N+1>()))==1 ? -1 : N+1>(); } 4fd800b3a8 2011-02-23 kinaba: template<> void Run_<-1>() {} 4fd800b3a8 2011-02-23 kinaba: int main() { Run_<0>(); } 4fd800b3a8 2011-02-23 kinaba: // END CUT HERE 4fd800b3a8 2011-02-23 kinaba: