#include using namespace std; typedef signed long long ll; #undef _P #define _P(...) (void)printf(__VA_ARGS__) #define FOR(x,to) for(x=0;x<(to);x++) #define FORR(x,arr) for(auto& x:arr) #define ITR(x,c) for(__typeof(c.begin()) x=c.begin();x!=c.end();x++) #define ALL(a) (a.begin()),(a.end()) #define ZERO(a) memset(a,0,sizeof(a)) #define MINUS(a) memset(a,0xff,sizeof(a)) //------------------------------------------------------- int N; int T; double P,W,V,GX,GY; pair pre(pair p) { double d=(1+V)*(1+V)+W*W; double x2=(1+V)*p.first+W*p.second; double y2=-W*p.first+(1+V)*p.second; return {x2/d,y2/d}; } pair nex(pair p) { double x2=(1+V)*p.first-W*p.second; double y2=W*p.first+(1+V)*p.second; return {x2,y2}; } void solve() { int i,j,k,l,r,x,y; string s; cin>>N; while(N--) { cin>>T>>P>>W>>V>>GX>>GY; pair A[16],B[16]; A[0]={1,0}; B[T]={GX,GY}; FOR(i,T) A[i+1]=nex(A[i]); for(i=T-1;i>=0;i--) B[i]=pre(B[i+1]); for(i=1;i<=T;i++) { double dx=B[i].first-A[i].first; double dy=B[i].second-A[i].second; if(dx*dx+dy*dy<=P) { for(j=1;j<=T;j++) { if(j!=i) _P("0 0\n"); else _P("%.12lf %.12lf\n",dx,dy); } break; } } if(i<=T) continue; double rs=hypot(A[0].first,A[0].second); double rt=hypot(GX,GY); double ds=atan2(A[0].first,A[0].second); double dt=atan2(GX,GY); double scale=sqrt((1+V)*(1+V)+W*W); double deg=atan2(1+V,W); rt=rt/rs/pow(scale,T-1); dt-=ds+deg*(T-1); // scale first for(x=1;x<=T;x++) for(y=1;y