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 <functional> 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: #include <cstring> 4fd800b3a8 2011-02-23 kinaba: using namespace std; 4fd800b3a8 2011-02-23 kinaba: typedef long long LL; 4fd800b3a8 2011-02-23 kinaba: typedef complex<double> CMP; 4fd800b3a8 2011-02-23 kinaba: 4fd800b3a8 2011-02-23 kinaba: class AfraidOfEven { public: 4fd800b3a8 2011-02-23 kinaba: vector <int> restoreProgression(vector <int> seq) 4fd800b3a8 2011-02-23 kinaba: { 4fd800b3a8 2011-02-23 kinaba: LL a = seq[0]; 4fd800b3a8 2011-02-23 kinaba: LL b = seq[1]; 4fd800b3a8 2011-02-23 kinaba: for(LL aa=a; aa<=0x7fffffff; aa*=2) 4fd800b3a8 2011-02-23 kinaba: for(LL bb=b; bb<=0x7fffffff; bb*=2) 4fd800b3a8 2011-02-23 kinaba: if( possible(seq,aa,bb) ) 4fd800b3a8 2011-02-23 kinaba: return seq; 4fd800b3a8 2011-02-23 kinaba: assert(false); 4fd800b3a8 2011-02-23 kinaba: } 4fd800b3a8 2011-02-23 kinaba: bool possible(vector<int>& seq, LL a, LL b) 4fd800b3a8 2011-02-23 kinaba: { 4fd800b3a8 2011-02-23 kinaba: LL d = b-a; 4fd800b3a8 2011-02-23 kinaba: vector<LL> q(seq.size()); 4fd800b3a8 2011-02-23 kinaba: q[0] = a; 4fd800b3a8 2011-02-23 kinaba: q[1] = b; 4fd800b3a8 2011-02-23 kinaba: for(int i=2; i<seq.size(); ++i) 4fd800b3a8 2011-02-23 kinaba: { 4fd800b3a8 2011-02-23 kinaba: q[i] = q[i-1] + d; 4fd800b3a8 2011-02-23 kinaba: if( q[i] <= 0 ) 4fd800b3a8 2011-02-23 kinaba: return false; 4fd800b3a8 2011-02-23 kinaba: if( q[i] % seq[i] != 0 ) 4fd800b3a8 2011-02-23 kinaba: return false; 4fd800b3a8 2011-02-23 kinaba: LL f = q[i] / seq[i]; 4fd800b3a8 2011-02-23 kinaba: if( f & (f-1) ) 4fd800b3a8 2011-02-23 kinaba: return false; 4fd800b3a8 2011-02-23 kinaba: } 4fd800b3a8 2011-02-23 kinaba: for(int i=0; i<seq.size(); ++i) 4fd800b3a8 2011-02-23 kinaba: seq[i] = int(q[i]); 4fd800b3a8 2011-02-23 kinaba: return true; 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() 4fd800b3a8 2011-02-23 kinaba: { ostringstream os; os << " (" << int((clock()-start_time)/CLOCKS_PER_SEC*1000) << " msec)"; return os.str(); } 4fd800b3a8 2011-02-23 kinaba: template<typename T> ostream& operator<<(ostream& os, const vector<T>& v) 4fd800b3a8 2011-02-23 kinaba: { os << "{ "; 4fd800b3a8 2011-02-23 kinaba: for(typename vector<T>::const_iterator it=v.begin(); it!=v.end(); ++it) 4fd800b3a8 2011-02-23 kinaba: os << '\"' << *it << '\"' << (it+1==v.end() ? "" : ", "); os << " }"; return os; } 4fd800b3a8 2011-02-23 kinaba: void verify_case(const vector <int>& Expected, const vector <int>& Received) { 4fd800b3a8 2011-02-23 kinaba: bool ok = (Expected == Received); 4fd800b3a8 2011-02-23 kinaba: if(ok) cerr << "PASSED" << timer() << endl; else { cerr << "FAILED" << timer() << endl; 4fd800b3a8 2011-02-23 kinaba: cerr << "\to: " << Expected << endl << "\tx: " << Received << endl; } } 4fd800b3a8 2011-02-23 kinaba: #define CASE(N) {cerr << "Test Case #" << N << "..." << flush; start_time=clock(); 4fd800b3a8 2011-02-23 kinaba: #define END verify_case(_, AfraidOfEven().restoreProgression(seq));} 4fd800b3a8 2011-02-23 kinaba: int main(){ 4fd800b3a8 2011-02-23 kinaba: 4fd800b3a8 2011-02-23 kinaba: CASE(0) 4fd800b3a8 2011-02-23 kinaba: int seq_[] = {1, 1, 3, 1, 5}; 4fd800b3a8 2011-02-23 kinaba: vector <int> seq(seq_, seq_+sizeof(seq_)/sizeof(*seq_)); 4fd800b3a8 2011-02-23 kinaba: int __[] = {1, 2, 3, 4, 5 }; 4fd800b3a8 2011-02-23 kinaba: vector <int> _(__, __+sizeof(__)/sizeof(*__)); 4fd800b3a8 2011-02-23 kinaba: END 4fd800b3a8 2011-02-23 kinaba: CASE(1) 4fd800b3a8 2011-02-23 kinaba: int seq_[] = {9, 7, 5, 3, 1}; 4fd800b3a8 2011-02-23 kinaba: vector <int> seq(seq_, seq_+sizeof(seq_)/sizeof(*seq_)); 4fd800b3a8 2011-02-23 kinaba: int __[] = {9, 7, 5, 3, 1 }; 4fd800b3a8 2011-02-23 kinaba: vector <int> _(__, __+sizeof(__)/sizeof(*__)); 4fd800b3a8 2011-02-23 kinaba: END 4fd800b3a8 2011-02-23 kinaba: CASE(2) 4fd800b3a8 2011-02-23 kinaba: int seq_[] = {999, 999, 999, 999}; 4fd800b3a8 2011-02-23 kinaba: vector <int> seq(seq_, seq_+sizeof(seq_)/sizeof(*seq_)); 4fd800b3a8 2011-02-23 kinaba: int __[] = {999, 999, 999, 999 }; 4fd800b3a8 2011-02-23 kinaba: vector <int> _(__, __+sizeof(__)/sizeof(*__)); 4fd800b3a8 2011-02-23 kinaba: END 4fd800b3a8 2011-02-23 kinaba: CASE(3) 4fd800b3a8 2011-02-23 kinaba: int seq_[] = {7, 47, 5, 113, 73, 179, 53}; 4fd800b3a8 2011-02-23 kinaba: vector <int> seq(seq_, seq_+sizeof(seq_)/sizeof(*seq_)); 4fd800b3a8 2011-02-23 kinaba: int __[] = {14, 47, 80, 113, 146, 179, 212 }; 4fd800b3a8 2011-02-23 kinaba: vector <int> _(__, __+sizeof(__)/sizeof(*__)); 4fd800b3a8 2011-02-23 kinaba: END 4fd800b3a8 2011-02-23 kinaba: CASE(4) 4fd800b3a8 2011-02-23 kinaba: int seq_[] = {749, 999, 125, 1}; 4fd800b3a8 2011-02-23 kinaba: vector <int> seq(seq_, seq_+sizeof(seq_)/sizeof(*seq_)); 4fd800b3a8 2011-02-23 kinaba: int __[] = {1498, 999, 500, 1 }; 4fd800b3a8 2011-02-23 kinaba: vector <int> _(__, __+sizeof(__)/sizeof(*__)); 4fd800b3a8 2011-02-23 kinaba: END 4fd800b3a8 2011-02-23 kinaba: CASE(5) 4fd800b3a8 2011-02-23 kinaba: int seq_[] = {3,7,1,9}; 4fd800b3a8 2011-02-23 kinaba: vector <int> seq(seq_, seq_+sizeof(seq_)/sizeof(*seq_)); 4fd800b3a8 2011-02-23 kinaba: int __[] = {6,7,8,9}; 4fd800b3a8 2011-02-23 kinaba: vector <int> _(__, __+sizeof(__)/sizeof(*__)); 4fd800b3a8 2011-02-23 kinaba: END 4fd800b3a8 2011-02-23 kinaba: /* 4fd800b3a8 2011-02-23 kinaba: CASE(6) 4fd800b3a8 2011-02-23 kinaba: int seq_[] = ; 4fd800b3a8 2011-02-23 kinaba: vector <int> seq(seq_, seq_+sizeof(seq_)/sizeof(*seq_)); 4fd800b3a8 2011-02-23 kinaba: int __[] = ; 4fd800b3a8 2011-02-23 kinaba: vector <int> _(__, __+sizeof(__)/sizeof(*__)); 4fd800b3a8 2011-02-23 kinaba: END 4fd800b3a8 2011-02-23 kinaba: */ 4fd800b3a8 2011-02-23 kinaba: } 4fd800b3a8 2011-02-23 kinaba: // END CUT HERE