結果
問題 | No.1036 Make One With GCD 2 |
ユーザー |
![]() |
提出日時 | 2020-04-24 22:32:00 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 334 ms / 2,000 ms |
コード長 | 2,031 bytes |
コンパイル時間 | 2,408 ms |
コンパイル使用メモリ | 198,956 KB |
最終ジャッジ日時 | 2025-01-10 00:06:37 |
ジャッジサーバーID (参考情報) |
judge3 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 41 |
ソースコード
#include<bits/stdc++.h>using namespace std;#define ALL(x) x.begin(),x.end()#define rep(i,n) for(int i=0;i<(n);i++)#define debug(v) cout<<#v<<":";for(auto x:v){cout<<x<<' ';}cout<<endl;#define INF 1000000000#define mod 1000000007using ll=long long;const ll LINF=1001002003004005006ll;int dx[]={1,0,-1,0};int dy[]={0,1,0,-1};// ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}template<class T>bool chmax(T &a,const T &b){if(a<b){a=b;return true;}return false;}template<class T>bool chmin(T &a,const T &b){if(b<a){a=b;return true;}return false;}template<typename Monoid>struct SegmentTree{using F=function<Monoid(Monoid,Monoid)>;int sz;vector<Monoid> seg;const F f;const Monoid gen;SegmentTree(int n,const F f,const Monoid &gen):f(f),gen(gen){sz=1;while(sz<n)sz<<=1;seg.assign(2*sz,gen);}void set(int k,const Monoid &x){seg[k+sz]=x;}void build(){for(int k=sz-1;k>0;k--) seg[k]=f(seg[2*k],seg[2*k+1]);}void update(int k,const Monoid &x){k+=sz;seg[k]=x;while(k>>=1) seg[k]=f(seg[2*k],seg[2*k+1]);}// [a,b)Monoid query(int a,int b){Monoid L=gen,R=gen;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 {return seg[k+sz];}};ll f(ll a,ll b){return gcd(a,b);}signed main(){cin.tie(0);ios::sync_with_stdio(0);ll n;cin>>n;vector<ll> a(n);SegmentTree<ll> seg(n,f,0);rep(i,n){cin>>a[i];seg.set(i,a[i]);}seg.build();ll r=0,res=0,x=0;for(ll l=0;l<n;l++){while(r<n and gcd(a[r],x)!=1){x=gcd(x,a[r]);r++;}res+=r-l;if(r==l){r++;x=0;}else{x=seg.query(l+1,r);}}cout<<n*(n+1)/2-res<<endl;return 0;}