#include using namespace std; using ll = long long; using uint = unsigned int; using ull = unsigned long long; #define rep(i,n) for(int i=0;i=0;i--) #define per1(i,n) for(int i=int(n);i>0;i--) #define all(c) c.begin(),c.end() #define si(x) int(x.size()) #define pb push_back #define eb emplace_back #define fs first #define sc second template using V = vector; template using VV = vector>; template bool chmax(T& x, U y){ if(x bool chmin(T& x, U y){ if(y void mkuni(V& v){sort(all(v));v.erase(unique(all(v)),v.end());} template int lwb(const V& v, const T& a){return lower_bound(all(v),a) - v.begin();} template V Vec(size_t a) { return V(a); } template auto Vec(size_t a, Ts... ts) { return V(ts...))>(a, Vec(ts...)); } template ostream& operator<<(ostream& o,const pair &p){ return o<<"("< ostream& operator<<(ostream& o,const vector &vc){ o<<"{"; for(const T& v:vc) o< D divFloor(D a, D b){ return a / b - (((a ^ b) < 0 && a % b != 0) ? 1 : 0); } template D divCeil(D a, D b) { return a / b + (((a ^ b) > 0 && a % b != 0) ? 1 : 0); } uint solve(){ int N,M; cin >> N >> M; V A(N); rep(i,N) cin >> A[i]; // x^2 mod 2^m = (x mod 2^{m-1})^2 // (m >= 1) if(N <= 40){ rep(t,N-1){ int n = N-1-t; rep(i,n) A[i] = (A[i]+A[i+1])*(A[i]+A[i+1]); } return A[0]%(1< B(32); rep(i,32){ // [i,i+N-32] int n = N-32+1; rep(j,n) if((n&j)==j) B[i] += A[i+j]; } rep(t,31){ int n = 32-1-t; rep(i,n) B[i] = (B[i]+B[i+1])*(B[i]+B[i+1]); } return B[0]%(1<> T; while(T--) cout << solve() << endl; }