#pragma GCC optimize("Ofast") #include using namespace std; using ll = long long int ; using ld = long double ; using P = pair; using Graph= vector>; struct edge{ll to ; ll cost ;} ; using graph =vector> ; #define rep(i,n) for (ll i=0; i < (n); ++i) #define rep2(i,n,m) for(ll i=n;i<=m;i++) #define rep3(i,n,m) for(ll i=n;i>=m;i--) #define pb push_back #define eb emplace_back #define ppb pop_back #define mpa make_pair #define fi first #define se second #define set20 cout<= mod) x -= mod; return *this; } mint& operator-=(const mint a) { if ((x += mod-a.x) >= mod) x -= mod; return *this; } mint& operator*=(const mint a) { (x *= a.x) %= mod; return *this;} mint operator+(const mint a) const { return mint(*this) += a;} mint operator-(const mint a) const { return mint(*this) -= a;} mint operator*(const mint a) const { return mint(*this) *= a;} mint pow(ll t) const { if (!t) return 1; mint a = pow(t>>1); a *= a; if (t&1) a *= *this; return a; } // for prime mod mint inv() const { return pow(mod-2);} mint& operator/=(const mint a) { return *this *= a.inv();} mint operator/(const mint a) const { return mint(*this) /= a;} }; istream& operator>>(istream& is, const mint& a) { return is >> a.x;} ostream& operator<<(ostream& os, const mint& a) { return os << a.x;} //昆布 struct combination { vector fact, ifact; combination(int n):fact(n+1),ifact(n+1) { assert(n < mod); //任意modではここ消すcombmain内に fact[0] = 1; for (int i = 1; i <= n; ++i) fact[i] = fact[i-1]*i; ifact[n] = fact[n].inv(); for (int i = n; i >= 1; --i) ifact[i-1] = ifact[i]*i; } mint operator()(int n, int k) { if (k < 0 || k > n) return 0; return fact[n]*ifact[k]*ifact[n-k]; } mint p(int n,int k){ return fact[n]*ifact[n-k] ; //kは個数 } } co(2000005); mint modpow(ll a,ll b){ if(b==0) return 1 ; mint c= modpow(a,b/2) ; if(b%2==1) return c*c*a ; else return c*c ; } mint mmodpow(mint a,ll b){ if(b==0) return 1ll ; mint c=mmodpow(a,(b/2)) ; if(b%2==1) return c*c*a ; else return c*c ; } mint komb(ll n,ll m){ mint x=1 ;mint y=1 ; rep(i,m){ x*= n-i ; y*= i+1 ; } return x/y ; } map factor(ll n){ //素因数とオーダーをマップで管理 map ord ; for(ll i=2;i*i<=n;i++){ if(n%i==0){ int res=0; while(n%i==0){ n/=i; res++; } ord[i]=res; } } if(n!=1) ord[n]++; return ord ; } vector alldiv(ll n){ vector ans; if(n<0) n*=-1; for(ll i=1;i*i<=n;++i){ if(n%i==0){ if(i*i==n) ans.pb(i); else{ans.pb(i);ans.pb(n/i);} } } return ans; } struct UnionFind { vector d; UnionFind(int n=0): d(n,-1) {} int find(int x) { if (d[x] < 0) return x; return d[x] = find(d[x]); } bool unite(int x, int y) { x = find(x); y = find(y); if (x == y) return false; if (d[x] > d[y]) swap(x,y); d[x] += d[y]; d[y] = x; return true; } bool same(int x, int y) { return find(x) == find(y);} int size(int x) { return -d[find(x)];} }; // sum(x) x以下の和 // sum(a,b) a以上b以下の和 template struct BIT { int n; vector d; BIT(int n=0):n(n),d(n+1) {} void add(int i, T x=1) { //x=1ならsumは個数のカウント for (i++; i <= n; i += i&-i) { d[i] += x; } } T sum(int i) { T x = 0; for (i++; i; i -= i&-i) { x += d[i]; } return x; } T sum(int i,int j) { if(i>0) return sum(j)-sum(i-1); else return sum(j); } }; // 返り値: a と b の最大公約数 // ax + by = gcd(a, b) を満たす (x, y) が格納される long long extGCD(long long a, long long b, long long &x, long long &y) { if (b == 0) { x = 1; y = 0; return a; } long long d = extGCD(b, a%b, y, x); y -= a/b * x; return d; } /* RMQ:[0,n-1] について、区間ごとの最大値を管理する構造体 update(i,x): i 番目の要素を x に更新。O(log(n)) query(a,b): [a,b) での最大の要素を取得。O(log(n)) */ template struct RMQ { const T INF = numeric_limits::max(); int n; // 葉の数 vector dat; // 完全二分木の配列 RMQ(int n_) : n(), dat(n_ * 4, -INF) { // 葉の数は 2^x の形 int x = 1; while (n_ > x) { x *= 2; } n = x; } void update(int i, T x) { i += n - 1; dat[i] = x; while (i > 0) { i = (i - 1) / 2; // parent dat[i] = max(dat[i * 2 + 1], dat[i * 2 + 2]); } } // the minimum element of [a,b) T query(int a, int b) { return query_sub(a, b, 0, 0, n); } T query_sub(int a, int b, int k, int l, int r) { if (r <= a || b <= l) { return -INF; } else if (a <= l && r <= b) { return dat[k]; } else { T vl = query_sub(a, b, k * 2 + 1, l, (l + r) / 2); T vr = query_sub(a, b, k * 2 + 2, (l + r) / 2, r); return max(vl, vr); } } }; ll p(ll x,ll y,ll m){ if(y==0) return 1; else{ ll ans=p(x,y/2,m); ans=(ans*ans)%m; if(y%2==0) return ans; else return (ans*x)%m ; } } pair dp[200005][31]; int main(){ ios::sync_with_stdio(false) ; cin.tie(nullptr) ; ll n;cin>>n; vector A(n); rep(i,n) cin>>A[i]; Graph g(n); rep(i,n){ g[i].pb(i); if(i==0) continue; vector x=g[i-1]; ll k=x.size(); ll now=A[i]; rep(j,k){ if((now|A[x[j]])!=now){ now|=A[x[j]]; g[i].pb(x[j]); } } } rep(i,n) g[i].pb(-1); rep(i,n) reverse(g[i].begin(),g[i].end()); dp[0][0]={0,1ll}; rep(i,n-1){ ll k=g[i+1].size(); vector ans(k-1); rep(j,31){ ll x=dp[i][j].fi; mint y=dp[i][j].se; rep(w,k-1){ if(x<=g[i+1][w+1]) ans[w]+=y; } } rep(j,k-1) dp[i+1][j]={g[i+1][j]+1,ans[j]}; } mint sum=0ll; rep(i,31) sum+=dp[n-1][i].se; cout<