#include #include namespace my{ using namespace std; #define done(...) return pp(__VA_ARGS__) #define LL(...) ll __VA_ARGS__;lin(__VA_ARGS__) #define VRD(T,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__);i=i##stop;i-=i##step) #define fe(a,i,...) for(auto&&__VA_OPT__([)i __VA_OPT__(,__VA_ARGS__]):a) #define single_testcase void solve();}int main(){my::io();my::solve();}namespace my{ void io(){cin.tie(nullptr)->sync_with_stdio(0);cout<auto range(bool s,A...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 nl=10; constexpr char sp=32; constexpr ll log2_floor(auto x){ll r=0;while(x>1)x>>=1,++r;return r;} auto ceil(auto x,auto y){if(y<0)x=-x,y=-y;return x<=0?x/y:(x-1)/y+1;} auto min(const auto&...a){return min(initializer_list>{a...});} templatestruct pair{ A a;B b; pair()=default; pair(A a,B b):a(a),b(b){} pair(const std::pair&p):a(p.first),b(p.second){} auto operator<=>(const pair&)const=default; friend ostream&operator<<(ostream&o,const pair&p){return o<ostream&operator<<(ostream&o,const std::pair&p){return 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<?nl:sp);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++(){fe(*this,e)++e;return*this;} vec&operator--(){fe(*this,e)--e;return*this;} auto scan(const auto&f)const{pair,bool>r{};fe(*this,e)if constexpr(!vectorial)r.b?f(r.a,e),r:r={e,1};else if(auto s=e.scan(f);s.b)r.b?f(r.a,s.a),r:r=s;return r;} auto min()const{return scan([](auto&a,const auto&b){a>b?a=b:0;;}).a;} }; 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<>block_size_exponent;} static constexpr size_t bit_index(size_t i){return i&(block_size-1);} size_t n; vecbits; size_t size()const{return n;} bitvec(size_t n):n(n),bits(ceil(n,block_size)){} bitvec(const string&s):n(s.size()),bits(ceil(s.size(),block_size)){fo(i,n)bits[block_index(i)]|=u64(s[i]-'0')<{x}.swap(bits);return*this;} struct proxy{ decltype(bits)&bits_proxy; size_t I; proxy(decltype(bits)&bits_proxy,size_t I):bits_proxy(bits_proxy),I(I){} auto&operator=(bool f){ bits_proxy[block_index(I)]&=~(u64(1)<auto&operator=(X x){return this->operator=(static_cast(x&1));} auto&operator=(const proxy&p){return this->operator=(static_cast(p));} operator bool()const{return bits_proxy[block_index(I)]>>bit_index(I)&1;} }; decltype(auto)operator[](size_t i){return proxy{bits,i};} decltype(auto)operator[](size_t i)const{return bits[block_index(i)]>>bit_index(i)&1;} auto&operator&=(const bitvec&b){fo(i,bits.size())bits[i]&=b.bits[i];return*this;} auto&operator|=(const bitvec&b){fo(i,bits.size())bits[i]|=b.bits[i];return*this;} auto&operator^=(const bitvec&b){fo(i,bits.size())bits[i]^=b.bits[i];return*this;} auto&operator<<=(int k){ int q=k/block_size,r=k-q*block_size; if(r==0){ of(i,bits.size(),q)bits[i]=bits[i-q]; fo(i,q)bits[i]=0; }else{ of(i,bits.size(),q+1)bits[i]=bits[i-q]<>(block_size-r); bits[q]=(bits[0]<>=(int k){ int q=k/block_size,r=k-q*block_size; if(r==0){ fo(i,bits.size()-q)bits[i]=bits[i+q]; fo(i,bits.size()-q,bits.size())bits[i]=0; }else{ fo(i,bits.size()-q-1)bits[i]=bits[i+q]>>r|bits[i+q+1]<<(block_size-r); bits[bits.size()-q-1]=bits[bits.size()-1]>>r; fo(i,bits.size()-q,bits.size())bits[i]=0; } return*this; } auto operator&(const bitvec&b)const{return bitvec{*this}&=b;} auto operator|(const bitvec&b)const{return bitvec{*this}|=b;} auto operator^(const bitvec&b)const{return bitvec{*this}^=b;} auto operator<<(int k)const{return bitvec{*this}<<=k;} auto operator>>(int k)const{return bitvec{*this}>>=k;} ll count()const{ll r=0;fe(bits,e)r+=std::popcount(e);return r;} string to_string()const{string s{};fo(i,size())s+=(*this)[i]+'0';return s;} friend istream&operator>>(istream&i,bitvec&b){string s;i>>s;b=s;return i;} friend ostream&operator<<(ostream&o,const bitvec&b){return o<mem(block_num,bitvec(K+1)); bitvec cu(K+1); cu[0]=1; fo(i,N){ if(i%block_size==0)mem[i/block_size]=cu; cu|=cu<dp(block_size,bitvec(K+1)); ll ans=0; cu=bitvec(K+1); cu[K]=1; of(i,block_num){ ll L=min(block_size,N-i*block_size); dp[0]=mem[i]; fo(j,L-1)dp[j+1]=dp[j]<>a[i*block_size+j]; } } pp(ans); }}