#include #include using namespace std; namespace my{ #define LL(...) ll __VA_ARGS__;lin(__VA_ARGS__) #define VL(n,...) vec__VA_ARGS__;setsize({n},__VA_ARGS__);lin(__VA_ARGS__) #define FO(n) for(ll ij=n;ij--;) #define FOR(i,...) for(auto[i,i##stop,i##step]=range(0,__VA_ARGS__);isync_with_stdio(0);cout<>(istream&i,ulll&x){ull t;i>>t;x=t;return i;} ostream&operator<<(ostream&o,const ulll&x){return(x<10?o:o<>(istream&i,lll&x){ll t;i>>t;x=t;return i;} ostream&operator<<(ostream&o,const lll&x){return o<0?x:-x);} constexpr auto range(bool s,auto...a){arrayr{0,0,1};ll I=0;((r[I++]=a),...);if(!s&&I==1)swap(r[0],r[1]);r[0]-=s;return r;} constexpr char newline=10; constexpr char space=32; lll pow(lll x,ll n){assert(n>=0);lll r=1;while(n)n&1?r*=x:r,x*=x,n>>=1;return r;} templateostream&operator<<(ostream&o,const map&m){fe(m,e)o<concept vectorial=is_base_of_v,V>; templatestruct vec_attr{using core_type=T;static constexpr int d=0;}; templatestruct vec_attr{using core_type=typename vec_attr::core_type;static constexpr int d=vec_attr::d+1;}; templateusing core_t=vec_attr::core_type; templateistream&operator>>(istream&i,vector&v){fe(v,e)i>>e;return i;} templateostream&operator<<(ostream&o,const vector&v){fe(v,e)o<?newline:space);return o;} templatestruct vec:vector{ using vector::vector; vec(const vector&v){vector::operator=(v);} vec&operator^=(const vec&u){this->insert(this->end(),u.begin(),u.end());return*this;} vec operator^(const vec&u)const{return vec{*this}^=u;} vec&operator+=(const vec&u){vec&v=*this;fo(i,v.size())v[i]+=u[i];return v;} vec&operator-=(const vec&u){vec&v=*this;fo(i,v.size())v[i]-=u[i];return v;} vec operator+(const vec&u)const{return vec{*this}+=u;} vec operator-(const vec&u)const{return vec{*this}-=u;} vec&operator++(){fe(*this,e)++e;return*this;} vec&operator--(){fe(*this,e)--e;return*this;} vec operator-()const{vec v=*this;fe(v,e)e=-e;return v;} auto upper_bound(const V&x)const{return std::upper_bound(this->begin(),this->end(),x);} }; templateauto make_vec(const ll(&s)[n],T x={}){if constexpr(n==i+1)return vec(s[i],x);else{auto X=make_vec(s,x);return vec(s[i],X);}} templatevoid setsize(const ll(&l)[n],A&...a){((a=make_vec(l,core_t())),...);} void lin(auto&...a){(cin>>...>>a);} templatevoid pp(const auto&...a){ll n=sizeof...(a);((cout<0,c)),...);cout<1); if(a>b)return 0LL; auto r=a; ll k=1; auto b_a=b/a; while(r<=b_a)r*=a,k++; return k; } templateauto zeta_on_map(mapm){ if(m.size()==0)return map{}; for(auto p=m.begin();next(p)!=m.end();p++)next(p)->second+=p->second; return m; } single_testcase void solve(){ LL(N,S); VL(N,a); mapdp; dp[0]=1; fo(i,N/2){ mapndp; fe(dp,s,x)fo(k,1,max_exponent(a[i],S-s)+1)ndp[s+pow(a[i],k)]+=x; swap(ndp,dp); } dp=zeta_on_map(dp); mapdp2; dp2[0]=1; fo(i,N/2,N){ mapndp; fe(dp2,s,x)fo(k,1,max_exponent(a[i],S-s)+1)ndp[s+pow(a[i],k)]+=x; swap(ndp,dp2); } ll ans=0; fe(dp2,s,x){ auto p=dp.upper_bound(S-s); if(p==dp.begin())continue; ans+=prev(p)->second*x; } pp(ans); }}