#pragma region include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include //#include #include #pragma endregion //#include ///////// #define REP(i, x, n) for(int i = x; i < n; ++i) #define rep(i,n) REP(i,0,n) #define ALL(X) X.begin(), X.end() ///////// #pragma region typedef typedef long long LL; typedef long double LD; typedef unsigned long long ULL; typedef std::pair PLL;// typedef std::pair PII;// #pragma endregion //typedef ////定数 const int INF = (int)1e9; const LL MOD = (LL)1e9+7; const LL LINF = (LL)4e18+20; const LD PI = acos(-1.0); const double EPS = 1e-9; ///////// using namespace::std; void solve(){ int N; cin >> N; vector< vector > base(N),tree(N); map > edgeUV; for(int i=0;i> a >> b; base[a].push_back( b ); base[b].push_back( a ); vector temp(2); temp[0] = min(a,b); temp[1] = max(a,b); edgeUV[i] = temp; } vector oya(N,-1);//i頂点の親 vector edgeC(N,-1);//i番目の辺の子方向の頂点番号 vector use(N,false); queue que; que.push( 0 ); map< vector,int > edgeNo;//(u,v)->No while( que.size() ){ int v = que.front(); que.pop(); if( use[v] ){ continue; } use[v] = true; int size = base[v].size(); for(int i=0;i dat(2); dat[0] = min(u,v); dat[1] = max(u,v); edgeNo[dat] = u; que.push( u ); } } vector< vector< vector > > X(N,vector< vector >(2,vector(2,0))); for(int i=0;i> Q; while(Q--){ cin >> type; if( type == 'x' ){ int pos,x00,x01,x10,x11; cin >> pos >> x00 >> x01 >> x10 >> x11; vector temp = edgeUV[pos]; if( oya[ temp[0] ] == temp[1] ){ pos = temp[0]; }else{ pos = temp[1]; } X[pos][0][0] = x00; X[pos][0][1] = x01; X[pos][1][0] = x10; X[pos][1][1] = x11; }else{ int to,from; cin >> to >> from; vector< vector > ret(2,vector(2,0)); vector< vector > temp(2,vector(2,0)); ret[0][0] = 1; ret[1][1] = 1; int pos = from; while(pos != to){ //ret = X[pos] * ret; if(true){ temp[0][0] = X[pos][0][0]*ret[0][0]+X[pos][0][1]*ret[1][0]; temp[0][1] = X[pos][0][0]*ret[0][1]+X[pos][0][1]*ret[1][1]; temp[1][0] = X[pos][1][0]*ret[0][0]+X[pos][1][1]*ret[1][0]; temp[1][1] = X[pos][1][0]*ret[0][1]+X[pos][1][1]*ret[1][1]; if(temp[0][0]>= MOD)temp[0][0] %= MOD; if(temp[0][1]>= MOD)temp[0][1] %= MOD; if(temp[1][0]>= MOD)temp[1][0] %= MOD; if(temp[1][1]>= MOD)temp[1][1] %= MOD; ret = temp; } if( pos == to ){ break; } pos = tree[pos][0]; } cout << ret[0][0] << " " << ret[0][1] << " "; cout << ret[1][0] << " " << ret[1][1] << '\n'; } } cout << flush; } #pragma region main signed main(void){ std::cin.tie(0); std::ios::sync_with_stdio(false); std::cout << std::fixed;//小数を10進数表示 cout << setprecision(16);//小数点以下の桁数を指定//coutとcerrで別 solve(); } #pragma endregion //main()