#line 1 "main.cpp" #define PROBLEM "https://yukicoder.me/problems/no/1036" #include #include #line 2 "/home/user/GitHub/competitive-programming-library/utils/macros.hpp" #define REP(i, n) for (int i = 0; (i) < (int)(n); ++ (i)) #define REP3(i, m, n) for (int i = (m); (i) < (int)(n); ++ (i)) #define REP_R(i, n) for (int i = (int)(n) - 1; (i) >= 0; -- (i)) #define REP3R(i, m, n) for (int i = (int)(n) - 1; (i) >= (int)(m); -- (i)) #define ALL(x) std::begin(x), std::end(x) #line 2 "/home/user/GitHub/competitive-programming-library/data_structure/sparse_table.hpp" #include #line 5 "/home/user/GitHub/competitive-programming-library/data_structure/sparse_table.hpp" /** * @brief Sparse Table (idempotent monoid) * @note the unit is required just for convenience * @note $O(N \log N)$ space */ template struct sparse_table { typedef typename IdempotentMonoid::value_type value_type; std::vector > table; IdempotentMonoid mon; sparse_table() = default; /** * @note $O(N \log N)$ time */ template sparse_table(InputIterator first, InputIterator last, const IdempotentMonoid & mon_ = IdempotentMonoid()) : mon(mon_) { table.emplace_back(first, last); int n = table[0].size(); int log_n = 32 - __builtin_clz(n); table.resize(log_n, std::vector(n)); REP (k, log_n - 1) { REP (i, n) { table[k + 1][i] = i + (1ll << k) < n ? mon.mult(table[k][i], table[k][i + (1ll << k)]) : table[k][i]; } } } /** * @note $O(1)$ */ value_type range_get(int l, int r) const { if (l == r) return mon.unit(); // if there is no unit, remove this line assert (0 <= l and l < r and r <= (int)table[0].size()); int k = 31 - __builtin_clz(r - l); // log2 return mon.mult(table[k][l], table[k][r - (1ll << k)]); } }; #line 2 "/home/user/GitHub/competitive-programming-library/number/gcd.hpp" #include /** * @note if arguments are negative, the result may be negative */ template T gcd(T a, T b) { while (a) { b %= a; std::swap(a, b); } return b; } template T lcm(T a, T b) { return a / gcd(a, b) * b; } #line 3 "/home/user/GitHub/competitive-programming-library/monoids/gcd.hpp" /** * @note a semilattice */ template struct gcd_monoid { typedef Integer value_type; Integer unit() const { return 0; } Integer mult(Integer a, Integer b) const { return gcd(a, b); } }; #line 7 "main.cpp" long long solve(int n, const std::vector & a) { sparse_table > table(ALL(a)); int r = 0; long long ans = 0; REP (l, n) { while (r < n and table.range_get(l, r) != 1) { ++ r; } if (r == n and table.range_get(l, r) != 1) { break; } ans += n - r + 1; } return ans; } int main() { int n; scanf("%d", &n); std::vector a(n); REP (i, n) { scanf("%lld", &a[i]); } printf("%lld\n", solve(n, a)); return 0; }