#include #define rep(i,a) for(int i=(int)0;i<(int)a;++i) #define rrep(i,a) for(int i=(int)a-1;i>=0;--i) #define REP(i,a,b) for(int i=(int)a;i<(int)b;++i) #define RREP(i,a,b) for(int i=(int)a-1;i>=b;--i) #define pb push_back #define eb emplace_back #define all(x) x.begin(),x.end() #define rall(x) x.rbegin(),x.rend() typedef std::vector vi; typedef std::vector> vvi; typedef std::vector vl; typedef std::vector> vvl; #define out(x) cout< inline bool chmin(T& a, T b) { if (a > b) { a = b; return true; } return false; } template inline bool chmax(T& a, T b) { if (a < b) { a = b; return true; } return false; } ll gcd(ll n, ll m) { ll tmp; while (m!=0) { tmp = n % m; n = m; m = tmp; } return n; } ll lcm(ll n, ll m) { return abs(n) / gcd(n, m)*abs(m);//gl=xy } using namespace std; template struct Segmenttree{//0-indexed using F=function; int sz; vectorseg; const F f;//2つの区間の要素をマージする二項演算 const Monoid unit;//単位元 Segmenttree(int n,const F f,const Monoid &unit):f(f),unit(unit){ sz=1; while(sz0;k--)seg[k]=f(seg[2*k+0],seg[2*k+1]); } void update(int k,const Monoid &x){//k番目の要素をxに変更する k+=sz; seg[k]=x; while(k>>=1){ seg[k]=f(seg[2*k+0],seg[2*k+1]); } } Monoid query(int a, int b) {//区間[a,b)に対して二項演算した結果を返す Monoid L = unit, R = unit; for(a += sz, b += sz; a < b; a >>= 1, b >>= 1) { if(a & 1) L = f(L, seg[a++]); if(b & 1) R = f(seg[--b], R); } return f(L, R); } Monoid operator[](const int &k) const {//k番目の要素を返す return seg[k + sz]; } }; void solve() { ll n; cin >> n; vector a(n); rep(i, n) cin >> a[i]; auto f = [](ll a, ll b) { return gcd(a, b); }; Segmenttree sg(n,f,0LL); rep(i,n)sg.update(i, a[i]); ll x = 0; ll ans = 0; ll pos = 0; for (ll i = 0; i < n; ++i) { x = sg.query(i, pos); for (ll j = pos; j < n; ++j) { x = gcd(x, a[j]); if (x == 1) { ans += n - pos; break; } ++pos; } pos = max(pos, i+1); } cout << ans << endl; } int main(){ ios::sync_with_stdio(false); cin.tie(0); cout<