#include using namespace std; using ll = long long; #define SPEED cin.tie(0);ios::sync_with_stdio(false); 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])); } }; int ctz(long long a) { return __builtin_ctzll(a); } template T gcd(T a, T b) { a = abs(a), b = abs(b); if (a == 0) return b; int k = ctz(a | b); a >>= ctz(a); while (b) { b >>= ctz(b); if (a > b) swap(a, b); b -= a; } return a << k; } template struct nodeGCD { typeNode unitNode = 0; typeNode funcNode(typeNode l,typeNode r){return gcd(l,r);} }; // solution by binary search in arbitary range on disjn sparse table int main() { SPEED ll N; cin >> N; assert(1 <= N && N <= 500000); vector A(N+1,1); for(int i = 0; i < N; ++i) cin >> A[i]; for(int i = 0; i < N; ++i) assert(1 <= A[i] && A[i] <= 1000000000000000000LL); DisjointSparseTable> dst(A); ll ans=0; int r=1; for(int i=0; i