#include using namespace std; using ll = long long; using ull = unsigned long long; #define ALL(obj) (obj).begin(),(obj).end() #define SPEED cin.tie(0);ios::sync_with_stdio(false); template using PQ = priority_queue; template using PQR = priority_queue,greater>; constexpr long long MOD = (long long)1e9 + 7; constexpr long long MOD2 = 998244353; constexpr long long HIGHINF = (long long)1e18; constexpr long long LOWINF = (long long)1e15; constexpr long double PI = 3.1415926535897932384626433; template vector multivector(size_t N,T init){return vector(N,init);} template auto multivector(size_t N,T... t){return vector(N,multivector(t...));} template void corner(bool flg, T hoge) {if (flg) {cout << hoge << endl; exit(0);}} template ostream &operator<<(ostream &o, const map&obj) {o << "{"; for (auto &x : obj) o << " {" << x.first << " : " << x.second << "}" << ","; o << " }"; return o;} template ostream &operator<<(ostream &o, const set&obj) {o << "{"; for (auto itr = obj.begin(); itr != obj.end(); ++itr) o << (itr != obj.begin() ? ", " : "") << *itr; o << "}"; return o;} template ostream &operator<<(ostream &o, const multiset&obj) {o << "{"; for (auto itr = obj.begin(); itr != obj.end(); ++itr) o << (itr != obj.begin() ? ", " : "") << *itr; o << "}"; return o;} template ostream &operator<<(ostream &o, const vector&obj) {o << "{"; for (int i = 0; i < (int)obj.size(); ++i)o << (i > 0 ? ", " : "") << obj[i]; o << "}"; return o;} template ostream &operator<<(ostream &o, const pair&obj) {o << "{" << obj.first << ", " << obj.second << "}"; return o;} void print(void) {cout << endl;} template void print(Head&& head) {cout << head;print();} template void print(Head&& head, Tail&&... tail) {cout << head << " ";print(forward(tail)...);} template void chmax(T& a, const T b){a=max(a,b);} template void chmin(T& a, const T b){a=min(a,b);} void YN(bool flg) {cout << (flg ? "YES" : "NO") << endl;} void Yn(bool flg) {cout << (flg ? "Yes" : "No") << endl;} void yn(bool flg) {cout << (flg ? "yes" : "no") << endl;} template class DisjointSparseTable{ public: Operator Op; using typeNode = decltype(Op.unitNode); size_t depth; size_t length; vector> node; vector msb; DisjointSparseTable(const vector& vec) { for(depth = 0;(1<>j) msb[i] = j; //init value node.resize(depth,vector(length,Op.unitNode)); for(int i = 0; i < vec.size(); ++i) node[0][i] = vec[i]; for(int i = 1; i < depth; ++i) { for(int r = (1<r||l<0||length<=r) ? Op.unitNode: (l==r ? node[0][l] : Op.funcNode(node[msb[l^r]][l],node[msb[l^r]][r])); } }; //一点更新 区間最小 template struct nodeGCD { typeNode unitNode = 0; typeNode funcNode(typeNode l,typeNode r){return ((r == 0) ? l : funcNode(r, l % r));} }; int main() { SPEED int N; cin >> N; vector A(N+1,1); for(int i = 0; i < N; ++i) cin >> A[i]; DisjointSparseTable> dst(A); ll ans = 0; for(int i = 0; i < N; ++i) { if(A[i] == 1){ ans += (N-i); } else{ int ok = N,ng = i,md; while(ok-ng>1){ md = (ok+ng)/2; (dst.get(i,md+1)==1?ok:ng)=md; } ans += (N-ok); } } cout << ans << endl; return 0; }