n=int(input())
a=list(map(int,input().split()))
for i in range(1,n):
  a[i]^=a[i-1]
M=10**9+7
q=[0]*(n+1)
d={}
for i in range(n-1):
  if a[i] in d:
    q[i]=q[i-1]-q[d[a[i]]-1]
  else:
    q[i]=q[i-1]+1
  q[i]%=M
  q[i]+=q[i-1]
  d[a[i]]=i
print((q[n-2]+1)%M)