#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include /* template start */ #define rep(i, a, b) for (long long i = (a); (i) < (b); (i)++) #define all(i) i.begin(), i.end() #define debug(...) std::cerr << "[" << #__VA_ARGS__ << "]:", debug_out(__VA_ARGS__) void debug_out(){std::cerr< void debug_out(Head h,Tail... t){ std::cerr<<" "<0)std::cout<<" :"; debug_out(t...); } template std::ostream& operator<<(std::ostream& os, std::pair pa) { return os << pa.first << " " << pa.second; } template std::ostream& operator<<(std::ostream& os, std::vector vec) { for (std::size_t i = 0; i < vec.size(); i++)os << vec[i] << (i + 1 == vec.size() ? "" : " "); return os; } template inline bool chmax(T1& a,T2 b){return a inline bool chmin(T1& a,T2 b){return a>b && (a=b,true);} long long pow_mod(long long a, long long b, long long mod=-1) { if(b==0)return 1; if ((a == 0)||(mod!=-1&&(a+mod)%mod==0))return 0; long long x = 1; while (b > 0) { if (b & 1)x = (mod!=-1)?(x * a) % mod:x*a; a = (mod!=-1)?(a * a) % mod:a*a; b >>= 1; } return x; } template class modint { using u64 = std::uint_fast64_t; public: u64 a; constexpr modint(const u64 x = 0) noexcept : a(x % Modulus) {} constexpr u64 &value() noexcept { return a; } constexpr const u64 &value() const noexcept { return a; } constexpr modint operator+(const modint rhs) const noexcept { return modint(*this) += rhs; } constexpr modint operator-(const modint rhs) const noexcept { return modint(*this) -= rhs; } constexpr modint operator*(const modint rhs) const noexcept { return modint(*this) *= rhs; } constexpr modint operator/(const modint rhs) const noexcept { return modint(*this) /= rhs; } constexpr modint &operator+=(const modint rhs) noexcept { a += rhs.a; if (a >= Modulus) { a -= Modulus; } return *this; } constexpr modint &operator-=(const modint rhs) noexcept { if (a < rhs.a) { a += Modulus; } a -= rhs.a; return *this; } constexpr modint &operator*=(const modint rhs) noexcept { a = a * rhs.a % Modulus; return *this; } constexpr modint &operator/=(modint rhs) noexcept { u64 exp = Modulus - 2; while (exp) { if (exp % 2) { *this *= rhs; } rhs *= rhs; exp /= 2; } return *this; } }; /* template end */ //SWAG スライド区間上の演算ができる template class SWAG{ private: using F=std::function; F fn; std::stack front,back; std::vector front_fold,back_fold; T id; public: SWAG(T id_,F fn_):id(id_),fn(fn_),front_fold(1,id_),back_fold(1,id_){} void push(T val){ back.push(val); back_fold.push_back(fn(back_fold.back(),val)); } void pop(){ if(front.empty()){ assert(front_fold.size()==1); while(!back.empty()){ T temp=back.top(); back.pop(); back_fold.pop_back(); front.push(temp); front_fold.push_back(fn(temp,front_fold.back())); } } assert(!front.empty()); front.pop(); front_fold.pop_back(); } T fold_all(){ return fn(front_fold.back(),back_fold.back()); } }; using ll = long long; using P=std::pair; int main() { std::cin.tie(nullptr); std::ios::sync_with_stdio(false); ll n; std::cin>>n; std::vector a(n); rep(i,0,n)std::cin>>a[i]; a.push_back(1); SWAG swag(0ll,[](ll a,ll b){return std::gcd(a,b);}); swag.push(a[0]); ll right=0; ll ans=0; rep(i,0,n){ while(swag.fold_all()!=1){ right++; swag.push(a[right]); } ans+=n-right; if(right==i){ right++; swag.push(a[right]); } swag.pop(); } std::cout<